大家好!我打算在R中用ape包进行简单的聚类分析,具体的操作代码如下:
library(ape)
trnwk <- "((((Homo:0.21,Pongo:0.21):0.28,Macaca:0.49):0.13,Ateles:0.62)"
trnwk[2] <- ":0.38,Galago:1.00);"
tr <- read.tree(text = trnwk)
plot(tr)
axisPhylo()
这个代码是在参考书上找到的,我的数据和参考书上不同,我的物种数是15个,我有每两个物种的基于cytb分析的遗传距离,不知道该怎么将这些成对的遗传距离读入R程序中,做出聚类图,麻烦各位高手指点一二,先谢谢了!