FractalParameterEstimation (version 1.1.2)

simMatrix: 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.

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

# NOT RUN {
simMatrix(p = 0.3, sierp = TRUE)
simMatrix(p = 0.8, sierp = FALSE)
# }