powered by
Simple transformation of a clustering or grouping to a normal mixture object (class "norMix", see, norMix.
"norMix"
norMix
clus2norMix(gr, x, name = deparse(sys.call()))
A call to norMix() with (mu, sig2, w) set to the empirical values of the groups (as defined by split(x,gr).
norMix()
(mu, sig2, w)
split(x,gr)
a grouping/clustering vector with values in \(\{1,..,K\}\); possibly a factor.
factor
numeric vector of (original) data (of the same length as gr).
gr
name for norMix() object; constructed from the call by default.
Martin Maechler, Dec. 2007
norMix; further pam() (or clara()) from package cluster for sensible clusterings.
pam()
clara()
x9 <- rnorMix(500, MW.nm9) require("cluster") pxc <- pam(x9, k=3) plot(pxc, which = 2)# silhouette (nm.p9 <- clus2norMix(pxc$clustering, x9)) plot(nm.p9, p.norm=FALSE) lines(MW.nm9, col="thistle")
Run the code above in your browser using DataLab