Estimating Neighbor Joining Tree using GTR model
1
0
Entering edit mode
5.0 years ago
mm8511 • 0

I am wondering if anyone can help me figure out how to build a neighbor joining tree in R, where the pairwise distances (from sequence alignments) are estimated under the GTR model. I have tried the dist.dna() method from R's phanghorn package, but it does not have support for the GTR model. I am new to phylogenetics, and any help would be great.

NJ R GTR phylogenetics • 2.1k views
ADD COMMENT
0
Entering edit mode
5.0 years ago
Brice Sarver ★ 3.6k

While you can estimate distances under a handful of models in R, to my knowledge you are only able to estimate trees under the GTR model when using likelihood-based approaches (as, for example, in phangorn's optim.pml() and others). Unfortunately, I think you'll have to do this outside of R. Most sophisticated phylogenetic estimation isn't done in R, anyway.

I don't have time to skim through the PAUP manual and provide detailed instructionsl, but you can calculate distances under the GTR model and estimate a tree generally following this tutorial.

ADD COMMENT
0
Entering edit mode

This is great. I really appreciate it

ADD REPLY

Login before adding your answer.

Traffic: 1691 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6