SpatialEpi: zones – R documentation

Become an expert in R — Interactive courses, Cheat Sheets, certificates and more!

zones

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

Methods and Data for Spatial Epidemiology

v1.2.3

GPL-2

Authors

Cici Chen [ctb], Albert Y. Kim [aut, cre], Michelle Ross [ctb], Jon Wakefield [aut]

Initial release