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

Description Usage Arguments Value Author(s) References Examples

View source: R/zones.R

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

1
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

1
2
3
4
5
6
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 May 20, 2017, 3:11 a.m.
Search within the SpatialEpi package
Search all R packages, documentation and source code

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.

Please suggest features or report bugs in the GitHub issue tracker.

All documentation is copyright its authors; we didn't write any of that.