[R-sig-Geo] test of CSR in 3D

[Hamid Ghorbani]

Dear list,
let L be the L function for a point pattern in 3D.
using spatstat I need a plot of \hat{L(r)} - r vs. r along with the
corresponding envelopes for a simulated pattern.

The K function have to be estimated for  r's  in [0,rmax] with "nrval"
equal to let say 1001 (nraval is the number of values of r for which
K3(r) will be estimated).

I use the following codes. using the ccommand env$r shows the rmax and
nrval and  has not taken in to account. Furthermore the y axis is
\hat{L(r)} not \hat{L(r)} - r.

Any correction please. many thanks in advance.

Hamid Gh.
=============
lambda<-470
r0<-73
X<-runifpoint3(lambda, domain = box3(c(150,350),c(150,350),c(300,400)))

env = envelope.pp3(X, fun = K3est,r=seq(0,r0, 1001), nsim = 99,
correction = "iso",transform=expression((3*./(4*pi))^(1/3)))
plot(env)
env$r