discretize.jointly: Discretize Multivariate Continuous Data by a...
In GridOnClusters: Cluster-Preserving Multivariate Joint Grid Discretization

Description Usage Arguments Details Value Author(s) References See Also Examples

View source: R/discretize_jointly.R

Discretize multivariate continuous data using a grid that captures the joint distribution via preserving clusters in the original data

discretize.jointly(
  data,
  k = c(2:10),
  min_level = 1,
  cluster_method = c("Ball+BIC", "kmeans+silhouette", "PAM"),
  grid_method = c("Sort+split", "MultiChannel.WUC"),
  cluster_label = NULL
)

`data`	a matrix containing two or more continuous variables. Columns are variables, rows are observations.
`k`	either an integer, a vector of integers, or `Inf`, specifying different ways to find clusters in `data`. The default is a vector containing integers from 2 to 10. If 'k' is a single number, `data` will be grouped into into exactly 'k' clusters. If 'k' is an integer vector, an optimal 'k' is chosen from among the integers, that maximizes the average silhouette width. If 'k' is set to `Inf`, an optimal k is chosen among 2 to `nrow(data)`. If `cluster_label` is specified, `k` is ignored.
`min_level`	integer or vector, signifying the minimum number of levels along each dimension. If a vector of size `ncol(data)`, then each element will be mapped 1:1 to each dimension in order. If an integer, then all dimensions will have the same minimum number of levels.
`cluster_method`	the clustering method to be used. Ignored if cluster labels are given "kmeans+silhouette" will use k-means to cluster `data` and the average Silhouette score to select the number of clusters k. "Ball+BIC" will use Mclust (modelNames = "VII") to cluster `data` and BIC score to select the number of cluster k.
`grid_method`	the discretization method to be used. "Sort+split" will sort the cluster by cluster mean in each dimension and then split consecutive pairs only if the sum of the error rate of each cluster is less than or equal to 50 in a certain dimension. The maximum number of lines is the number of clusters minus one. "MultiChannel.WUC" will split each dimension by weighted with-in cluster sum of squared distances by "Ckmeans.1d.dp::MultiChannel.WUC". Applied in each projection on each dimension. The channel of each point is defined by its multivariate cluster label.
`cluster_label`	a vector of user-specified cluster labels for each observation in `data`. The user is free to choose any clustering. If unspecified, k-means clustering is used by default.

The function implements algorithms described in \insertCiteJwang2020BCBGridOnClusters.

A list that contains four items:

`D`	a matrix that contains the discretized version of the original `data`. Discretized values are one(1)-based.
`grid`	a list of vectors containing decision boundaries for each variable/dimension.
`clabels`	a vector containing cluster labels for each observation in `data`.
`csimilarity`	a similarity score between clusters from joint discretization `D` and cluster labels `clabels`. The score is the adjusted Rand index.

Jiandong Wang, Sajal Kumar and Mingzhou Song

\insertAllCited

See Ckmeans.1d.dp for discretizing univariate continuous data.

# using a specified k
x = rnorm(100)
y = sin(x)
z = cos(x)
data = cbind(x, y, z)
discretized_data = discretize.jointly(data, k=5)$D

# using a range of k
x = rnorm(100)
y = log1p(abs(x))
z = tan(x)
data = cbind(x, y, z)
discretized_data = discretize.jointly(data, k=c(3:10))$D

# using k = Inf
x = c()
y = c()
mns = seq(0,1200,100)
for(i in 1:12){
  x = c(x,runif(n=20, min=mns[i], max=mns[i]+20))
  y = c(y,runif(n=20, min=mns[i], max=mns[i]+20))
}
data = cbind(x, y)
discretized_data = discretize.jointly(data, k=Inf)$D

# using an alternate clustering method to k-means
library(cluster)
x = rnorm(100)
y = log1p(abs(x))
z = sin(x)
data = cbind(x, y, z)

# pre-cluster the data using partition around medoids (PAM)
cluster_label = pam(x=data, diss = FALSE, metric = "euclidean", k = 5)$clustering
discretized_data = discretize.jointly(data, cluster_label = cluster_label)$D