zones: Create set of all single zones and output geographical...

In SpatialEpi: Methods and Data for Spatial Epidemiology

Create set of all single zones and output geographical information

Description

Based on the population counts and centroid coordinates of each of n areas, output the set of n.zones single zones as defined by Kulldorff and other geographical information.

Usage

zones(geo, population, pop.upper.bound)

Arguments

geo	n x 2 table of the (x,y)-coordinates of the area centroids
population	a vector of population counts of each area
pop.upper.bound	maximum proportion of study region each zone can contain

Value

A list containing

nearest.neighbors	list of n elements, where each element is a vector of the nearest neighbors in order of distance up until pop.upper.bound of the total population is attained
cluster.coords	n.zones x 2 table of the center and the radial area for each zone
dist	n x n inter-point distance matrix of the centroids

Author(s)

Albert Y. Kim

References

Kulldorff, M. (1997) A spatial scan statistic. Communications in Statistics: Theory and Methods, 26, 1481–1496. Kulldorff M. and Nagarwalla N. (1995) Spatial disease clusters: Detection and Inference. Statistics in Medicine, 14, 799–810.

Examples

data(pennLC)
geo <- pennLC$geo[,2:3]
geo <- latlong2grid(geo)
population <- tapply(pennLC$data$population, pennLC$data$county, sum)
pop.upper.bound <- 0.5
geo.info <- zones(geo, population, pop.upper.bound)

SpatialEpi documentation built on March 7, 2023, 8 p.m.