Last chance! 50% off unlimited learning
Sale ends in
For a multitype point pattern,
estimate the inhomogeneous version of the dot
Ldot.inhom(X, i, ..., correction)
An object of class "fv"
(see fv.object
).
Essentially a data frame containing numeric columns
the values of the argument
the theoretical value of
together with a column or columns named
"border"
, "bord.modif"
,
"iso"
and/or "trans"
,
according to the selected edge corrections. These columns contain
estimates of the function
obtained by the edge corrections named.
The observed point pattern,
from which an estimate of the inhomogeneous cross type
The type (mark value)
of the points in X
from which distances are measured.
A character string (or something that will be converted to a
character string).
Defaults to the first level of marks(X)
.
Other arguments passed to Kdot.inhom
.
The argument i
is interpreted as
a level of the factor X$marks
. It is converted to a character
string if it is not already a character string.
The value i=1
does not
refer to the first level of the factor.
Adrian Baddeley Adrian.Baddeley@curtin.edu.au, Rolf Turner rolfturner@posteo.net and Ege Rubak rubak@math.aau.dk
This a generalisation of the function Ldot
to include an adjustment for spatially inhomogeneous intensity,
in a manner similar to the function Linhom
.
All the arguments are passed to Kdot.inhom
, which
estimates the inhomogeneous multitype K function
Moller, J. and Waagepetersen, R. Statistical Inference and Simulation for Spatial Point Processes Chapman and Hall/CRC Boca Raton, 2003.
Ldot
,
Linhom
,
Kdot.inhom
,
Lcross.inhom
.
# Lansing Woods data
lan <- lansing
lan <- lan[seq(1,npoints(lan), by=10)]
ma <- split(lan)$maple
lg <- unmark(lan)
# Estimate intensities by nonparametric smoothing
lambdaM <- density(ma, sigma=0.15, at="points")
lambdadot <- density(lg, sigma=0.15, at="points")
L <- Ldot.inhom(lan, "maple", lambdaI=lambdaM,
lambdadot=lambdadot)
# synthetic example: type A points have intensity 50,
# type B points have intensity 50 + 100 * x
lamB <- as.im(function(x,y){50 + 100 * x}, owin())
lamdot <- as.im(function(x,y) { 100 + 100 * x}, owin())
X <- superimpose(A=runifpoispp(50), B=rpoispp(lamB))
L <- Ldot.inhom(X, "B", lambdaI=lamB, lambdadot=lamdot)
Run the code above in your browser using DataLab