powered by
Two more functions with different initialization method.
starts.via.svd(x, nclass = 1, method = c("em", "kmeans"), EMC = .EMC) emgroup(x, nclass = 1, EMC = .EMC)
The starts.via.svd returns an object with class svd, and the emgroup returns and object emobj with class
starts.via.svd
svd
emgroup
emobj
emret.
emret
the data matrix, dimension \(n\times p\).
the desired number of clusters, \(K\).
method with the svd initializations.
the control for the EM iterations.
Wei-Chen Chen wccsnow@gmail.com and Ranjan Maitra.
The starts.via.svd utilizes SVD to initial parameters, and the emgroup runs the EM algorithm starting from the initial.
https://www.stat.iastate.edu/people/ranjan-maitra
init.EM, .EMControl.
init.EM
.EMControl
library(EMCluster, quietly = TRUE) set.seed(1234) x1 <- da1$da emobj <- emgroup(x1, nclass = 10) summary(emobj) ret.0 <- starts.via.svd(x1, nclass = 10, method = "kmeans") summary(ret.0)
Run the code above in your browser using DataLab