plot.locmincon {spatstat.local}

Plot a Locally Fitted Cluster or Cox Point Process Model

Description

Plot an object of class "locmincon" representing a locally-fitted cluster or Cox point process model.

Usage

## S3 method for class 'locmincon'
plot(x, ...,
               how = c("exact", "smoothed"),
               which = NULL, sigma = NULL, do.points = TRUE)

Arguments

x

The model to be plotted. A locally-fitted cluster or Cox point process model (object of class "locmincon" or "loccit").

...

Arguments passed to plot.ppp or plot.im to control the plot.

how

Character string determining whether to display the fitted parameter values at the data points (how="exact") or the smoothed fitted parameters as pixel images (how="smoothed").

which

Optional. Which component(s) of the vector-valued quantity to display. An index or index vector. Default is to plot all components.

sigma

Numeric. Smoothing bandwidth to be used if how="smoothed".

do.points

Logical. Whether to display the original point data as well.

Details

This is a method for the generic command plot for the class "locmincon".

The argument which, if present, specifies which fitted parameters are displayed. It may be any kind of index for a numeric vector.

The quantities are computed at irregularly-placed points. If how="exact" the exact computed values will be displayed as circles centred at the locations where they were computed. If how="smoothed" these values will be kernel-smoothed using Smooth.ppp and displayed as a pixel image.

Value

NULL.

Author(s)

Adrian Baddeley [e-mail].

References

Baddeley, A. (2017) Local composite likelihood for spatial point patterns. Spatial Statistics 22, 261–295. DOI: 10.1016/j.spasta.2017.03.001

Baddeley, A., Rubak, E. and Turner, R. (2015) Spatial Point Patterns: Methodology and Applications with R. Chapman and Hall/CRC Press.

See Also

locmincon, methods.locmincon, plot, plot.default

Examples

   X <- redwood[owin(c(0,1), c(-1,-1/2))]
   fitm <- locmincon(X, ~1, "Thomas", sigma=0.07)
   plot(fitm, how="smoothed")  
   plot(fitm, how="exact")