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.

178 Marked <strong>po<strong>in</strong>t</strong> <strong>patterns</strong><br />

29 Marked <strong>po<strong>in</strong>t</strong> <strong>patterns</strong><br />

29.1 Marked <strong>po<strong>in</strong>t</strong> <strong>patterns</strong><br />

Each <strong>po<strong>in</strong>t</strong> <strong>in</strong> a <strong>spatial</strong> <strong>po<strong>in</strong>t</strong> pattern may carry additional <strong>in</strong>formation called a ‘mark’. For example,<br />

<strong>po<strong>in</strong>t</strong>s which are classified <strong>in</strong>to two or more different types (on/off, case/control, species,<br />

colour, etc) may be regarded as marked <strong>po<strong>in</strong>t</strong>s, with a mark which identifies which type they<br />

are. Data record<strong>in</strong>g the locations and heights of trees <strong>in</strong> a forest can be regarded as a marked<br />

<strong>po<strong>in</strong>t</strong> pattern where the mark attached to a tree’s location is the tree height.<br />

Many of the functions <strong>in</strong>spatstat handle marked <strong>po<strong>in</strong>t</strong> <strong>patterns</strong> <strong>in</strong> which the mark attached<br />

to each <strong>po<strong>in</strong>t</strong> is either<br />

acont<strong>in</strong>uous variate or “real number”. An example is the Longleaf P<strong>in</strong>es dataset<br />

(longleaf) <strong>in</strong> which each tree is marked with its diameter at breast height. The marks<br />

component must be a numeric vector such that marks[i] is the mark value associated<br />

with the ith <strong>po<strong>in</strong>t</strong>. We say the <strong>po<strong>in</strong>t</strong> pattern has cont<strong>in</strong>uous marks.<br />

a categorical variate. An example is the Amacr<strong>in</strong>e Cells dataset (amacr<strong>in</strong>e) <strong>in</strong> which<br />

each cell is identified as either “on” or “off”. Such <strong>po<strong>in</strong>t</strong> <strong>patterns</strong> may be regarded as<br />

consist<strong>in</strong>g of <strong>po<strong>in</strong>t</strong>s of different “types”. The marks component must be a factor such<br />

that marks[i] is the label or type of the ith <strong>po<strong>in</strong>t</strong>. We call this a multitype <strong>po<strong>in</strong>t</strong> pattern<br />

and the levels of the factor are the possible types.<br />

longleaf<br />

amacr<strong>in</strong>e<br />

Note that, <strong>in</strong> some other packages, a <strong>po<strong>in</strong>t</strong> pattern dataset consist<strong>in</strong>g of <strong>po<strong>in</strong>t</strong>s of two different<br />

types (A and B say) is represented by two datasets, one represent<strong>in</strong>g the <strong>po<strong>in</strong>t</strong>s of type A and<br />

another conta<strong>in</strong><strong>in</strong>g the <strong>po<strong>in</strong>t</strong>s of type B. In spatstat we take a different approach, <strong>in</strong> which<br />

all the <strong>po<strong>in</strong>t</strong>s are collected together <strong>in</strong> one <strong>po<strong>in</strong>t</strong> pattern, and the <strong>po<strong>in</strong>t</strong>s are then labelled by<br />

the type to which they belong. An advantage of this approach is that it is easy to deal with<br />

multitype <strong>po<strong>in</strong>t</strong> <strong>patterns</strong> with more than 2 types. For example the classic Lans<strong>in</strong>g Woods dataset<br />

represents the positions of trees of 6 different species. This is available <strong>in</strong> spatstat as a s<strong>in</strong>gle<br />

dataset, a marked <strong>po<strong>in</strong>t</strong> pattern, with the marks hav<strong>in</strong>g 6 levels.<br />

29.2 Formulation<br />

A mark variable may be <strong>in</strong>terpreted as an additional coord<strong>in</strong>ate for the <strong>po<strong>in</strong>t</strong>: for example<br />

a <strong>po<strong>in</strong>t</strong> process of earthquake epicentre locations (longitude, latitude), with marks giv<strong>in</strong>g the<br />

occurrence time of each earthquake, can alternatively be viewed as a <strong>po<strong>in</strong>t</strong> process <strong>in</strong> space-time<br />

with coord<strong>in</strong>ates (longitude, latitude, time).<br />

Copyright<strong>CSIRO</strong> 2010

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

Saved successfully!

Ooh no, something went wrong!