
Gibbs Mo dels 537
Now we compute the envelopes of the L-function based on simulations from the fitted model:
> swedenv <- envelope(swedfit, Lest, nsim=39, global=TRUE,
savepatterns=TRUE)
where we set savepatterns=TRUE to save the simulated point patterns for re-use in other calcu-
lations. The centred L-function is plotted in Figure 13.29 using the plo t formula.
~
. - r. This
does not suggest significant departure from the fitted mod el.
0 5 10 15 20
−4 −2 0 1 2
r (one unit = 0.1 metres)
L
(
r
)
− r
Figure 13.29. Centered L- function for the copper data (solid black line) and global envelopes of
the same function from 39 simulations of the fitted inhomogeneo u