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

Input data from file "683-tomodel-phylip/chr3_281.phylip" (7 taxa, 1135 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	-3589.8375	7211.6749	0.4043	0.4043
GTRG	20	-3586.6387	7213.2773	0.1814	0.5857
HKYIG	17	-3589.8293	7213.6586	0.1499	0.7356
HKYI	16	-3591.0504	7214.1008	0.1202	0.8558
GTRI	20	-3587.5228	7215.0455	0.0749	0.9308
GTRIG	21	-3586.6024	7215.2048	0.0692	1.0000
HKY	15	-3620.3608	7270.7215	0.0000	1.0000
GTR	19	-3617.2422	7272.4843	0.0000	1.0000
K2PG	13	-3656.9745	7339.9489	0.0000	1.0000
K2PIG	14	-3656.9366	7341.8733	0.0000	1.0000
SYMG	17	-3654.2785	7342.5570	0.0000	1.0000
K2PI	13	-3658.5756	7343.1511	0.0000	1.0000
SYMIG	18	-3654.2083	7344.4166	0.0000	1.0000
SYMI	17	-3655.7346	7345.4693	0.0000	1.0000
K2P	12	-3688.3924	7400.7847	0.0000	1.0000
SYM	16	-3686.0159	7404.0318	0.0000	1.0000
F81G	15	-3712.2203	7454.4405	0.0000	1.0000
F81IG	16	-3712.2038	7456.4077	0.0000	1.0000
F81I	15	-3713.4230	7456.8460	0.0000	1.0000
F81	14	-3739.2537	7506.5075	0.0000	1.0000
JC69G	12	-3764.8923	7553.7845	0.0000	1.0000
JC69IG	13	-3764.8627	7555.7254	0.0000	1.0000
JC69I	12	-3766.3896	7556.7793	0.0000	1.0000
JC69	11	-3793.1306	7608.2613	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-3589.8375	7292.2251	0.7534	0.7534
HKYI	16	-3591.0504	7294.6510	0.2240	0.9774
HKYIG	17	-3589.8293	7299.2432	0.0225	1.0000
GTRG	20	-3586.6387	7313.9651	0.0000	1.0000
GTRI	20	-3587.5228	7315.7333	0.0000	1.0000
GTRIG	21	-3586.6024	7320.9269	0.0000	1.0000
HKY	15	-3620.3608	7346.2374	0.0000	1.0000
GTR	19	-3617.2422	7368.1377	0.0000	1.0000
K2PG	13	-3656.9745	7405.3960	0.0000	1.0000
K2PI	13	-3658.5756	7408.5981	0.0000	1.0000
K2PIG	14	-3656.9366	7412.3547	0.0000	1.0000
SYMG	17	-3654.2785	7428.1416	0.0000	1.0000
SYMI	17	-3655.7346	7431.0539	0.0000	1.0000
SYMIG	18	-3654.2083	7435.0356	0.0000	1.0000
K2P	12	-3688.3924	7461.1974	0.0000	1.0000
SYM	16	-3686.0159	7484.5820	0.0000	1.0000
F81G	15	-3712.2203	7529.9563	0.0000	1.0000
F81I	15	-3713.4230	7532.3619	0.0000	1.0000
F81IG	16	-3712.2038	7536.9579	0.0000	1.0000
F81	14	-3739.2537	7576.9889	0.0000	1.0000
JC69G	12	-3764.8923	7614.1972	0.0000	1.0000
JC69I	12	-3766.3896	7617.1919	0.0000	1.0000
JC69IG	13	-3764.8627	7621.1724	0.0000	1.0000
JC69	11	-3793.1306	7663.6395	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 683-tomodel-phylip/chr3_281.phylip = 1 - 1135;
 Partition Dummy = 1:683-tomodel-phylip/chr3_281.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 683-tomodel-phylip/chr3_281.phylip = 1 - 1135;
 Partition Dummy = 1:683-tomodel-phylip/chr3_281.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 683-tomodel-phylip/chr3_281.phylip = 1 - 1135;
 Partition Dummy = 1:683-tomodel-phylip/chr3_281.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

