plot.localdiag: Plot of the local diagnostics' result on a spatio-temporal point process model

Description

This function plots the result of the local diagnostics performed with localdiag on either a stp or stlp object. It highlights the points of the analysed spatio-temporal point pattern X which are identified as outlying by the previously performed local diagnostics; the remaining points of X are also represented.

It also shows the underlying linear network, if the local diagnostics has been applied to point patterns occurring on the same linear network, that is, if localdiag has been applied to a stlp object.

Usage

# S3 method for localdiag
plot(x, marg = TRUE, col = "grey", col2 = "red", cols = "lightgrey", ...)

Arguments

x: A localdiag object
marg: Default to TRUE. If marg = F, only the spatio-temporal point pattern is plotted
col: Color of the outlying points
col2: Color of the network (if applicable)
cols: Color of the non-outlying points
...: additional unused argument

Author

Nicoletta D'Angelo and Giada Adelfio

References

Adelfio, G., Siino, M., Mateu, J., and Rodríguez-Cortés, F. J. (2020). Some properties of local weighted second-order statistics for spatio-temporal point processes. Stochastic Environmental Research and Risk Assessment, 34(1), 149-168.

D’Angelo, N., Adelfio, G. and Mateu, J. (2022) Local inhomogeneous second-order characteristics for spatio-temporal point processes on linear networks. Stat Papers. https://doi.org/10.1007/s00362-022-01338-4

Examples

Run this code


set.seed(2)
inh <- rstpp(lambda = function(x, y, t, a) {exp(a[1] + a[2]*x)}, 
             par = c(.3, 6))

mod1 <- stppm(inh, formula = ~ 1)

resmod1 <- localdiag(inh, mod1$l, p = .9)

plot(resmod1)
plot(resmod1, marg = FALSE)