nimbleSCR (version 0.2.1)

integrateIntensityLocal_exp: Integrate the multivariate exponential intensity with local evaluation

Description

Calculate the integral of the intensity function with an isotropic multivariate exponential kernel over a set of windows. The local evaluation technique is implemented.

Usage

integrateIntensityLocal_exp(
  lowerCoords,
  upperCoords,
  s,
  baseIntensities,
  lambda,
  numLocalWindows,
  localWindows
)

Value

A vector of integrated intensities over all local windows.

Arguments

lowerCoords, upperCoords

Matrices of lower and upper x- and y-coordinates of a set of windows. One row for each window.

s

Vector of x- and y-coordinates of the isotropic multivariate exponential distribution mean (AC location).

baseIntensities

Vector of baseline intensities for all windows.

lambda

Rate parameter of the isotropic multivariate exponential distribution.

numLocalWindows

Number of windows that are close to the activity center

localWindows

Vector of indices of the windows that are close to the activity center.

Author

Cyril Milleret and Wei Zhang

References

W. Zhang, J. D. Chipperfield, J. B. Illian, P. Dupont, C. Milleret, P. de Valpine and R. Bischof. 2020. A hierarchical point process model for spatial capture-recapture data. bioRxiv. DOI 10.1101/2020.10.06.325035

Examples

lowerCoords <- matrix(c(0, 0, 1, 0, 0, 1, 1, 1), nrow = 4, byrow = TRUE)
upperCoords <- matrix(c(1, 1, 2, 1, 1, 2, 2, 2), nrow = 4, byrow = TRUE)
s <- c(0.1, 0.9)
lambda <- 0.1
baseIntensities <- c(1:4)
numLocalWindows <- 2
localWindows <- c(1, 3)
integrateIntensityLocal_exp(lowerCoords, upperCoords,
                            s, baseIntensities, lambda, 
                            numLocalWindows, localWindows)