Hello.
I'm trying to see clusters into a certain number of classes as seeing they are
too many of them see entirely on dendrogram. So I did the following:
> d <- dist(as.matrix(medData))> hc <- hclust(d, method="average")> plot(hc)>
> plot(hc, hang = -1)
But when it came to this:> rect.hclust(hc, 3)Error in rect(m[which[n]] + 0.66,
par("usr")[3L], m[which[n] + 1] + 0.33, : plot.new has not been called yet
Here's my dput for hc. It's too long so I'm just showing you the essential
information:
structure(list(merge = structure(c(-716L, -698L, -758L, -257L, -186L, -519L,
-38L, -532L, -797L, -70L, -955L, -389L, -869L,
...........................................1058L, 1065L, 1064L, 1060L, 1067L,
1068L, 1062L, 1037L, 1071L, 1073L, 1074L), .Dim = c(1075L, 2L)), height = c(10,
25.1793566240283, 31.6631867410376, 38.2753184180093, 40.7308237088326,
42.8835632847831, ..................................398.00063581242,
421.396275641104, 428.384890986889), order = c(407L, 42L, 106L, 601L, 576L,
247L, 249L, 237L, 643L, 226L, 62L, 543L,
....................................985L, 458L, 1063L), labels = NULL, method =
"average", call = hclust(d = d, method = "average"), dist.method =
"euclidean"), .Names = c("merge", "height", "order", "labels", "method",
"call", "dist.method"), class = "hclust")
Please help,Suhaila
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list
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.
