聚类分析-R语言
2018-07-28 本文已影响14人
育种数据分析之放飞自我
1,原始数据,是矩阵,有行头和列名:
![](https://img.haomeiwen.com/i4277952/5d027fdba6f1c2f1.png)
2,amap包的聚类分析:
library(amap)
clu <- hclusterpar(matx)
plot(clu,sub="",hang = -1,xlab = NA,ylab = NA,main = NA)
![](https://img.haomeiwen.com/i4277952/9f01cf5cacf9b9b8.png)
用法:
hcluster(x, method = "euclidean", diag = FALSE, upper = FALSE,
link = "complete", members = NULL, nbproc = 2,
doubleprecision = TRUE)
3,cluster包
library(cluster)
agnx <- agnes(matx,method = "complete")
pltree(agnx)
agnes(x, diss = inherits(x, "dist"), metric = "euclidean",
stand = FALSE, method = "average", par.method,
keep.diss = n < 100, keep.data = !diss, trace.lev = 0)
默认方法不一样,这里改为complete
![](https://img.haomeiwen.com/i4277952/40ed8e89aecf8efc.png)
转化为横着排放的格式:
dagn <- as.dendrogram(as.hclust(agnx))
plot(dagn, horiz = TRUE, center = TRUE,
nodePar = list(lab.cex = 0.6, lab.col = "forest green", pch = NA),
main = deparse(agn$call))
![](https://img.haomeiwen.com/i4277952/17f1803158ce23b5.png)