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

Input data from file "591-tomodel-phylip/chr2_22024.phylip" (7 taxa, 759 characters)

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: GTRG


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
GTRG	20	-1944.1230	3929.3841	0.5535	0.5535
GTRIG	21	-1943.8601	3930.9740	0.2500	0.8035
GTRI	20	-1945.1586	3931.4554	0.1965	1.0000
HKYG	16	-1960.4972	3953.7276	0.0000	1.0000
HKYIG	17	-1959.9383	3954.7026	0.0000	1.0000
HKYI	16	-1961.7059	3956.1450	0.0000	1.0000
GTR	19	-1970.1936	3979.4156	0.0000	1.0000
SYMG	17	-1974.4527	3983.7313	0.0000	1.0000
SYMIG	18	-1974.1218	3985.1679	0.0000	1.0000
SYMI	17	-1975.3670	3985.5599	0.0000	1.0000
K2PG	13	-1985.0175	3996.5235	0.0000	1.0000
K2PIG	14	-1984.4618	3997.4881	0.0000	1.0000
K2PI	13	-1986.1691	3998.8269	0.0000	1.0000
HKY	15	-1985.2343	4001.1146	0.0000	1.0000
SYM	16	-1998.8279	4030.3889	0.0000	1.0000
F81G	15	-2000.9852	4032.6164	0.0000	1.0000
F81IG	16	-2000.5844	4033.9020	0.0000	1.0000
F81I	15	-2002.1093	4034.8647	0.0000	1.0000
K2P	12	-2009.4713	4043.3609	0.0000	1.0000
JC69G	12	-2021.3286	4067.0755	0.0000	1.0000
JC69IG	13	-2020.9341	4068.3568	0.0000	1.0000
JC69I	12	-2022.4129	4069.2440	0.0000	1.0000
F81	14	-2023.7614	4076.0874	0.0000	1.0000
JC69	11	-2044.1471	4110.6477	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-1944.1230	3928.2459	0.5454	0.5454
GTRIG	21	-1943.8601	3929.7202	0.2610	0.8064
GTRI	20	-1945.1586	3930.3171	0.1936	1.0000
HKYG	16	-1960.4972	3952.9945	0.0000	1.0000
HKYIG	17	-1959.9383	3953.8767	0.0000	1.0000
HKYI	16	-1961.7059	3955.4118	0.0000	1.0000
GTR	19	-1970.1936	3978.3872	0.0000	1.0000
SYMG	17	-1974.4527	3982.9053	0.0000	1.0000
SYMIG	18	-1974.1218	3984.2436	0.0000	1.0000
SYMI	17	-1975.3670	3984.7340	0.0000	1.0000
K2PG	13	-1985.0175	3996.0349	0.0000	1.0000
K2PIG	14	-1984.4618	3996.9236	0.0000	1.0000
K2PI	13	-1986.1691	3998.3383	0.0000	1.0000
HKY	15	-1985.2343	4000.4686	0.0000	1.0000
SYM	16	-1998.8279	4029.6558	0.0000	1.0000
F81G	15	-2000.9852	4031.9703	0.0000	1.0000
F81IG	16	-2000.5844	4033.1688	0.0000	1.0000
F81I	15	-2002.1093	4034.2187	0.0000	1.0000
K2P	12	-2009.4713	4042.9427	0.0000	1.0000
JC69G	12	-2021.3286	4066.6572	0.0000	1.0000
JC69IG	13	-2020.9341	4067.8682	0.0000	1.0000
JC69I	12	-2022.4129	4068.8257	0.0000	1.0000
F81	14	-2023.7614	4075.5229	0.0000	1.0000
JC69	11	-2044.1471	4110.2942	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	20	-1944.1230	4020.8860	0.6835	0.6835
GTRI	20	-1945.1586	4022.9572	0.2427	0.9262
GTRIG	21	-1943.8601	4026.9923	0.0323	0.9585
HKYG	16	-1960.4972	4027.1065	0.0305	0.9890
HKYI	16	-1961.7059	4029.5238	0.0091	0.9981
HKYIG	17	-1959.9383	4032.6207	0.0019	1.0000
K2PG	13	-1985.0175	4056.2509	0.0000	1.0000
K2PI	13	-1986.1691	4058.5543	0.0000	1.0000
SYMG	17	-1974.4527	4061.6494	0.0000	1.0000
K2PIG	14	-1984.4618	4061.7716	0.0000	1.0000
SYMI	17	-1975.3670	4063.4780	0.0000	1.0000
GTR	19	-1970.1936	4066.3952	0.0000	1.0000
SYMIG	18	-1974.1218	4067.6196	0.0000	1.0000
HKY	15	-1985.2343	4069.9486	0.0000	1.0000
K2P	12	-2009.4713	4098.5267	0.0000	1.0000
F81G	15	-2000.9852	4101.4504	0.0000	1.0000
F81I	15	-2002.1093	4103.6987	0.0000	1.0000
SYM	16	-1998.8279	4103.7678	0.0000	1.0000
F81IG	16	-2000.5844	4107.2809	0.0000	1.0000
JC69G	12	-2021.3286	4122.2413	0.0000	1.0000
JC69I	12	-2022.4129	4124.4097	0.0000	1.0000
JC69IG	13	-2020.9341	4128.0842	0.0000	1.0000
F81	14	-2023.7614	4140.3709	0.0000	1.0000
JC69	11	-2044.1471	4161.2463	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 591-tomodel-phylip/chr2_22024.phylip = 1 - 759;
 Partition Dummy = 1:591-tomodel-phylip/chr2_22024.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 (GTRG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr2_22024.phylip = 1 - 759;
 Partition Dummy = 1:591-tomodel-phylip/chr2_22024.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 BIC model (GTRG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr2_22024.phylip = 1 - 759;
 Partition Dummy = 1:591-tomodel-phylip/chr2_22024.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;


-------------------------------------------------------------
End of Output

