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

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

Minimum AIC  model: GTRIG
Minimum AICc model: GTRIG
Minimum BIC  model: GTRG


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	21	-2687.1715	5416.3430	0.5151	0.5151
GTRG	20	-2688.2418	5416.4837	0.4801	0.9952
GTRI	20	-2692.8418	5425.6836	0.0048	1.0000
HKYIG	17	-2709.6787	5453.3574	0.0000	1.0000
HKYG	16	-2711.0706	5454.1412	0.0000	1.0000
HKYI	16	-2715.0065	5462.0130	0.0000	1.0000
SYMIG	18	-2713.6650	5463.3301	0.0000	1.0000
SYMG	17	-2714.8845	5463.7689	0.0000	1.0000
SYMI	17	-2718.1334	5470.2667	0.0000	1.0000
K2PIG	14	-2726.5368	5481.0737	0.0000	1.0000
K2PG	13	-2727.9064	5481.8128	0.0000	1.0000
K2PI	13	-2731.4008	5488.8016	0.0000	1.0000
GTR	19	-2745.9030	5529.8059	0.0000	1.0000
F81IG	16	-2765.5833	5563.1666	0.0000	1.0000
F81G	15	-2766.6066	5563.2133	0.0000	1.0000
F81I	15	-2770.1262	5570.2525	0.0000	1.0000
HKY	15	-2770.6806	5571.3611	0.0000	1.0000
SYM	16	-2772.1147	5576.2295	0.0000	1.0000
JC69IG	13	-2778.7178	5583.4356	0.0000	1.0000
JC69G	12	-2779.7498	5583.4997	0.0000	1.0000
JC69I	12	-2783.1244	5590.2488	0.0000	1.0000
K2P	12	-2788.2102	5600.4204	0.0000	1.0000
F81	14	-2822.9271	5673.8541	0.0000	1.0000
JC69	11	-2836.6557	5695.3115	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	20	-2688.2418	5512.1048	0.9022	0.9022
GTRIG	21	-2687.1715	5516.7452	0.0886	0.9908
GTRI	20	-2692.8418	5521.3048	0.0091	0.9999
HKYG	16	-2711.0706	5530.6381	0.0001	1.0000
HKYIG	17	-2709.6787	5534.6354	0.0000	1.0000
HKYI	16	-2715.0065	5538.5099	0.0000	1.0000
K2PG	13	-2727.9064	5543.9665	0.0000	1.0000
SYMG	17	-2714.8845	5545.0469	0.0000	1.0000
K2PIG	14	-2726.5368	5548.0085	0.0000	1.0000
SYMIG	18	-2713.6650	5549.3891	0.0000	1.0000
K2PI	13	-2731.4008	5550.9553	0.0000	1.0000
SYMI	17	-2718.1334	5551.5447	0.0000	1.0000
GTR	19	-2745.9030	5620.6460	0.0000	1.0000
F81G	15	-2766.6066	5634.9291	0.0000	1.0000
F81IG	16	-2765.5833	5639.6635	0.0000	1.0000
JC69G	12	-2779.7498	5640.8724	0.0000	1.0000
F81I	15	-2770.1262	5641.9683	0.0000	1.0000
HKY	15	-2770.6806	5643.0770	0.0000	1.0000
JC69IG	13	-2778.7178	5645.5894	0.0000	1.0000
JC69I	12	-2783.1244	5647.6215	0.0000	1.0000
SYM	16	-2772.1147	5652.7264	0.0000	1.0000
K2P	12	-2788.2102	5657.7931	0.0000	1.0000
F81	14	-2822.9271	5740.7889	0.0000	1.0000
JC69	11	-2836.6557	5747.9031	0.0000	1.0000



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

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

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


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

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


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr6_8822.phylip = 1 - 881;
 Partition Dummy = 1:591-tomodel-phylip/chr6_8822.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;


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

