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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
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
GTRG	20	-779.5491	1602.8651	0.5712	0.5712
GTRI	20	-780.6571	1605.0810	0.1886	0.7599
GTRIG	21	-779.5492	1605.2605	0.1724	0.9323
SYMG	17	-786.3956	1609.4992	0.0207	0.9530
K2PG	13	-791.6686	1610.9197	0.0102	0.9632
SYMI	17	-787.4030	1611.5140	0.0076	0.9708
GTR	19	-785.2131	1611.8191	0.0065	0.9773
SYMIG	18	-786.4027	1611.8454	0.0064	0.9837
HKYG	16	-788.9231	1612.2427	0.0053	0.9889
K2PI	13	-792.6528	1612.8882	0.0038	0.9927
K2PIG	14	-791.6694	1613.1729	0.0033	0.9960
HKYI	16	-789.8934	1614.1833	0.0020	0.9980
HKYIG	17	-788.9239	1614.5557	0.0017	0.9997
SYM	16	-792.0293	1618.4551	0.0002	0.9999
K2P	12	-798.0101	1621.3708	0.0001	1.0000
HKY	15	-795.2661	1622.6375	0.0000	1.0000
JC69G	12	-807.5581	1640.4668	0.0000	1.0000
F81G	15	-804.3692	1640.8437	0.0000	1.0000
JC69I	12	-807.8458	1641.0423	0.0000	1.0000
F81I	15	-804.5759	1641.2570	0.0000	1.0000
JC69IG	13	-807.5607	1642.7039	0.0000	1.0000
F81IG	16	-804.3723	1643.1411	0.0000	1.0000
F81	14	-809.2301	1648.2942	0.0000	1.0000
JC69	11	-812.8166	1648.7712	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-779.5491	1599.0983	0.5671	0.5671
GTRIG	21	-779.5492	1601.0984	0.2086	0.7758
GTRI	20	-780.6571	1601.3142	0.1873	0.9631
SYMG	17	-786.3956	1606.7912	0.0121	0.9752
GTR	19	-785.2131	1608.4262	0.0053	0.9805
SYMIG	18	-786.4027	1608.8054	0.0044	0.9849
SYMI	17	-787.4030	1608.8061	0.0044	0.9894
K2PG	13	-791.6686	1609.3371	0.0034	0.9927
HKYG	16	-788.9231	1609.8462	0.0026	0.9954
K2PI	13	-792.6528	1611.3056	0.0013	0.9966
K2PIG	14	-791.6694	1611.3389	0.0012	0.9979
HKYI	16	-789.8934	1611.7868	0.0010	0.9989
HKYIG	17	-788.9239	1611.8478	0.0010	0.9999
SYM	16	-792.0293	1616.0586	0.0001	1.0000
K2P	12	-798.0101	1620.0201	0.0000	1.0000
HKY	15	-795.2661	1620.5323	0.0000	1.0000
F81G	15	-804.3692	1638.7384	0.0000	1.0000
JC69G	12	-807.5581	1639.1162	0.0000	1.0000
F81I	15	-804.5759	1639.1517	0.0000	1.0000
JC69I	12	-807.8458	1639.6917	0.0000	1.0000
F81IG	16	-804.3723	1640.7446	0.0000	1.0000
JC69IG	13	-807.5607	1641.1213	0.0000	1.0000
F81	14	-809.2301	1646.4601	0.0000	1.0000
JC69	11	-812.8166	1647.6333	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	13	-791.6686	1654.8003	0.6769	0.6769
K2PI	13	-792.6528	1656.7688	0.2530	0.9299
K2PIG	14	-791.6694	1660.2992	0.0433	0.9732
K2P	12	-798.0101	1661.9862	0.0186	0.9918
HKYG	16	-788.9231	1665.8009	0.0028	0.9946
SYMG	17	-786.3956	1666.2431	0.0022	0.9968
HKYI	16	-789.8934	1667.7415	0.0010	0.9979
SYMI	17	-787.4030	1668.2579	0.0008	0.9987
GTRG	20	-779.5491	1669.0416	0.0005	0.9992
GTRI	20	-780.6571	1671.2576	0.0002	0.9994
HKYIG	17	-788.9239	1671.2996	0.0002	0.9996
SYMIG	18	-786.4027	1671.7544	0.0001	0.9997
SYM	16	-792.0293	1672.0133	0.0001	0.9999
HKY	15	-795.2661	1672.9898	0.0001	0.9999
GTRIG	21	-779.5492	1674.5389	0.0000	1.0000
GTR	19	-785.2131	1674.8724	0.0000	1.0000
JC69G	12	-807.5581	1681.0822	0.0000	1.0000
JC69I	12	-807.8458	1681.6577	0.0000	1.0000
JC69	11	-812.8166	1686.1022	0.0000	1.0000
JC69IG	13	-807.5607	1686.5845	0.0000	1.0000
F81G	15	-804.3692	1691.1959	0.0000	1.0000
F81I	15	-804.5759	1691.6092	0.0000	1.0000
F81	14	-809.2301	1695.4205	0.0000	1.0000
F81IG	16	-804.3723	1696.6993	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/chr2_7915.phylip = 1 - 244;
 Partition Dummy = 1:591-tomodel-phylip/chr2_7915.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 (GTRG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr2_7915.phylip = 1 - 244;
 Partition Dummy = 1:591-tomodel-phylip/chr2_7915.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 BIC model (K2PG)]

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

