
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
21:33 04/04/2011

Input data from file "591-tomodel-phylip/chr13_2870.phylip" (7 taxa, 970 characters)

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: HKYG


(Output sorted by AIC. w: Akaike weight, cw: cumulative w.)

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-3391.3718	6822.7436	0.5238	0.5238
GTRIG	21	-3391.2783	6824.5566	0.2116	0.7353
HKYG	16	-3396.9491	6825.8983	0.1082	0.8435
GTRI	20	-3393.3519	6826.7039	0.0723	0.9158
HKYIG	17	-3396.7096	6827.4192	0.0506	0.9664
HKYI	16	-3398.1214	6828.2427	0.0335	0.9999
K2PG	13	-3407.6402	6841.2805	0.0000	0.9999
K2PIG	14	-3407.3094	6842.6188	0.0000	0.9999
K2PI	13	-3408.4402	6842.8803	0.0000	1.0000
SYMG	17	-3404.7432	6843.4864	0.0000	1.0000
SYMIG	18	-3404.4827	6844.9654	0.0000	1.0000
SYMI	17	-3405.7189	6845.4378	0.0000	1.0000
GTR	19	-3437.8414	6913.6828	0.0000	1.0000
HKY	15	-3444.1358	6918.2716	0.0000	1.0000
SYM	16	-3452.5846	6937.1691	0.0000	1.0000
K2P	12	-3456.6423	6937.2845	0.0000	1.0000
F81G	15	-3492.8490	7015.6980	0.0000	1.0000
F81I	15	-3493.3397	7016.6795	0.0000	1.0000
F81IG	16	-3492.6526	7017.3052	0.0000	1.0000
JC69G	12	-3501.5981	7027.1961	0.0000	1.0000
JC69I	12	-3502.2234	7028.4468	0.0000	1.0000
JC69IG	13	-3501.3567	7028.7134	0.0000	1.0000
F81	14	-3536.6224	7101.2449	0.0000	1.0000
JC69	11	-3546.4692	7114.9384	0.0000	1.0000



(Output sorted by BIC. w: Akaike weight, cw: cumulative w.)

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-3396.9491	6903.9350	0.4206	0.4206
K2PG	13	-3407.6402	6904.6853	0.2890	0.7096
HKYI	16	-3398.1214	6906.2795	0.1303	0.8399
K2PI	13	-3408.4402	6906.2852	0.1299	0.9698
HKYIG	17	-3396.7096	6910.3332	0.0172	0.9869
K2PIG	14	-3407.3094	6910.9009	0.0129	0.9999
GTRG	20	-3391.3718	6920.2895	0.0001	1.0000
GTRI	20	-3393.3519	6924.2498	0.0000	1.0000
SYMG	17	-3404.7432	6926.4004	0.0000	1.0000
GTRIG	21	-3391.2783	6926.9799	0.0000	1.0000
SYMI	17	-3405.7189	6928.3519	0.0000	1.0000
SYMIG	18	-3404.4827	6932.7567	0.0000	1.0000
HKY	15	-3444.1358	6991.4311	0.0000	1.0000
K2P	12	-3456.6423	6995.8121	0.0000	1.0000
GTR	19	-3437.8414	7006.3514	0.0000	1.0000
SYM	16	-3452.5846	7015.2059	0.0000	1.0000
JC69G	12	-3501.5981	7085.7237	0.0000	1.0000
JC69I	12	-3502.2234	7086.9743	0.0000	1.0000
F81G	15	-3492.8490	7088.8574	0.0000	1.0000
F81I	15	-3493.3397	7089.8389	0.0000	1.0000
JC69IG	13	-3501.3567	7092.1183	0.0000	1.0000
F81IG	16	-3492.6526	7095.3420	0.0000	1.0000
JC69	11	-3546.4692	7168.5886	0.0000	1.0000
F81	14	-3536.6224	7169.5270	0.0000	1.0000



To specify default values for best AIC, AICc, or BIC models in MrBayes:

[Mrbayes block for the best AIC model (GTRG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr13_2870.phylip = 1 - 970;
 Partition Dummy = 1:591-tomodel-phylip/chr13_2870.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 Prset applyto=(1) revmatpr=Dirichlet(1.0,1.0,1.0,1.0,1.0,1.0) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


[Mrbayes block for the best AICc model (GTRG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr13_2870.phylip = 1 - 970;
 Partition Dummy = 1:591-tomodel-phylip/chr13_2870.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 Prset applyto=(1) revmatpr=Dirichlet(1.0,1.0,1.0,1.0,1.0,1.0) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


[Mrbayes block for the best BIC model (HKYG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr13_2870.phylip = 1 - 970;
 Partition Dummy = 1:591-tomodel-phylip/chr13_2870.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


-------------------------------------------------------------
End of Output

