markstat (version 0.1.5)

lfmcorrinhom: Local mark correlation functions for inhomogeneous point patterns with function-valued marks.

Description

Local mark correlation functions for inhomogeneous point patterns with function-valued marks.

Usage

lfmcorrinhom(X,
ftype = c("variogram", "stoyan", "rcorr", "shimatani", "beisbart",
 "isham", "stoyancov", "schlather"),
r = NULL,
lambda = NULL,
method_lambda = c("kernel", "Voronoi"),
bw = NULL,
f = NULL,
method = c("density", "loess"),
correction = c("Ripley", "translate", "none"),
normalise = TRUE,
tol = 0.01,
...)

Value

a data.frame which gives the estimated overall local mark correlation function and the distance vector \(r\) at which the local mark correlation function is estimated. The outputs of the local mark correlation functions for each time point are stored as an attribute, which can be extracted as attr(., "ests.time"). The outputs of the local mark correlation functions for each data point are stored as an attribute, which can be extracted as attr(., "ests.points").

Arguments

X

An object of class ppp or lpp.

ftype

Type of the test function \(t_f\). Currently any selection of "variogram", "stoyan", "rcorr", "shimatani", "beisbart", "isham", "stoyancov", "schlather".

r

Optional. The values of the argument r at which the mark correlation function should be evaluated.

lambda

Estimated intensity at data points. If not given, it will be estimated internally. See details.

method_lambda

The method to be used for estimating intensity at data points, if lambda = NULL.

bw

Bandwidth method to be used for estimating intensity at data points if lambda = NULL and method_lambda = "kernel".

f

Optional. Test function \(t_f\) used in the definition of the mark correlation function. If ftype is given, \(t_f\) should be NULL.

method

Type of smoothing, either density or loess.

correction

Type of edge correction to be applied, either of "Ripley", "translate", "none". This is used only if X is of class ppp.

normalise

If normalise=FALSE, only the numerator of the expression for the mark correlation function will be computed.

tol

Tolerance used in the calculation of the conditional mean of marks. This is used only if ftype is schlather.

...

Arguments passed to unnormdensity or loess.

Author

Mehdi Moradi and Matthias Eckardt

Details

This function computes local mark correlation functions for an inhomogeneous point pattern with a function-valued mark. See the details of test functions used in fmcorrinhom. Technical details are given in Eckardt and Moradi (2025).

References

Eckardt, M., Mateu, J., & Moradi, M. (2024). Function‐Valued Marked Spatial Point Processes on Linear Networks: Application to Urban Cycling Profiles. Stat, 13(4), e70013.

Eckardt, M., & Moradi, M. (2025). Local indicators of mark association for marked spatial point processes.

Moradi, M., & Eckardt, M. (2025). Inhomogeneous mark correlation functions for general marked point processes. arXiv e-prints, arXiv-2505.

See Also

fmcorr, mcorrinhom.ppp, mcorrinhom.lpp.

Examples

 library(spatstat.random)
 library(spatstat.geom)
 library(spatstat.explore)
 X <- rpoispp(100)
 marks(X) <- data.frame(
 t1 = runif(npoints(X),1,10),
 t2 = runif(npoints(X),1,10),
 t3 = runif(npoints(X),1,10),
 t4 = runif(npoints(X),1,10),
 t5 = runif(npoints(X),1,10))
 lfmcorrinhom(X, ftype = "stoyan", method = "density", method_lambda = "kernel", bw = bw.scott)