mergeknclusters {DClusterm}

Merges clusters so that they are identifed as levels of a factor.

Description

Given a data frame with clusters that do not overlap this function merges the clusters and construct a factor. The levels of the factor are "NCL" if the polygon of the map is not in any cluster, and "CL" if the polygon i is in cluster i.

Usage

mergeknclusters(datamap, knresults, indClustersPlot)

Arguments

datamap

data of the SpatialPolygonsDataFrame with the polygons of the map.

knresults

Data frame with information of the detected clusters. Each row represents the information of one of the clusters. It contains the coordinates of the center, the size, the start and end dates, the log-likelihood ratio, a boolean indicating if it is a cluster (TRUE in all cases), and the p-value of the cluster.

indClustersPlot

rows of knresults that denote the clusters to be plotted.

Value

factor with levels that represent the clusters.

Examples

library("DClusterm")
library("RColorBrewer")

data("NY8")
data("NY8_clusters")

stcl <- mergeknclusters(NY8, ny.cl1, 1:2)
#Get first cluster
NY8$CLUSTER <- stcl

#Plot cluster
spplot(NY8, "CLUSTER", col.regions = c("white", "lightgray", "gray"))