
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
11:53 08/01/2011

Input data from file "261-aggressive-phylip/7733_hg18.phylip" (29 taxa, 342 characters)

Minimum AIC  model: GTRG
Minimum AICc model: HKYG
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
HKYG	60	-898.7049	1943.4596	0.4261	0.4261
HKYI	60	-899.6224	1945.2947	0.1702	0.5963
HKYIG	61	-898.2722	1945.5587	0.1492	0.7455
GTRG	64	-893.7843	1945.6047	0.1458	0.8913
GTRI	64	-894.7723	1947.5806	0.0543	0.9456
GTRIG	65	-893.4287	1947.9444	0.0453	0.9908
K2PG	57	-907.7195	1952.7206	0.0042	0.9950
K2PIG	58	-907.2896	1954.7630	0.0015	0.9965
K2PI	57	-908.7564	1954.7946	0.0015	0.9979
SYMG	61	-903.1967	1955.4077	0.0011	0.9990
SYMIG	62	-902.8360	1957.6720	0.0003	0.9994
SYMI	61	-904.3848	1957.7839	0.0003	0.9997
HKY	59	-907.9327	1958.9719	0.0002	0.9999
GTR	63	-902.5597	1960.1266	0.0001	1.0000
F81G	59	-911.8621	1966.8305	0.0000	1.0000
K2P	56	-916.5743	1967.5487	0.0000	1.0000
F81I	59	-912.8683	1968.8431	0.0000	1.0000
F81IG	60	-911.4526	1968.9550	0.0000	1.0000
SYM	60	-911.5573	1969.1643	0.0000	1.0000
JC69G	56	-919.5804	1973.5609	0.0000	1.0000
JC69IG	57	-919.0038	1975.2893	0.0000	1.0000
JC69I	56	-920.6451	1975.6902	0.0000	1.0000
F81	58	-920.6204	1981.4246	0.0000	1.0000
JC69	55	-928.2490	1988.0364	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-893.7843	1915.5686	0.3723	0.3723
GTRIG	65	-893.4287	1916.8574	0.1954	0.5678
HKYG	60	-898.7049	1917.4098	0.1483	0.7160
GTRI	64	-894.7723	1917.5445	0.1386	0.8547
HKYIG	61	-898.2722	1918.5444	0.0841	0.9387
HKYI	60	-899.6224	1919.2449	0.0592	0.9980
SYMG	61	-903.1967	1928.3934	0.0006	0.9986
K2PG	57	-907.7195	1929.4389	0.0004	0.9990
SYMIG	62	-902.8360	1929.6720	0.0003	0.9993
K2PIG	58	-907.2896	1930.5792	0.0002	0.9995
SYMI	61	-904.3848	1930.7696	0.0002	0.9997
GTR	63	-902.5597	1931.1194	0.0002	0.9998
K2PI	57	-908.7564	1931.5129	0.0001	1.0000
HKY	59	-907.9327	1933.8655	0.0000	1.0000
F81G	59	-911.8621	1941.7242	0.0000	1.0000
F81IG	60	-911.4526	1942.9052	0.0000	1.0000
SYM	60	-911.5573	1943.1145	0.0000	1.0000
F81I	59	-912.8683	1943.7367	0.0000	1.0000
K2P	56	-916.5743	1945.1487	0.0000	1.0000
JC69G	56	-919.5804	1951.1609	0.0000	1.0000
JC69IG	57	-919.0038	1952.0076	0.0000	1.0000
JC69I	56	-920.6451	1953.2902	0.0000	1.0000
F81	58	-920.6204	1957.2408	0.0000	1.0000
JC69	55	-928.2490	1966.4979	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-898.7049	2147.4985	0.3853	0.3853
K2PG	57	-907.7195	2148.0231	0.2964	0.6817
HKYI	60	-899.6224	2149.3335	0.1539	0.8356
K2PI	57	-908.7564	2150.0971	0.1051	0.9407
HKYIG	61	-898.2722	2152.4679	0.0321	0.9728
K2PIG	58	-907.2896	2152.9983	0.0246	0.9975
K2P	56	-916.5743	2159.8981	0.0008	0.9982
HKY	59	-907.9327	2160.1193	0.0007	0.9989
GTRG	64	-893.7843	2160.9964	0.0005	0.9994
SYMG	61	-903.1967	2162.3168	0.0002	0.9996
GTRI	64	-894.7723	2162.9724	0.0002	0.9998
SYMI	61	-904.3848	2164.6931	0.0001	0.9999
JC69G	56	-919.5804	2165.9103	0.0000	0.9999
GTRIG	65	-893.4287	2166.1201	0.0000	0.9999
SYMIG	62	-902.8360	2167.4302	0.0000	1.0000
F81G	59	-911.8621	2167.9780	0.0000	1.0000
JC69I	56	-920.6451	2168.0396	0.0000	1.0000
F81I	59	-912.8683	2169.9905	0.0000	1.0000
JC69IG	57	-919.0038	2170.5919	0.0000	1.0000
GTR	63	-902.5597	2172.7124	0.0000	1.0000
F81IG	60	-911.4526	2172.9939	0.0000	1.0000
SYM	60	-911.5573	2173.2032	0.0000	1.0000
JC69	55	-928.2490	2177.4125	0.0000	1.0000
F81	58	-920.6204	2179.6599	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 261-aggressive-phylip/7733_hg18.phylip = 1 - 342;
 Partition Dummy = 1:261-aggressive-phylip/7733_hg18.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 (HKYG)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/7733_hg18.phylip = 1 - 342;
 Partition Dummy = 1:261-aggressive-phylip/7733_hg18.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 BIC model (HKYG)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/7733_hg18.phylip = 1 - 342;
 Partition Dummy = 1:261-aggressive-phylip/7733_hg18.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

