[R-sig-Geo] Kest (spatstat) and "r" values

Rolf Turner r.turner at auckland.ac.nz
Thu Oct 23 03:07:14 CEST 2014


On 23/10/14 08:59, Marcelino de la Cruz wrote:
> Maybe the information related to r you are looking for is in the Warning
> section:
>
> "Warnings
>
> The estimator of K(r) is approximately unbiased for each fixed r. Bias
> increases with r and depends on the window geometry. For a rectangular
> window it is prudent to restrict the r values to a maximum of 1/4 of the
> smaller side length of the rectangle".
>
> Maybe I'm wrong but I think that  Kest computest K(r) by default for r
> up to 1/4 of the shorter side of the *enclosing rectangle* of your window.

That is indeed correct, modulo a possible adjustment (which I believe 
rarely gets applied in practice) in terms of the estimated intensity of
the process.  See the code for rmax.rule() if you are interested in the 
details.

cheers,

Rolf

-- 
Rolf Turner
Technical Editor ANZJS



More information about the R-sig-Geo mailing list