#从excel中导入花一下午输入的TM-score的矩阵;这里有一个问题:理论上由一半的数据,另一半需要吗?
#安装Rtool
writeLines('PATH="${RTOOLS40_HOME}\\usr\\bin;${PATH}"', con = "~/.Renviron")
Sys.which("make")
#使用stats包中的hclust函数,其中可直接进行HPGMA聚类;导入的excel中的metrix需要先用dist函数处理
library('stats');help(hclust)
data1 <- dist(data)
data2 <- hclust(data1, "average")
METTL.name <- c("1","2","2B","3","4","5","6","7A","7B","8","9","10","11A","11B",
"12","13","14","15","16","17","18","19","20","21A","21B","21C",
"21D","22","23","24","25","25B","26","27")
plot(data2,hang=1,cex=0.5,labels = METTL.name)
利用TM-score矩阵绘制进化树图
最新推荐文章于 2025-05-21 10:13:51 发布