
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
21:56 04/04/2011

Input data from file "591-tomodel-phylip/chrZ_495.phylip" (7 taxa, 797 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	16	-1946.5659	3925.8292	0.3975	0.3975
HKYI	16	-1946.6380	3925.9733	0.3698	0.7673
HKYIG	17	-1946.5771	3927.9398	0.1383	0.9056
GTRG	20	-1944.6403	3930.3630	0.0412	0.9468
GTRI	20	-1944.6969	3930.4762	0.0389	0.9857
GTRIG	21	-1944.6465	3932.4854	0.0143	1.0000
HKY	15	-1968.8179	3968.2503	0.0000	1.0000
GTR	19	-1967.2108	3973.3996	0.0000	1.0000
F81G	15	-1985.3748	4001.3641	0.0000	1.0000
F81I	15	-1985.4846	4001.5838	0.0000	1.0000
F81IG	16	-1985.3741	4003.4457	0.0000	1.0000
K2PG	13	-1990.3573	4007.1794	0.0000	1.0000
K2PI	13	-1990.4516	4007.3681	0.0000	1.0000
SYMG	17	-1986.6731	4008.1318	0.0000	1.0000
SYMI	17	-1986.6963	4008.1781	0.0000	1.0000
K2PIG	14	-1990.3525	4009.2422	0.0000	1.0000
SYMIG	18	-1986.7199	4010.3189	0.0000	1.0000
F81	14	-2005.8468	4040.2307	0.0000	1.0000
K2P	12	-2012.3039	4049.0058	0.0000	1.0000
SYM	16	-2008.4211	4049.5396	0.0000	1.0000
JC69G	12	-2026.3850	4077.1679	0.0000	1.0000
JC69I	12	-2026.4543	4077.3067	0.0000	1.0000
JC69IG	13	-2026.3721	4079.2090	0.0000	1.0000
JC69	11	-2046.4984	4115.3331	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-1946.5659	3925.1317	0.3869	0.3869
HKYI	16	-1946.6380	3925.2759	0.3600	0.7469
HKYIG	17	-1946.5771	3927.1542	0.1408	0.8877
GTRG	20	-1944.6403	3929.2806	0.0486	0.9363
GTRI	20	-1944.6969	3929.3937	0.0459	0.9822
GTRIG	21	-1944.6465	3931.2931	0.0178	1.0000
HKY	15	-1968.8179	3967.6357	0.0000	1.0000
GTR	19	-1967.2108	3972.4215	0.0000	1.0000
F81G	15	-1985.3748	4000.7495	0.0000	1.0000
F81I	15	-1985.4846	4000.9692	0.0000	1.0000
F81IG	16	-1985.3741	4002.7482	0.0000	1.0000
K2PG	13	-1990.3573	4006.7146	0.0000	1.0000
K2PI	13	-1990.4516	4006.9032	0.0000	1.0000
SYMG	17	-1986.6731	4007.3462	0.0000	1.0000
SYMI	17	-1986.6963	4007.3925	0.0000	1.0000
K2PIG	14	-1990.3525	4008.7051	0.0000	1.0000
SYMIG	18	-1986.7199	4009.4398	0.0000	1.0000
F81	14	-2005.8468	4039.6936	0.0000	1.0000
K2P	12	-2012.3039	4048.6079	0.0000	1.0000
SYM	16	-2008.4211	4048.8421	0.0000	1.0000
JC69G	12	-2026.3850	4076.7699	0.0000	1.0000
JC69I	12	-2026.4543	4076.9087	0.0000	1.0000
JC69IG	13	-2026.3721	4078.7442	0.0000	1.0000
JC69	11	-2046.4984	4114.9968	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1946.5659	4000.0254	0.5088	0.5088
HKYI	16	-1946.6380	4000.1696	0.4734	0.9822
HKYIG	17	-1946.5771	4006.7287	0.0178	1.0000
GTRG	20	-1944.6403	4022.8977	0.0000	1.0000
GTRI	20	-1944.6969	4023.0108	0.0000	1.0000
GTRIG	21	-1944.6465	4029.5910	0.0000	1.0000
HKY	15	-1968.8179	4037.8485	0.0000	1.0000
GTR	19	-1967.2108	4061.3578	0.0000	1.0000
K2PG	13	-1990.3573	4067.5657	0.0000	1.0000
K2PI	13	-1990.4516	4067.7543	0.0000	1.0000
F81G	15	-1985.3748	4070.9623	0.0000	1.0000
F81I	15	-1985.4846	4071.1820	0.0000	1.0000
K2PIG	14	-1990.3525	4074.2370	0.0000	1.0000
F81IG	16	-1985.3741	4077.6419	0.0000	1.0000
SYMG	17	-1986.6731	4086.9207	0.0000	1.0000
SYMI	17	-1986.6963	4086.9670	0.0000	1.0000
SYMIG	18	-1986.7199	4093.6951	0.0000	1.0000
K2P	12	-2012.3039	4104.7781	0.0000	1.0000
F81	14	-2005.8468	4105.2256	0.0000	1.0000
SYM	16	-2008.4211	4123.7358	0.0000	1.0000
JC69G	12	-2026.3850	4132.9402	0.0000	1.0000
JC69I	12	-2026.4543	4133.0790	0.0000	1.0000
JC69IG	13	-2026.3721	4139.5953	0.0000	1.0000
JC69	11	-2046.4984	4166.4862	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 591-tomodel-phylip/chrZ_495.phylip = 1 - 797;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_495.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 591-tomodel-phylip/chrZ_495.phylip = 1 - 797;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_495.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 591-tomodel-phylip/chrZ_495.phylip = 1 - 797;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_495.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

