
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
01:00 02/17/2011

Input data from file "183-phylip/chr8_3279.phylip" (29 taxa, 287 characters)

Minimum AIC  model: HKYG
Minimum AICc model: K2P
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
K2P	56	-893.3375	1926.4315	0.2826	0.2826
K2PG	57	-891.8120	1926.4973	0.2734	0.5560
K2PI	57	-891.8723	1926.6179	0.2574	0.8134
HKY	59	-890.3189	1929.8273	0.0517	0.8652
HKYG	60	-888.7632	1929.9158	0.0495	0.9147
HKYI	60	-888.8103	1930.0100	0.0472	0.9619
HKYIG	61	-888.7660	1933.1498	0.0098	0.9717
K2PIG	58	-893.6135	1933.2446	0.0094	0.9811
JC69	55	-898.9884	1934.6434	0.0047	0.9857
JC69G	56	-897.4965	1934.7496	0.0044	0.9901
JC69I	56	-897.5418	1934.8402	0.0042	0.9944
SYM	60	-892.7601	1937.9097	0.0009	0.9953
SYMG	61	-891.1879	1937.9936	0.0009	0.9961
F81	58	-896.1305	1938.2785	0.0008	0.9969
F81G	59	-894.6374	1938.4643	0.0007	0.9976
F81I	59	-894.6783	1938.5461	0.0007	0.9983
GTRG	64	-887.0569	1939.5912	0.0004	0.9986
GTR	63	-888.7882	1939.7378	0.0004	0.9990
GTRI	64	-887.1323	1939.7422	0.0004	0.9994
SYMIG	62	-891.1912	1941.2573	0.0002	0.9995
SYMI	61	-893.0213	1941.6604	0.0001	0.9997
F81IG	60	-894.6367	1941.6629	0.0001	0.9998
JC69IG	57	-899.6820	1942.2374	0.0001	0.9999
GTRIG	65	-887.0578	1942.9391	0.0001	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-888.7632	1897.5264	0.1737	0.1737
HKYI	60	-888.8103	1897.6206	0.1657	0.3393
K2PG	57	-891.8120	1897.6239	0.1654	0.5048
K2PI	57	-891.8723	1897.7446	0.1557	0.6605
HKY	59	-890.3189	1898.6379	0.0996	0.7601
K2P	56	-893.3375	1898.6749	0.0978	0.8579
HKYIG	61	-888.7660	1899.5320	0.0637	0.9216
GTRG	64	-887.0569	1902.1138	0.0175	0.9391
GTRI	64	-887.1323	1902.2647	0.0162	0.9554
K2PIG	58	-893.6135	1903.2271	0.0100	0.9654
GTR	63	-888.7882	1903.5763	0.0084	0.9739
GTRIG	65	-887.0578	1904.1156	0.0064	0.9803
SYMG	61	-891.1879	1904.3758	0.0057	0.9859
SYM	60	-892.7601	1905.5203	0.0032	0.9891
SYMIG	62	-891.1912	1906.3823	0.0021	0.9912
JC69G	56	-897.4965	1906.9931	0.0015	0.9927
JC69I	56	-897.5418	1907.0837	0.0015	0.9942
F81G	59	-894.6374	1907.2748	0.0013	0.9955
F81I	59	-894.6783	1907.3566	0.0013	0.9968
JC69	55	-898.9884	1907.9768	0.0009	0.9977
SYMI	61	-893.0213	1908.0426	0.0009	0.9986
F81	58	-896.1305	1908.2610	0.0008	0.9994
F81IG	60	-894.6367	1909.2735	0.0005	0.9999
JC69IG	57	-899.6820	1913.3640	0.0001	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2P	56	-893.3375	2103.6059	0.6148	0.6148
K2PG	57	-891.8120	2106.2144	0.1668	0.7816
K2PI	57	-891.8723	2106.3351	0.1571	0.9387
JC69	55	-898.9884	2109.2483	0.0366	0.9753
JC69G	56	-897.4965	2111.9241	0.0096	0.9849
JC69I	56	-897.5418	2112.0147	0.0092	0.9940
HKY	59	-890.3189	2114.5473	0.0026	0.9966
K2PIG	58	-893.6135	2115.4770	0.0016	0.9983
HKYG	60	-888.7632	2117.0954	0.0007	0.9990
HKYI	60	-888.8103	2117.1896	0.0007	0.9997
F81	58	-896.1305	2120.5110	0.0001	0.9998
JC69IG	57	-899.6820	2121.9545	0.0001	0.9999
HKYIG	61	-888.7660	2122.7605	0.0000	0.9999
F81G	59	-894.6374	2123.1843	0.0000	0.9999
F81I	59	-894.6783	2123.2661	0.0000	1.0000
SYM	60	-892.7601	2125.0892	0.0000	1.0000
SYMG	61	-891.1879	2127.6042	0.0000	1.0000
F81IG	60	-894.6367	2128.8424	0.0000	1.0000
SYMI	61	-893.0213	2131.2710	0.0000	1.0000
SYMIG	62	-891.1912	2133.2702	0.0000	1.0000
GTR	63	-888.7882	2134.1237	0.0000	1.0000
GTRG	64	-887.0569	2136.3206	0.0000	1.0000
GTRI	64	-887.1323	2136.4716	0.0000	1.0000
GTRIG	65	-887.0578	2141.9819	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 183-phylip/chr8_3279.phylip = 1 - 287;
 Partition Dummy = 1:183-phylip/chr8_3279.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 AICc model (K2P)]

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


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

BEGIN MRBAYES;
 Charset 183-phylip/chr8_3279.phylip = 1 - 287;
 Partition Dummy = 1:183-phylip/chr8_3279.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

