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

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

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


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-1995.0880	4030.1759	0.3698	0.3698
GTRI	20	-1995.1976	4030.3952	0.3314	0.7012
GTRIG	21	-1995.0164	4032.0328	0.1461	0.8473
HKYG	16	-2000.8271	4033.6541	0.0650	0.9123
HKYI	16	-2000.8750	4033.7499	0.0619	0.9742
HKYIG	17	-2000.8063	4035.6125	0.0244	0.9986
GTR	19	-2001.8377	4041.6755	0.0012	0.9998
HKY	15	-2007.6525	4045.3051	0.0002	1.0000
SYMG	17	-2017.8507	4069.7013	0.0000	1.0000
SYMI	17	-2017.9497	4069.8995	0.0000	1.0000
SYMIG	18	-2017.7904	4071.5808	0.0000	1.0000
K2PG	13	-2023.1294	4072.2587	0.0000	1.0000
K2PI	13	-2023.1742	4072.3484	0.0000	1.0000
K2PIG	14	-2023.1074	4074.2148	0.0000	1.0000
SYM	16	-2024.5856	4081.1712	0.0000	1.0000
K2P	12	-2030.0750	4084.1499	0.0000	1.0000
F81G	15	-2031.6726	4093.3453	0.0000	1.0000
F81I	15	-2031.6856	4093.3713	0.0000	1.0000
F81IG	16	-2031.6746	4095.3491	0.0000	1.0000
F81	14	-2037.9493	4103.8986	0.0000	1.0000
JC69G	12	-2052.3692	4128.7383	0.0000	1.0000
JC69I	12	-2052.3841	4128.7681	0.0000	1.0000
JC69IG	13	-2052.3693	4130.7385	0.0000	1.0000
JC69	11	-2058.7473	4139.4946	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2000.8271	4112.0174	0.4949	0.4949
HKYI	16	-2000.8750	4112.1132	0.4718	0.9667
HKY	15	-2007.6525	4118.7706	0.0169	0.9836
HKYIG	17	-2000.8063	4118.8735	0.0161	0.9997
GTRG	20	-1995.0880	4128.1300	0.0002	0.9998
GTRI	20	-1995.1976	4128.3493	0.0001	1.0000
GTR	19	-2001.8377	4134.7319	0.0000	1.0000
GTRIG	21	-1995.0164	4134.8846	0.0000	1.0000
K2PG	13	-2023.1294	4135.9289	0.0000	1.0000
K2PI	13	-2023.1742	4136.0185	0.0000	1.0000
K2PIG	14	-2023.1074	4142.7827	0.0000	1.0000
K2P	12	-2030.0750	4142.9224	0.0000	1.0000
SYMG	17	-2017.8507	4152.9623	0.0000	1.0000
SYMI	17	-2017.9497	4153.1604	0.0000	1.0000
SYM	16	-2024.5856	4159.5344	0.0000	1.0000
SYMIG	18	-2017.7904	4159.7395	0.0000	1.0000
F81G	15	-2031.6726	4166.8109	0.0000	1.0000
F81I	15	-2031.6856	4166.8369	0.0000	1.0000
F81	14	-2037.9493	4172.4665	0.0000	1.0000
F81IG	16	-2031.6746	4173.7124	0.0000	1.0000
JC69G	12	-2052.3692	4187.5108	0.0000	1.0000
JC69I	12	-2052.3841	4187.5406	0.0000	1.0000
JC69	11	-2058.7473	4193.3693	0.0000	1.0000
JC69IG	13	-2052.3693	4194.4087	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/chr1_15412.phylip = 1 - 990;
 Partition Dummy = 1:591-tomodel-phylip/chr1_15412.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/chr1_15412.phylip = 1 - 990;
 Partition Dummy = 1:591-tomodel-phylip/chr1_15412.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 (HKYG)]

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

