R语言画聚类分析树形图

以iris数据集为例：

data(iris)

newiris<-iris

newiris$Species<-NULL

library(colorspace)

require(amap, quietly=TRUE)

require(ggplot2, quietly=TRUE)

require(ggdendro, quietly=TRUE)

hclust <- hclusterpar(na.omit(newiris), method="euclidean", link="ward", nbproc=2)#nbproc为处理器数,系统聚类分析

ddata <- dendro_data(hclust, type="rectangle")
g <- ggplot(segment(ddata))
g <- g + geom_segment(aes(x = y, y = x, xend = yend, yend = xend))#将x和y互换，即将树的形状转90度
g <- g + scale_y_discrete(labels = ddata$label$label)#y轴的显示，即x的值，如公司名称等
g <- g + labs(x="高度", y="公司")#x和y轴的label
print(g)

2、枝叶在横轴上的树形图

源码：