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

Input data from file "591-tomodel-phylip/chr6_3100.phylip" (7 taxa, 758 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	-1760.0895	3561.3187	0.4687	0.4687
GTRIG	21	-1759.5835	3562.4224	0.2699	0.7386
GTRI	20	-1760.6841	3562.5079	0.2586	0.9972
HKYG	16	-1770.1479	3573.0300	0.0013	0.9985
HKYI	16	-1770.6760	3574.0861	0.0008	0.9993
HKYIG	17	-1769.7837	3574.3944	0.0007	1.0000
GTR	19	-1774.4115	3587.8527	0.0000	1.0000
HKY	15	-1782.4702	3595.5872	0.0000	1.0000
SYMG	17	-1786.3683	3607.5636	0.0000	1.0000
SYMI	17	-1786.8246	3608.4762	0.0000	1.0000
SYMIG	18	-1786.0225	3608.9706	0.0000	1.0000
F81G	15	-1789.6434	3609.9336	0.0000	1.0000
K2PG	13	-1791.7507	3609.9907	0.0000	1.0000
F81I	15	-1790.0532	3610.7533	0.0000	1.0000
K2PI	13	-1792.1885	3610.8663	0.0000	1.0000
K2PIG	14	-1791.4608	3611.4868	0.0000	1.0000
F81IG	16	-1789.4073	3611.5487	0.0000	1.0000
F81	14	-1800.9955	3630.5562	0.0000	1.0000
SYM	16	-1798.9123	3630.5587	0.0000	1.0000
K2P	12	-1803.7729	3631.9646	0.0000	1.0000
JC69G	12	-1808.9794	3642.3777	0.0000	1.0000
JC69I	12	-1809.3292	3643.0773	0.0000	1.0000
JC69IG	13	-1808.7868	3644.0629	0.0000	1.0000
JC69	11	-1820.2784	3662.9106	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-1760.0895	3560.1790	0.4615	0.4615
GTRIG	21	-1759.5835	3561.1670	0.2816	0.7431
GTRI	20	-1760.6841	3561.3682	0.2546	0.9977
HKYG	16	-1770.1479	3572.2959	0.0011	0.9988
HKYI	16	-1770.6760	3573.3519	0.0006	0.9994
HKYIG	17	-1769.7837	3573.5674	0.0006	1.0000
GTR	19	-1774.4115	3586.8229	0.0000	1.0000
HKY	15	-1782.4702	3594.9403	0.0000	1.0000
SYMG	17	-1786.3683	3606.7366	0.0000	1.0000
SYMI	17	-1786.8246	3607.6492	0.0000	1.0000
SYMIG	18	-1786.0225	3608.0450	0.0000	1.0000
F81G	15	-1789.6434	3609.2867	0.0000	1.0000
K2PG	13	-1791.7507	3609.5014	0.0000	1.0000
F81I	15	-1790.0532	3610.1064	0.0000	1.0000
K2PI	13	-1792.1885	3610.3771	0.0000	1.0000
F81IG	16	-1789.4073	3610.8146	0.0000	1.0000
K2PIG	14	-1791.4608	3610.9215	0.0000	1.0000
SYM	16	-1798.9123	3629.8246	0.0000	1.0000
F81	14	-1800.9955	3629.9910	0.0000	1.0000
K2P	12	-1803.7729	3631.5458	0.0000	1.0000
JC69G	12	-1808.9794	3641.9589	0.0000	1.0000
JC69I	12	-1809.3292	3642.6585	0.0000	1.0000
JC69IG	13	-1808.7868	3643.5736	0.0000	1.0000
JC69	11	-1820.2784	3662.5567	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1770.1479	3646.3868	0.5856	0.5856
HKYI	16	-1770.6760	3647.4429	0.3454	0.9309
HKYIG	17	-1769.7837	3652.2890	0.0306	0.9616
GTRG	20	-1760.0895	3652.7926	0.0238	0.9854
GTRI	20	-1760.6841	3653.9818	0.0131	0.9985
GTRIG	21	-1759.5835	3658.4113	0.0014	0.9999
HKY	15	-1782.4702	3664.4006	0.0001	1.0000
K2PG	13	-1791.7507	3669.7003	0.0000	1.0000
K2PI	13	-1792.1885	3670.5759	0.0000	1.0000
GTR	19	-1774.4115	3674.8059	0.0000	1.0000
K2PIG	14	-1791.4608	3675.7511	0.0000	1.0000
F81G	15	-1789.6434	3678.7470	0.0000	1.0000
F81I	15	-1790.0532	3679.5666	0.0000	1.0000
F81IG	16	-1789.4073	3684.9055	0.0000	1.0000
SYMG	17	-1786.3683	3685.4582	0.0000	1.0000
SYMI	17	-1786.8246	3686.3708	0.0000	1.0000
K2P	12	-1803.7729	3687.1140	0.0000	1.0000
SYMIG	18	-1786.0225	3691.3973	0.0000	1.0000
F81	14	-1800.9955	3694.8205	0.0000	1.0000
JC69G	12	-1808.9794	3697.5271	0.0000	1.0000
JC69I	12	-1809.3292	3698.2267	0.0000	1.0000
JC69IG	13	-1808.7868	3703.7725	0.0000	1.0000
SYM	16	-1798.9123	3703.9155	0.0000	1.0000
JC69	11	-1820.2784	3713.4942	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/chr6_3100.phylip = 1 - 758;
 Partition Dummy = 1:591-tomodel-phylip/chr6_3100.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/chr6_3100.phylip = 1 - 758;
 Partition Dummy = 1:591-tomodel-phylip/chr6_3100.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/chr6_3100.phylip = 1 - 758;
 Partition Dummy = 1:591-tomodel-phylip/chr6_3100.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

