This is a test version of Biostars. For the public version, visit https://www.biostars.org.
phylogenetic analyses in R program

I compared two phylogenetic trees based on two different genes via Shimodaira–Hasegawa tests in R program. Here is the code:

library(ape)
library(phangorn)
tr1 <- read.tree("ND1_besttree.out")
PCGstr <- read.tree("PCGs_besttree.out")
phyDat1 <- read.phyDat("ND1.phy")
PCGsphyDat <- read.phyDat("PCGs.phy")
fit1 <- pml(tr1, phyDat1)
fit2 <- pml(PCGstr, PCGsphyDat)
X <- SH.test(fit1, fit2, B=10000)
X

And here are the results:

     Trees      ln L      Diff ln L    p-value
[1,]   1     -24837.17    195.1926     0.0000
[2,]   2     -24641.98     0.0000      0.5079

I am a beginner in R and phylogeny. I wonder how to analyse the result and how to calculate ΔLn L? And I want to know if any significant differences between the two trees. Thanks very much for your help!

r

thank you i have beeen lookin for something similar. can you tell me what is that phyDat command for? every time i give my tree for the command "pml it says that the tree should be of class phylo" and when i convrt my tree to class phylo it diplaysError in " pml(tree2, seq1) : tip labels are not in data" please do help me any reply regarding myqueri would be very helpful

This is not an answer, please use ADD COMMENT.

Hello slx865004802!

We believe that this post does not fit the main topic of this site.

In fact this is not a question, either.

For this reason we have closed your question. This allows us to keep the site focused on the topics that the community can help with.

If you disagree please tell us why in a reply below, we'll be happy to talk about it.

Cheers!

0 answers

No answers yet.

Log in to answer this question.