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

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

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


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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-2719.2789	5470.5579	0.5244	0.5244
HKYIG	17	-2719.1812	5472.3625	0.2127	0.7371
HKYI	16	-2720.4328	5472.8656	0.1654	0.9025
GTRG	20	-2717.5011	5475.0022	0.0568	0.9594
GTRIG	21	-2717.4117	5476.8234	0.0229	0.9822
GTRI	20	-2718.6641	5477.3282	0.0178	1.0000
HKY	15	-2742.2977	5514.5954	0.0000	1.0000
GTR	19	-2741.0178	5520.0355	0.0000	1.0000
K2PG	13	-2755.0601	5536.1201	0.0000	1.0000
K2PIG	14	-2754.9875	5537.9751	0.0000	1.0000
K2PI	13	-2756.0609	5538.1218	0.0000	1.0000
SYMG	17	-2753.8272	5541.6544	0.0000	1.0000
SYMI	17	-2754.7087	5543.4174	0.0000	1.0000
SYMIG	18	-2753.7595	5543.5190	0.0000	1.0000
K2P	12	-2778.4049	5580.8097	0.0000	1.0000
F81G	15	-2776.9758	5583.9515	0.0000	1.0000
F81I	15	-2777.6411	5585.2822	0.0000	1.0000
F81IG	16	-2776.9622	5585.9243	0.0000	1.0000
SYM	16	-2777.5677	5587.1354	0.0000	1.0000
F81	14	-2797.8367	5623.6735	0.0000	1.0000
JC69G	12	-2807.1956	5638.3912	0.0000	1.0000
JC69I	12	-2807.7843	5639.5687	0.0000	1.0000
JC69IG	13	-2807.1882	5640.3764	0.0000	1.0000
JC69	11	-2828.3399	5678.6799	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2719.2789	5548.9373	0.7405	0.7405
HKYI	16	-2720.4328	5551.2451	0.2336	0.9741
HKYIG	17	-2719.1812	5555.6406	0.0259	1.0000
GTRG	20	-2717.5011	5572.9765	0.0000	1.0000
GTRI	20	-2718.6641	5575.3025	0.0000	1.0000
GTRIG	21	-2717.4117	5579.6964	0.0000	1.0000
HKY	15	-2742.2977	5588.0761	0.0000	1.0000
K2PG	13	-2755.0601	5599.8034	0.0000	1.0000
K2PI	13	-2756.0609	5601.8050	0.0000	1.0000
K2PIG	14	-2754.9875	5606.5571	0.0000	1.0000
GTR	19	-2741.0178	5613.1111	0.0000	1.0000
SYMG	17	-2753.8272	5624.9326	0.0000	1.0000
SYMI	17	-2754.7087	5626.6956	0.0000	1.0000
SYMIG	18	-2753.7595	5631.6959	0.0000	1.0000
K2P	12	-2778.4049	5639.5943	0.0000	1.0000
F81G	15	-2776.9758	5657.4323	0.0000	1.0000
F81I	15	-2777.6411	5658.7629	0.0000	1.0000
F81IG	16	-2776.9622	5664.3038	0.0000	1.0000
SYM	16	-2777.5677	5665.5149	0.0000	1.0000
F81	14	-2797.8367	5692.2555	0.0000	1.0000
JC69G	12	-2807.1956	5697.1758	0.0000	1.0000
JC69I	12	-2807.7843	5698.3533	0.0000	1.0000
JC69IG	13	-2807.1882	5704.0597	0.0000	1.0000
JC69	11	-2828.3399	5732.5657	0.0000	1.0000



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

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

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


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

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


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

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

