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

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

Minimum AIC  model: GTRG
Minimum AICc model: K2PG
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
K2PG	57	-639.7815	1443.6538	0.2620	0.2620
HKYG	60	-633.8067	1444.3576	0.1843	0.4463
K2PI	57	-640.3703	1444.8316	0.1454	0.5917
GTRG	64	-625.2059	1444.9719	0.1356	0.7273
HKYI	60	-634.5517	1445.8476	0.0875	0.8148
GTRI	64	-625.9132	1446.3863	0.0668	0.8816
K2PIG	58	-639.4927	1447.2296	0.0438	0.9254
HKYIG	61	-633.4284	1447.9505	0.0306	0.9560
GTRIG	65	-624.7170	1448.6275	0.0218	0.9778
SYMG	61	-634.4169	1449.9276	0.0114	0.9892
SYMI	61	-634.6266	1450.3469	0.0092	0.9984
SYMIG	62	-634.1756	1453.8631	0.0016	1.0000
K2P	56	-655.2666	1470.5332	0.0000	1.0000
HKY	59	-649.6216	1471.7047	0.0000	1.0000
GTR	63	-643.1671	1476.3343	0.0000	1.0000
SYM	60	-650.8753	1478.4947	0.0000	1.0000
JC69G	56	-662.8474	1485.6948	0.0000	1.0000
JC69I	56	-663.5576	1487.1152	0.0000	1.0000
JC69IG	57	-662.3465	1488.7839	0.0000	1.0000
F81G	59	-658.8499	1490.1613	0.0000	1.0000
F81I	59	-659.4995	1491.4604	0.0000	1.0000
F81IG	60	-658.8478	1494.4399	0.0000	1.0000
JC69	55	-678.2842	1512.5385	0.0000	1.0000
F81	58	-674.3292	1516.9027	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-625.2059	1378.4119	0.4720	0.4720
GTRIG	65	-624.7170	1379.4340	0.2831	0.7551
GTRI	64	-625.9132	1379.8263	0.2327	0.9878
HKYG	60	-633.8067	1387.6134	0.0047	0.9925
HKYIG	61	-633.4284	1388.8567	0.0025	0.9951
HKYI	60	-634.5517	1389.1034	0.0023	0.9973
SYMG	61	-634.4169	1390.8338	0.0009	0.9983
SYMI	61	-634.6266	1391.2532	0.0008	0.9991
SYMIG	62	-634.1756	1392.3513	0.0004	0.9995
K2PG	57	-639.7815	1393.5629	0.0002	0.9997
K2PI	57	-640.3703	1394.7407	0.0001	0.9999
K2PIG	58	-639.4927	1394.9853	0.0001	1.0000
GTR	63	-643.1671	1412.3343	0.0000	1.0000
HKY	59	-649.6216	1417.2431	0.0000	1.0000
SYM	60	-650.8753	1421.7505	0.0000	1.0000
K2P	56	-655.2666	1422.5332	0.0000	1.0000
F81G	59	-658.8499	1435.6998	0.0000	1.0000
F81I	59	-659.4995	1436.9989	0.0000	1.0000
JC69G	56	-662.8474	1437.6948	0.0000	1.0000
F81IG	60	-658.8478	1437.6957	0.0000	1.0000
JC69IG	57	-662.3465	1438.6930	0.0000	1.0000
JC69I	56	-663.5576	1439.1152	0.0000	1.0000
F81	58	-674.3292	1464.6585	0.0000	1.0000
JC69	55	-678.2842	1466.5684	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-639.7815	1578.6433	0.5163	0.5163
K2PI	57	-640.3703	1579.8210	0.2865	0.8028
HKYG	60	-633.8067	1582.4349	0.0775	0.8803
K2PIG	58	-639.4927	1583.3127	0.0500	0.9303
HKYI	60	-634.5517	1583.9248	0.0368	0.9672
GTRG	64	-625.2059	1586.2214	0.0117	0.9788
HKYIG	61	-633.4284	1586.9252	0.0082	0.9870
GTRI	64	-625.9132	1587.6359	0.0058	0.9928
SYMG	61	-634.4169	1588.9023	0.0031	0.9959
SYMI	61	-634.6266	1589.3216	0.0025	0.9983
GTRIG	65	-624.7170	1590.4906	0.0014	0.9997
SYMIG	62	-634.1756	1593.6668	0.0003	1.0000
K2P	56	-655.2666	1604.3665	0.0000	1.0000
HKY	59	-649.6216	1608.8176	0.0000	1.0000
SYM	60	-650.8753	1616.5720	0.0000	1.0000
GTR	63	-643.1671	1616.8968	0.0000	1.0000
JC69G	56	-662.8474	1619.5282	0.0000	1.0000
JC69I	56	-663.5576	1620.9485	0.0000	1.0000
JC69IG	57	-662.3465	1623.7734	0.0000	1.0000
F81G	59	-658.8499	1627.2742	0.0000	1.0000
F81I	59	-659.4995	1628.5733	0.0000	1.0000
F81IG	60	-658.8478	1632.5171	0.0000	1.0000
JC69	55	-678.2842	1645.1547	0.0000	1.0000
F81	58	-674.3292	1652.9859	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/chr10_4411.phylip = 1 - 190;
 Partition Dummy = 1:183-phylip/chr10_4411.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 (K2PG)]

BEGIN MRBAYES;
 Charset 183-phylip/chr10_4411.phylip = 1 - 190;
 Partition Dummy = 1:183-phylip/chr10_4411.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;


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

BEGIN MRBAYES;
 Charset 183-phylip/chr10_4411.phylip = 1 - 190;
 Partition Dummy = 1:183-phylip/chr10_4411.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

