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);
发布评论