plot.localdiag {stopp}

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 class '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(s)

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

See Also

infl, print.localdiag, summary.localdiag

Examples

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)