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

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

Minimum AIC  model: GTRG
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	-751.4805	1536.5610	0.2202	0.2202
HKYI	16	-751.5544	1536.7088	0.2045	0.4247
GTRG	20	-747.1552	1536.8104	0.1944	0.6191
GTRI	20	-747.2153	1536.9306	0.1830	0.8021
HKYIG	17	-751.4517	1538.7088	0.0752	0.8774
GTRIG	21	-747.1486	1539.0553	0.0633	0.9407
GTR	19	-750.0726	1540.4004	0.0323	0.9729
HKY	15	-754.6746	1540.7568	0.0270	1.0000
K2PG	13	-764.6540	1556.3693	0.0000	1.0000
K2PI	13	-764.7252	1556.5117	0.0000	1.0000
K2PIG	14	-764.6010	1558.4301	0.0000	1.0000
SYMG	17	-761.6880	1559.1813	0.0000	1.0000
SYMI	17	-761.7557	1559.3167	0.0000	1.0000
K2P	12	-768.2085	1561.3240	0.0000	1.0000
SYMIG	18	-761.6679	1561.3594	0.0000	1.0000
SYM	16	-764.8808	1563.3615	0.0000	1.0000
F81I	15	-769.0798	1569.5672	0.0000	1.0000
F81G	15	-769.1038	1569.6151	0.0000	1.0000
F81IG	16	-769.0634	1571.7269	0.0000	1.0000
F81	14	-771.6622	1572.5524	0.0000	1.0000
JC69I	12	-781.2437	1587.3943	0.0000	1.0000
JC69G	12	-781.2640	1587.4350	0.0000	1.0000
JC69IG	13	-781.2737	1589.6086	0.0000	1.0000
JC69	11	-784.0320	1590.8292	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-747.1552	1534.3104	0.2377	0.2377
GTRI	20	-747.2153	1534.4306	0.2239	0.4616
HKYG	16	-751.4805	1534.9610	0.1717	0.6333
HKYI	16	-751.5544	1535.1088	0.1595	0.7928
GTRIG	21	-747.1486	1536.2971	0.0880	0.8809
HKYIG	17	-751.4517	1536.9035	0.0650	0.9459
GTR	19	-750.0726	1538.1452	0.0349	0.9808
HKY	15	-754.6746	1539.3492	0.0191	1.0000
K2PG	13	-764.6540	1555.3081	0.0000	1.0000
K2PI	13	-764.7252	1555.4505	0.0000	1.0000
K2PIG	14	-764.6010	1557.2021	0.0000	1.0000
SYMG	17	-761.6880	1557.3760	0.0000	1.0000
SYMI	17	-761.7557	1557.5114	0.0000	1.0000
SYMIG	18	-761.6679	1559.3357	0.0000	1.0000
K2P	12	-768.2085	1560.4170	0.0000	1.0000
SYM	16	-764.8808	1561.7615	0.0000	1.0000
F81I	15	-769.0798	1568.1596	0.0000	1.0000
F81G	15	-769.1038	1568.2075	0.0000	1.0000
F81IG	16	-769.0634	1570.1269	0.0000	1.0000
F81	14	-771.6622	1571.3243	0.0000	1.0000
JC69I	12	-781.2437	1586.4873	0.0000	1.0000
JC69G	12	-781.2640	1586.5281	0.0000	1.0000
JC69IG	13	-781.2737	1588.5474	0.0000	1.0000
JC69	11	-784.0320	1590.0640	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-751.4805	1597.0048	0.3581	0.3581
HKYI	16	-751.5544	1597.1526	0.3326	0.6908
HKY	15	-754.6746	1597.5152	0.2775	0.9682
HKYIG	17	-751.4517	1602.8250	0.0195	0.9878
K2PG	13	-764.6540	1605.7186	0.0046	0.9923
K2PI	13	-764.7252	1605.8611	0.0043	0.9966
K2P	12	-768.2085	1606.9498	0.0025	0.9991
K2PIG	14	-764.6010	1611.4904	0.0003	0.9994
GTR	19	-750.0726	1611.8222	0.0002	0.9996
GTRG	20	-747.1552	1611.8651	0.0002	0.9998
GTRI	20	-747.2153	1611.9853	0.0002	1.0000
GTRIG	21	-747.1486	1617.7296	0.0000	1.0000
SYMG	17	-761.6880	1623.2975	0.0000	1.0000
SYMI	17	-761.7557	1623.4329	0.0000	1.0000
SYM	16	-764.8808	1623.8053	0.0000	1.0000
F81	14	-771.6622	1625.6126	0.0000	1.0000
F81I	15	-769.0798	1626.3257	0.0000	1.0000
F81G	15	-769.1038	1626.3736	0.0000	1.0000
SYMIG	18	-761.6679	1629.1350	0.0000	1.0000
F81IG	16	-769.0634	1632.1706	0.0000	1.0000
JC69	11	-784.0320	1632.7191	0.0000	1.0000
JC69I	12	-781.2437	1633.0201	0.0000	1.0000
JC69G	12	-781.2640	1633.0609	0.0000	1.0000
JC69IG	13	-781.2737	1638.9580	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 591-tomodel-phylip/chr4_10567.phylip = 1 - 357;
 Partition Dummy = 1:591-tomodel-phylip/chr4_10567.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 (HKYG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr4_10567.phylip = 1 - 357;
 Partition Dummy = 1:591-tomodel-phylip/chr4_10567.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/chr4_10567.phylip = 1 - 357;
 Partition Dummy = 1:591-tomodel-phylip/chr4_10567.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

