integrateIntensityLocal_normal {nimbleSCR}

Integrate the multivariate normal intensity with local evaluation

Description

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

Usage

integrateIntensityLocal_normal(
  lowerCoords,
  upperCoords,
  s,
  baseIntensities,
  sd,
  numLocalWindows,
  localWindows
)

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.

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.

Value

A vector of integrated intensities over all local windows.

Author(s)

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)
 sd <- 0.1
 baseIntensities <- c(1:4)
 numLocalWindows <- 2
 localWindows <- c(1, 3)
 integrateIntensityLocal_normal(lowerCoords, upperCoords, s,
                                baseIntensities, sd,
                                numLocalWindows, localWindows)