compute: Compute exceedance probability in a cylinder.

In mcavallaro/rancovr: Cluster detection in R with RAndom Neighbourhood COVeRing

Compute exceedance probability in a cylinder.

Description

A cylinder is defined by the circle coordinated (say, x,y, and radius) and lower and upper height limits (aay, t.low and t.upp, respectively). For a given cylinder, this function computes the number of observed events (n_cases) in the cylinder according to observation.matrix, the expected number mu of events according the Poisson point model (with intensity defined in baseline.matrix), and the probability that . The function returns c(n_cases, mu, p.val).

Usage

compute(cylinder, observation.matrix, baseline.matrix, postcode.locations)

Arguments

cylinder
observation.matrix	A sparseMatrix object encoding the events.
baseline.matrix	A Matrix object encoding the baseline.
postcode.locations	A data.frame that maps the rows of observation.matrix to geographical coordinates.

Value

A numeric vector of dimension 3.

Examples

exceedance=compute(c(x,y,rho,t.low,t.upp), observation.matrix, baseline.matrix, postcode.locations)

mcavallaro/rancovr documentation built on April 17, 2025, 7:21 p.m.