Dear All
I am constructing maximum parsimony trees:
library(phangorn)
primates = read.phyDat("primates.dna", format="phylip", type="DNA")
dm = dist.dna(as.DNAbin(primates))
treeNJ = NJ(dm)
treePars = optim.parsimony(treeNJ, primates)
##and
treeRatchet = pratchet(primates, trace = 0)
, and I would like to know how to do parsimony bootstrap analysis. Anyone know
how this can be done?
Regards,
Brian.
[[alternative HTML version deleted]]
_______________________________________________
R-sig-phylo mailing list
R-sig-phylo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-phylo
