R help - Dec 2005 - O-ring statistic

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

Thanks a lot for your reply

I'll look into that and let you know if I have further questions

Rainer


Adrian Baddeley wrote:
> 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
> 
> 
-- 
--
Rainer M. Krug, Dipl. Phys. (Germany), MSc Conservation Biology (UCT)

Department of Conservation Ecology
University of Stellenbosch
Matieland 7602
South Africa

Tel:        +27 - (0)72 808 2975 (w)
Fax:        +27 - (0)21 808 3304
Cell:        +27 - (0)83 9479 042

email:    RKrug at sun.ac.za
           Rainer at krugs.de