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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
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
GTRG	20	-1883.1699	3807.4643	0.2815	0.2815
HKYG	16	-1887.4037	3807.5317	0.2721	0.5536
HKYI	16	-1888.1930	3809.1103	0.1236	0.6772
GTRIG	21	-1882.9547	3809.1481	0.1213	0.7985
HKYIG	17	-1887.2276	3809.2712	0.1140	0.9125
GTRI	20	-1884.3819	3809.8883	0.0838	0.9963
F81G	15	-1893.3569	3817.3520	0.0020	0.9983
F81I	15	-1894.1478	3818.9338	0.0009	0.9992
F81IG	16	-1893.2143	3819.1530	0.0008	1.0000
HKY	15	-1901.9269	3834.4920	0.0000	1.0000
GTR	19	-1897.9931	3835.0022	0.0000	1.0000
F81	14	-1907.6079	3843.7735	0.0000	1.0000
SYMG	17	-1918.1033	3871.0226	0.0000	1.0000
K2PG	13	-1922.4230	3871.3288	0.0000	1.0000
K2PI	13	-1923.0717	3872.6261	0.0000	1.0000
SYMIG	18	-1917.9235	3872.7603	0.0000	1.0000
K2PIG	14	-1922.2796	3873.1170	0.0000	1.0000
SYMI	17	-1919.4143	3873.6446	0.0000	1.0000
JC69G	12	-1929.5429	3883.4991	0.0000	1.0000
JC69I	12	-1930.1615	3884.7362	0.0000	1.0000
JC69IG	13	-1929.4394	3885.3616	0.0000	1.0000
K2P	12	-1936.7339	3897.8811	0.0000	1.0000
SYM	16	-1932.8341	3898.3926	0.0000	1.0000
JC69	11	-1943.6826	3909.7143	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-1883.1699	3806.3398	0.3065	0.3065
HKYG	16	-1887.4037	3806.8073	0.2426	0.5491
GTRIG	21	-1882.9547	3807.9095	0.1398	0.6889
HKYI	16	-1888.1930	3808.3860	0.1102	0.7991
HKYIG	17	-1887.2276	3808.4552	0.1064	0.9056
GTRI	20	-1884.3819	3808.7638	0.0912	0.9968
F81G	15	-1893.3569	3816.7137	0.0017	0.9985
F81I	15	-1894.1478	3818.2955	0.0008	0.9993
F81IG	16	-1893.2143	3818.4287	0.0007	1.0000
HKY	15	-1901.9269	3833.8537	0.0000	1.0000
GTR	19	-1897.9931	3833.9862	0.0000	1.0000
F81	14	-1907.6079	3843.2157	0.0000	1.0000
SYMG	17	-1918.1033	3870.2066	0.0000	1.0000
K2PG	13	-1922.4230	3870.8460	0.0000	1.0000
SYMIG	18	-1917.9235	3871.8471	0.0000	1.0000
K2PI	13	-1923.0717	3872.1433	0.0000	1.0000
K2PIG	14	-1922.2796	3872.5593	0.0000	1.0000
SYMI	17	-1919.4143	3872.8286	0.0000	1.0000
JC69G	12	-1929.5429	3883.0859	0.0000	1.0000
JC69I	12	-1930.1615	3884.3229	0.0000	1.0000
JC69IG	13	-1929.4394	3884.8789	0.0000	1.0000
K2P	12	-1936.7339	3897.4679	0.0000	1.0000
SYM	16	-1932.8341	3897.6683	0.0000	1.0000
JC69	11	-1943.6826	3909.3651	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1887.4037	3881.1080	0.6231	0.6231
HKYI	16	-1888.1930	3882.6866	0.2830	0.9060
F81G	15	-1893.3569	3886.3706	0.0449	0.9509
HKYIG	17	-1887.2276	3887.3996	0.0268	0.9777
F81I	15	-1894.1478	3887.9524	0.0203	0.9980
F81IG	16	-1893.2143	3892.7293	0.0019	0.9999
GTRG	20	-1883.1699	3899.2156	0.0001	1.0000
GTRI	20	-1884.3819	3901.6396	0.0000	1.0000
HKY	15	-1901.9269	3903.5106	0.0000	1.0000
GTRIG	21	-1882.9547	3905.4290	0.0000	1.0000
F81	14	-1907.6079	3908.2288	0.0000	1.0000
GTR	19	-1897.9931	3922.2182	0.0000	1.0000
K2PG	13	-1922.4230	3931.2153	0.0000	1.0000
K2PI	13	-1923.0717	3932.5126	0.0000	1.0000
K2PIG	14	-1922.2796	3937.5723	0.0000	1.0000
JC69G	12	-1929.5429	3938.8113	0.0000	1.0000
JC69I	12	-1930.1615	3940.0484	0.0000	1.0000
JC69IG	13	-1929.4394	3945.2481	0.0000	1.0000
SYMG	17	-1918.1033	3949.1510	0.0000	1.0000
SYMI	17	-1919.4143	3951.7730	0.0000	1.0000
K2P	12	-1936.7339	3953.1933	0.0000	1.0000
SYMIG	18	-1917.9235	3955.4353	0.0000	1.0000
JC69	11	-1943.6826	3960.4468	0.0000	1.0000
SYM	16	-1932.8341	3971.9689	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/chr5_2303.phylip = 1 - 768;
 Partition Dummy = 1:591-tomodel-phylip/chr5_2303.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/chr5_2303.phylip = 1 - 768;
 Partition Dummy = 1:591-tomodel-phylip/chr5_2303.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 (HKYG)]

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

