cholera (version 0.8.0)

neighborhoodVoronoi: Compute Voronoi pump neighborhoods.

Description

Group cases into neighborhoods using Voronoi tessellation.

Usage

neighborhoodVoronoi(pump.select = NULL, vestry = FALSE,
  case.location = "address", pump.location = "nominal",
  polygon.vertices = FALSE)

Value

An R list with 12 objects.

• pump.id: vector of selected pumps

• voronoi: output from deldir::deldir().

• snow.colors: neighborhood color based on snowColors().

• x.rng: range of x for plot.

• y.rng: range of y for plot.

• select.string: description of "pump.select" for plot title.

• expected.data: expected neighborhood fatality counts, based on Voronoi cell area.

• coordinates: polygon vertices of Voronoi cells.

• statistic.data: observed neighborhood fatality counts.

• pump.select: "pump.select" from neighborhoodVoronoi().

• statistic: "statistic" from neighborhoodVoronoi().

• vestry: "vestry" from neighborhoodVoronoi().

Arguments

pump.select

Numeric. Vector of numeric pump IDs to define pump neighborhoods (i.e., the "population"). Negative selection possible. NULL selects all pumps.

vestry

Logical. TRUE uses the 14 pumps from the Vestry report. FALSE uses the 13 in the original map.

case.location

Character. "address" or "nominal". "address" uses the x-y coordinates of ortho.proj. "nominal" uses the x-y coordinates of fatalities.

pump.location

Character. "address" or "nominal". "address" uses the x-y coordinates of ortho.proj.pump or ortho.proj.pump.vestry. "nominal" uses the x-y coordinates of pumps or pumps.vestry.

polygon.vertices

Logical. TRUE returns a list of x-y coordinates of the vertices of Voronoi cells. Useful for sp::point.in.polygon() as used in print.voronoi() method.

Examples

neighborhoodVoronoi()
neighborhoodVoronoi(vestry = TRUE)
neighborhoodVoronoi(pump.select = 6:7)
neighborhoodVoronoi(pump.select = -6)
neighborhoodVoronoi(pump.select = -6, polygon.vertices = TRUE)

# coordinates for vertices also available in the returned object.
dat <- neighborhoodVoronoi(pump.select = -6)
dat$coordinates