qualitySet: Goodness of classifications for a set of k clusters.

In evaluomeR: Evaluation of Bioinformatics Metrics

Description

The goodness of the classifications are assessed by validating the clusters generated for a range of k values. For this purpose, we use the Silhouette width as validity index. This index computes and compares the quality of the clustering outputs found by the different metrics, thus enabling to measure the goodness of the classification for both instances and metrics. More precisely, this measurement provides an assessment of how similar an instance is to other instances from the same cluster and dissimilar to the rest of clusters. The average on all the instances quantifies how the instances appropriately are clustered. Kaufman and Rousseeuw suggested the interpretation of the global Silhouette width score as the effectiveness of the clustering structure. The values are in the range [0,1], having the following meaning:

• There is no substantial clustering structure: [-1, 0.25].

• The clustering structure is weak and could be artificial: ]0.25, 0.50].

• There is a reasonable clustering structure: ]0.50, 0.70].

• A strong clustering structure has been found: ]0.70, 1].

Usage

qualitySet(data, k.set = c(2, 4), cbi = "kmeans", getImages = TRUE,
  seed = NULL)

Arguments

data	A SummarizedExperiment. The SummarizedExperiment must contain an assay with the following structure: A valid header with names. The first column of the header is the ID or name of the instance of the dataset (e.g., ontology, pathway, etc.) on which the metrics are measured. The other columns of the header contains the names of the metrics. The rows contains the measurements of the metrics for each instance in the dataset.
k.set	A list of integer values of k, as in c(2,4,8). The values must be contained in [2,15] range.
cbi	Clusterboot interface name (default: "kmeans"): "kmeans", "clara", "clara_pam", "hclust", "pamk", "pamk_pam", "pamk". Any CBI appended with '_pam' makes use of pam. The method used in 'hclust' CBI is "ward.D2".
getImages	Boolean. If true, a plot is displayed.
seed	Positive integer. A seed for internal bootstrap.

Value

A list of SummarizedExperiment containing the silhouette width measurements and cluster sizes from k.set.

References

\insertRef

kaufman2009findingevaluomeR

Examples

# Using example data from our package
data("rnaMetrics")
# Without plotting
dataFrameList = qualitySet(rnaMetrics, k.set=c(2,3), getImages = FALSE)

evaluomeR documentation built on March 15, 2021, 6 p.m.