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

Input data from file "183-phylip/chr5_2978.phylip" (29 taxa, 370 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	-915.0423	1973.7739	0.3140	0.3140
HKYI	60	-915.3270	1974.3432	0.2362	0.5501
GTRG	64	-909.8380	1974.9547	0.1740	0.7241
HKYIG	61	-914.6559	1975.8702	0.1101	0.8342
GTRI	64	-910.3584	1975.9955	0.1034	0.9375
GTRIG	65	-909.3908	1977.0053	0.0624	0.9999
GTR	63	-920.1158	1992.5845	0.0000	1.0000
HKY	59	-926.0223	1992.8834	0.0000	1.0000
F81G	59	-927.5893	1996.0172	0.0000	1.0000
F81I	59	-927.7591	1996.3568	0.0000	1.0000
F81IG	60	-927.3285	1998.3462	0.0000	1.0000
K2PG	57	-935.8876	2006.9675	0.0000	1.0000
K2PI	57	-936.0276	2007.2475	0.0000	1.0000
K2PIG	58	-935.5914	2009.1892	0.0000	1.0000
SYMG	61	-932.6580	2011.8744	0.0000	1.0000
SYMI	61	-932.9043	2012.3670	0.0000	1.0000
SYMIG	62	-932.3509	2014.1481	0.0000	1.0000
F81	58	-938.9707	2015.9479	0.0000	1.0000
JC69G	56	-946.7143	2025.8248	0.0000	1.0000
JC69I	56	-946.8311	2026.0583	0.0000	1.0000
K2P	56	-947.7872	2027.9706	0.0000	1.0000
JC69IG	57	-946.4881	2028.1685	0.0000	1.0000
SYM	60	-943.2898	2030.2688	0.0000	1.0000
JC69	55	-958.4859	2046.5897	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-909.8380	1947.6760	0.3499	0.3499
GTRI	64	-910.3584	1948.7168	0.2080	0.5579
GTRIG	65	-909.3908	1948.7816	0.2013	0.7592
HKYG	60	-915.0423	1950.0846	0.1050	0.8642
HKYI	60	-915.3270	1950.6539	0.0790	0.9431
HKYIG	61	-914.6559	1951.3117	0.0568	1.0000
GTR	63	-920.1158	1966.2316	0.0000	1.0000
HKY	59	-926.0223	1970.0447	0.0000	1.0000
F81G	59	-927.5893	1973.1785	0.0000	1.0000
F81I	59	-927.7591	1973.5181	0.0000	1.0000
F81IG	60	-927.3285	1974.6569	0.0000	1.0000
K2PG	57	-935.8876	1985.7752	0.0000	1.0000
K2PI	57	-936.0276	1986.0552	0.0000	1.0000
K2PIG	58	-935.5914	1987.1827	0.0000	1.0000
SYMG	61	-932.6580	1987.3160	0.0000	1.0000
SYMI	61	-932.9043	1987.8085	0.0000	1.0000
SYMIG	62	-932.3509	1988.7018	0.0000	1.0000
F81	58	-938.9707	1993.9415	0.0000	1.0000
JC69G	56	-946.7143	2005.4286	0.0000	1.0000
JC69I	56	-946.8311	2005.6621	0.0000	1.0000
SYM	60	-943.2898	2006.5795	0.0000	1.0000
JC69IG	57	-946.4881	2006.9762	0.0000	1.0000
K2P	56	-947.7872	2007.5744	0.0000	1.0000
JC69	55	-958.4859	2026.9719	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-915.0423	2184.8948	0.5460	0.5460
HKYI	60	-915.3270	2185.4641	0.4107	0.9568
HKYIG	61	-914.6559	2190.0354	0.0418	0.9985
GTRG	64	-909.8380	2198.1402	0.0007	0.9993
GTRI	64	-910.3584	2199.1810	0.0004	0.9997
HKY	59	-926.0223	2200.9414	0.0002	0.9999
GTRIG	65	-909.3908	2203.1593	0.0001	0.9999
F81G	59	-927.5893	2204.0752	0.0000	1.0000
F81I	59	-927.7591	2204.4148	0.0000	1.0000
K2PG	57	-935.8876	2208.8449	0.0000	1.0000
K2PI	57	-936.0276	2209.1249	0.0000	1.0000
F81IG	60	-927.3285	2209.4671	0.0000	1.0000
GTR	63	-920.1158	2212.7823	0.0000	1.0000
K2PIG	58	-935.5914	2214.1659	0.0000	1.0000
F81	58	-938.9707	2220.9247	0.0000	1.0000
JC69G	56	-946.7143	2224.5848	0.0000	1.0000
JC69I	56	-946.8311	2224.8183	0.0000	1.0000
SYMG	61	-932.6580	2226.0397	0.0000	1.0000
SYMI	61	-932.9043	2226.5322	0.0000	1.0000
K2P	56	-947.7872	2226.7306	0.0000	1.0000
JC69IG	57	-946.4881	2230.0459	0.0000	1.0000
SYMIG	62	-932.3509	2231.3390	0.0000	1.0000
SYM	60	-943.2898	2241.3897	0.0000	1.0000
JC69	55	-958.4859	2242.2145	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/chr5_2978.phylip = 1 - 370;
 Partition Dummy = 1:183-phylip/chr5_2978.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/chr5_2978.phylip = 1 - 370;
 Partition Dummy = 1:183-phylip/chr5_2978.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/chr5_2978.phylip = 1 - 370;
 Partition Dummy = 1:183-phylip/chr5_2978.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

