spatstat.model (version 3.7-2)

vcov.dppm: Variance-Covariance Matrix for a Fitted Determinantal Point Process Model

Description

Returns the variance-covariance matrix of the estimates of the parameters of a fitted determinantal point process model.

Usage

# S3 method for dppm
vcov(object, ...,
          what=c("vcov", "corr", "fisher"),
          fast = NULL, rmax = NULL, eps.rmax = 0.01,
          verbose = TRUE)

Value

A square matrix.

Arguments

object

A fitted determinantal point process model (an object of class "dppm".)

...

Ignored.

what

Character string (partially-matched) that specifies what matrix is returned. Options are "vcov" for the variance-covariance matrix, "corr" for the correlation matrix, and "fisher" for the Fisher information matrix.

fast

Logical specifying whether tapering (using sparse matrices from Matrix) should be used to speed up calculations. Warning: This is expected to underestimate the true asymptotic variances/covariances.

rmax

Optional. The dependence range. Not usually specified by the user. Only used when fast=TRUE.

eps.rmax

Numeric. A small positive number which is used to determine rmax from the tail behaviour of the pair correlation function when fast option (fast=TRUE) is used. Namely rmax is the smallest value of \(r\) at which \((1 - g(r))/(1 - g(0))\) falls below eps.rmax. Only used when fast=TRUE. Ignored if rmax is provided.

verbose

Logical value indicating whether to print progress reports during very long calculations.

Author

Abdollah Jalilian and Rasmus Waagepetersen. Ported to spatstat by Adrian Baddeley and Ege Rubak . Adapted for dppm by Adrian Baddeley .

Details

This function computes the asymptotic variance-covariance matrix of the estimates of the canonical (regression) parameters in the cluster point process model object. It is a method for the generic function vcov.

The result is an n * n matrix where n = length(coef(model)).

To calculate a confidence interval for a regression parameter, use confint as shown in the examples.

References

Waagepetersen, R. (2007) Estimating functions for inhomogeneous spatial point processes with incomplete covariate data. Biometrika 95, 351--363.

See Also

dppm, vcov, vcov.ppm

Examples

   fit <- dppm(japanesepines ~ x, dppGauss)
   vcov(fit)
   vcov(fit, what="corr")

   # confidence interval (95 percent)
   confint(fit)
   # cross-check the confidence interval by hand:
   sd <- sqrt(diag(vcov(fit)))
   coef(fit) + qnorm(0.975) * outer(sd, c(lower=-1, upper=1))