Displaying 1 result from an estimated 1 matches for "srndcents20000".
2011 Sep 12
1
hclust and cutree: identifying branches as classes
...t(dist(srndpc$x[1:1000,1:3]),method="ward") #creation of hierarchical tree
plclust(srndpchc,hmin=20000) #visualisation
srndpchc20000 = cutree(srndpchc,h=20000) #returns 4 classes
table(srndpchc20000 )
srndclass20000 = cbind(srnd@data[1:1000,],srndpchc20000) #assigning classes to objects
srndcents20000 <- aggregate(srndclass20000, by=list(srndpchc20000), FUN=median)
matplot(1:36,t(srndcents20000[,-c(1,38)]))
But how can we make sure that, for example, class 1 is the first branch in the tree plottedĀ by plclust() ?
Regards
Laurent
[[alternative HTML version deleted]]