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

Input data from file "183-phylip/chr2_24825.phylip" (29 taxa, 315 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	-925.2683	1999.3555	0.5044	0.5044
GTRG	64	-920.1198	2001.5196	0.1709	0.6754
HKYI	60	-926.6383	2002.0955	0.1282	0.8035
GTRI	64	-920.8756	2003.0312	0.0803	0.8838
HKYIG	61	-925.5835	2003.0642	0.0790	0.9628
GTRIG	65	-920.0626	2004.5831	0.0370	0.9997
K2PG	57	-938.1210	2015.9696	0.0001	0.9999
HKY	59	-935.9484	2017.6615	0.0001	0.9999
K2PI	57	-939.1751	2018.0779	0.0000	0.9999
K2PIG	58	-937.8881	2018.5106	0.0000	1.0000
GTR	63	-931.2753	2020.6782	0.0000	1.0000
SYMG	61	-935.6955	2023.2883	0.0000	1.0000
SYMI	61	-936.5624	2025.0220	0.0000	1.0000
SYMIG	62	-935.5072	2026.0145	0.0000	1.0000
F81G	59	-945.2776	2036.3199	0.0000	1.0000
K2P	56	-949.9868	2036.7178	0.0000	1.0000
F81I	59	-946.4521	2038.6688	0.0000	1.0000
F81IG	60	-945.0570	2038.9330	0.0000	1.0000
SYM	60	-947.5819	2043.9828	0.0000	1.0000
JC69G	56	-957.0899	2050.9240	0.0000	1.0000
JC69IG	57	-956.7886	2053.3049	0.0000	1.0000
JC69I	56	-958.3378	2053.4198	0.0000	1.0000
F81	58	-956.5291	2055.7926	0.0000	1.0000
JC69	55	-969.0677	2071.9192	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-920.1198	1968.2396	0.4270	0.4270
GTRI	64	-920.8756	1969.7512	0.2005	0.6275
GTRIG	65	-920.0626	1970.1252	0.1663	0.7938
HKYG	60	-925.2683	1970.5366	0.1354	0.9292
HKYIG	61	-925.5835	1973.1670	0.0363	0.9656
HKYI	60	-926.6383	1973.2766	0.0344	1.0000
GTR	63	-931.2753	1988.5507	0.0000	1.0000
HKY	59	-935.9484	1989.8968	0.0000	1.0000
K2PG	57	-938.1210	1990.2420	0.0000	1.0000
K2PIG	58	-937.8881	1991.7762	0.0000	1.0000
K2PI	57	-939.1751	1992.3503	0.0000	1.0000
SYMG	61	-935.6955	1993.3911	0.0000	1.0000
SYMIG	62	-935.5072	1995.0145	0.0000	1.0000
SYMI	61	-936.5624	1995.1248	0.0000	1.0000
F81G	59	-945.2776	2008.5552	0.0000	1.0000
F81IG	60	-945.0570	2010.1141	0.0000	1.0000
F81I	59	-946.4521	2010.9041	0.0000	1.0000
K2P	56	-949.9868	2011.9736	0.0000	1.0000
SYM	60	-947.5819	2015.1639	0.0000	1.0000
JC69G	56	-957.0899	2026.1798	0.0000	1.0000
JC69IG	57	-956.7886	2027.5773	0.0000	1.0000
JC69I	56	-958.3378	2028.6756	0.0000	1.0000
F81	58	-956.5291	2029.0582	0.0000	1.0000
JC69	55	-969.0677	2048.1354	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-925.2683	2195.6909	0.7581	0.7581
HKYI	60	-926.6383	2198.4310	0.1926	0.9507
HKYIG	61	-925.5835	2202.0739	0.0312	0.9819
K2PG	57	-938.1210	2204.1386	0.0111	0.9930
K2PI	57	-939.1751	2206.2469	0.0039	0.9969
GTRG	64	-920.1198	2208.4043	0.0013	0.9982
K2PIG	58	-937.8881	2209.4254	0.0008	0.9990
GTRI	64	-920.8756	2209.9159	0.0006	0.9996
HKY	59	-935.9484	2211.2986	0.0003	0.9999
GTRIG	65	-920.0626	2214.0425	0.0001	1.0000
K2P	56	-949.9868	2222.1177	0.0000	1.0000
SYMG	61	-935.6955	2222.2980	0.0000	1.0000
SYMI	61	-936.5624	2224.0317	0.0000	1.0000
GTR	63	-931.2753	2224.9627	0.0000	1.0000
SYMIG	62	-935.5072	2227.6740	0.0000	1.0000
F81G	59	-945.2776	2229.9570	0.0000	1.0000
F81I	59	-946.4521	2232.3059	0.0000	1.0000
F81IG	60	-945.0570	2235.2684	0.0000	1.0000
JC69G	56	-957.0899	2236.3238	0.0000	1.0000
JC69I	56	-958.3378	2238.8197	0.0000	1.0000
SYM	60	-947.5819	2240.3182	0.0000	1.0000
JC69IG	57	-956.7886	2241.4739	0.0000	1.0000
F81	58	-956.5291	2246.7074	0.0000	1.0000
JC69	55	-969.0677	2254.5269	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_24825.phylip = 1 - 315;
 Partition Dummy = 1:183-phylip/chr2_24825.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 183-phylip/chr2_24825.phylip = 1 - 315;
 Partition Dummy = 1:183-phylip/chr2_24825.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 183-phylip/chr2_24825.phylip = 1 - 315;
 Partition Dummy = 1:183-phylip/chr2_24825.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

