
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
16:17 04/04/2011

Input data from file "683-tomodel-phylip/chr5_11648.phylip" (7 taxa, 462 characters)

Minimum AIC  model: HKYI
Minimum AICc model: HKYI
Minimum BIC  model: K2P


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	-920.1892	1873.6008	0.3173	0.3173
HKYG	16	-920.2361	1873.6946	0.3028	0.6201
HKY	15	-921.5714	1874.2190	0.2330	0.8531
HKYIG	17	-920.2465	1875.8714	0.1020	0.9551
GTRI	20	-919.3078	1880.5203	0.0100	0.9650
GTRG	20	-919.3682	1880.6411	0.0094	0.9744
GTR	19	-920.7324	1881.1842	0.0072	0.9816
F81I	15	-925.7008	1882.4779	0.0037	0.9853
F81G	15	-925.7405	1882.5573	0.0036	0.9889
GTRIG	21	-919.3822	1882.8645	0.0031	0.9920
F81	14	-926.9939	1882.9274	0.0030	0.9950
F81IG	16	-925.7480	1884.7184	0.0012	0.9963
K2PI	13	-929.0190	1884.8504	0.0011	0.9974
K2PG	13	-929.0609	1884.9342	0.0011	0.9985
K2P	12	-930.2715	1885.2379	0.0009	0.9994
K2PIG	14	-929.0680	1887.0757	0.0004	0.9998
SYMI	17	-928.0897	1891.5577	0.0000	0.9999
SYMG	17	-928.1417	1891.6618	0.0000	0.9999
SYM	16	-929.2787	1891.7799	0.0000	0.9999
JC69I	12	-934.2500	1893.1948	0.0000	0.9999
JC69G	12	-934.2888	1893.2725	0.0000	1.0000
JC69	11	-935.4491	1893.4848	0.0000	1.0000
SYMIG	18	-928.1490	1893.8420	0.0000	1.0000
JC69IG	13	-934.2946	1895.4018	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	16	-920.1892	1872.3783	0.3168	0.3168
HKYG	16	-920.2361	1872.4722	0.3022	0.6190
HKY	15	-921.5714	1873.1427	0.2161	0.8351
HKYIG	17	-920.2465	1874.4930	0.1100	0.9452
GTRI	20	-919.3078	1878.6155	0.0140	0.9592
GTRG	20	-919.3682	1878.7363	0.0132	0.9724
GTR	19	-920.7324	1879.4647	0.0092	0.9815
GTRIG	21	-919.3822	1880.7645	0.0048	0.9863
F81I	15	-925.7008	1881.4016	0.0035	0.9898
F81G	15	-925.7405	1881.4810	0.0033	0.9931
F81	14	-926.9939	1881.9878	0.0026	0.9957
F81IG	16	-925.7480	1883.4960	0.0012	0.9970
K2PI	13	-929.0190	1884.0379	0.0009	0.9979
K2PG	13	-929.0609	1884.1217	0.0009	0.9988
K2P	12	-930.2715	1884.5430	0.0007	0.9995
K2PIG	14	-929.0680	1886.1361	0.0003	0.9998
SYMI	17	-928.0897	1890.1793	0.0000	0.9999
SYMG	17	-928.1417	1890.2834	0.0000	0.9999
SYM	16	-929.2787	1890.5574	0.0000	0.9999
SYMIG	18	-928.1490	1892.2980	0.0000	1.0000
JC69I	12	-934.2500	1892.4999	0.0000	1.0000
JC69G	12	-934.2888	1892.5776	0.0000	1.0000
JC69	11	-935.4491	1892.8981	0.0000	1.0000
JC69IG	13	-934.2946	1894.5893	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2P	12	-930.2715	1934.1698	0.4184	0.4184
HKY	15	-921.5714	1935.1762	0.2529	0.6713
K2PI	13	-929.0190	1937.8003	0.0681	0.7394
K2PG	13	-929.0609	1937.8841	0.0653	0.8047
JC69	11	-935.4491	1938.3894	0.0507	0.8554
HKYI	16	-920.1892	1938.5474	0.0469	0.9023
HKYG	16	-920.2361	1938.6412	0.0447	0.9470
F81	14	-926.9939	1939.8857	0.0240	0.9710
JC69I	12	-934.2500	1942.1267	0.0078	0.9789
JC69G	12	-934.2888	1942.2044	0.0075	0.9864
F81I	15	-925.7008	1943.4351	0.0041	0.9905
F81G	15	-925.7405	1943.5145	0.0039	0.9944
K2PIG	14	-929.0680	1944.0340	0.0030	0.9974
HKYIG	17	-920.2465	1944.7976	0.0021	0.9995
JC69IG	13	-934.2946	1948.3516	0.0003	0.9998
F81IG	16	-925.7480	1949.6650	0.0002	1.0000
SYM	16	-929.2787	1956.7265	0.0000	1.0000
GTR	19	-920.7324	1958.0405	0.0000	1.0000
SYMI	17	-928.0897	1960.4839	0.0000	1.0000
SYMG	17	-928.1417	1960.5880	0.0000	1.0000
GTRI	20	-919.3078	1961.3268	0.0000	1.0000
GTRG	20	-919.3682	1961.4476	0.0000	1.0000
SYMIG	18	-928.1490	1966.7382	0.0000	1.0000
GTRIG	21	-919.3822	1967.6114	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 683-tomodel-phylip/chr5_11648.phylip = 1 - 462;
 Partition Dummy = 1:683-tomodel-phylip/chr5_11648.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 683-tomodel-phylip/chr5_11648.phylip = 1 - 462;
 Partition Dummy = 1:683-tomodel-phylip/chr5_11648.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 (K2P)]

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr5_11648.phylip = 1 - 462;
 Partition Dummy = 1:683-tomodel-phylip/chr5_11648.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0);
END;


-------------------------------------------------------------
End of Output

