使用R语言构建简单的NJ树

熊荣川

xiong rongchuan

六盘水师范学院生物信息学实验室

xiongrongchuan@126.com

http://blog.sciencenet.cn/u/Bearjazz

path = "D:/ziliao/优先阅读/R语言学习/R语言与系统发育/使用R语言构建nj树"

setwd(path)

library(ape)

library(seqinr)

Aln = read.alignment(file ="Chou.fas", format = "fasta") 

Dnj = dist.alignment(Aln, matrix =c("similarity", "identity")) #计算距离矩阵

Otr = nj(Dnj) #使用 Saitou andNei (1987)方法构建NJ树

plot(Otr)

Otr2 = bionj(Dnj) #Gascuel (1997)构建NJ树

plot(Otr2)