integrateIntensity_exp {nimbleSCR}

Integrate the multivariate exponential intensity

Description

Calculate the integral of the intensity function with an isotropic multivariate exponential kernel over a set of windows.

Usage

integrateIntensity_exp(
  lowerCoords,
  upperCoords,
  s,
  baseIntensities,
  lambda,
  numWindows
)

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 AC location.

baseIntensities

Vector of baseline intensities for all windows.

lambda

Rate parameter of the isotropic multivariate exponential distribution.

numWindows

Total number of windows. This value (positive integer) is used to truncate lowerCoords and upperCoords so that extra rows beyond numWindows are ignored.

Value

A vector of integrated intensities over all windows.

Author(s)

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(1, 1)
lambda <- 1.0
baseIntensities <- c(1:4)
numWindows <- 4
integrateIntensity_exp(lowerCoords, upperCoords, s, baseIntensities, lambda, numWindows)