smerc (version 1.3)

plot.smerc_cluster: Plot object of class smerc_cluster.

Description

Plot clusters (the centroids of the regions in each cluster) in different colors. The most likely cluster is plotted with solid red circles by default. Points not in a cluster are black open circles. The other cluster points are plotted with different symbols and colors.

Usage

# S3 method for smerc_cluster
plot(
  x,
  ...,
  nclusters = length(x$clusters),
  ccol = NULL,
  cpch = NULL,
  add = FALSE,
  usemap = FALSE,
  mapargs = list()
)

Arguments

x

An object of class scan to be plotted.

...

Additional graphical parameters passed to the plot function.

nclusters

Number of clusters to plot.

ccol

Fill color of the plotted points. Default is grDevices::hcl.colors(nclusters, palette = "viridis").

cpch

Plotting character to use for points in each cluster. Default is NULL, indicating pch = 20 for the most likely cluster and then pch = 2, 3, .., up to the remaining number of clusters.

add

A logical indicating whether results should be drawn on existing map.

usemap

Logical indicating whether the maps::map function should be used to create a plot background for the coordinates. Default is FALSE. Use TRUE if you have longitude/latitude coordinates.

mapargs

A list of arguments for the map function.

Examples

# NOT RUN {
data(nydf)
coords = with(nydf, cbind(longitude, latitude))
out = scan.test(coords = coords, cases = floor(nydf$cases),
                pop = nydf$pop, nsim = 0,
                longlat = TRUE, alpha = 1)
plot(out, nclusters = 3)
## plot output for new york state
# specify desired argument values
mapargs = list(database = "county", region = "new york",
               xlim = range(out$coords[,1]),
               ylim = range(out$coords[,2]))
# needed for "county" database (unless you execute library(maps))
data(countyMapEnv, package = "maps")
plot(out, nclusters = 3, usemap = TRUE, mapargs = mapargs)
# }