kmeansFit¶
Purpose¶
Partitions data into k clusters, using the kmeans algorithm.
Format¶
-
mdl =
kmeansFit
(X, clusters[, ctl])¶ - Parameters
X (NxP matrix) – The training data.
clusters (Scalar) – The number of clusters, or a matrix containing the initial centroids.
ctl –
Optional input, an instance of a
kmeansControl
structure.ctl.initMethod
Scalar specifying the algorithm used to create the initial centroids. Options include:
0
kmeans++ (default).
1
parallel k-means++
2
\(k\) randomly-selected observations.
ctl.nStarts
Scalar, the number of times to run the kmeans algorithm with new starting centroids. Note: this input will be ignored if the clusters input is a starting centroid.
ctl.seed
Seed for the random number generator which creates the initial centroids. Note: this input will be ignored if the clusters input is a starting centroid.
ctl.tolerance
Scalar, the convergence tolerance for the kmeans algorithm.
ctl.maxIters
Scalar, the maximum number of iterations to allow each of the nStarts to run before forcing convergence.
- Returns
mdl (struct) –
An instance of a
kmeansModel
structure.mdl.centroids
kxP matrix, containing the centroids with the lowest intra-cluster sum of squares.
mdl.assignments
Nx1 matrix, containing the centroid assignment for the corresponding observation of the input matrix.
mdl.clusterSS
Scalar, sum of squared differences between each observation and its assigned centroid.
mdl.elapsedIters
Scalar, the number of iterations taken by the start with the lowest clusterSS.
Examples¶
new;
library gml;
// Get dataset with full name
fname = getGAUSSHome() $+ "pkgs/gml/examples/iris.csv";
// Load data
X = loadd(fname, ". -species");
// For repeatable sample
rndseed 234234;
// Split data into x_train and x_test
{ x_train, x_test } = splitData(X, 0.70);
// Number of clusters
n_clusters = 3;
// Declare kmeansModel struct
struct kmeansModel mdl;
// Fit kmeans model
mdl = kmeansFit(x_train , n_clusters);
// Print the k centroids
print mdl.centroids;
The above code will print the following:
5.824 2.710 4.350 1.421
5.009 3.406 1.475 0.250
6.887 3.061 5.732 2.061
References¶
Parallel Kmeans++ initialization. B. Bahmani, B. Moseley, A. Vattani, R. Kumar, S. Vassilvitskii. Scalable K-means++. Proceedings of the VLDB Endowment, 2012.
See also
kmeansFit()
, kmeansControlCreate()