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

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

Minimum AIC  model: SYMI
Minimum AICc model: SYMI
Minimum BIC  model: SYMI


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

Model	df	lnL		AIC		wAIC	cwAIC
SYMI	17	-3754.9825	7543.9649	0.5095	0.5095
SYMG	17	-3755.9601	7545.9202	0.1917	0.7012
SYMIG	18	-3754.9807	7545.9614	0.1878	0.8890
GTRI	20	-3754.0645	7548.1291	0.0635	0.9525
GTRG	20	-3755.0401	7550.0803	0.0239	0.9765
GTRIG	21	-3754.0585	7550.1170	0.0235	1.0000
K2PI	13	-3773.1272	7572.2545	0.0000	1.0000
K2PG	13	-3773.3954	7572.7907	0.0000	1.0000
K2PIG	14	-3773.0009	7574.0018	0.0000	1.0000
HKYI	16	-3772.0294	7576.0587	0.0000	1.0000
HKYG	16	-3772.1797	7576.3595	0.0000	1.0000
HKYIG	17	-3771.8586	7577.7173	0.0000	1.0000
SYM	16	-3777.1366	7586.2733	0.0000	1.0000
GTR	19	-3776.4225	7590.8451	0.0000	1.0000
K2P	12	-3794.5159	7613.0318	0.0000	1.0000
HKY	15	-3793.4599	7616.9199	0.0000	1.0000
JC69I	12	-3828.4317	7680.8633	0.0000	1.0000
JC69G	12	-3829.0852	7682.1704	0.0000	1.0000
JC69IG	13	-3828.4282	7682.8564	0.0000	1.0000
F81I	15	-3826.6800	7683.3600	0.0000	1.0000
F81G	15	-3827.2331	7684.4662	0.0000	1.0000
F81IG	16	-3826.6627	7685.3255	0.0000	1.0000
JC69	11	-3846.7762	7715.5524	0.0000	1.0000
F81	14	-3845.0628	7718.1256	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
SYMI	17	-3754.9825	7628.9706	0.6968	0.6968
SYMG	17	-3755.9601	7630.9259	0.2621	0.9590
SYMIG	18	-3754.9807	7635.9674	0.0211	0.9801
K2PI	13	-3773.1272	7637.2588	0.0111	0.9911
K2PG	13	-3773.3954	7637.7951	0.0085	0.9996
K2PIG	14	-3773.0009	7644.0065	0.0004	0.9999
GTRI	20	-3754.0645	7648.1357	0.0000	1.0000
GTRG	20	-3755.0401	7650.0870	0.0000	1.0000
GTRIG	21	-3754.0585	7655.1240	0.0000	1.0000
HKYI	16	-3772.0294	7656.0641	0.0000	1.0000
HKYG	16	-3772.1797	7656.3648	0.0000	1.0000
HKYIG	17	-3771.8586	7662.7230	0.0000	1.0000
SYM	16	-3777.1366	7666.2786	0.0000	1.0000
K2P	12	-3794.5159	7673.0358	0.0000	1.0000
GTR	19	-3776.4225	7685.8515	0.0000	1.0000
HKY	15	-3793.4599	7691.9249	0.0000	1.0000
JC69I	12	-3828.4317	7740.8673	0.0000	1.0000
JC69G	12	-3829.0852	7742.1744	0.0000	1.0000
JC69IG	13	-3828.4282	7747.8607	0.0000	1.0000
F81I	15	-3826.6800	7758.3651	0.0000	1.0000
F81G	15	-3827.2331	7759.4712	0.0000	1.0000
F81IG	16	-3826.6627	7765.3308	0.0000	1.0000
JC69	11	-3846.7762	7770.5561	0.0000	1.0000
F81	14	-3845.0628	7788.1303	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chrUn_random_3551.phylip = 1 - 1097;
 Partition Dummy = 1:591-tomodel-phylip/chrUn_random_3551.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=Fixed(Equal) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chrUn_random_3551.phylip = 1 - 1097;
 Partition Dummy = 1:591-tomodel-phylip/chrUn_random_3551.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=Fixed(Equal) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chrUn_random_3551.phylip = 1 - 1097;
 Partition Dummy = 1:591-tomodel-phylip/chrUn_random_3551.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=Fixed(Equal) pinvarpr=Uniform(0.0,1.0);
END;


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

