Skip to content

Instantly share code, notes, and snippets.

@brantfaircloth
Created March 30, 2011 20:55
Show Gist options
  • Save brantfaircloth/895282 to your computer and use it in GitHub Desktop.
Save brantfaircloth/895282 to your computer and use it in GitHub Desktop.
Substitution models in mrbayes

##GTR lset applyto=() nst=6 # GTR lset applyto=() nst=6 rates=propinv # GTR + I lset applyto=() nst=6 rates=gamma # GTR + gamma lset applyto=() nst=6 rates=invgamma # GTR + I + gamma

##SYM

lset applyto=() nst=6                           # SYM
prset applyto=() statefreqpr=fixed(equal)

lset applyto=() nst=6 rates=propinv             # SYM + I
prset applyto=() statefreqpr=fixed(equal)

lset applyto=() nst=6 rates=gamma               # SYM + gamma
prset applyto=() statefreqpr=fixed(equal)

lset applyto=() nst=6 rates=invgamma            # SYM + I + gamma
prset applyto=() statefreqpr=fixed(equal)

##HKY

lset applyto=() nst=2                           # HKY
lset applyto=() nst=2 rates=propinv             # HKY + I
lset applyto=() nst=2 rates=gamma               # HKY + gamma
lset applyto=() nst=2 rates=invgamma            # HKY + I + gamma

##K2P

lset applyto=() nst=2                           # K2P
prset applyto=() statefreqpr=fixed(equal)

lset applyto=() nst=2 rates=propinv             # K2P + I
prset applyto=() statefreqpr=fixed(equal)

lset applyto=() nst=2 rates=gamma               # K2P + gamma
prset applyto=() statefreqpr=fixed(equal)

lset applyto=() nst=2 rates=invgamma            # K2P + I + gamma
prset applyto=() statefreqpr=fixed(equal)

##F81

lset applyto=() nst=1                           # F81
lset applyto=() nst=1 rates=propinv             # F81 + I
lset applyto=() nst=1 rates=gamma               # F81 + gamma
lset applyto=() nst=1 rates=invgamma            # F81 + I + gamma

Jukes Cantor

lset applyto=() nst=1                           # JC
prset applyto=() statefreqpr=fixed(equal)

lset applyto=() nst=1 rates=propinv             # JC + I
prset applyto=() statefreqpr=fixed(equal)

lset applyto=() nst=1 rates=gamma               # JC + gamma
prset applyto=() statefreqpr=fixed(equal)

lset applyto=() nst=1 rates=incgamma            # JC + I + gamma
prset applyto=() statefreqpr=fixed(equal)
@oleon12
Copy link

oleon12 commented May 15, 2015

Excellent !! This is so usefull, thank you !!

@weatherwatcher42
Copy link

Is it possible to make the TIM2, TIM3 and TIM3ef models in MrBayes 3? I keep getting them as my lowest AICc score but can't figure out how to implement them.
Here's a reference list with more info about the TIM models and others:
http://www.molecularevolution.org/molevolfiles/models/submodels_final.pdf

@Shooroosh
Copy link

Thanks

@NwaBiafra-debug
Copy link

How do I set TVMef+G in MrBayes?

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment