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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: K2PG


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
GTRG	64	-640.1134	1463.6935	0.5601	0.5601
GTRI	64	-641.0040	1465.4746	0.2299	0.7899
GTRIG	65	-639.8028	1467.1895	0.0975	0.8874
HKYG	60	-650.5298	1468.5920	0.0484	0.9358
HKYI	60	-651.1475	1469.8274	0.0261	0.9619
K2PG	57	-657.7179	1471.5504	0.0110	0.9729
HKYIG	61	-650.0726	1471.5831	0.0108	0.9838
SYMG	61	-650.5287	1472.4952	0.0069	0.9906
SYMI	61	-650.9796	1473.3972	0.0044	0.9950
K2PIG	58	-657.1198	1474.1115	0.0031	0.9981
SYMIG	62	-650.1502	1475.6950	0.0014	0.9994
K2PI	57	-660.7156	1477.5459	0.0005	1.0000
GTR	63	-655.2319	1489.8677	0.0000	1.0000
HKY	59	-665.5539	1494.7852	0.0000	1.0000
SYM	60	-666.5929	1500.7183	0.0000	1.0000
K2P	56	-674.1713	1500.7476	0.0000	1.0000
F81G	59	-670.0617	1503.8009	0.0000	1.0000
F81I	59	-670.4293	1504.5360	0.0000	1.0000
JC69G	56	-676.3112	1505.0275	0.0000	1.0000
F81IG	60	-669.3910	1506.3146	0.0000	1.0000
JC69I	56	-676.9606	1506.3262	0.0000	1.0000
JC69IG	57	-675.5275	1507.1695	0.0000	1.0000
F81	58	-685.8504	1531.5725	0.0000	1.0000
JC69	55	-692.7273	1534.1967	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-640.1134	1408.2269	0.5216	0.5216
GTRIG	65	-639.8028	1409.6056	0.2618	0.7834
GTRI	64	-641.0040	1410.0079	0.2141	0.9975
HKYG	60	-650.5298	1421.0595	0.0009	0.9983
HKYIG	61	-650.0726	1422.1452	0.0005	0.9988
HKYI	60	-651.1475	1422.2950	0.0005	0.9993
SYMG	61	-650.5287	1423.0573	0.0003	0.9996
SYMI	61	-650.9796	1423.9593	0.0002	0.9998
SYMIG	62	-650.1502	1424.3003	0.0002	1.0000
K2PG	57	-657.7179	1429.4357	0.0000	1.0000
K2PIG	58	-657.1198	1430.2397	0.0000	1.0000
K2PI	57	-660.7156	1435.4312	0.0000	1.0000
GTR	63	-655.2319	1436.4637	0.0000	1.0000
HKY	59	-665.5539	1449.1078	0.0000	1.0000
SYM	60	-666.5929	1453.1859	0.0000	1.0000
F81G	59	-670.0617	1458.1235	0.0000	1.0000
F81IG	60	-669.3910	1458.7821	0.0000	1.0000
F81I	59	-670.4293	1458.8586	0.0000	1.0000
K2P	56	-674.1713	1460.3425	0.0000	1.0000
JC69G	56	-676.3112	1464.6224	0.0000	1.0000
JC69IG	57	-675.5275	1465.0549	0.0000	1.0000
JC69I	56	-676.9606	1465.9211	0.0000	1.0000
F81	58	-685.8504	1487.7007	0.0000	1.0000
JC69	55	-692.7273	1495.4546	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-657.7179	1621.5621	0.4217	0.4217
HKYG	60	-650.5298	1623.2978	0.1771	0.5988
GTRG	64	-640.1134	1623.9477	0.1279	0.7267
HKYI	60	-651.1475	1624.5332	0.0955	0.8222
GTRI	64	-641.0040	1625.7288	0.0525	0.8747
K2PIG	58	-657.1198	1625.7367	0.0523	0.9270
K2PI	57	-660.7156	1627.5576	0.0210	0.9480
HKYIG	61	-650.0726	1627.7541	0.0191	0.9671
SYMG	61	-650.5287	1628.6662	0.0121	0.9792
GTRIG	65	-639.8028	1628.6971	0.0119	0.9911
SYMI	61	-650.9796	1629.5682	0.0077	0.9988
SYMIG	62	-650.1502	1633.2799	0.0012	1.0000
HKY	59	-665.5539	1647.9754	0.0000	1.0000
GTR	63	-655.2319	1648.8139	0.0000	1.0000
K2P	56	-674.1713	1649.0982	0.0000	1.0000
JC69G	56	-676.3112	1653.3782	0.0000	1.0000
JC69I	56	-676.9606	1654.6768	0.0000	1.0000
SYM	60	-666.5929	1655.4241	0.0000	1.0000
F81G	59	-670.0617	1656.9911	0.0000	1.0000
JC69IG	57	-675.5275	1657.1813	0.0000	1.0000
F81I	59	-670.4293	1657.7262	0.0000	1.0000
F81IG	60	-669.3910	1661.0204	0.0000	1.0000
JC69	55	-692.7273	1680.8397	0.0000	1.0000
F81	58	-685.8504	1683.1977	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/chr20_5193.phylip = 1 - 215;
 Partition Dummy = 1:183-phylip/chr20_5193.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 183-phylip/chr20_5193.phylip = 1 - 215;
 Partition Dummy = 1:183-phylip/chr20_5193.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 (K2PG)]

BEGIN MRBAYES;
 Charset 183-phylip/chr20_5193.phylip = 1 - 215;
 Partition Dummy = 1:183-phylip/chr20_5193.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

