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

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

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


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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-1670.9488	3373.8976	0.3517	0.3517
HKYI	16	-1670.9497	3373.8995	0.3514	0.7031
HKYIG	17	-1670.9545	3375.9090	0.1286	0.8317
GTRG	20	-1668.5495	3377.0991	0.0710	0.9027
GTRI	20	-1668.5699	3377.1398	0.0695	0.9722
GTRIG	21	-1668.4936	3378.9872	0.0276	0.9998
HKY	15	-1679.7203	3389.4407	0.0001	1.0000
GTR	19	-1677.0257	3392.0513	0.0000	1.0000
F81G	15	-1693.4328	3416.8655	0.0000	1.0000
F81I	15	-1693.4331	3416.8662	0.0000	1.0000
F81IG	16	-1693.4407	3418.8814	0.0000	1.0000
F81	14	-1702.2292	3432.4583	0.0000	1.0000
K2PI	13	-1733.7847	3493.5695	0.0000	1.0000
K2PG	13	-1733.7852	3493.5703	0.0000	1.0000
K2PIG	14	-1733.7919	3495.5839	0.0000	1.0000
SYMG	17	-1731.4357	3496.8715	0.0000	1.0000
SYMI	17	-1731.4371	3496.8742	0.0000	1.0000
SYMIG	18	-1731.4277	3498.8554	0.0000	1.0000
K2P	12	-1742.6684	3509.3368	0.0000	1.0000
SYM	16	-1740.1674	3512.3349	0.0000	1.0000
JC69I	12	-1752.3935	3528.7870	0.0000	1.0000
JC69G	12	-1752.3940	3528.7880	0.0000	1.0000
JC69IG	13	-1752.4021	3530.8043	0.0000	1.0000
JC69	11	-1761.2744	3544.5487	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1670.9488	3450.4127	0.4909	0.4909
HKYI	16	-1670.9497	3450.4146	0.4904	0.9813
HKYIG	17	-1670.9545	3457.2063	0.0164	0.9977
HKY	15	-1679.7203	3461.1736	0.0023	1.0000
GTRG	20	-1668.5495	3472.7429	0.0000	1.0000
GTRI	20	-1668.5699	3472.7836	0.0000	1.0000
GTRIG	21	-1668.4936	3479.4132	0.0000	1.0000
GTR	19	-1677.0257	3482.9130	0.0000	1.0000
F81G	15	-1693.4328	3488.5984	0.0000	1.0000
F81I	15	-1693.4331	3488.5991	0.0000	1.0000
F81IG	16	-1693.4407	3495.3965	0.0000	1.0000
F81	14	-1702.2292	3499.4090	0.0000	1.0000
K2PI	13	-1733.7847	3555.7380	0.0000	1.0000
K2PG	13	-1733.7852	3555.7388	0.0000	1.0000
K2PIG	14	-1733.7919	3562.5345	0.0000	1.0000
K2P	12	-1742.6684	3566.7231	0.0000	1.0000
SYMG	17	-1731.4357	3578.1687	0.0000	1.0000
SYMI	17	-1731.4371	3578.1714	0.0000	1.0000
SYMIG	18	-1731.4277	3584.9349	0.0000	1.0000
JC69I	12	-1752.3935	3586.1733	0.0000	1.0000
JC69G	12	-1752.3940	3586.1743	0.0000	1.0000
SYM	16	-1740.1674	3588.8499	0.0000	1.0000
JC69IG	13	-1752.4021	3592.9728	0.0000	1.0000
JC69	11	-1761.2744	3597.1529	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr13_727.phylip = 1 - 882;
 Partition Dummy = 1:591-tomodel-phylip/chr13_727.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;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr13_727.phylip = 1 - 882;
 Partition Dummy = 1:591-tomodel-phylip/chr13_727.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;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr13_727.phylip = 1 - 882;
 Partition Dummy = 1:591-tomodel-phylip/chr13_727.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

