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

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

Minimum AIC  model: GTRI
Minimum AICc model: GTRI
Minimum BIC  model: HKYI


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-3661.4089	7362.8178	0.5886	0.5886
GTRG	20	-3662.0804	7364.1607	0.3007	0.8893
GTRIG	21	-3662.0805	7366.1611	0.1106	1.0000
HKYI	16	-3675.5464	7383.0927	0.0000	1.0000
HKYG	16	-3676.2728	7384.5456	0.0000	1.0000
HKYIG	17	-3676.2730	7386.5461	0.0000	1.0000
SYMI	17	-3721.4409	7476.8817	0.0000	1.0000
SYMG	17	-3722.0414	7478.0828	0.0000	1.0000
SYMIG	18	-3722.0415	7480.0831	0.0000	1.0000
K2PI	13	-3727.4533	7480.9065	0.0000	1.0000
K2PG	13	-3727.6640	7481.3280	0.0000	1.0000
K2PIG	14	-3727.6641	7483.3282	0.0000	1.0000
GTR	19	-3729.3215	7496.6431	0.0000	1.0000
F81I	15	-3741.9004	7513.8007	0.0000	1.0000
F81G	15	-3743.1608	7516.3215	0.0000	1.0000
F81IG	16	-3743.2974	7518.5948	0.0000	1.0000
HKY	15	-3745.8532	7521.7064	0.0000	1.0000
JC69I	12	-3784.3679	7592.7358	0.0000	1.0000
JC69G	12	-3785.0366	7594.0732	0.0000	1.0000
JC69IG	13	-3785.0367	7596.0735	0.0000	1.0000
SYM	16	-3788.7470	7609.4940	0.0000	1.0000
K2P	12	-3797.6464	7619.2928	0.0000	1.0000
F81	14	-3806.4658	7640.9316	0.0000	1.0000
JC69	11	-3849.7926	7721.5853	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-3675.5464	7464.5872	0.3390	0.3390
GTRI	20	-3661.4089	7464.6859	0.3227	0.6617
GTRG	20	-3662.0804	7466.0288	0.1649	0.8266
HKYG	16	-3676.2728	7466.0401	0.1640	0.9905
GTRIG	21	-3662.0805	7473.1226	0.0048	0.9953
HKYIG	17	-3676.2730	7473.1339	0.0047	1.0000
K2PI	13	-3727.4533	7547.1208	0.0000	1.0000
K2PG	13	-3727.6640	7547.5423	0.0000	1.0000
K2PIG	14	-3727.6641	7554.6359	0.0000	1.0000
SYMI	17	-3721.4409	7563.4696	0.0000	1.0000
SYMG	17	-3722.0414	7564.6706	0.0000	1.0000
SYMIG	18	-3722.0415	7571.7643	0.0000	1.0000
F81I	15	-3741.9004	7590.2018	0.0000	1.0000
F81G	15	-3743.1608	7592.7226	0.0000	1.0000
GTR	19	-3729.3215	7593.4178	0.0000	1.0000
HKY	15	-3745.8532	7598.1075	0.0000	1.0000
F81IG	16	-3743.2974	7600.0893	0.0000	1.0000
JC69I	12	-3784.3679	7653.8567	0.0000	1.0000
JC69G	12	-3785.0366	7655.1940	0.0000	1.0000
JC69IG	13	-3785.0367	7662.2877	0.0000	1.0000
K2P	12	-3797.6464	7680.4136	0.0000	1.0000
SYM	16	-3788.7470	7690.9885	0.0000	1.0000
F81	14	-3806.4658	7712.2392	0.0000	1.0000
JC69	11	-3849.7926	7777.6127	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr1_5213.phylip = 1 - 1204;
 Partition Dummy = 1:591-tomodel-phylip/chr1_5213.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr1_5213.phylip = 1 - 1204;
 Partition Dummy = 1:591-tomodel-phylip/chr1_5213.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr1_5213.phylip = 1 - 1204;
 Partition Dummy = 1:591-tomodel-phylip/chr1_5213.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

