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

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

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


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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-3550.4318	7132.8636	0.6152	0.6152
HKYIG	17	-3550.4193	7134.8386	0.2292	0.8444
HKYI	16	-3552.5134	7137.0269	0.0767	0.9211
GTRG	20	-3548.8976	7137.7952	0.0523	0.9734
GTRIG	21	-3548.8781	7139.7562	0.0196	0.9930
GTRI	20	-3550.9053	7141.8106	0.0070	1.0000
HKY	15	-3572.1753	7174.3507	0.0000	1.0000
GTR	19	-3570.7899	7179.5797	0.0000	1.0000
K2PG	13	-3585.0417	7196.0834	0.0000	1.0000
K2PIG	14	-3584.8938	7197.7876	0.0000	1.0000
K2PI	13	-3586.7389	7199.4777	0.0000	1.0000
SYMG	17	-3582.9614	7199.9229	0.0000	1.0000
SYMIG	18	-3582.7647	7201.5293	0.0000	1.0000
SYMI	17	-3583.9947	7201.9894	0.0000	1.0000
K2P	12	-3607.7678	7239.5355	0.0000	1.0000
SYM	16	-3605.6856	7243.3711	0.0000	1.0000
F81G	15	-3641.7037	7313.4075	0.0000	1.0000
F81I	15	-3641.9229	7313.8457	0.0000	1.0000
F81IG	16	-3641.4063	7314.8125	0.0000	1.0000
F81	14	-3658.0048	7344.0095	0.0000	1.0000
JC69I	12	-3671.0396	7366.0792	0.0000	1.0000
JC69G	12	-3671.1472	7366.2945	0.0000	1.0000
JC69IG	13	-3670.6908	7367.3816	0.0000	1.0000
JC69	11	-3687.8462	7397.6924	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-3550.4318	7211.1133	0.8643	0.8643
HKYI	16	-3552.5134	7215.2766	0.1078	0.9721
HKYIG	17	-3550.4193	7217.9790	0.0279	1.0000
GTRG	20	-3548.8976	7235.6074	0.0000	1.0000
GTRI	20	-3550.9053	7239.6227	0.0000	1.0000
GTRIG	21	-3548.8781	7242.4590	0.0000	1.0000
HKY	15	-3572.1753	7247.7098	0.0000	1.0000
K2PG	13	-3585.0417	7259.6614	0.0000	1.0000
K2PI	13	-3586.7389	7263.0556	0.0000	1.0000
K2PIG	14	-3584.8938	7266.2561	0.0000	1.0000
GTR	19	-3570.7899	7272.5013	0.0000	1.0000
SYMG	17	-3582.9614	7283.0633	0.0000	1.0000
SYMI	17	-3583.9947	7285.1298	0.0000	1.0000
SYMIG	18	-3582.7647	7289.5603	0.0000	1.0000
K2P	12	-3607.7678	7298.2228	0.0000	1.0000
SYM	16	-3605.6856	7321.6209	0.0000	1.0000
F81G	15	-3641.7037	7386.7666	0.0000	1.0000
F81I	15	-3641.9229	7387.2049	0.0000	1.0000
F81IG	16	-3641.4063	7393.0623	0.0000	1.0000
F81	14	-3658.0048	7412.4781	0.0000	1.0000
JC69I	12	-3671.0396	7424.7665	0.0000	1.0000
JC69G	12	-3671.1472	7424.9818	0.0000	1.0000
JC69IG	13	-3670.6908	7430.9596	0.0000	1.0000
JC69	11	-3687.8462	7451.4891	0.0000	1.0000



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

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

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


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

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


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

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

