linearKcross.inhom(X, i, j, lambdaI, lambdaJ,
r=NULL, ..., correction="Ang", normalise=TRUE)
"lpp"
which
must be a multitype point pattern (a marked point pattern
whose marks are a fX
from which distances are measured.
Defaults to the first level of marks(X)
.X
to which distances are measured.
Defaults to the second level of marks(X)
.i
. Either a numeric vector,
a function
, a pixel image
(object of class "im"
or "linim"
) or
a fitted point process model (object of class "ppm"
j
. Either a numeric vector,
a function
, a pixel image
(object of class "im"
or "linim"
) or
a fitted point process model (object of class "ppm"
"none"
or "Ang"
. See Details.lambdaI
and lambdaJ
if
they are functions.TRUE
(the default), the denominator of the estimator is
data-dependent (equal to the sum of the reciprocal intensities at
the points of type i
), which reduces the sampling variability.
If FALSE
"fv"
(see fv.object
).i
and j
are interpreted as
levels of the factor marks(X)
. Beware of the usual
trap with factors: numerical values are not
interpreted in the same way as character values.Kcross.inhom
for a point pattern on a linear network (object of class "lpp"
). The arguments i
and j
will be interpreted as
levels of the factor marks(X)
.
If i
and j
are missing, they default to the first
and second level of the marks factor, respectively.
The argument r
is the vector of values for the
distance $r$ at which $K_{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 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
.
linearKdot
,
linearK
.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)) }
K <- linearKcross.inhom(chicago, "assault", "robbery", lamI, lamJ)
fit <- lppm(chicago, ~marks + x)
linearKcross.inhom(chicago, "assault", "robbery", fit, fit)
Run the code above in your browser using DataLab