linearpcfcross.inhom {spatstat.linnet} | R Documentation |
Inhomogeneous Multitype Pair Correlation Function (Cross-type) for Linear Point Pattern
Description
For a multitype point pattern on a linear network,
estimate the inhomogeneous multitype pair correlation function
from points of type i
to points of type j
.
Usage
linearpcfcross.inhom(X, i, j, lambdaI, lambdaJ, r=NULL, ...,
correction="Ang", normalise=TRUE,
sigma=NULL, adjust.sigma=1,
bw="nrd0", adjust.bw=1)
Arguments
X |
The observed point pattern,
from which an estimate of the |
i |
Number or character string identifying the type (mark value)
of the points in |
j |
Number or character string identifying the type (mark value)
of the points in |
lambdaI |
Intensity values for the points of type |
lambdaJ |
Intensity values for the points of type |
r |
numeric vector. The values of the argument |
correction |
Geometry correction.
Either |
... |
Arguments passed to |
normalise |
Logical. If |
sigma |
Smoothing bandwidth passed to |
adjust.sigma |
Numeric value. |
bw |
Smoothing bandwidth (passed to |
adjust.bw |
Numeric value. |
Details
This is a counterpart of the function pcfcross.inhom
for a point pattern on a linear network (object of class "lpp"
).
The argument i
will be interpreted as
levels of the factor marks(X)
.
If i
is missing, it defaults to the first
level of the marks factor.
The argument r
is the vector of values for the
distance r
at which g_{ij}(r)
should be evaluated.
The values of r
must be increasing nonnegative numbers
and the maximum r
value must not exceed the radius of the
largest disc contained in the window.
If lambdaI
or lambdaJ
is missing or NULL
, it will
be estimated by kernel smoothing using density.lpp
.
If lambdaI
or lambdaJ
is a fitted point process model,
the default behaviour is to update the model by re-fitting it to
the data, before computing the fitted intensity.
This can be disabled by setting update=FALSE
.
Value
An object of class "fv"
(see fv.object
).
Warnings
The argument i
is interpreted as a
level of the factor marks(X)
. Beware of the usual
trap with factors: numerical values are not
interpreted in the same way as character values.
Author(s)
Adrian Baddeley Adrian.Baddeley@curtin.edu.au
References
Baddeley, A, Jammalamadaka, A. and Nair, G. (2014) Multitype point process analysis of spines on the dendrite network of a neuron. Applied Statistics (Journal of the Royal Statistical Society, Series C), 63, 673–694.
See Also
linearpcfdot
,
linearpcf
,
pcfcross.inhom
.
Examples
lam <- table(marks(chicago))/(summary(chicago)$totlength)
lamI <- function(x,y,const=lam[["assault"]]){ rep(const, length(x)) }
lamJ <- function(x,y,const=lam[["robbery"]]){ rep(const, length(x)) }
g <- linearpcfcross.inhom(chicago, "assault", "robbery", lamI, lamJ)
# using fitted models for intensity
# fit <- lppm(chicago ~marks + x)
# linearpcfcross.inhom(chicago, "assault", "robbery", fit, fit)