Extract Quadrature Scheme Used to Fit a Point Process Model
Given a fitted point process model, this function extracts the quadrature scheme used to fit the model.
- fitted point process model (an object of class
- Logical value determining whether to delete quadrature points that were not used to fit the model.
An object of class
"ppm" represents a point process model
that has been fitted to data. It is typically produced by
the model-fitting algorithm
The maximum pseudolikelihood algorithm in
approximates the pseudolikelihood
integral by a sum over a finite set of quadrature points,
which is constructed by augmenting
the original data point pattern by a set of ``dummy'' points.
The fitted model object returned by
contains complete information about this quadrature scheme.
ppm.object for further
quad.ppm extracts the quadrature scheme.
A typical use of this function would be to inspect the quadrature scheme
(points and weights) to gauge the accuracy of the approximation to the
It may happen that some quadrature points are not actually used in
fitting the model (typically because the value of a covariate is
at these points). The argument
drop specifies whether these
unused quadrature points shall be deleted (
drop=FALSE) in the return value.
This function can also be applied to objects of class
- A quadrature scheme (object of class
data(cells) fit <- ppm(cells, ~1, Strauss(r=0.1)) Q <- quad.ppm(fit) plot(Q) Q$data$n Q$dummy$n