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

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

Minimum AIC  model: GTRI
Minimum AICc model: GTRI
Minimum BIC  model: HKYI


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-1998.5867	4037.1733	0.2812	0.2812
GTRG	20	-1998.8811	4037.7622	0.2095	0.4907
HKYI	16	-2002.9490	4037.8979	0.1958	0.6865
HKYG	16	-2003.2126	4038.4253	0.1504	0.8369
GTRIG	21	-1998.8885	4039.7770	0.0765	0.9134
HKYIG	17	-2003.2101	4040.4202	0.0555	0.9688
HKY	15	-2006.4639	4042.9277	0.0158	0.9847
GTR	19	-2002.4965	4042.9931	0.0153	1.0000
SYMI	17	-2014.7994	4063.5987	0.0000	1.0000
K2PI	13	-2018.9952	4063.9904	0.0000	1.0000
SYMG	17	-2015.1213	4064.2427	0.0000	1.0000
K2PG	13	-2019.2805	4064.5610	0.0000	1.0000
SYMIG	18	-2014.8098	4065.6197	0.0000	1.0000
K2PIG	14	-2019.0049	4066.0098	0.0000	1.0000
SYM	16	-2017.7215	4067.4430	0.0000	1.0000
K2P	12	-2021.8982	4067.7964	0.0000	1.0000
F81I	15	-2026.7177	4083.4353	0.0000	1.0000
F81G	15	-2027.0011	4084.0022	0.0000	1.0000
F81IG	16	-2026.7271	4085.4542	0.0000	1.0000
F81	14	-2029.6717	4087.3435	0.0000	1.0000
JC69I	12	-2041.1138	4106.2277	0.0000	1.0000
JC69G	12	-2041.4071	4106.8142	0.0000	1.0000
JC69IG	13	-2041.1238	4108.2475	0.0000	1.0000
JC69	11	-2043.7811	4109.5621	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-2002.9490	4114.0831	0.3734	0.3734
HKY	15	-2006.4639	4114.3513	0.3265	0.6998
HKYG	16	-2003.2126	4114.6104	0.2868	0.9867
HKYIG	17	-2003.2101	4121.3669	0.0098	0.9964
K2P	12	-2021.8982	4124.9352	0.0016	0.9981
K2PI	13	-2018.9952	4125.8909	0.0010	0.9991
K2PG	13	-2019.2805	4126.4614	0.0008	0.9999
GTRI	20	-1998.5867	4132.4048	0.0000	0.9999
K2PIG	14	-2019.0049	4132.6719	0.0000	0.9999
GTRG	20	-1998.8811	4132.9936	0.0000	1.0000
GTR	19	-2002.4965	4133.4629	0.0000	1.0000
GTRIG	21	-1998.8885	4139.7701	0.0000	1.0000
SYM	16	-2017.7215	4143.6282	0.0000	1.0000
SYMI	17	-2014.7994	4144.5455	0.0000	1.0000
SYMG	17	-2015.1213	4145.1894	0.0000	1.0000
SYMIG	18	-2014.8098	4151.3280	0.0000	1.0000
F81	14	-2029.6717	4154.0055	0.0000	1.0000
F81I	15	-2026.7177	4154.8589	0.0000	1.0000
F81G	15	-2027.0011	4155.4258	0.0000	1.0000
F81IG	16	-2026.7271	4161.6394	0.0000	1.0000
JC69	11	-2043.7811	4161.9394	0.0000	1.0000
JC69I	12	-2041.1138	4163.3666	0.0000	1.0000
JC69G	12	-2041.4071	4163.9531	0.0000	1.0000
JC69IG	13	-2041.1238	4170.1480	0.0000	1.0000



To specify default values for best AIC, AICc, or BIC models in MrBayes:

[Mrbayes block for the best AIC model (GTRI)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr7_6238.phylip = 1 - 864;
 Partition Dummy = 1:591-tomodel-phylip/chr7_6238.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


[Mrbayes block for the best AICc model (GTRI)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr7_6238.phylip = 1 - 864;
 Partition Dummy = 1:591-tomodel-phylip/chr7_6238.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


[Mrbayes block for the best BIC model (HKYI)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr7_6238.phylip = 1 - 864;
 Partition Dummy = 1:591-tomodel-phylip/chr7_6238.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

