unfoldr (version 0.2)

simSphereSystem: Simulation of sphere system

Description

The function simulates a Poisson sphere system of intensity lam where each sphere center is uniformly distributed in a box. The function returns a list of spheres with elements id, center and radius r.

Usage

simSphereSystem(theta, rdist, box = list(c(0, 1)), pl = 0)

Arguments

theta

list of named arguments lam,radii

rdist

string, radii random generating function name

box

simualtion box

pl

print level

Value

• list of class spheres if pl>100 or empty list

Details

Any random generating function, passed as a name, for the radii distribution is accepted as long as the formal function parameter names match the actual parameter names exactly as defined in the parameter list theta$radii.

The simulation box is of type list. The vector arguments correspond to the lower and upper points in x,y and z direction. If box has only one element, i.e. list(c(0,1), the same extent is used for the other dimensions. The argument pl denotes the print level of information during simulation. Currently, only pl=0 for no output and pl>100 is implemented. Argument cond$rdist is of type string naming the (user defined) radii random generating function.

Examples

theta <- list("lam"=100,"radii"=list("meanlog"=-2.5,"sdlog"=0.2))
 S <- simSphereSystem(theta=theta,rdist="rlnorm",pl=101)