For a marked point pattern, 
estimate the multitype pair correlation function
using kernel methods.

pcfmulti

Functionality for exploratory data analysis and nonparametric analysis of
spatial data, mainly spatial point patterns,
in the 'spatstat' family of packages.
(Excludes analysis of spatial data on a linear network,
which is covered by the separate package 'spatstat.linnet'.)
Methods include quadrat counts, K-functions and their simulation envelopes, nearest neighbour distance and empty space statistics, Fry plots, pair correlation function, kernel smoothed intensity, relative risk estimation with cross-validated bandwidth selection, mark correlation functions, segregation indices, mark dependence diagnostics, and kernel estimates of covariate effects. Formal hypothesis tests of random pattern (chi-squared, Kolmogorov-Smirnov, Monte Carlo, Diggle-Cressie-Loosmore-Ford, Dao-Genton, two-stage Monte Carlo) and tests for covariate effects (Cox-Berman-Waller-Lawson, Kolmogorov-Smirnov, ANOVA) are also supported.

Adrian Baddeley

spatstat.explore

Exploratory Data Analysis for the 'spatstat' Family

Rolf Turner

Ege Rubak

Kasper Klitgaard Berthelsen

Warick Brown

Achmad Choiruddin

Ya-Mei Chang

Jean-Francois Coeurjolly

Lucia Cobo Sanchez

Ottmar Cronie

Tilman Davies

Julian Gilbey

Jonatan Gonzalez

Yongtao Guan

Ute Hahn

Martin Hazelton

Kassel Hingee

Abdollah Jalilian

Frederic Lavancier

Marie-Colette van Lieshout

Greg McSwiggan

Robin K Milne

Tuomas Rajala

Suman Rakshit

Dominic Schuhmacher

Rasmus Plenge Waagepetersen

Hangsheng Wang

Tingting Zhan

pcfmulti function

Marked pair correlation function — pcfmulti

<dl>

 <dt>X</dt>
<dd>The observed point pattern, 
 from which an estimate of the multitype pair correlation function
 \(g_{IJ}(r)\) will be computed.
 It must be a multitype point pattern (a marked point pattern
 whose marks are a factor).</dd>

 <dt>I</dt>
<dd>Subset index specifying the points of <code>X</code>
 from which distances are measured.</dd>

 <dt>J</dt>
<dd>Subset index specifying the points in <code>X</code>
 to which distances are measured.</dd>

 <dt>...</dt>
<dd>Ignored.</dd>

 <dt>r</dt>
<dd>Vector of values for the argument \(r\) at which \(g(r)\) 
 should be evaluated. There is a sensible default.</dd>

 <dt>rmax</dt>
<dd>Optional. Maximum desired value of the argument \(r\).
 A single numeric value. There is a sensible default.</dd>

 <dt>adaptive</dt>
<dd>Logical value specifying whether to use adaptive kernel smoothing
 (<code>adaptive=TRUE</code>) or fixed-bandwidth kernel smoothing
 (<code>adaptive=FALSE</code>, the default).</dd>

 <dt>kernel</dt>
<dd>Choice of smoothing kernel,
 passed to <code><a href='https://rdrr.io/r/stats/density.html'>density.default</a></code>.</dd>

 <dt>bw</dt>
<dd>Bandwidth for smoothing kernel.
 Either a single numeric value giving the standard deviation of the
 kernel, or a character string specifying a bandwidth selection rule,
 or a function that computes the selected bandwidth.
 See Details.</dd>

 <dt>h</dt>
<dd>Kernel halfwidth \(h\) (incompatible with argument <code>bw</code>).
 A numerical value.
 The parameter <code>h</code> is defined as the
 half-width of the support of the kernel, except for the Gaussian
 kernel where <code>h</code> is the standard deviation.</dd>

 <dt>bw.args</dt>
<dd>Optional. List of additional arguments to be passed to <code>bw</code>
 when <code>bw</code> is a function. Alternatively, <code>bw</code> may be a
 function that should be applied to <code>X</code> to produce a list of
 additional arguments.</dd>

 <dt>stoyan</dt>
<dd>Coefficient for default bandwidth rule.</dd>

 <dt>adjust</dt>
<dd>Numerical adjustment factor for the bandwidth.
 The bandwidth actually used is <code>adjust * bw</code>.
 This makes it easy to specify choices like &#8216;half the
 selected bandwidth&#8217;.</dd>

 <dt>correction</dt>
<dd>String (partially matched) specifying the choice or choices
 of spatial edge correction. Options include <code>"translate"</code> for
 the translation correction, <code>"isotropic"</code> or <code>"Ripley"</code>
 for Ripley's isotropic correction, and <code>"none"</code> for no edge correction.</dd>

 <dt>divisor</dt>
<dd>String specifying the choice of estimator.
 See <code>pcf.ppp</code>.</dd>

 <dt>zerocor</dt>
<dd>String (partially matched) specifying a correction for the boundary effect
 bias at \(r=0\). Possible values are
 <code>"none"</code>, <code>"weighted"</code>, <code>"convolution"</code>,
 <code>"reflection"</code>, <code>"bdrykern"</code> and <code>"JonesFoster"</code>.
 See <code>pcf.ppp</code>.</dd>

 <dt>nsmall</dt>
<dd>Optional. Integer. The maximum number of data points
 for which the default value of <code>zerocor</code> will be
 <code>"JonesFoster"</code>.</dd>

 <dt>gref</dt>
<dd>Optional. A pair correlation function that will be used as the
 reference for the transformation to uniformity, when
 <code>divisor="t"</code>. Either a <code>function</code> in the R language
 giving the pair correlation function, or a fitted model
 (object of class <code>"kppm"</code>, <code>"dppm"</code>, <code>"ppm"</code>
 or <code>"slrm"</code>) or a theoretical point process model
 (object of class <code>"zclustermodel"</code> or <code>"detpointprocfamily"</code>)
 for which the pair correlation function
 can be computed.</dd>

 <dt>tau</dt>
<dd>Optional shrinkage coefficient. A single numeric value.</dd>

 <dt>Iname,Jname</dt>
<dd>Optional. Character strings describing the members of
 the subsets <code>I</code> and <code>J</code>.</dd>

 <dt>IJexclusive</dt>
<dd>Logical value indicating whether the subsets <code>I</code> and <code>J</code>
 are guaranteed to be mutually exclusive.</dd>

 <dt>ratio</dt>
<dd>Logical. 
 If <code>TRUE</code>, the numerator and denominator of
 each edge-corrected estimate will also be saved,
 for use in analysing replicated point patterns.</dd>

 <dt>close</dt>
<dd>Advanced use only. Precomputed data obtained from <code>crosspairs</code>.</dd>

</dl>

Arguments

Adrian Baddeley <a href='mailto:Adrian.Baddeley@curtin.edu.au'>Adrian.Baddeley@curtin.edu.au</a>, Rolf Turner <a href='mailto:rolfturner@posteo.net'>rolfturner@posteo.net</a>, Ege Rubak <a href='mailto:rubak@math.aau.dk'>rubak@math.aau.dk</a>, Tilman Davies <a href='mailto:Tilman.Davies@otago.ac.nz'>Tilman.Davies@otago.ac.nz</a> and Martin Hazelton <a href='mailto:Martin.Hazelton@otago.ac.nz'>Martin.Hazelton@otago.ac.nz</a>.

Author

Marked pair correlation function

pcfmulti: Marked pair correlation function

Description

Usage

Arguments

Value

Details

See Also

Examples