Running LG+I+G in phylobayes

Hello,

Prottest reveals LG+I+G (shocker) to be the best model for my amino acid alignment.

How do I specify this in phylobayes? Is this the correct way?

mpirun -np 8 pb_mpi -d gappyout.phy -cat -lg -dgam 4 chain1

From what I understand this is the +G part, but how do we specify +I?

Adrian

phylogenetics phylobayes alignment

• 2,016 views • link

updated 3.9 years ago by Ram

• written 10.9 years ago by Adrian Pelin

0 answers

No answers yet.

More posts like this

PAML/codeml Poxvirus dN/dS
written 2.9 years ago by Adrian Pelin

I have orthologous clusters of viral genes between related species. I am working with Vaccinia virus and want to compare Vaccinia genes from the Copenhagen …
Bar and line graphs using Bokeh Module in python
written 4.3 years ago by anasjamshed

I have a gene sequence. I translate it and calculate its amino acid properties. Now I want to take x values as amino acid residues …
Bi-gram feature extraction from a protein sequence.
written 8.7 years ago by arijeetsat •

**Can someone explain me , how to extract bi-gram feature distribution from a given protein sequence? For eg: Consider the amino acid sequence as MNNPQMNPQRS. …
ProtTest and PhyML
written 11.2 years ago by Michel Edwar

I am drawing a phylogenetic tree. To select the most appropriate model, I ran ProtTest with the aligned amino acid sequences. The ProtTest suggested that …
MrBayes, how to judge the quality of a run?
written 12.3 years ago by Adrian Pelin

Hello, So I got the development version of MrBayes to try the now available LG model. I also have the analysis on the same analysis …
MrBayes: Implementing the LG model
written 12.3 years ago by Adrian Pelin

Hello, When I run prottest on my amino acid alignment I am told LG+G to be the best model, which is common among amino acid …
How To Set Tim1 Model In Mrbayes
written 12.5 years ago by Adrian Pelin

Hello, jModelTest suggested using the TIM1 model for my analysis. I got this as output: <pre><code>Model selected: Model = TIM1+I+G partition = 012230 -lnL = …
Deltaaic Phyml/Mrbayes
written 12.5 years ago by Adrian Pelin

Hello, jModelTest tells me my best model is TIM1+I+G with an AIC of 58750.6, the third best model is the famous GTR+I+G has a deltaAIC …
Mr. Bayes Protein Models OR How to implement prottest output in mrbayes?
written 13.0 years ago by kanwalashraf07 •

Hello, I am using Mr. Bayes for my amino acid data set, The model predicted by prottest is rtRev + G + I (rtRev+G+I). I …
How To Set The Substitution Model In Phyml?
written 14.7 years ago by Dejian

I am drawing two phylogenetic trees. To select the most approapriate model, I ran ProtTest with the aligned amino acid sequences. The ProtTest suggested that …