
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
16:12 04/04/2011

Input data from file "683-tomodel-phylip/chr21_1559.phylip" (7 taxa, 913 characters)

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


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-2671.9203	5383.8405	0.6337	0.6337
GTRIG	21	-2671.6584	5385.3169	0.3029	0.9366
GTRI	20	-2674.2533	5388.5067	0.0615	0.9981
HKYG	16	-2682.2165	5396.4331	0.0012	0.9993
HKYIG	17	-2681.8674	5397.7348	0.0006	0.9999
HKYI	16	-2684.4759	5400.9517	0.0001	1.0000
SYMG	17	-2691.0574	5416.1147	0.0000	1.0000
SYMIG	18	-2690.7697	5417.5394	0.0000	1.0000
SYMI	17	-2693.3499	5420.6999	0.0000	1.0000
K2PG	13	-2698.2499	5422.4998	0.0000	1.0000
K2PIG	14	-2697.9247	5423.8495	0.0000	1.0000
K2PI	13	-2700.1379	5426.2759	0.0000	1.0000
GTR	19	-2713.5093	5465.0187	0.0000	1.0000
HKY	15	-2726.4136	5482.8271	0.0000	1.0000
SYM	16	-2733.7901	5499.5802	0.0000	1.0000
K2P	12	-2743.4868	5510.9736	0.0000	1.0000
F81G	15	-2770.6932	5571.3863	0.0000	1.0000
F81IG	16	-2770.2959	5572.5918	0.0000	1.0000
F81I	15	-2771.4410	5572.8819	0.0000	1.0000
JC69G	12	-2792.9337	5609.8675	0.0000	1.0000
JC69IG	13	-2792.8019	5611.6038	0.0000	1.0000
JC69I	12	-2794.2798	5612.5595	0.0000	1.0000
F81	14	-2809.9241	5647.8482	0.0000	1.0000
JC69	11	-2832.9530	5687.9060	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2682.2165	5473.5009	0.8365	0.8365
HKYI	16	-2684.4759	5478.0195	0.0873	0.9238
HKYIG	17	-2681.8674	5479.6193	0.0393	0.9631
GTRG	20	-2671.9203	5480.1752	0.0297	0.9928
GTRI	20	-2674.2533	5484.8414	0.0029	0.9957
K2PG	13	-2698.2499	5485.1174	0.0025	0.9982
GTRIG	21	-2671.6584	5486.4683	0.0013	0.9995
K2PI	13	-2700.1379	5488.8935	0.0004	0.9999
K2PIG	14	-2697.9247	5491.2838	0.0001	1.0000
SYMG	17	-2691.0574	5497.9992	0.0000	1.0000
SYMI	17	-2693.3499	5502.5844	0.0000	1.0000
SYMIG	18	-2690.7697	5504.2406	0.0000	1.0000
HKY	15	-2726.4136	5555.0782	0.0000	1.0000
GTR	19	-2713.5093	5556.5367	0.0000	1.0000
K2P	12	-2743.4868	5568.7745	0.0000	1.0000
SYM	16	-2733.7901	5576.6480	0.0000	1.0000
F81G	15	-2770.6932	5643.6374	0.0000	1.0000
F81I	15	-2771.4410	5645.1330	0.0000	1.0000
F81IG	16	-2770.2959	5649.6596	0.0000	1.0000
JC69G	12	-2792.9337	5667.6683	0.0000	1.0000
JC69I	12	-2794.2798	5670.3603	0.0000	1.0000
JC69IG	13	-2792.8019	5674.2213	0.0000	1.0000
F81	14	-2809.9241	5715.2825	0.0000	1.0000
JC69	11	-2832.9530	5740.8901	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr21_1559.phylip = 1 - 913;
 Partition Dummy = 1:683-tomodel-phylip/chr21_1559.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;


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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr21_1559.phylip = 1 - 913;
 Partition Dummy = 1:683-tomodel-phylip/chr21_1559.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;


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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr21_1559.phylip = 1 - 913;
 Partition Dummy = 1:683-tomodel-phylip/chr21_1559.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

