2. http://cran.r-project.org/web/packages/phangorn/
data(yeast)
tree <- NJ(dist.logDet(yeast))
fit <- pml(tree, yeast, k=4, inv= .2)
fit <- optim.pml(fit, optNni=TRUE, optGamma=TRUE, optInv=TRUE, model=“GTR”)
BS <- bootstrap.pml(fit, optNni=TRUE)
plotBS(fit$tree, BS, type = “phylogram”)
没有评论:
发表评论