integrateIntensity_exp: Integrate the multivariate exponential intensity

In nimbleSCR: Spatial Capture-Recapture (SCR) Methods Using 'nimble'

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)

nimbleSCR documentation built on Dec. 1, 2022, 1:17 a.m.