Computes the smoothed partial residuals, a diagnostic for transformation of a covariate in a Poisson point process model on a linear network.
# S3 method for lppm
parres(model, covariate, ...,
smooth.effect=FALSE, subregion=NULL,
bw = "nrd0", adjust=1, from = NULL, to = NULL, n = 512,
bw.input = c("points", "quad"), bw.restrict=FALSE, covname)A function value table (object of class "fv")
containing the values of the smoothed partial residual,
the estimated variance, and the fitted effect of the covariate.
Also belongs to the class "parres"
which has methods for print and plot.
This command computes the smoothed partial residual diagnostic (Baddeley, Chang, Song and Turner, 2012) for the transformation of a covariate in a Poisson point process model.
The function parres is generic, with methods for different
classes of point process models. This page documents the method
parres.lppm. The argument model must be a fitted
point process model on a linear network.
The diagnostic works in two different ways:
The argument covariate may be a character string
which is the name of one of the canonical covariates in the
model.
The canonical covariates are the
functions \(Z_j\) that appear
in the expression for the Poisson point process intensity
$$
\lambda(u) = \exp(\beta_1 Z_1(u) + \ldots + \beta_p Z_p(u))
$$
at spatial location \(u\).
Type names(coef(model)) to see the names of the
canonical covariates in model.
If the selected covariate is \(Z_j\), then
the diagnostic plot concerns the model term
\(\beta_j Z_j(u)\). The plot shows a smooth
estimate of a function \(h(z)\) that should replace this linear
term, that is, \(\beta_j Z_j(u)\) should be
replaced by \(h(Z_j(u))\). The linear function is
also plotted as a dotted line.
If the argument covariate is a pixel image
(object of class "im") or a function(x,y),
it is assumed to provide the values of a covariate that is
not present in the model.
Alternatively covariate can be the name of a
covariate that was supplied when the model was fitted
(i.e. in the call to ppm)
but which does not feature in the model formula.
In either case we speak of a new covariate \(Z(u)\).
If the fitted model intensity is \(\lambda(u)\)
then we consider modifying this to
\(\lambda(u) \exp(h(Z(u)))\)
where \(h(z)\) is some function. The diagnostic plot shows
an estimate of \(h(z)\).
Warning: in this case the diagnostic is not theoretically
justified. This option is provided for research purposes.
Alternatively covariate can be one of the character strings
"x" or "y" signifying the Cartesian coordinates.
The behaviour here depends on whether the coordinate was one of the
canonical covariates in the model.
If there is more than one canonical covariate in the model
that depends on the specified covariate, then
the covariate effect is computed using all these canonical covariates.
For example in a log-quadratic model which includes the terms x and
I(x^2), the quadratic effect involving both these terms
will be computed.
There are two choices for the algorithm.
If smooth.effect=TRUE, the fitted covariate effect (according
to model) is added to the point process residuals, then
smoothing is applied to these values. If smooth.effect=FALSE,
the point process residuals are smoothed first, and then the fitted
covariate effect is added to the result.
The smoothing bandwidth is controlled by the arguments
bw, adjust, bw.input and bw.restrict.
If bw is a numeric value, then
the bandwidth is taken to be adjust * bw.
If bw is a string representing a bandwidth selection rule
(recognised by density.default)
then the bandwidth is selected by this rule.
The data used for automatic bandwidth selection are
specified by bw.input and bw.restrict.
If bw.input="points" (the default) then bandwidth selection is
based on the covariate values at the points of the original point
pattern dataset to which the model was fitted.
If bw.input="quad" then bandwidth selection is
based on the covariate values at every quadrature point used to
fit the model.
If bw.restrict=TRUE then the bandwidth selection is performed
using only data from inside the subregion.
Baddeley, A., Chang, Y.-M., Song, Y. and Turner, R. (2013) Residual diagnostics for covariate effects in spatial point process models. Journal of Computational and Graphical Statistics, 22, 886--905.
addvar,
rhohat,
rho2hat
X <- rpoispp(function(x,y){exp(3+x+2*x^2)})
model <- ppm(X ~x+y)
tra <- parres(model, "x")
plot(tra)
tra
plot(parres(model, "x", subregion=square(0.5)))
model2 <- ppm(X ~x+I(x^2)+y)
plot(parres(model2, "x"))
Z <- setcov(owin())
plot(parres(model2, Z))
#' when the model involves only one covariate
modelb <- ppm(bei ~ elev + I(elev^2), data=bei.extra)
plot(parres(modelb))
Run the code above in your browser using DataLab