SCRIPT FOR 28S LONGIDORUS TABERNENSIS SP. NOV.  CAI ET AL. 2020 JZS
-------------------------------------------------------------------
begin mrbayes;
lset nst=6 rates=gamma;
prset  statefreqpr=fixed(0.2305,0.2274,0.2821,0.2600);
prset revmatpr=fixed(0.9002,2.6499,1.5039,0.4583,4.9910,1.0000);
prset shapepr=fixed(0.6870);
prset Pinvarpr=fixed(0.2580);
end;
mcmc ngen=2000000 printfreq=100 samplefreq=100 nchains=4 savebrlens=yes filename=BI;
sumt filename=BI burnin=2001 contype=halfcompat;
sump filename=BI burnin=2001;
end;


SCRIPT FOR 18S LONGIDORUS TABERNENSIS SP. NOV.  CAI ET AL. 2020 JZS
-------------------------------------------------------------------
begin mrbayes;
lset nst=6 rates=gamma;
prset  statefreqpr=fixed(0.2701,0.2176,0.2623,0.2500);
prset revmatpr=fixed(1.1859,3.3754,1.8179,0.5470,6.2308,1.0000);
prset shapepr=fixed(0.6050);
prset Pinvarpr=fixed(0.7450);
end;
mcmc ngen=2000000 printfreq=100 samplefreq=100 nchains=4 savebrlens=yes filename=BI;
sumt filename=BI burnin=2001 contype=halfcompat;
sump filename=BI burnin=2001;
end;

SCRIPT FOR COXI LONGIDORUS TABERNENSIS SP. NOV.  CAI ET AL. 2020 JZS
---------------------------------------------------------------------
begin mrbayes;
lset nst=6 rates=gamma;
prset  statefreqpr=fixed(0.2942,0.1777,0.1758,0.3523);
prset revmatpr=fixed(0.7189,7.7638,0.7189,1.0000,18.2723,1.0000);
prset shapepr=fixed(0.4080);
prset Pinvarpr=fixed(0.3900);
end;
mcmc ngen=2000000 printfreq=100 samplefreq=100 nchains=4 savebrlens=yes filename=BI;
sumt filename=BI burnin=2001 contype=halfcompat;
sump filename=BI burnin=2001;
end;
-------------------------------------------------------------------------------------------------------------
-------------------------------------------------------------------------------------------------------------

SCRIPT FOR 28S XIPHINEMA SUBBAETENSE SP. NOV.  CAI ET AL. 2020 JZS
-------------------------------------------------------------------
begin mrbayes;
lset nst=6 rates=gamma;
prset  statefreqpr=fixed(0.2279,0.2329,0.3060,0.2333);
prset revmatpr=fixed(1.3863,3.8696,3.3174,0.6756,5.2302,1.0000);
prset shapepr=fixed(0.6750);
prset Pinvarpr=fixed(0.3870);
end;
mcmc ngen=2000000 printfreq=100 samplefreq=100 nchains=4 savebrlens=yes filename=BI;
sumt filename=BI burnin=2001 contype=halfcompat;
sump filename=BI burnin=2001;
end;


SCRIPT FOR 18S XIPHINEMA SUBBAETENSE SP. NOV.  CAI ET AL. 2020 JZS
-------------------------------------------------------------------
begin mrbayes;
lset nst=6 rates=gamma;
prset statefreqpr=fixed(0.2744,0.2161,0.2583,0.2513);
prset revmatpr=fixed(0.8996,2.1358,1.2539,0.3543,4.8968,1.0000);
prset shapepr=fixed(0.6900);
prset Pinvarpr=fixed(0.7590);

mcmc ngen=2000000 printfreq=100 samplefreq=100 nchains=4 savebrlens=yes filename=BI;
sumt filename=BI burnin=2001 contype=halfcompat;
sump filename=BI burnin=2001;
end;


SCRIPT FOR COXI XIPHINEMA SUBBAETENSE SP. NOV.  CAI ET AL. 2020 JZS
-------------------------------------------------------------------
begin mrbayes;
lset nst=6 rates=gamma;
prset  statefreqpr=fixed(0.3639,0.1336,0.1385,0.3641);
prset revmatpr=fixed(1.4799,13.5288,1.2934,4.8253,27.8124,1.0000);
prset shapepr=fixed(0.4200);
prset Pinvarpr=fixed(0.2760);
end;
mcmc ngen=2000000 printfreq=100 samplefreq=100 nchains=4 savebrlens=yes filename=BI;
sumt filename=BI burnin=2001 contype=halfcompat;
sump filename=BI burnin=2001;
end;