
Cluster a nd Cox Models 469
> mu <- as.im(function(x,y){ exp(2 * x + 1) }, owin())
> X <- rMatClust(10, 0.05, mu)
The result is sh own in Figure 12.19.
Figure 12.19. Example of a simulated inhomogeneous Matérn cluster process with reference inten-
sity
µ
(x,y) = exp(2x + 1 ) in the unit square.
The function rNeymanScott can simulate a general Neyman-Scott process. The parents are
Poisson with inten sity kappa, where kappa is given as a numeric value, a function, or a pixel im-
age. The cluster mechanism is specified by th e a rgument rcluster which can be either a function
that generates ran dom point patterns, or a list(mu, f) where mu is the mean number ...