利用TM-score矩阵绘制进化树图

#从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)

评论 4
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值