29.01.2015 Views

Analysing spatial point patterns in R - CSIRO

Analysing spatial point patterns in R - CSIRO

Analysing spatial point patterns in R - CSIRO

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

31.3 Distance methods and summary functions 193<br />

array of G functions for lans<strong>in</strong>g.<br />

hickory<br />

maple<br />

Gcrosshickory, hickory(r)<br />

0.0 0.2 0.4 0.6 0.8<br />

Gcrosshickory, maple(r)<br />

0.0 0.2 0.4 0.6 0.8<br />

0.000 0.005 0.010 0.015 0.020 0.025 0.030<br />

r (one unit = 924 feet)<br />

0.000 0.005 0.010 0.015 0.020 0.025 0.030<br />

r (one unit = 924 feet)<br />

maple<br />

Gcrossmaple, hickory(r)<br />

0.0 0.2 0.4 0.6 0.8<br />

Gcrossmaple, maple(r)<br />

0.0 0.2 0.4 0.6 0.8<br />

hickory<br />

0.000 0.005 0.010 0.015 0.020 0.025 0.030<br />

r (one unit = 924 feet)<br />

0.000 0.005 0.010 0.015 0.020 0.025 0.030<br />

r (one unit = 924 feet)<br />

31.3.3 One type to any type<br />

Also def<strong>in</strong>ed are the “i-to-any” summaries<br />

G i• (r), the distribution function of the distance from a <strong>po<strong>in</strong>t</strong> of type i to the nearest other<br />

<strong>po<strong>in</strong>t</strong> of any type;<br />

K i• (r) is 1/λ times the expected number of <strong>po<strong>in</strong>t</strong>s of any type with<strong>in</strong> a distance r of a<br />

typical <strong>po<strong>in</strong>t</strong> of type i. Here λ = ∑ j λ j is the <strong>in</strong>tensity of the entire process X.<br />

L i• (r) is the correspond<strong>in</strong>g L-function<br />

√<br />

Ki• (r)<br />

L i• (r) = .<br />

π<br />

g i• (r) is the correspond<strong>in</strong>g analogue of the pair correlation function<br />

g i• (r) = K′ i• (r)<br />

2πr .<br />

J i• is def<strong>in</strong>ed by<br />

J i• (r) = 1 − G i•<br />

1 − F(r)<br />

These are comput<strong>in</strong>g by Gdot,Kdot, Ldot,pcfdot andJdot respectively, or us<strong>in</strong>g alltypes.<br />

> plot(Gdot(amacr<strong>in</strong>e, "on"))<br />

Copyright<strong>CSIRO</strong> 2010

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!