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

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

Minimum AIC  model: HKYI
Minimum AICc model: HKYI
Minimum BIC  model: HKYI


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
HKYI	16	-667.7190	1369.3740	0.4004	0.4004
HKYG	16	-667.7228	1369.3816	0.3989	0.7994
HKYIG	17	-667.7229	1371.6316	0.1295	0.9289
GTRI	20	-665.7315	1374.4955	0.0309	0.9598
GTRG	20	-665.7998	1374.6320	0.0289	0.9887
GTRIG	21	-665.7996	1376.9470	0.0091	0.9978
HKY	15	-674.0831	1379.8683	0.0021	0.9999
GTR	19	-672.4156	1385.5650	0.0001	1.0000
F81G	15	-682.6776	1397.0572	0.0000	1.0000
F81I	15	-682.7334	1397.1690	0.0000	1.0000
F81IG	16	-682.5845	1399.1050	0.0000	1.0000
K2PI	13	-688.4541	1404.1899	0.0000	1.0000
K2PG	13	-689.4447	1406.1710	0.0000	1.0000
K2PIG	14	-688.4673	1406.4188	0.0000	1.0000
F81	14	-688.8577	1407.1994	0.0000	1.0000
SYMI	17	-686.4985	1409.1826	0.0000	1.0000
SYMG	17	-687.5567	1411.2990	0.0000	1.0000
SYMIG	18	-687.4547	1413.3610	0.0000	1.0000
K2P	12	-696.4501	1417.9950	0.0000	1.0000
SYM	16	-692.9074	1419.7507	0.0000	1.0000
JC69I	12	-702.6058	1430.3063	0.0000	1.0000
JC69G	12	-702.6161	1430.3269	0.0000	1.0000
JC69IG	13	-702.6232	1432.5282	0.0000	1.0000
JC69	11	-708.5840	1440.0910	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	16	-667.7190	1367.4381	0.3742	0.3742
HKYG	16	-667.7228	1367.4456	0.3728	0.7470
HKYIG	17	-667.7229	1369.4458	0.1371	0.8842
GTRI	20	-665.7315	1371.4630	0.0500	0.9342
GTRG	20	-665.7998	1371.5995	0.0467	0.9809
GTRIG	21	-665.7996	1373.5991	0.0172	0.9981
HKY	15	-674.0831	1378.1662	0.0018	0.9998
GTR	19	-672.4156	1382.8311	0.0002	1.0000
F81G	15	-682.6776	1395.3551	0.0000	1.0000
F81I	15	-682.7334	1395.4669	0.0000	1.0000
F81IG	16	-682.5845	1397.1691	0.0000	1.0000
K2PI	13	-688.4541	1402.9082	0.0000	1.0000
K2PG	13	-689.4447	1404.8893	0.0000	1.0000
K2PIG	14	-688.4673	1404.9347	0.0000	1.0000
F81	14	-688.8577	1405.7153	0.0000	1.0000
SYMI	17	-686.4985	1406.9969	0.0000	1.0000
SYMG	17	-687.5567	1409.1133	0.0000	1.0000
SYMIG	18	-687.4547	1410.9094	0.0000	1.0000
K2P	12	-696.4501	1416.9002	0.0000	1.0000
SYM	16	-692.9074	1417.8148	0.0000	1.0000
JC69I	12	-702.6058	1429.2116	0.0000	1.0000
JC69G	12	-702.6161	1429.2322	0.0000	1.0000
JC69IG	13	-702.6232	1431.2465	0.0000	1.0000
JC69	11	-708.5840	1439.1680	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-667.7190	1426.5916	0.4799	0.4799
HKYG	16	-667.7228	1426.5991	0.4781	0.9580
HKYIG	17	-667.7229	1432.2964	0.0277	0.9856
HKY	15	-674.0831	1433.6226	0.0143	0.9999
GTRI	20	-665.7315	1445.4049	0.0000	1.0000
GTRG	20	-665.7998	1445.5414	0.0000	1.0000
F81G	15	-682.6776	1450.8115	0.0000	1.0000
F81I	15	-682.7334	1450.9233	0.0000	1.0000
K2PI	13	-688.4541	1450.9704	0.0000	1.0000
GTRIG	21	-665.7996	1451.2381	0.0000	1.0000
K2PG	13	-689.4447	1452.9515	0.0000	1.0000
GTR	19	-672.4156	1453.0759	0.0000	1.0000
F81IG	16	-682.5845	1456.3226	0.0000	1.0000
K2PIG	14	-688.4673	1456.6940	0.0000	1.0000
F81	14	-688.8577	1457.4746	0.0000	1.0000
K2P	12	-696.4501	1461.2653	0.0000	1.0000
SYMI	17	-686.4985	1469.8475	0.0000	1.0000
SYMG	17	-687.5567	1471.9639	0.0000	1.0000
JC69I	12	-702.6058	1473.5767	0.0000	1.0000
JC69G	12	-702.6161	1473.5973	0.0000	1.0000
SYM	16	-692.9074	1476.9683	0.0000	1.0000
SYMIG	18	-687.4547	1477.4571	0.0000	1.0000
JC69IG	13	-702.6232	1479.3087	0.0000	1.0000
JC69	11	-708.5840	1479.8360	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr2_25261.phylip = 1 - 298;
 Partition Dummy = 1:591-tomodel-phylip/chr2_25261.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr2_25261.phylip = 1 - 298;
 Partition Dummy = 1:591-tomodel-phylip/chr2_25261.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr2_25261.phylip = 1 - 298;
 Partition Dummy = 1:591-tomodel-phylip/chr2_25261.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

