Jeff Fang
2011-Feb-27 07:00 UTC
[R] Problem in simulating envelopes with R package "Spatstat"
Hi all,
While useing the R package "Spatstat" to detect the spatial point
pattern of
my data, I met a problem. When I computes simulation envelopes with a
fitted point process model: Poisson cluster (Thomas) process
(kappa=2.751010e-05; sigma=5.634274e+01; mu=4.943639e+02), I cannot get the
high values of the envelope, which show "NA" , but the
"mmean"
values computed by averaging simulated values are available, this result
confused me. Is it a bug in this package or only caused by my data? How can
I solve this issue?
Thanks a lot
#my code is as follow:
>dist = seq(0, 12.5, 0.1)
>species.c = kppm(species, ~ 1, "Thomas")
>species.cenv = envelope(species.c, Kest, r=dist, nsim=999, nrank=5)
Jeff
[[alternative HTML version deleted]]