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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: K2PG


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-3733.7225	7507.4450	0.4494	0.4494
HKYG	16	-3738.2898	7508.5796	0.2549	0.7043
GTRIG	21	-3733.7226	7509.4453	0.1653	0.8696
HKYIG	17	-3738.2899	7510.5799	0.0937	0.9634
GTRI	20	-3736.9415	7513.8829	0.0180	0.9813
HKYI	16	-3741.3021	7514.6042	0.0125	0.9939
K2PG	13	-3745.7527	7517.5054	0.0029	0.9968
SYMG	17	-3742.5085	7519.0170	0.0014	0.9982
K2PIG	14	-3745.7528	7519.5056	0.0011	0.9993
SYMIG	18	-3742.5087	7521.0173	0.0005	0.9998
K2PI	13	-3748.6483	7523.2965	0.0002	0.9999
SYMI	17	-3745.6673	7525.3346	0.0001	1.0000
GTR	19	-3752.3954	7542.7908	0.0000	1.0000
HKY	15	-3757.7303	7545.4605	0.0000	1.0000
SYM	16	-3761.4166	7554.8332	0.0000	1.0000
K2P	12	-3765.5425	7555.0851	0.0000	1.0000
F81G	15	-3825.2515	7680.5030	0.0000	1.0000
F81IG	16	-3825.2515	7682.5031	0.0000	1.0000
F81I	15	-3827.5469	7685.0938	0.0000	1.0000
JC69G	12	-3836.3887	7696.7773	0.0000	1.0000
JC69IG	13	-3836.3887	7698.7775	0.0000	1.0000
JC69I	12	-3838.6843	7701.3686	0.0000	1.0000
F81	14	-3843.1907	7714.3815	0.0000	1.0000
JC69	11	-3853.8659	7729.7319	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	13	-3745.7527	7581.7534	0.8750	0.8750
K2PI	13	-3748.6483	7587.5446	0.0484	0.9233
HKYG	16	-3738.2898	7587.6541	0.0458	0.9691
K2PIG	14	-3745.7528	7588.6958	0.0272	0.9963
HKYI	16	-3741.3021	7593.6787	0.0023	0.9986
HKYIG	17	-3738.2899	7594.5965	0.0014	1.0000
SYMG	17	-3742.5085	7603.0336	0.0000	1.0000
GTRG	20	-3733.7225	7606.2881	0.0000	1.0000
SYMI	17	-3745.6673	7609.3513	0.0000	1.0000
SYMIG	18	-3742.5087	7609.9761	0.0000	1.0000
GTRI	20	-3736.9415	7612.7261	0.0000	1.0000
GTRIG	21	-3733.7226	7613.2306	0.0000	1.0000
K2P	12	-3765.5425	7614.3909	0.0000	1.0000
HKY	15	-3757.7303	7619.5929	0.0000	1.0000
SYM	16	-3761.4166	7633.9077	0.0000	1.0000
GTR	19	-3752.3954	7636.6918	0.0000	1.0000
F81G	15	-3825.2515	7754.6354	0.0000	1.0000
JC69G	12	-3836.3887	7756.0832	0.0000	1.0000
F81I	15	-3827.5469	7759.2261	0.0000	1.0000
JC69I	12	-3838.6843	7760.6744	0.0000	1.0000
F81IG	16	-3825.2515	7761.5776	0.0000	1.0000
JC69IG	13	-3836.3887	7763.0255	0.0000	1.0000
F81	14	-3843.1907	7783.5717	0.0000	1.0000
JC69	11	-3853.8659	7784.0956	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 683-tomodel-phylip/chr5_3203.phylip = 1 - 1035;
 Partition Dummy = 1:683-tomodel-phylip/chr5_3203.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 683-tomodel-phylip/chr5_3203.phylip = 1 - 1035;
 Partition Dummy = 1:683-tomodel-phylip/chr5_3203.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 (K2PG)]

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr5_3203.phylip = 1 - 1035;
 Partition Dummy = 1:683-tomodel-phylip/chr5_3203.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

