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

Input data from file "591-tomodel-phylip/chr5_11269.phylip" (7 taxa, 1025 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	-2613.8632	5259.7265	0.2220	0.2220
GTRI	20	-2609.9606	5259.9212	0.2014	0.4235
GTRG	20	-2609.9683	5259.9367	0.1999	0.6234
HKYI	16	-2613.9882	5259.9764	0.1960	0.8193
HKYIG	17	-2613.7603	5261.5205	0.0905	0.9099
GTRIG	21	-2609.7648	5261.5296	0.0901	1.0000
GTR	19	-2637.6735	5313.3469	0.0000	1.0000
HKY	15	-2641.7529	5313.5057	0.0000	1.0000
K2PG	13	-2650.1356	5326.2711	0.0000	1.0000
K2PI	13	-2650.5060	5327.0120	0.0000	1.0000
K2PIG	14	-2650.0591	5328.1181	0.0000	1.0000
SYMG	17	-2648.3805	5330.7610	0.0000	1.0000
SYMI	17	-2648.6245	5331.2490	0.0000	1.0000
SYMIG	18	-2648.3515	5332.7029	0.0000	1.0000
F81I	15	-2673.1929	5376.3857	0.0000	1.0000
F81G	15	-2673.2517	5376.5035	0.0000	1.0000
K2P	12	-2676.6907	5377.3813	0.0000	1.0000
F81IG	16	-2673.3174	5378.6348	0.0000	1.0000
SYM	16	-2674.1743	5380.3487	0.0000	1.0000
F81	14	-2697.6071	5423.2141	0.0000	1.0000
JC69G	12	-2703.3897	5430.7794	0.0000	1.0000
JC69I	12	-2703.4218	5430.8435	0.0000	1.0000
JC69IG	13	-2703.3983	5432.7966	0.0000	1.0000
JC69	11	-2727.4030	5476.8060	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2613.8632	5338.6456	0.5216	0.5216
HKYI	16	-2613.9882	5338.8955	0.4603	0.9819
HKYIG	17	-2613.7603	5345.3721	0.0181	1.0000
GTRI	20	-2609.9606	5358.5701	0.0000	1.0000
GTRG	20	-2609.9683	5358.5856	0.0000	1.0000
GTRIG	21	-2609.7648	5365.1110	0.0000	1.0000
HKY	15	-2641.7529	5387.4924	0.0000	1.0000
K2PG	13	-2650.1356	5390.3929	0.0000	1.0000
K2PI	13	-2650.5060	5391.1339	0.0000	1.0000
K2PIG	14	-2650.0591	5397.1724	0.0000	1.0000
GTR	19	-2637.6735	5407.0634	0.0000	1.0000
SYMG	17	-2648.3805	5414.6126	0.0000	1.0000
SYMI	17	-2648.6245	5415.1006	0.0000	1.0000
SYMIG	18	-2648.3515	5421.4870	0.0000	1.0000
K2P	12	-2676.6907	5436.5707	0.0000	1.0000
F81I	15	-2673.1929	5450.3725	0.0000	1.0000
F81G	15	-2673.2517	5450.4902	0.0000	1.0000
F81IG	16	-2673.3174	5457.5540	0.0000	1.0000
SYM	16	-2674.1743	5459.2678	0.0000	1.0000
JC69G	12	-2703.3897	5489.9688	0.0000	1.0000
JC69I	12	-2703.4218	5490.0329	0.0000	1.0000
F81	14	-2697.6071	5492.2684	0.0000	1.0000
JC69IG	13	-2703.3983	5496.9185	0.0000	1.0000
JC69	11	-2727.4030	5531.0629	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/chr5_11269.phylip = 1 - 1025;
 Partition Dummy = 1:591-tomodel-phylip/chr5_11269.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/chr5_11269.phylip = 1 - 1025;
 Partition Dummy = 1:591-tomodel-phylip/chr5_11269.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/chr5_11269.phylip = 1 - 1025;
 Partition Dummy = 1:591-tomodel-phylip/chr5_11269.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

