
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
00:59 02/17/2011

Input data from file "183-phylip/chr7_11467.phylip" (29 taxa, 231 characters)

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


Note: number of parameters compared to sample size
is low (Nchar/Nparams < 40). Consider the use of AICc.


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

Model	df	lnL		AICc		wAICc	cwAICc
HKYG	60	-590.7983	1344.6554	0.3144	0.3144
HKY	59	-592.6513	1344.7062	0.3066	0.6210
HKYI	60	-590.8519	1344.7627	0.2980	0.9190
HKYIG	61	-590.7963	1348.3501	0.0496	0.9686
GTR	63	-588.4530	1351.1935	0.0120	0.9806
GTRG	64	-586.7600	1351.6405	0.0096	0.9901
GTRI	64	-586.8793	1351.8790	0.0085	0.9986
GTRIG	65	-586.7605	1355.5209	0.0014	1.0000
K2P	56	-614.0880	1376.8657	0.0000	1.0000
K2PI	57	-612.8158	1377.8512	0.0000	1.0000
K2PG	57	-612.8320	1377.8837	0.0000	1.0000
K2PIG	58	-612.8328	1381.4563	0.0000	1.0000
F81I	59	-612.5385	1384.4806	0.0000	1.0000
F81G	59	-612.5748	1384.5531	0.0000	1.0000
F81	58	-614.7138	1385.2184	0.0000	1.0000
SYMG	61	-609.5256	1385.8085	0.0000	1.0000
SYM	60	-611.8289	1386.7167	0.0000	1.0000
F81IG	60	-612.5763	1388.2115	0.0000	1.0000
SYMI	61	-610.7642	1388.2858	0.0000	1.0000
SYMIG	62	-610.7961	1392.0921	0.0000	1.0000
JC69G	56	-629.8197	1408.3291	0.0000	1.0000
JC69	55	-632.3552	1409.9104	0.0000	1.0000
JC69I	56	-631.1453	1410.9803	0.0000	1.0000
JC69IG	57	-631.2014	1414.6224	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-586.7600	1301.5200	0.1854	0.1854
HKYG	60	-590.7983	1301.5966	0.1784	0.3638
HKYI	60	-590.8519	1301.7039	0.1691	0.5328
GTRI	64	-586.8793	1301.7585	0.1645	0.6974
GTR	63	-588.4530	1302.9061	0.0927	0.7901
HKY	59	-592.6513	1303.3027	0.0760	0.8661
GTRIG	65	-586.7605	1303.5209	0.0682	0.9342
HKYIG	61	-590.7963	1303.5927	0.0658	1.0000
K2PI	57	-612.8158	1339.6316	0.0000	1.0000
K2PG	57	-612.8320	1339.6640	0.0000	1.0000
K2P	56	-614.0880	1340.1761	0.0000	1.0000
SYMG	61	-609.5256	1341.0511	0.0000	1.0000
K2PIG	58	-612.8328	1341.6656	0.0000	1.0000
F81I	59	-612.5385	1343.0771	0.0000	1.0000
F81G	59	-612.5748	1343.1496	0.0000	1.0000
SYMI	61	-610.7642	1343.5284	0.0000	1.0000
SYM	60	-611.8289	1343.6579	0.0000	1.0000
F81IG	60	-612.5763	1345.1527	0.0000	1.0000
F81	58	-614.7138	1345.4277	0.0000	1.0000
SYMIG	62	-610.7961	1345.5921	0.0000	1.0000
JC69G	56	-629.8197	1371.6395	0.0000	1.0000
JC69I	56	-631.1453	1374.2906	0.0000	1.0000
JC69	55	-632.3552	1374.7104	0.0000	1.0000
JC69IG	57	-631.2014	1376.4028	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKY	59	-592.6513	1506.4053	0.5413	0.5413
HKYG	60	-590.7983	1508.1417	0.2272	0.7685
HKYI	60	-590.8519	1508.2490	0.2153	0.9839
HKYIG	61	-590.7963	1513.5801	0.0150	0.9989
GTR	63	-588.4530	1519.7784	0.0007	0.9995
GTRG	64	-586.7600	1521.8347	0.0002	0.9998
GTRI	64	-586.8793	1522.0733	0.0002	1.0000
GTRIG	65	-586.7605	1527.2781	0.0000	1.0000
K2P	56	-614.0880	1532.9515	0.0000	1.0000
K2PI	57	-612.8158	1535.8494	0.0000	1.0000
K2PG	57	-612.8320	1535.8818	0.0000	1.0000
K2PIG	58	-612.8328	1541.3258	0.0000	1.0000
F81	58	-614.7138	1545.0879	0.0000	1.0000
F81I	59	-612.5385	1546.1797	0.0000	1.0000
F81G	59	-612.5748	1546.2522	0.0000	1.0000
SYM	60	-611.8289	1550.2029	0.0000	1.0000
SYMG	61	-609.5256	1551.0386	0.0000	1.0000
F81IG	60	-612.5763	1551.6977	0.0000	1.0000
SYMI	61	-610.7642	1553.5159	0.0000	1.0000
SYMIG	62	-610.7961	1559.0220	0.0000	1.0000
JC69	55	-632.3552	1564.0434	0.0000	1.0000
JC69G	56	-629.8197	1564.4149	0.0000	1.0000
JC69I	56	-631.1453	1567.0660	0.0000	1.0000
JC69IG	57	-631.2014	1572.6206	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 183-phylip/chr7_11467.phylip = 1 - 231;
 Partition Dummy = 1:183-phylip/chr7_11467.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 (HKYG)]

BEGIN MRBAYES;
 Charset 183-phylip/chr7_11467.phylip = 1 - 231;
 Partition Dummy = 1:183-phylip/chr7_11467.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 (HKY)]

BEGIN MRBAYES;
 Charset 183-phylip/chr7_11467.phylip = 1 - 231;
 Partition Dummy = 1:183-phylip/chr7_11467.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0);
END;


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

