Analysing spatial point patterns in R - CSIRO
Analysing spatial point patterns in R - CSIRO
Analysing spatial point patterns in R - CSIRO
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