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.

158 Gibbs models<br />

Hard core process<br />

26.3.2 Strauss process<br />

Generalis<strong>in</strong>g the hard core process, suppose we take b(u) ≡ β and<br />

c(u,v) =<br />

{ 1 if ||u − v|| > r<br />

γ if ||u − v|| ≤ r<br />

(36)<br />

where γ is a parameter. Then the density becomes<br />

f(x) = αβ n(x) γ s(x) (37)<br />

where s(x) is the number of pairs of dist<strong>in</strong>ct <strong>po<strong>in</strong>t</strong>s <strong>in</strong> x that lie closer than r units apart.<br />

The parameter γ controls the ‘strength’ of <strong>in</strong>teraction between <strong>po<strong>in</strong>t</strong>s. If γ = 1 the model<br />

reduces to a Poisson process with <strong>in</strong>tensity β. If γ = 0 the model is a hard core process. For<br />

values 0 < γ < 1, the process exhibits <strong>in</strong>hibition (negative association) between <strong>po<strong>in</strong>t</strong>s.<br />

Strauss(γ = 0.2) Strauss(γ = 0.7)<br />

For γ > 1, the density (37) is not <strong>in</strong>tegrable. Hence the Strauss process is def<strong>in</strong>ed only for<br />

0 ≤ γ ≤ 1 and is a model for <strong>in</strong>hibition between <strong>po<strong>in</strong>t</strong>s. This is typical of most Gibbs models.<br />

Copyright<strong>CSIRO</strong> 2010

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

Saved successfully!

Ooh no, something went wrong!