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

Input data from file "683-tomodel-phylip/chr8_8088.phylip" (7 taxa, 1055 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	-3523.3878	7086.7755	0.6603	0.6603
GTRIG	21	-3523.3820	7088.7639	0.2443	0.9047
GTRI	20	-3525.3261	7090.6522	0.0951	0.9997
HKYG	16	-3535.6366	7103.2731	0.0002	0.9999
HKYIG	17	-3535.6568	7105.3137	0.0001	1.0000
HKYI	16	-3537.0029	7106.0057	0.0000	1.0000
SYMG	17	-3555.5507	7145.1015	0.0000	1.0000
SYMIG	18	-3555.5799	7147.1599	0.0000	1.0000
SYMI	17	-3556.6395	7147.2790	0.0000	1.0000
K2PG	13	-3561.2881	7148.5762	0.0000	1.0000
K2PIG	14	-3561.3230	7150.6461	0.0000	1.0000
K2PI	13	-3562.4428	7150.8855	0.0000	1.0000
GTR	19	-3566.6512	7171.3024	0.0000	1.0000
HKY	15	-3579.4636	7188.9272	0.0000	1.0000
SYM	16	-3599.6034	7231.2067	0.0000	1.0000
K2P	12	-3605.4597	7234.9193	0.0000	1.0000
F81G	15	-3642.4753	7314.9506	0.0000	1.0000
F81I	15	-3643.1166	7316.2333	0.0000	1.0000
F81IG	16	-3642.5332	7317.0664	0.0000	1.0000
JC69G	12	-3661.9574	7347.9147	0.0000	1.0000
JC69I	12	-3662.5523	7349.1045	0.0000	1.0000
JC69IG	13	-3661.9854	7349.9709	0.0000	1.0000
F81	14	-3681.7963	7391.5925	0.0000	1.0000
JC69	11	-3701.7592	7425.5183	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-3535.6366	7182.6539	0.6642	0.6642
HKYI	16	-3537.0029	7185.3865	0.1694	0.8336
GTRG	20	-3523.3878	7186.0015	0.1246	0.9582
HKYIG	17	-3535.6568	7189.6557	0.0200	0.9782
GTRI	20	-3525.3261	7189.8781	0.0179	0.9961
GTRIG	21	-3523.3820	7192.9511	0.0039	1.0000
K2PG	13	-3561.2881	7213.0730	0.0000	1.0000
K2PI	13	-3562.4428	7215.3823	0.0000	1.0000
K2PIG	14	-3561.3230	7220.1042	0.0000	1.0000
SYMG	17	-3555.5507	7229.4435	0.0000	1.0000
SYMI	17	-3556.6395	7231.6210	0.0000	1.0000
SYMIG	18	-3555.5799	7236.4632	0.0000	1.0000
HKY	15	-3579.4636	7263.3467	0.0000	1.0000
GTR	19	-3566.6512	7265.5671	0.0000	1.0000
K2P	12	-3605.4597	7294.4549	0.0000	1.0000
SYM	16	-3599.6034	7310.5875	0.0000	1.0000
F81G	15	-3642.4753	7389.3700	0.0000	1.0000
F81I	15	-3643.1166	7390.6527	0.0000	1.0000
F81IG	16	-3642.5332	7396.4471	0.0000	1.0000
JC69G	12	-3661.9574	7407.4503	0.0000	1.0000
JC69I	12	-3662.5523	7408.6401	0.0000	1.0000
JC69IG	13	-3661.9854	7414.4677	0.0000	1.0000
F81	14	-3681.7963	7461.0507	0.0000	1.0000
JC69	11	-3701.7592	7480.0926	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/chr8_8088.phylip = 1 - 1055;
 Partition Dummy = 1:683-tomodel-phylip/chr8_8088.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/chr8_8088.phylip = 1 - 1055;
 Partition Dummy = 1:683-tomodel-phylip/chr8_8088.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/chr8_8088.phylip = 1 - 1055;
 Partition Dummy = 1:683-tomodel-phylip/chr8_8088.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

