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

Input data from file "591-tomodel-phylip/chr1_5365.phylip" (7 taxa, 347 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	-764.2110	1562.0705	0.3371	0.3371
HKYG	16	-764.2909	1562.2302	0.3112	0.6483
GTRI	20	-760.8133	1564.2033	0.1160	0.7644
HKYIG	17	-764.2910	1564.4421	0.1030	0.8674
GTRG	20	-760.9649	1564.5064	0.0997	0.9671
GTRIG	21	-760.9650	1566.7730	0.0321	0.9992
HKY	15	-771.7150	1574.8801	0.0006	0.9997
GTR	19	-768.1089	1576.5420	0.0002	1.0000
K2PI	13	-779.1741	1585.4414	0.0000	1.0000
K2PG	13	-779.1759	1585.4448	0.0000	1.0000
K2PIG	14	-779.1759	1587.6169	0.0000	1.0000
F81I	15	-778.9961	1589.4424	0.0000	1.0000
F81G	15	-779.0509	1589.5519	0.0000	1.0000
SYMI	17	-777.5050	1590.8701	0.0000	1.0000
SYMG	17	-777.5272	1590.9146	0.0000	1.0000
F81IG	16	-779.0510	1591.7504	0.0000	1.0000
SYMIG	18	-777.5273	1593.1399	0.0000	1.0000
K2P	12	-787.6643	1600.2627	0.0000	1.0000
F81	14	-786.6263	1602.5176	0.0000	1.0000
SYM	16	-785.5429	1604.7343	0.0000	1.0000
JC69I	12	-792.8382	1610.6105	0.0000	1.0000
JC69G	12	-792.8543	1610.6427	0.0000	1.0000
JC69IG	13	-792.8544	1612.8018	0.0000	1.0000
JC69	11	-800.9636	1624.7153	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	16	-764.2110	1560.4220	0.2893	0.2893
HKYG	16	-764.2909	1560.5817	0.2671	0.5564
GTRI	20	-760.8133	1561.6266	0.1584	0.7148
GTRG	20	-760.9649	1561.9298	0.1361	0.8509
HKYIG	17	-764.2910	1562.5819	0.0982	0.9492
GTRIG	21	-760.9650	1563.9300	0.0501	0.9993
HKY	15	-771.7150	1573.4299	0.0004	0.9997
GTR	19	-768.1089	1574.2178	0.0003	1.0000
K2PI	13	-779.1741	1584.3483	0.0000	1.0000
K2PG	13	-779.1759	1584.3517	0.0000	1.0000
K2PIG	14	-779.1759	1586.3518	0.0000	1.0000
F81I	15	-778.9961	1587.9922	0.0000	1.0000
F81G	15	-779.0509	1588.1017	0.0000	1.0000
SYMI	17	-777.5050	1589.0099	0.0000	1.0000
SYMG	17	-777.5272	1589.0544	0.0000	1.0000
F81IG	16	-779.0510	1590.1019	0.0000	1.0000
SYMIG	18	-777.5273	1591.0546	0.0000	1.0000
K2P	12	-787.6643	1599.3286	0.0000	1.0000
F81	14	-786.6263	1601.2525	0.0000	1.0000
SYM	16	-785.5429	1603.0858	0.0000	1.0000
JC69I	12	-792.8382	1609.6764	0.0000	1.0000
JC69G	12	-792.8543	1609.7086	0.0000	1.0000
JC69IG	13	-792.8544	1611.7087	0.0000	1.0000
JC69	11	-800.9636	1623.9272	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-764.2110	1622.0112	0.5031	0.5031
HKYG	16	-764.2909	1622.1709	0.4645	0.9675
HKYIG	17	-764.2910	1628.0205	0.0249	0.9925
HKY	15	-771.7150	1631.1698	0.0052	0.9976
K2PI	13	-779.1741	1634.3895	0.0010	0.9987
K2PG	13	-779.1759	1634.3930	0.0010	0.9997
GTRI	20	-760.8133	1638.6131	0.0001	0.9998
GTRG	20	-760.9649	1638.9163	0.0001	0.9999
K2PIG	14	-779.1759	1640.2423	0.0001	1.0000
GTRIG	21	-760.9650	1644.7658	0.0000	1.0000
K2P	12	-787.6643	1645.5205	0.0000	1.0000
F81I	15	-778.9961	1645.7321	0.0000	1.0000
F81G	15	-779.0509	1645.8416	0.0000	1.0000
GTR	19	-768.1089	1647.3550	0.0000	1.0000
F81IG	16	-779.0510	1651.6911	0.0000	1.0000
SYMI	17	-777.5050	1654.4484	0.0000	1.0000
SYMG	17	-777.5272	1654.4930	0.0000	1.0000
F81	14	-786.6263	1655.1430	0.0000	1.0000
JC69I	12	-792.8382	1655.8683	0.0000	1.0000
JC69G	12	-792.8543	1655.9005	0.0000	1.0000
SYMIG	18	-777.5273	1660.3424	0.0000	1.0000
JC69IG	13	-792.8544	1661.7500	0.0000	1.0000
SYM	16	-785.5429	1664.6750	0.0000	1.0000
JC69	11	-800.9636	1666.2698	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/chr1_5365.phylip = 1 - 347;
 Partition Dummy = 1:591-tomodel-phylip/chr1_5365.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/chr1_5365.phylip = 1 - 347;
 Partition Dummy = 1:591-tomodel-phylip/chr1_5365.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/chr1_5365.phylip = 1 - 347;
 Partition Dummy = 1:591-tomodel-phylip/chr1_5365.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

