stabilitySet: Stability index for a set of k clusters.

In neobernad/evaluomeR: Evaluation of Bioinformatics Metrics

Stability index for a set of k clusters.

Description

This analysis permits to estimate whether the clustering is meaningfully affected by small variations in the sample. For a set of k values (k.set), a clustering using the k-means algorithm is carried out. Then, the stability index is the mean of the Jaccard coefficient values of a number of bs bootstrap replicates. The values are in the range [0,1], having the following meaning:

• Unstable: [0, 0.60[.

• Doubtful: [0.60, 0.75].

• Stable: ]0.75, 0.85].

• Highly Stable: ]0.85, 1].

Usage

stabilitySet(
  data,
  k.set = c(2, 3),
  bs = 100,
  cbi = "kmeans",
  getImages = FALSE,
  all_metrics = FALSE,
  seed = NULL,
  gold_standard = 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.
bs	Positive integer. Bootstrap value to perform the resampling.
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.
all_metrics	Boolean. If true, clustering is performed upon all the dataset.
seed	Positive integer. A seed for internal bootstrap.

Value

A ExperimentList containing the stability and cluster measurements of the list of k clusters.

References

\insertRef

milligan1996measuringevaluomeR

\insertRef

jaccard1901distributionevaluomeR

Examples

# Using example data from our package
data("rnaMetrics")
result <- stabilitySet(rnaMetrics, k.set=c(2,3))

neobernad/evaluomeR documentation built on Nov. 12, 2024, 11:29 p.m.