Hi Kevin

If you look carefully at the figure on the web site, you should see that both trees have the same topology: they are plotted in different ways. You can also see that the screenshot was taken some time ago: the code of nj() has been improved since then to avoid this (among other things).

Best,

Emmanuel

Le 18/03/2014 12:13, Kevin Livingstone a écrit :
I'm trying to recreate the screenshot that illustrates the different
trees obtained by the K80 and TN93 models
<http://ape-package.ird.fr/image/sc3.png> for a class exercise, but both
of my output trees are the same. Here's the code I'm entering:

 > library(ape)
 > data(woodmouse)
 > M.K80 <- dist.dna(woodmouse,model="K80")
 > M.TN93 <- dist.dna(woodmouse, model="TN93")
 > tr.K80 <- nj(M.K80)
 > tr.TN93 <- nj(M.TN93)
 > layout(matrix(1:2,1))
 > plot(tr.K80)
 > plot(tr.TN93)

I have attached the output figure, and both trees appear to match the
TN93 tree from the original figure. I'm concerned this means I have
missed or mis-installed something big.

thanks!
-Kevin


_______________________________________________
R-sig-phylo mailing list - [email protected]
https://stat.ethz.ch/mailman/listinfo/r-sig-phylo
Searchable archive at http://www.mail-archive.com/[email protected]/


_______________________________________________
R-sig-phylo mailing list - [email protected]
https://stat.ethz.ch/mailman/listinfo/r-sig-phylo
Searchable archive at http://www.mail-archive.com/[email protected]/

Reply via email to