2024年4月30日发(作者:)

核苷酸替换模型参数NST修正:

JC nst=1 base=equal rates=equal pinv=0

JC+I nst=1 base=equal rates=equal pinv=est

JC+G nst=1 base=equal rates=gamma shape=est pinv=0

JC+I+G nst=1 base=equal rates=gamma shape=est pinv=est

F81 nst=1 base=est rates=equal pinv=0

F81+I nst=1 base=est rates=equal pinv=est

F81+G nst=1 base=est rates=gamma shape=est pinv=0

F81+I+G nst=1 base=est rates=invgamma shape=est pinv=est

K80 nst=2 base=equal tratio=est rates=equal pinv=0

K80+I nst=2 base=equal tratio=est rates=equal pin=est

K80+G nst=2 base=equal tratio=est rates=gamma shape=est pinv=0

K80+I+G nst=2 base=equal tratio=est rates=invgamma shape=est pinv=est

HKY nst=2 base=est tratio=est rates=equal pinv=0

HKY+I nst=2 base=est tratio=est rates=equal pinv=est

HKY+G nst=2 base=est tratio=est rates=gamma shape=est pinv=0

HKY+I+G nst=2 base=est tratio=est rates=invgamma shape=est pinv=est

TrNef nst=6 base=equal rmat=est rclass=(a b a a e a) rates=equal pinv=0

TrNef+I nst=6 base=equal rmat=est rclass=(a b a a e a) rates=equal pinv=est

TrNef+G nst=6 base=equal rmat=est rclass=(a b a a e a) rates=gamma shape=est pinv=0

TrNef+I+G nst=6 base=equal rmat=est rclass=(a b a a e a) rates=invgamma shape=est pinv=est

TrN nst=6 base=est rmat=est rclass=(a b a a e a) rates=equal pinv=0

TrN+I nst=6 base=est rmat=est rclass=(a b a a e a) rates=equal pinv=est

TrN+G nst=6 base=est rmat=est rclass=(a b a a e a) rates=gamma shape=est pinv=0

TrN+I+G nst=6 base=est rmat=est rclass=(a b a a e a) rates=invgamma shape=est pinv=est

K81 nst=6 base=equal rmat=est rclass=(a b c c b a) rates=equal pinv=0

K81+I nst=6 base=equal rmat=est rclass=(a b c c b a) rates=equal pinv=est

K81+G nst=6 base=equal rmat=est rclass=(a b c c b a) rates=gamma shape=est pinv=0

K81+I+G nst=6 base=equal rmat=est rclass=(a b c c b a) rates=invgamma shape=est pinv=est

K81uf nst=6 base=est rmat=est rclass=(a b c c b a) rates=equal pinv=0

K81uf+I nst=6 base=est rmat=est rclass=(a b c c b a) rates=equal pinv=est

K81uf+G nst=6 base=est rmat=est rclass=(a b c c b a) rates=gamma shape=est pinv=0

K81uf+I+G nst=6 base=est rmat=est rclass=(a b c c b a) rates=invgamma shape=est pinv=est

TIMef nst=6 base=equal rmat=est rclass=(a b c c e a) rates=equal pinv=0

TIMef+I nst=6 base=equal rmat=est rclass=(a b c c e a) rates=equal pinv=est

TIMef+G nst=6 base=equal rmat=est rclass=(a b c c e a) rates=gamma shape=est pinv=0

TIMef+I+G nst=6 base=equal rmat=est rclass=(a b c c e a) rates=invgamma shape=est pinv=est

TIM nst=6 base=est rmat=est rclass=(a b c c e a) rates=equal pinv=0

TIM+I nst=6 base=est rmat=est rclass=(a b c c e a) rates=equal pinv=est

TIM+G nst=6 base=est rmat=est rclass=(a b c c e a) rates=gamma shape=est pinv=0

TIM+I+G nst=6 base=est rmatrix=est rclass=(a b c c e a) rates=invgamma shape=est

pinv=est

TVMef nst=6 base=equal rmat=est rclass=(a b c d b e) rates=equal pinv=0

TVMef+I nst=6 base=equal rmat=est rclass=(a b c d b e) rates=equal pinv=est

TVMef+G nst=6 base=equal rmat=est rclass=(a b c d b e) rates=gamma shape=est pinv=0

TVMef+I+G nst=6 base=equal rmat=est rclass=(a b c d b e) rates=invgamma shape=est pinv=est

TVM nst=6 base=est rmat=est rclass=(a b c d b e) rates=equal pinv=0

TVM+I nst=6 base=est rmat=est rclass=(a b c d b e) rates=equal pinv=est

TVM+G nst=6 base=est rmat=est rclass=(a b c d b e) rates=gamma shape=est pinv=0

TVM+I+G nst=6 base=est rmat=est rclass=(a b c d b e) rates=invgamma shape=est pinv=est

SYM nst=6 base=equal rmat=est rclass= (a b c d e f) rates=equal pinv=0

SYM+I nst=6 base=equal rmat=est rclass= (a b c d e f) rates=equal pinv=est

SYM+G nst=6 base=equal rmat=est rclass= (a b c d e f) rates=gamma shape=est pinv=0

SYM+I+G nst=6 base=equal rmat=est rclass= (a b c d e f) rates=invgamma shape=est pinv=est

GTR nst=6 base=est rmat=est rclass= (a b c d e f) rates=equal pinv=0

GTR+I nst=6 base=est rmat=est rclass= (a b c d e f) rates=equal pinv=est

GTR+G nst=6 base=est rmat=est rclass= (a b c d e f) rates=gamma shape=est pinv=0

GTR+I+G nst=6 base=est rmat=est rclass= (a b c d e f) rates=invgamma shape=est pinv=est

注意:其中参数base和rmat是paup模块里的,mrbayes里不需要设置,rclas也不是mrbayes里的

参数。

在mrbayes里的GTR模型示例如下:

lset nst=6 rates=equal ;

Prset statefreqpr=dirichlet(1,1,1,1) pinv=fix(0);