integrateIntensity_normal {nimbleSCR}

Integrate the multivariate normal intensity

Description

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

Usage

integrateIntensity_normal(
  lowerCoords,
  upperCoords,
  s,
  baseIntensities,
  sd,
  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 isotropic multivariate normal distribution mean (AC location).

baseIntensities

Vector of baseline intensities for all windows.

sd

Standard deviation of the isotropic multivariate normal 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)
sd <- 0.1
baseIntensities <- c(1:4)
numWindows <- 4
integrateIntensity_normal(lowerCoords, upperCoords, s, baseIntensities, sd, numWindows)