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

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

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


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-3510.5850	7061.1701	0.3257	0.3257
GTRI	20	-3511.0605	7062.1210	0.2024	0.5281
HKYG	16	-3515.1913	7062.3826	0.1776	0.7058
GTRIG	21	-3510.5794	7063.1587	0.1205	0.8263
HKYI	16	-3515.6828	7063.3656	0.1087	0.9349
HKYIG	17	-3515.1953	7064.3906	0.0651	1.0000
GTR	19	-3524.5322	7087.0643	0.0000	1.0000
HKY	15	-3528.9826	7087.9653	0.0000	1.0000
K2PG	13	-3533.1209	7092.2418	0.0000	1.0000
SYMG	17	-3529.4304	7092.8607	0.0000	1.0000
K2PI	13	-3533.6793	7093.3585	0.0000	1.0000
SYMI	17	-3529.8995	7093.7990	0.0000	1.0000
K2PIG	14	-3533.1369	7094.2738	0.0000	1.0000
SYMIG	18	-3529.4293	7094.8586	0.0000	1.0000
K2P	12	-3546.8208	7117.6416	0.0000	1.0000
SYM	16	-3543.2319	7118.4639	0.0000	1.0000
F81G	15	-3621.6756	7273.3512	0.0000	1.0000
F81I	15	-3622.2596	7274.5192	0.0000	1.0000
F81IG	16	-3621.6791	7275.3582	0.0000	1.0000
F81	14	-3633.4181	7294.8361	0.0000	1.0000
JC69G	12	-3636.1429	7296.2859	0.0000	1.0000
JC69I	12	-3636.6822	7297.3643	0.0000	1.0000
JC69IG	13	-3636.1609	7298.3218	0.0000	1.0000
JC69	11	-3647.8574	7317.7147	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-3515.1913	7142.7342	0.6088	0.6088
HKYI	16	-3515.6828	7143.7173	0.3724	0.9812
HKYIG	17	-3515.1953	7149.7642	0.0181	0.9993
K2PG	13	-3533.1209	7157.5275	0.0004	0.9997
K2PI	13	-3533.6793	7158.6442	0.0002	0.9999
GTRG	20	-3510.5850	7161.6096	0.0000	0.9999
GTRI	20	-3511.0605	7162.5606	0.0000	1.0000
HKY	15	-3528.9826	7163.2949	0.0000	1.0000
K2PIG	14	-3533.1369	7164.5815	0.0000	1.0000
GTRIG	21	-3510.5794	7168.6202	0.0000	1.0000
K2P	12	-3546.8208	7177.9054	0.0000	1.0000
SYMG	17	-3529.4304	7178.2343	0.0000	1.0000
SYMI	17	-3529.8995	7179.1726	0.0000	1.0000
GTR	19	-3524.5322	7182.4819	0.0000	1.0000
SYMIG	18	-3529.4293	7185.2542	0.0000	1.0000
SYM	16	-3543.2319	7198.8155	0.0000	1.0000
F81G	15	-3621.6756	7348.6808	0.0000	1.0000
F81I	15	-3622.2596	7349.8488	0.0000	1.0000
F81IG	16	-3621.6791	7355.7098	0.0000	1.0000
JC69G	12	-3636.1429	7356.5496	0.0000	1.0000
JC69I	12	-3636.6822	7357.6280	0.0000	1.0000
JC69IG	13	-3636.1609	7363.6075	0.0000	1.0000
F81	14	-3633.4181	7365.1438	0.0000	1.0000
JC69	11	-3647.8574	7372.9565	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_10375.phylip = 1 - 1121;
 Partition Dummy = 1:683-tomodel-phylip/chr5_10375.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_10375.phylip = 1 - 1121;
 Partition Dummy = 1:683-tomodel-phylip/chr5_10375.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 683-tomodel-phylip/chr5_10375.phylip = 1 - 1121;
 Partition Dummy = 1:683-tomodel-phylip/chr5_10375.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

