以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、枝叶在横轴上的树形图
源码: