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

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

Minimum AIC  model: GTRIG
Minimum AICc model: GTRIG
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
GTRIG	21	-1689.2021	3422.1153	0.4212	0.4212
GTRG	20	-1690.4069	3422.3665	0.3715	0.7927
HKYIG	17	-1695.2872	3425.6995	0.0702	0.8629
HKYG	16	-1696.5688	3426.1358	0.0564	0.9193
GTRI	20	-1692.4935	3426.5398	0.0461	0.9655
K2PIG	14	-1700.4936	3429.7551	0.0092	0.9747
K2PG	13	-1701.6943	3430.0528	0.0080	0.9827
SYMIG	18	-1696.6951	3430.6499	0.0059	0.9886
SYMG	17	-1697.8421	3430.8092	0.0055	0.9940
HKYI	16	-1699.0484	3431.0950	0.0047	0.9987
K2PI	13	-1704.1554	3434.9751	0.0007	0.9994
SYMI	17	-1700.0786	3435.2823	0.0006	1.0000
F81IG	16	-1721.7406	3476.4795	0.0000	1.0000
F81G	15	-1722.8958	3476.6708	0.0000	1.0000
JC69IG	13	-1726.5138	3479.6919	0.0000	1.0000
JC69G	12	-1727.5924	3479.7531	0.0000	1.0000
F81I	15	-1725.0978	3481.0748	0.0000	1.0000
JC69I	12	-1729.6988	3483.9659	0.0000	1.0000
HKY	15	-1733.6430	3498.1651	0.0000	1.0000
K2P	12	-1737.3094	3499.1871	0.0000	1.0000
GTR	19	-1729.9021	3499.2065	0.0000	1.0000
SYM	16	-1734.5100	3502.0182	0.0000	1.0000
F81	14	-1757.4245	3543.6168	0.0000	1.0000
JC69	11	-1761.0863	3544.6527	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	21	-1689.2021	3420.4041	0.4567	0.4567
GTRG	20	-1690.4069	3420.8139	0.3721	0.8288
HKYIG	17	-1695.2872	3424.5745	0.0568	0.8856
GTRI	20	-1692.4935	3424.9871	0.0462	0.9318
HKYG	16	-1696.5688	3425.1376	0.0428	0.9746
K2PIG	14	-1700.4936	3428.9873	0.0062	0.9809
K2PG	13	-1701.6943	3429.3885	0.0051	0.9860
SYMIG	18	-1696.6951	3429.3902	0.0051	0.9911
SYMG	17	-1697.8421	3429.6842	0.0044	0.9955
HKYI	16	-1699.0484	3430.0968	0.0036	0.9991
SYMI	17	-1700.0786	3434.1573	0.0005	0.9996
K2PI	13	-1704.1554	3434.3108	0.0004	1.0000
F81IG	16	-1721.7406	3475.4813	0.0000	1.0000
F81G	15	-1722.8958	3475.7917	0.0000	1.0000
JC69IG	13	-1726.5138	3479.0277	0.0000	1.0000
JC69G	12	-1727.5924	3479.1848	0.0000	1.0000
F81I	15	-1725.0978	3480.1957	0.0000	1.0000
JC69I	12	-1729.6988	3483.3976	0.0000	1.0000
HKY	15	-1733.6430	3497.2860	0.0000	1.0000
GTR	19	-1729.9021	3497.8043	0.0000	1.0000
K2P	12	-1737.3094	3498.6188	0.0000	1.0000
SYM	16	-1734.5100	3501.0201	0.0000	1.0000
F81	14	-1757.4245	3542.8490	0.0000	1.0000
JC69	11	-1761.0863	3544.1727	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	13	-1701.6943	3485.6980	0.8056	0.8056
K2PIG	14	-1700.4936	3489.6283	0.1129	0.9185
K2PI	13	-1704.1554	3490.6204	0.0687	0.9873
HKYG	16	-1696.5688	3494.4417	0.0102	0.9974
HKYIG	17	-1695.2872	3498.2100	0.0015	0.9990
HKYI	16	-1699.0484	3499.4008	0.0009	0.9998
SYMG	17	-1697.8421	3503.3198	0.0001	1.0000
SYMIG	18	-1696.6951	3507.3573	0.0000	1.0000
GTRG	20	-1690.4069	3507.4439	0.0000	1.0000
SYMI	17	-1700.0786	3507.7928	0.0000	1.0000
GTRIG	21	-1689.2021	3511.3657	0.0000	1.0000
GTRI	20	-1692.4935	3511.6171	0.0000	1.0000
JC69G	12	-1727.5924	3531.1628	0.0000	1.0000
JC69IG	13	-1726.5138	3535.3372	0.0000	1.0000
JC69I	12	-1729.6988	3535.3756	0.0000	1.0000
F81G	15	-1722.8958	3540.7642	0.0000	1.0000
F81IG	16	-1721.7406	3544.7853	0.0000	1.0000
F81I	15	-1725.0978	3545.1682	0.0000	1.0000
K2P	12	-1737.3094	3550.5968	0.0000	1.0000
HKY	15	-1733.6430	3562.2585	0.0000	1.0000
SYM	16	-1734.5100	3570.3241	0.0000	1.0000
GTR	19	-1729.9021	3580.1028	0.0000	1.0000
JC69	11	-1761.0863	3591.8192	0.0000	1.0000
F81	14	-1757.4245	3603.4900	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr9_4928.phylip = 1 - 562;
 Partition Dummy = 1:591-tomodel-phylip/chr9_4928.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr9_4928.phylip = 1 - 562;
 Partition Dummy = 1:591-tomodel-phylip/chr9_4928.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

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

