Dear all, I would like to draw a dendrogram and mark some parts/branches (by using "segments") including their labels. If I draw it without specifying the length of x axix, I am able to do that (as in My dendrogram 1 of the following codes). However, if I want to specify the x axix, I am not able to draw marking line (by using "segments") including labels (as in My dendrogram 2 of the following codes). Is there a way that I can include my labels into as well? Any help is deeply appreciated.
Kind Regards Seyit Ali Codes: set.seed(201) x<-matrix(rnorm(100, 0, 1), ncol=20, byrow=TRUE) myclust<- hclust(dist(cor(x), method='euclidian'), method="ave") mydend<-as.dendrogram(myclust) par(mfrow=c(1,2), mar=c(5,4,4,3), oma=c(1, 3, 1, 3)) plot(mydend, xlim=c(4, -0.2), horiz = TRUE, main="My Dendrogram 1") x0<- -0.35 y0<- 0.5 x1<- 1.8 y1<- 0.5 segments(x0, y0, x1, y1, col="red", lty=2) x0<- -0.35 y0<- 4.4 x1<- 1.8 y1<- 4.4 segments(x0, y0, x1, y1, col="red", lty=2) x0<- -0.35 y0<- 0.5 x1<- -0.35 y1<- 4.4 segments(x0, y0, x1, y1, col="red", lty=2) x0<- 1.8 y0<- 0.5 x1<- 1.8 y1<- 4.4 segments(x0, y0, x1, y1, col="red", lty=2) ################################################################################ plot(mydend, xlim=c(4, 0.5), horiz = TRUE, main="My Dendrogram 2", dLeaf=0.35)#, yaxt="n") x0<- 0.3 y0<- 0.5 x1<- 1.8 y1<- 0.5 segments(x0, y0, x1, y1, col="red", lty=2) x0<- 0.3 y0<- 4.4 x1<- 1.8 y1<- 4.4 segments(x0, y0, x1, y1, col="red", lty=2) x0<- 0.3 y0<- 0.5 x1<- 0.3 y1<- 4.4 segments(x0, y0, x1, y1, col="red", lty=2) x0<- 1.8 y0<- 0.5 x1<- 1.8 y1<- 4.4 segments(x0, y0, x1, y1, col="red", lty=2) ------------------------------------------------------------------------------------------------------------------ Dr. Seyit Ali KAYIS Selcuk University Faculty of Agriculture Kampus, Konya, TURKEY s_a_ka...@yahoo.com, s_a_ka...@hotmail.com Tell: +90 332 223 2830 Mobile: +90 535 587 1139 Fax: +90 332 241 0108 Greetings from Konya, TURKEY http://www.ziraat.selcuk.edu.tr/skayis/ ---------------------------------------------------------------------------------------------------------------------- _________________________________________________________________ [[elided Hotmail spam]] [[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.