Description Usage Arguments Details Value References See Also Examples

Runs joint dimension and clustering algorithms repeatedly for different numbers of clusters on bootstrap replica of the original data and returns corresponding cluster assignments, and cluster agreement indices comparing pairs of partitions.

1 2 3 4 |

`data` |
Continuous, Categorical ot Mixed data set |

`nclusrange` |
An integer or an integer vector with the number of clusters or a range of numbers of clusters (should be greater than one) |

`ndim` |
Dimensionality of the solution; if |

`method` |
Specifies the method. Options are RKM for Reduced K-means, FKM for Factorial K-means, mixedRKM for Mixed Reduced K-means, mixedFKM for Mixed Factorial K-means, MCAk for MCA K-means, iFCB for Iterative Factorial Clustering of Binary variables and clusCA for Cluster Correspondence Analysis. |

`nboot` |
Number of bootstrap pairs of partitions |

`alpha` |
Adjusts for the relative importance of (mixed) RKM and FKM in the objective function; |

`alphak` |
Non-negative scalar to adjust for the relative importance of MCA ( |

`center` |
A logical value indicating whether the metric variables should be shifted to be zero centered (default = |

`scale` |
A logical value indicating whether the metric variables should be scaled to have unit variance before the analysis takes place (default = |

`nstart` |
Number of random starts (default = 100) |

`smartStart` |
If |

`seed` |
An integer that is used as argument by |

`...` |
Not used |

The algorithm for assessing global cluster stability is similar to that in Dolnicar and Leisch (2010) and can be summarized in three steps:

*Step 1. Resampling:* Draw bootstrap samples S_i and T_i of size *n* from the data and use the original data, X, as evaluation set E_i = X. Apply the clustering method of choice to S_i and T_i and obtain C^S_i and C^T_i.

*Step 2. Mapping:* Assign each observation x_i to the closest centers of C^S_i and C^T_i using Euclidean distance, resulting in partitions C^XS_i and C^XT_i, where C^XS_i is the partition of the original data, X, predicted from clustering bootstrap sample S_i (same for T_i and C^XT_i).

*Step 3. Evaluation:* Use the Adjusted Rand Index (ARI, Hubert & Arabie, 1985) or the Measure of Concordance (MOC, Pfitzner 2008) as measure of agreement and stability.

Inspect the distributions of ARI/MOC to assess the global reproducibility of the clustering solutions.

While nboot = 100 is recommended, smaller run numbers could give quite informative results as well, if computation times become too high.

Note that the stability of a clustering solution is assessed, but stability is not the only important validity criterion - clustering solutions obtained by very inflexible clustering methods may be stable but not valid, as discussed in Hennig (2007).

`nclusrange` |
An integer or an integer vector with the number of clusters or a range of numbers of clusters |

`clust1` |
Partitions, C^XS_i of the original data, X, predicted from clustering bootstrap sample S_i (see Details) |

`clust2` |
Partitions, C^XT_i of the original data, X, predicted from clustering bootstrap sample T_i (see Details) |

`index1` |
Indices of the original data rows in bootstrap sample S_i |

`index2` |
Indices of the original data rows in bootstrap sample T_i |

`rand` |
Adjusted Rand Index values |

`moc` |
Measure of Concordance values |

Hennig, C. (2007). Cluster-wise assessment of cluster stability. *Computational Statistics and Data Analysis*, *52*, 258-271.

Pfitzner, D., Leibbrandt, R., & Powers, D. (2009). Characterization and evaluation of similarity measures for pairs of clusterings. *Knowledge and Information Systems*, *19*(3), 361-394.

Dolnicar, S., & Leisch, F. (2010). Evaluation of structure and reproducibility of cluster solutions using the bootstrap. *Marketing Letters*, *21*(1), 83-101.

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 | ```
## 5 bootstrap replicates and nstart = 1 for speed in example,
## use more for real applications
data(diamond)
boot_mixedRKM = global_bootclus(diamond[,-7], nclusrange = 3:5,
method = "mixedRKM", nboot = 5, nstart = 1, seed = 1234)
boxplot(boot_mixedRKM$rand, xlab = "Number of clusters", ylab =
"adjusted Rand Index")
## 5 bootstrap replicates and nstart = 10 for speed in example,
## use more for real applications
#data(macro)
#boot_RKM = global_bootclus(macro, nclusrange = 2:5,
#method = "RKM", nboot = 5, nstart = 10, seed = 1234)
#boxplot(boot_RKM$rand, xlab = "Number of clusters", ylab =
#"adjusted Rand Index")
## 5 bootstrap replicates and nstart = 1 for speed in example,
## use more for real applications
#data(bribery)
#boot_cluCA = global_bootclus(bribery, nclusrange = 2:5,
#method = "clusCA", nboot = 5, nstart = 1, seed = 1234)
#boxplot(boot_cluCA$rand, xlab = "Number of clusters", ylab =
#"adjusted Rand Index")
``` |

clustrd documentation built on May 8, 2019, 5:03 p.m.

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.