simMatrix {FractalParameterEstimation}

Simulating 3x3 matrix with binomial random variables

Description

This function simulates a matrix of size 3x3 and fills these with bernoulli random variables under a probability of p. An optional parameter sierp is added in order to set the center of the matrix to 0 as for the general Sierpinski-Carpet.

Usage

simMatrix(p, sierp)

Arguments

p

A numeric value between 0 and 1 indicating the probability of success (0 or 1) for the Bernoulli random variables of the matrix.

sierp

An optional logical parameter: if TRUE then the center of the matrix is automatically set to 0 as for the general Sierpinski-Carpet, else also a Bernoulli random variable is simulated.

Value

The function creates a 3x3 matrix containing simulated zeros and ones from Bernoulli distribution under given probability p.

Note

Shall not be called or changed by user.

Author(s)

Philipp Hermann; Jozef Kiselak; Milan Stehlik\ [e-mail]; [e-mail]; [e-mail]

References

Hermann, P., Mrkvicka, T., Mattfeldt, T., Minarova, M., Helisova, K., Nicolis, O., Wartner, F., and Stehlik, M. (2015). Fractal and stochastic Geometry Inference for Breast Cancer: a Case Study with Random Fractal Models and Quermass-Interaction Process. Statistics in Medicine, 34(18), 2636-2661. doi: 10.1002/sim.6497.

Examples

simMatrix(p = 0.3, sierp = TRUE)
simMatrix(p = 0.8, sierp = FALSE)