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

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

126 Distance methods for <strong>po<strong>in</strong>t</strong> <strong>patterns</strong><br />

> par(pty = "s")<br />

> plot(Kest(cells))<br />

Kest(cells)<br />

K(r)<br />

0.00 0.05 0.10 0.15 0.20<br />

iso<br />

trans<br />

border<br />

theo<br />

0.00 0.05 0.10 0.15 0.20 0.25<br />

r<br />

In this case, the <strong>in</strong>terpretation of all three summary statistics F, G and K is the same:<br />

emphatic evidence of a regular pattern. It is not always the case that these three summaries<br />

give equivalent messages.<br />

A commonly-used transformation of K is the L-function<br />

L(r) =<br />

√<br />

K(r)<br />

which transforms the Poisson K function to the straight l<strong>in</strong>e L pois (r) = r, mak<strong>in</strong>g visual assessment<br />

of the graph much easier. The square root transformation also approximately stabilises<br />

the variance of the estimator, mak<strong>in</strong>g it easier to assess deviations.<br />

To compute the estimated L function, use Lest.<br />

> L plot(L, ma<strong>in</strong> = "L function")<br />

π<br />

L function<br />

L(r)<br />

0.00 0.05 0.10 0.15 0.20 0.25<br />

iso<br />

trans<br />

border<br />

theo<br />

0.00 0.05 0.10 0.15 0.20 0.25<br />

r<br />

Another related summary function is the pair correlation function<br />

Copyright<strong>CSIRO</strong> 2010<br />

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

2πr

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

Saved successfully!

Ooh no, something went wrong!