Dear All,
Sorry for the novice question. I have gone through the PAML documentation but it is still very confusing for me.
I have two lineages in my phylogenetic tree. I suspect positive selection in both these lineages, one of them being under a greater influence of positive selection than the other. I have run site models which prove this. Now, I want to present a phylogenetic tree specifying separate omega values for every branch that form these two clades and find out which are the species that are contributing the most to the over all omega value of the lineage.
To achieve this, I used the branch-site model (model=2 NSsite=0). But it prints out a tree where all the foreground and background branches have a single omega. So, what I want to know is ...
- Is there any model that prints a tree with separate omega value for each branch in the tree (like the GA Branch test on the datamonkey server, which I can not use as it has a limitation of 20 sequences)?
- Is this done using using the branch site model A (model=2 and NSsite=2) or the clade models (C and D)? I did try the branch site model A but it only prints out the probabilities for each amino-acid site to be under a site class and not the individual dn/ds ratio for each branch I am interested in.
I would really appreciate your help with this.
Thank you,
Kartik
model
paml
•
written
by
Kartik Sunagar
