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

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

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


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	16	-797.9958	1629.7185	0.4323	0.4323
HKYIG	17	-797.3616	1630.6722	0.2683	0.7006
HKYI	16	-798.8289	1631.3847	0.1879	0.8885
GTRG	20	-795.5508	1633.8027	0.0561	0.9446
GTRIG	21	-794.9137	1634.8080	0.0339	0.9785
GTRI	20	-796.5098	1635.7206	0.0215	1.0000
HKY	15	-809.2799	1650.0787	0.0000	1.0000
GTR	19	-806.5231	1653.4820	0.0000	1.0000
K2PG	13	-819.5507	1666.2460	0.0000	1.0000
K2PIG	14	-818.8995	1667.1240	0.0000	1.0000
K2PI	13	-820.0576	1667.2598	0.0000	1.0000
SYMG	17	-817.3280	1670.6050	0.0000	1.0000
SYMIG	18	-816.7300	1671.6453	0.0000	1.0000
SYMI	17	-817.8833	1671.7157	0.0000	1.0000
F81G	15	-820.2850	1672.0889	0.0000	1.0000
F81IG	16	-819.6880	1673.1030	0.0000	1.0000
F81I	15	-820.9891	1673.4971	0.0000	1.0000
K2P	12	-832.0139	1689.0058	0.0000	1.0000
F81	14	-831.7449	1692.8146	0.0000	1.0000
SYM	16	-830.0674	1693.8617	0.0000	1.0000
JC69G	12	-840.2857	1705.5495	0.0000	1.0000
JC69IG	13	-839.7001	1706.5449	0.0000	1.0000
JC69I	12	-840.8793	1706.7367	0.0000	1.0000
JC69	11	-852.3977	1727.6205	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-797.9958	1627.9915	0.3895	0.3895
HKYIG	17	-797.3616	1628.7232	0.2702	0.6597
HKYI	16	-798.8289	1629.6578	0.1693	0.8290
GTRG	20	-795.5508	1631.1017	0.0823	0.9112
GTRIG	21	-794.9137	1631.8274	0.0572	0.9685
GTRI	20	-796.5098	1633.0197	0.0315	1.0000
HKY	15	-809.2799	1648.5597	0.0000	1.0000
GTR	19	-806.5231	1651.0461	0.0000	1.0000
K2PG	13	-819.5507	1665.1013	0.0000	1.0000
K2PIG	14	-818.8995	1665.7991	0.0000	1.0000
K2PI	13	-820.0576	1666.1151	0.0000	1.0000
SYMG	17	-817.3280	1668.6559	0.0000	1.0000
SYMIG	18	-816.7300	1669.4600	0.0000	1.0000
SYMI	17	-817.8833	1669.7667	0.0000	1.0000
F81G	15	-820.2850	1670.5700	0.0000	1.0000
F81IG	16	-819.6880	1671.3760	0.0000	1.0000
F81I	15	-820.9891	1671.9781	0.0000	1.0000
K2P	12	-832.0139	1688.0277	0.0000	1.0000
F81	14	-831.7449	1691.4897	0.0000	1.0000
SYM	16	-830.0674	1692.1347	0.0000	1.0000
JC69G	12	-840.2857	1704.5714	0.0000	1.0000
JC69IG	13	-839.7001	1705.4002	0.0000	1.0000
JC69I	12	-840.8793	1705.7586	0.0000	1.0000
JC69	11	-852.3977	1726.7955	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-797.9958	1688.8737	0.6500	0.6500
HKYI	16	-798.8289	1690.5399	0.2825	0.9325
HKYIG	17	-797.3616	1693.4105	0.0673	0.9997
HKY	15	-809.2799	1705.6367	0.0001	0.9999
GTRG	20	-795.5508	1707.2044	0.0001	1.0000
GTRI	20	-796.5098	1709.1224	0.0000	1.0000
GTRIG	21	-794.9137	1711.7352	0.0000	1.0000
K2PG	13	-819.5507	1714.5681	0.0000	1.0000
K2PI	13	-820.0576	1715.5819	0.0000	1.0000
K2PIG	14	-818.8995	1719.0710	0.0000	1.0000
GTR	19	-806.5231	1723.3437	0.0000	1.0000
F81G	15	-820.2850	1727.6470	0.0000	1.0000
F81I	15	-820.9891	1729.0551	0.0000	1.0000
F81IG	16	-819.6880	1732.2582	0.0000	1.0000
SYMG	17	-817.3280	1733.3432	0.0000	1.0000
K2P	12	-832.0139	1733.6893	0.0000	1.0000
SYMI	17	-817.8833	1734.4540	0.0000	1.0000
SYMIG	18	-816.7300	1737.9524	0.0000	1.0000
F81	14	-831.7449	1744.7616	0.0000	1.0000
JC69G	12	-840.2857	1750.2331	0.0000	1.0000
JC69I	12	-840.8793	1751.4203	0.0000	1.0000
SYM	16	-830.0674	1753.0169	0.0000	1.0000
JC69IG	13	-839.7001	1754.8670	0.0000	1.0000
JC69	11	-852.3977	1768.6519	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/chr3_226.phylip = 1 - 332;
 Partition Dummy = 1:591-tomodel-phylip/chr3_226.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/chr3_226.phylip = 1 - 332;
 Partition Dummy = 1:591-tomodel-phylip/chr3_226.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/chr3_226.phylip = 1 - 332;
 Partition Dummy = 1:591-tomodel-phylip/chr3_226.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

