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

Input data from file "183-phylip/chr2_23670.phylip" (29 taxa, 314 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	64	-878.7404	1918.8945	0.6506	0.6506
GTRI	64	-879.9879	1921.3895	0.1869	0.8375
GTRIG	65	-878.7427	1922.0821	0.1322	0.9697
HKYG	60	-888.6239	1926.1805	0.0170	0.9867
HKYIG	61	-888.2733	1928.5625	0.0052	0.9919
HKYI	60	-889.8348	1928.6024	0.0051	0.9970
GTR	63	-886.2761	1930.8081	0.0017	0.9986
SYMG	61	-890.1448	1932.3054	0.0008	0.9994
SYMIG	62	-889.7643	1934.6522	0.0002	0.9997
SYMI	61	-891.6105	1935.2369	0.0002	0.9999
F81G	59	-896.0168	1937.9076	0.0000	0.9999
HKY	59	-896.6941	1939.2623	0.0000	0.9999
F81IG	60	-895.6755	1940.2838	0.0000	1.0000
F81I	59	-897.2473	1940.3686	0.0000	1.0000
K2PG	57	-900.2982	1940.4246	0.0000	1.0000
K2PIG	58	-899.8941	1942.6275	0.0000	1.0000
K2PI	57	-901.5070	1942.8422	0.0000	1.0000
SYM	60	-898.3162	1945.5652	0.0000	1.0000
F81	58	-903.6642	1950.1675	0.0000	1.0000
JC69G	56	-907.1420	1951.1244	0.0000	1.0000
JC69IG	57	-906.7740	1953.3762	0.0000	1.0000
JC69I	56	-908.3744	1953.5892	0.0000	1.0000
K2P	56	-908.8287	1954.4978	0.0000	1.0000
JC69	55	-915.0872	1964.0504	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-878.7404	1885.4809	0.6020	0.6020
GTRIG	65	-878.7427	1887.4853	0.2210	0.8229
GTRI	64	-879.9879	1887.9759	0.1729	0.9958
HKYG	60	-888.6239	1897.2477	0.0017	0.9975
HKYIG	61	-888.2733	1898.5467	0.0009	0.9984
GTR	63	-886.2761	1898.5521	0.0009	0.9993
HKYI	60	-889.8348	1899.6696	0.0005	0.9998
SYMG	61	-890.1448	1902.2895	0.0001	0.9999
SYMIG	62	-889.7643	1903.5287	0.0001	1.0000
SYMI	61	-891.6105	1905.2210	0.0000	1.0000
F81G	59	-896.0168	1910.0335	0.0000	1.0000
F81IG	60	-895.6755	1911.3510	0.0000	1.0000
HKY	59	-896.6941	1911.3883	0.0000	1.0000
F81I	59	-897.2473	1912.4946	0.0000	1.0000
K2PG	57	-900.2982	1914.5964	0.0000	1.0000
K2PIG	58	-899.8941	1915.7883	0.0000	1.0000
SYM	60	-898.3162	1916.6324	0.0000	1.0000
K2PI	57	-901.5070	1917.0141	0.0000	1.0000
F81	58	-903.6642	1923.3283	0.0000	1.0000
JC69G	56	-907.1420	1926.2839	0.0000	1.0000
JC69IG	57	-906.7740	1927.5481	0.0000	1.0000
JC69I	56	-908.3744	1928.7487	0.0000	1.0000
K2P	56	-908.8287	1929.6573	0.0000	1.0000
JC69	55	-915.0872	1940.1745	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-888.6239	2122.2113	0.5712	0.5712
HKYI	60	-889.8348	2124.6332	0.1702	0.7413
GTRG	64	-878.7404	2125.4420	0.1136	0.8549
HKYIG	61	-888.2733	2127.2596	0.0458	0.9007
GTRI	64	-879.9879	2127.9370	0.0326	0.9333
K2PG	57	-900.2982	2128.3118	0.0270	0.9603
K2PI	57	-901.5070	2130.7295	0.0081	0.9684
SYMG	61	-890.1448	2131.0025	0.0070	0.9754
GTRIG	65	-878.7427	2131.1958	0.0064	0.9818
F81G	59	-896.0168	2131.2477	0.0062	0.9881
HKY	59	-896.6941	2132.6024	0.0032	0.9912
K2PIG	58	-899.8941	2133.2531	0.0023	0.9935
F81I	59	-897.2473	2133.7088	0.0018	0.9953
SYMI	61	-891.6105	2133.9340	0.0016	0.9970
GTR	63	-886.2761	2134.7639	0.0011	0.9980
SYMIG	62	-889.7643	2135.9910	0.0006	0.9986
JC69G	56	-907.1420	2136.2499	0.0005	0.9991
F81IG	60	-895.6755	2136.3145	0.0005	0.9996
JC69I	56	-908.3744	2138.7147	0.0001	0.9998
K2P	56	-908.8287	2139.6233	0.0001	0.9999
F81	58	-903.6642	2140.7931	0.0001	0.9999
JC69IG	57	-906.7740	2141.2635	0.0000	1.0000
SYM	60	-898.3162	2141.5960	0.0000	1.0000
JC69	55	-915.0872	2146.3911	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 183-phylip/chr2_23670.phylip = 1 - 314;
 Partition Dummy = 1:183-phylip/chr2_23670.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 183-phylip/chr2_23670.phylip = 1 - 314;
 Partition Dummy = 1:183-phylip/chr2_23670.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 183-phylip/chr2_23670.phylip = 1 - 314;
 Partition Dummy = 1:183-phylip/chr2_23670.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

