Description Usage Arguments Details Value Author(s) References See Also Examples

View source: R/residuals.mppm.R

Given a point process model fitted to multiple point patterns, compute residuals for each pattern.

1 2 3 | ```
## S3 method for class 'mppm'
residuals(object, type = "raw", ...,
fittedvalues = fitted.mppm(object))
``` |

`object` |
Fitted point process model (object of class |

`...` |
Ignored. |

`type` |
Type of residuals: either |

`fittedvalues` |
Advanced use only. Fitted values of the model to be used in the calculation. |

Baddeley et al (2005) defined residuals for the fit of
a point process model to spatial point pattern data.
For an explanation of these residuals, see the help file for
`residuals.ppm`

.

This function computes the residuals
for a point process model fitted to *multiple* point patterns.
The `object`

should be an object of class `"mppm"`

obtained
from `mppm`

.

The return value is a list.
The number of entries in the list equals the
number of point patterns in the original data. Each entry in the list
has the same format as the output of
`residuals.ppm`

.
That is, each entry in the list is a signed measure (object of
class `"msr"`

) giving the residual measure for the corresponding
point pattern.

A list of signed measures (objects of class `"msr"`

)
giving the residual measure for each of the original point patterns.
See Details.

, Ida-Maria Sintorn and Leanne Bischoff. Implemented by \adrian

\rolfand \ege

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., Rubak, E. and Turner, R. (2015)
*Spatial Point Patterns: Methodology and Applications with R*.
London: Chapman and Hall/CRC Press.

1 2 3 4 5 6 7 | ```
fit <- mppm(Bugs ~ x, hyperframe(Bugs=waterstriders))
r <- residuals(fit)
# compute total residual for each point pattern
rtot <- sapply(r, integral.msr)
# standardise the total residuals
areas <- sapply(windows.mppm(fit), area.owin)
rtot/sqrt(areas)
``` |

spatstat documentation built on Jan. 30, 2018, 1:05 a.m.

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.