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

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

Minimum AIC  model: HKYG
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	-879.1313	1903.4173	0.4546	0.4546
HKYI	60	-879.1848	1903.5241	0.4310	0.8856
HKYIG	61	-879.0791	1906.2410	0.1108	0.9964
GTRG	64	-878.8574	1914.7044	0.0016	0.9980
GTRI	64	-878.9268	1914.8432	0.0015	0.9995
GTRIG	65	-878.7912	1917.5824	0.0004	0.9999
HKY	59	-888.9192	1920.0850	0.0001	1.0000
GTR	63	-888.6619	1931.3239	0.0000	1.0000
F81G	59	-904.0736	1950.3938	0.0000	1.0000
F81I	59	-904.1309	1950.5084	0.0000	1.0000
K2PG	57	-908.0763	1952.6423	0.0000	1.0000
K2PI	57	-908.1043	1952.6984	0.0000	1.0000
F81IG	60	-904.0334	1953.2214	0.0000	1.0000
K2PIG	58	-908.0618	1955.4820	0.0000	1.0000
SYMG	61	-907.0783	1962.2393	0.0000	1.0000
SYMI	61	-907.1217	1962.3262	0.0000	1.0000
SYMIG	62	-907.0394	1965.1099	0.0000	1.0000
F81	58	-913.4024	1966.1632	0.0000	1.0000
K2P	56	-918.0195	1969.6798	0.0000	1.0000
SYM	60	-917.1343	1979.4232	0.0000	1.0000
JC69G	56	-929.4436	1992.5279	0.0000	1.0000
JC69I	56	-929.4904	1992.6215	0.0000	1.0000
JC69IG	57	-929.4091	1995.3079	0.0000	1.0000
JC69	55	-938.9760	2008.7627	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-879.1313	1878.2627	0.4181	0.4181
HKYI	60	-879.1848	1878.3695	0.3964	0.8145
HKYIG	61	-879.0791	1880.1583	0.1621	0.9765
GTRG	64	-878.8574	1885.7149	0.0101	0.9866
GTRI	64	-878.9268	1885.8536	0.0094	0.9960
GTRIG	65	-878.7912	1887.5824	0.0040	0.9999
HKY	59	-888.9192	1895.8384	0.0001	1.0000
GTR	63	-888.6619	1903.3239	0.0000	1.0000
F81G	59	-904.0736	1926.1472	0.0000	1.0000
F81I	59	-904.1309	1926.2618	0.0000	1.0000
F81IG	60	-904.0334	1928.0668	0.0000	1.0000
K2PG	57	-908.0763	1930.1525	0.0000	1.0000
K2PI	57	-908.1043	1930.2086	0.0000	1.0000
K2PIG	58	-908.0618	1932.1236	0.0000	1.0000
SYMG	61	-907.0783	1936.1566	0.0000	1.0000
SYMI	61	-907.1217	1936.2434	0.0000	1.0000
SYMIG	62	-907.0394	1938.0788	0.0000	1.0000
F81	58	-913.4024	1942.8049	0.0000	1.0000
K2P	56	-918.0195	1948.0391	0.0000	1.0000
SYM	60	-917.1343	1954.2686	0.0000	1.0000
JC69G	56	-929.4436	1970.8872	0.0000	1.0000
JC69I	56	-929.4904	1970.9809	0.0000	1.0000
JC69IG	57	-929.4091	1972.8181	0.0000	1.0000
JC69	55	-938.9760	1987.9519	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-879.1313	2110.0806	0.4987	0.4987
HKYI	60	-879.1848	2110.1874	0.4728	0.9715
HKYIG	61	-879.0791	2115.8398	0.0280	0.9995
HKY	59	-888.9192	2123.7926	0.0005	1.0000
GTRG	64	-878.8574	2132.9873	0.0000	1.0000
GTRI	64	-878.9268	2133.1260	0.0000	1.0000
GTRIG	65	-878.7912	2138.7184	0.0000	1.0000
GTR	63	-888.6619	2146.7326	0.0000	1.0000
K2PG	57	-908.0763	2150.3795	0.0000	1.0000
K2PI	57	-908.1043	2150.4356	0.0000	1.0000
F81G	59	-904.0736	2154.1014	0.0000	1.0000
F81I	59	-904.1309	2154.2161	0.0000	1.0000
K2PIG	58	-908.0618	2156.2142	0.0000	1.0000
F81IG	60	-904.0334	2159.8847	0.0000	1.0000
K2P	56	-918.0195	2164.4024	0.0000	1.0000
F81	58	-913.4024	2166.8955	0.0000	1.0000
SYMG	61	-907.0783	2171.8381	0.0000	1.0000
SYMI	61	-907.1217	2171.9249	0.0000	1.0000
SYMIG	62	-907.0394	2177.6239	0.0000	1.0000
SYM	60	-917.1343	2186.0865	0.0000	1.0000
JC69G	56	-929.4436	2187.2506	0.0000	1.0000
JC69I	56	-929.4904	2187.3442	0.0000	1.0000
JC69IG	57	-929.4091	2193.0451	0.0000	1.0000
JC69	55	-938.9760	2200.4516	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/chr10_4414.phylip = 1 - 352;
 Partition Dummy = 1:183-phylip/chr10_4414.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 (HKYG)]

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

