spatstat.model (version 3.3-5)

residualMeasure: Residual Measure for an Observed Point Pattern and a Fitted Intensity

Description

Given a point pattern and an estimate of its intensity function obtained in any fashion, compute the residual measure.

Usage

residualMeasure(Q, lambda,
               type = c("raw", "inverse", "Pearson", "pearson"),
               ...)

Value

A measure (object of class "msr").

Arguments

Q

A point pattern (object of class "ppp") or quadrature scheme (object of class "quad").

lambda

Predicted intensity. An image (object of class "im") or a list of images.

type

Character string (partially matched) specifying the type of residuals.

...

Arguments passed to quadscheme if Q is a point pattern.

Author

Adrian Baddeley .

Details

This command constructs the residual measure for the model in which Q is the observed point pattern or quadrature scheme, and lambda is the estimated intensity function obtained in some fashion.

References

Baddeley, A., Turner, R., Moller, J. and Hazelton, M. (2005) Residual analysis for spatial point processes. Journal of the Royal Statistical Society, Series B 67, 617--666.

Baddeley, A., Moller, J. and Pakes, A.G. (2008) Properties of residuals for spatial point processes. Annals of the Institute of Statistical Mathematics 60, 627--649.

See Also

residuals.ppm

Examples

  ## nonparametric regression estimate of intensity
  ## as a function of terrain elevation
  f <- rhohat(bei, bei.extra$elev)
  ## predicted intensity as a function of location
  lam <- predict(f)
  ## residuals
  res <- residualMeasure(bei, lam)
  res
  plot(res)