Michael Artz
2016-Apr-12 01:30 UTC
[R] Dissimilarity matrix and number clusters determination
Hi,
I already have a dissimilarity matrix and I am submitting the results to
the elbow.obj method to get an optimal number of clusters. Am I reading
the below output correctly that I should have 17 clusters?
code:
top150 <- sampleset[1:150,]
{cluster1 <- daisy(top150
, metric = c("gower")
, stand = TRUE
, type = list(symm = 1))
}
dist.obj <- dist(cluster1)
hclust.obj <- hclust(dist.obj)
css.obj <- css.hclust(dist.obj,hclust.obj)
elbow.obj <- elbow.batch(css.obj)
[1] "A \"good\" k=17 (EV=0.80) is detected when the EV is no less
than
0.8\nand the increment of EV is no more than 0.01 for a bigger k.\n"
attr(,"class")
[[alternative HTML version deleted]]
Luisfo Chiroque
2016-Apr-12 21:37 UTC
[R] Dissimilarity matrix and number clusters determination
Dear Michael, Yes, AFAIK you are correctly reading the results. You can print elbow.obj$k to obtain the optimal number of clusters, and ?visually? you can check it plotting the variance vs #clusters plot(css.obj$k, css.obj$ev) HTH Best, Luisfo Chiroque PhD Student IMDEA Networks Institute http://fourier.networks.imdea.org/people/~luis_nunez/ <http://fourier.networks.imdea.org/people/~luis_nunez/>> El 12 abr 2016, a las 4:30, Michael Artz <michaeleartz at gmail.com> escribi?: > > Hi, > I already have a dissimilarity matrix and I am submitting the results to > the elbow.obj method to get an optimal number of clusters. Am I reading > the below output correctly that I should have 17 clusters? > > code: > top150 <- sampleset[1:150,] > {cluster1 <- daisy(top150 > , metric = c("gower") > , stand = TRUE > , type = list(symm = 1)) > } > > dist.obj <- dist(cluster1) > hclust.obj <- hclust(dist.obj) > css.obj <- css.hclust(dist.obj,hclust.obj) > elbow.obj <- elbow.batch(css.obj) > > [1] "A \"good\" k=17 (EV=0.80) is detected when the EV is no less than > 0.8\nand the increment of EV is no more than 0.01 for a bigger k.\n" > attr(,"class") > > [[alternative HTML version deleted]] > > ______________________________________________ > R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code.[[alternative HTML version deleted]]