[R] O-ring statistic
Adrian Baddeley
adrian at maths.uwa.edu.au
Fri Dec 9 03:39:44 CET 2005
Rainer M Krug writes:
> Thorsten Wiegand used in his paper Wiegand T., and K. A. Moloney 2004.
> Rings, circles and null-models for point pattern analysis in ecology.
> Oikos 104: 209-229 a statistic he called O-Ring statistic which is
> similar to Ripley's K, only that it uses rings instead of circles.
>
> http://www.oesa.ufz.de/towi/towi_programita.html#ring
>
> Is this statistic included in one of the packages in R?
This kind of functionality is available in the R package `spatstat'
(available on CRAN or from www.spatstat.org)
According to the cited website, the O-ring statistic is a rescaled
version of the pair correlation function between two types of points:
O_12(r) = lambda_2 g_12(r)
In spatstat, pair correlation functions are computed by the function
'pcf'. To estimate the cross-type pair correlation function,
you do something like
pcf(Kcross(X, 1, 2))
where X is a marked point pattern containing points of types 1 and 2.
To estimate the intensity of type 2 points you use summary(X).
Here's an example for the bivariate point pattern dataset 'amacrine'
provided in the spatstat package. The dataset has points of two types
labelled "on" and "off".
data(amacrine)
K12 <- Kcross(amacrine, "on", "off")
g12 <- pcf(K12, method="d", spar=0.7)
lambda2 <- summary(amacrine)$marks["off","intensity"]
Oring <- eval.fv(lambda2 * g12)
plot(Oring, ylab="Oring(r)")
regards
Adrian Baddeley
More information about the R-help
mailing list