color.clusters: Color clusters

In smerc: Statistical Methods for Regional Counts

Color clusters

Description

color.clusters is a helper function to color clusters of regions produced by an appropriate method, e.g., scan.test or uls.test. Regions that are not part of any cluster have no color.

Usage

color.clusters(
  x,
  idx = seq_along(x$clusters),
  col = grDevices::hcl.colors(length(idx))
)

Arguments

x	An object of class scan produced by a function such as scan.test.
idx	An index vector indicating the elements of object$clusters to print information for. The default is all clusters.
col	A vector of colors to color the clusters in x. Should have same length as the number of clusters in x.

Value

Returns a vector with colors for each region/centroid for the data set used to construct x.

Author(s)

Joshua French

Examples

data(nydf)
coords <- with(nydf, cbind(longitude, latitude))
out <- scan.test(
  coords = coords, cases = floor(nydf$cases),
  pop = nydf$pop, alpha = 0.2, longlat = TRUE,
  nsim = 9
)
data(nypoly)
library(sp)
# plot all clusters
plot(nypoly, col = color.clusters(out), axes = TRUE)
# zoom in on small cluster
plot(nypoly,
  col = color.clusters(out),
  xlim = c(400000, 450000),
  ylim = c(4750000, 4800000)
)
# plot only clusters 1 and 3
plot(nypoly, col = color.clusters(out, idx = c(1, 3)))

smerc documentation built on Oct. 13, 2022, 9:07 a.m.