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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
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
GTRG	20	-1754.6319	3550.6296	0.5407	0.5407
GTRIG	21	-1754.2985	3552.1019	0.2590	0.7997
GTRI	20	-1755.9135	3553.1929	0.1501	0.9498
HKYG	16	-1761.7992	3556.4772	0.0291	0.9789
HKYIG	17	-1761.5459	3558.0821	0.0130	0.9919
HKYI	16	-1763.0753	3559.0295	0.0081	1.0000
K2PG	13	-1779.0062	3584.5976	0.0000	1.0000
GTR	19	-1773.1505	3585.5348	0.0000	1.0000
K2PIG	14	-1778.8076	3586.2914	0.0000	1.0000
K2PI	13	-1780.0383	3586.6617	0.0000	1.0000
SYMG	17	-1775.9847	3586.9596	0.0000	1.0000
SYMIG	18	-1775.8022	3588.7130	0.0000	1.0000
SYMI	17	-1776.9259	3588.8420	0.0000	1.0000
HKY	15	-1779.1587	3589.0917	0.0000	1.0000
F81G	15	-1780.9783	3592.7308	0.0000	1.0000
F81IG	16	-1780.8242	3594.5273	0.0000	1.0000
F81I	15	-1782.0896	3594.9534	0.0000	1.0000
K2P	12	-1795.4099	3615.3205	0.0000	1.0000
JC69G	12	-1795.4909	3615.4826	0.0000	1.0000
JC69IG	13	-1795.3554	3617.2960	0.0000	1.0000
JC69I	12	-1796.4807	3617.4621	0.0000	1.0000
SYM	16	-1792.3289	3617.5367	0.0000	1.0000
F81	14	-1796.5625	3621.8014	0.0000	1.0000
JC69	11	-1810.7952	3644.0134	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-1754.6319	3549.2638	0.5362	0.5362
GTRIG	21	-1754.2985	3550.5970	0.2753	0.8116
GTRI	20	-1755.9135	3551.8270	0.1489	0.9604
HKYG	16	-1761.7992	3555.5984	0.0226	0.9830
HKYIG	17	-1761.5459	3557.0918	0.0107	0.9937
HKYI	16	-1763.0753	3558.1507	0.0063	1.0000
K2PG	13	-1779.0062	3584.0124	0.0000	1.0000
GTR	19	-1773.1505	3584.3010	0.0000	1.0000
K2PIG	14	-1778.8076	3585.6151	0.0000	1.0000
SYMG	17	-1775.9847	3585.9693	0.0000	1.0000
K2PI	13	-1780.0383	3586.0765	0.0000	1.0000
SYMIG	18	-1775.8022	3587.6044	0.0000	1.0000
SYMI	17	-1776.9259	3587.8517	0.0000	1.0000
HKY	15	-1779.1587	3588.3175	0.0000	1.0000
F81G	15	-1780.9783	3591.9566	0.0000	1.0000
F81IG	16	-1780.8242	3593.6485	0.0000	1.0000
F81I	15	-1782.0896	3594.1792	0.0000	1.0000
K2P	12	-1795.4099	3614.8197	0.0000	1.0000
JC69G	12	-1795.4909	3614.9818	0.0000	1.0000
SYM	16	-1792.3289	3616.6578	0.0000	1.0000
JC69IG	13	-1795.3554	3616.7108	0.0000	1.0000
JC69I	12	-1796.4807	3616.9613	0.0000	1.0000
F81	14	-1796.5625	3621.1250	0.0000	1.0000
JC69	11	-1810.7952	3643.5904	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1761.7992	3626.8816	0.7489	0.7489
HKYI	16	-1763.0753	3629.4338	0.2090	0.9580
HKYIG	17	-1761.5459	3632.8302	0.0383	0.9962
GTRG	20	-1754.6319	3638.3677	0.0024	0.9986
GTRI	20	-1755.9135	3640.9310	0.0007	0.9993
K2PG	13	-1779.0062	3641.9300	0.0004	0.9997
K2PI	13	-1780.0383	3643.9941	0.0001	0.9998
GTRIG	21	-1754.2985	3644.1561	0.0001	1.0000
K2PIG	14	-1778.8076	3647.9879	0.0000	1.0000
HKY	15	-1779.1587	3655.1454	0.0000	1.0000
F81G	15	-1780.9783	3658.7846	0.0000	1.0000
F81I	15	-1782.0896	3661.0072	0.0000	1.0000
SYMG	17	-1775.9847	3661.7077	0.0000	1.0000
SYMI	17	-1776.9259	3663.5901	0.0000	1.0000
F81IG	16	-1780.8242	3664.9317	0.0000	1.0000
SYMIG	18	-1775.8022	3667.7980	0.0000	1.0000
K2P	12	-1795.4099	3668.2821	0.0000	1.0000
JC69G	12	-1795.4909	3668.4442	0.0000	1.0000
GTR	19	-1773.1505	3668.9498	0.0000	1.0000
JC69I	12	-1796.4807	3670.4237	0.0000	1.0000
JC69IG	13	-1795.3554	3674.6283	0.0000	1.0000
F81	14	-1796.5625	3683.4978	0.0000	1.0000
SYM	16	-1792.3289	3687.9410	0.0000	1.0000
JC69	11	-1810.7952	3692.5975	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/chr11_3253.phylip = 1 - 636;
 Partition Dummy = 1:591-tomodel-phylip/chr11_3253.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/chr11_3253.phylip = 1 - 636;
 Partition Dummy = 1:591-tomodel-phylip/chr11_3253.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/chr11_3253.phylip = 1 - 636;
 Partition Dummy = 1:591-tomodel-phylip/chr11_3253.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

