
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
16:18 04/04/2011

Input data from file "683-tomodel-phylip/chr7_10502.phylip" (7 taxa, 883 characters)

Minimum AIC  model: HKYI
Minimum AICc model: HKYI
Minimum BIC  model: HKY


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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	16	-1590.9658	3213.9316	0.3706	0.3706
HKYG	16	-1590.9674	3213.9347	0.3701	0.7407
HKYIG	17	-1590.9712	3215.9425	0.1356	0.8763
HKY	15	-1593.8914	3217.7827	0.0540	0.9304
GTRG	20	-1589.4820	3218.9640	0.0299	0.9603
GTRI	20	-1589.4830	3218.9660	0.0299	0.9902
GTRIG	21	-1590.0167	3222.0334	0.0065	0.9967
GTR	19	-1592.6745	3223.3491	0.0033	1.0000
F81G	15	-1607.7359	3245.4719	0.0000	1.0000
F81I	15	-1608.1099	3246.2197	0.0000	1.0000
F81IG	16	-1607.9311	3247.8622	0.0000	1.0000
F81	14	-1610.2964	3248.5927	0.0000	1.0000
K2PG	13	-1627.9963	3281.9927	0.0000	1.0000
K2PI	13	-1628.4822	3282.9645	0.0000	1.0000
K2PIG	14	-1628.3017	3284.6034	0.0000	1.0000
K2P	12	-1630.6972	3285.3944	0.0000	1.0000
SYMG	17	-1627.6586	3289.3172	0.0000	1.0000
SYMI	17	-1627.9306	3289.8612	0.0000	1.0000
SYMIG	18	-1627.7602	3291.5205	0.0000	1.0000
SYM	16	-1630.4237	3292.8475	0.0000	1.0000
JC69I	12	-1642.6289	3309.2579	0.0000	1.0000
JC69G	12	-1642.7887	3309.5773	0.0000	1.0000
JC69IG	13	-1642.4617	3310.9234	0.0000	1.0000
JC69	11	-1645.3299	3312.6598	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKY	15	-1593.8914	3289.5326	0.4396	0.4396
HKYI	16	-1590.9658	3290.4648	0.2758	0.7154
HKYG	16	-1590.9674	3290.4679	0.2754	0.9908
HKYIG	17	-1590.9712	3297.2590	0.0092	1.0000
GTR	19	-1592.6745	3314.2323	0.0000	1.0000
GTRG	20	-1589.4820	3314.6305	0.0000	1.0000
GTRI	20	-1589.4830	3314.6325	0.0000	1.0000
F81	14	-1610.2964	3315.5593	0.0000	1.0000
F81G	15	-1607.7359	3317.2217	0.0000	1.0000
F81I	15	-1608.1099	3317.9696	0.0000	1.0000
GTRIG	21	-1590.0167	3322.4832	0.0000	1.0000
F81IG	16	-1607.9311	3324.3954	0.0000	1.0000
K2P	12	-1630.6972	3342.7943	0.0000	1.0000
K2PG	13	-1627.9963	3344.1759	0.0000	1.0000
K2PI	13	-1628.4822	3345.1477	0.0000	1.0000
K2PIG	14	-1628.3017	3351.5699	0.0000	1.0000
JC69	11	-1645.3299	3365.2764	0.0000	1.0000
JC69I	12	-1642.6289	3366.6578	0.0000	1.0000
JC69G	12	-1642.7887	3366.9772	0.0000	1.0000
SYM	16	-1630.4237	3369.3807	0.0000	1.0000
SYMG	17	-1627.6586	3370.6338	0.0000	1.0000
SYMI	17	-1627.9306	3371.1777	0.0000	1.0000
JC69IG	13	-1642.4617	3373.1066	0.0000	1.0000
SYMIG	18	-1627.7602	3377.6203	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr7_10502.phylip = 1 - 883;
 Partition Dummy = 1:683-tomodel-phylip/chr7_10502.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;


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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr7_10502.phylip = 1 - 883;
 Partition Dummy = 1:683-tomodel-phylip/chr7_10502.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;


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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr7_10502.phylip = 1 - 883;
 Partition Dummy = 1:683-tomodel-phylip/chr7_10502.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0);
END;


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

