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

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

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


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-4728.8002	9497.6005	0.7319	0.7319
GTRIG	21	-4728.8098	9499.6197	0.2667	0.9985
GTRI	20	-4735.0166	9510.0332	0.0015	1.0000
HKYG	16	-4749.1580	9530.3161	0.0000	1.0000
HKYIG	17	-4749.1317	9532.2635	0.0000	1.0000
HKYI	16	-4752.4840	9536.9679	0.0000	1.0000
GTR	19	-4783.9973	9605.9947	0.0000	1.0000
SYMG	17	-4789.0519	9612.1038	0.0000	1.0000
SYMIG	18	-4789.0325	9614.0649	0.0000	1.0000
SYMI	17	-4792.3136	9618.6271	0.0000	1.0000
K2PG	13	-4796.8290	9619.6580	0.0000	1.0000
K2PIG	14	-4796.7258	9621.4515	0.0000	1.0000
K2PI	13	-4799.1895	9624.3791	0.0000	1.0000
HKY	15	-4801.6585	9633.3170	0.0000	1.0000
SYM	16	-4844.5073	9721.0145	0.0000	1.0000
K2P	12	-4851.3578	9726.7157	0.0000	1.0000
F81G	15	-4877.6448	9785.2896	0.0000	1.0000
F81IG	16	-4877.4336	9786.8672	0.0000	1.0000
F81I	15	-4879.0531	9788.1061	0.0000	1.0000
JC69G	12	-4911.8774	9847.7548	0.0000	1.0000
JC69IG	13	-4911.6152	9849.2305	0.0000	1.0000
JC69I	12	-4913.1714	9850.3427	0.0000	1.0000
F81	14	-4922.4417	9872.8835	0.0000	1.0000
JC69	11	-4958.5383	9939.0766	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	20	-4728.8002	9600.6145	0.9688	0.9688
GTRIG	21	-4728.8098	9607.7844	0.0269	0.9956
HKYG	16	-4749.1580	9612.7273	0.0023	0.9979
GTRI	20	-4735.0166	9613.0473	0.0019	0.9999
HKYI	16	-4752.4840	9619.3792	0.0001	0.9999
HKYIG	17	-4749.1317	9619.8254	0.0001	1.0000
K2PG	13	-4796.8290	9686.6171	0.0000	1.0000
K2PI	13	-4799.1895	9691.3382	0.0000	1.0000
K2PIG	14	-4796.7258	9693.5613	0.0000	1.0000
SYMG	17	-4789.0519	9699.6657	0.0000	1.0000
GTR	19	-4783.9973	9703.8580	0.0000	1.0000
SYMI	17	-4792.3136	9706.1891	0.0000	1.0000
SYMIG	18	-4789.0325	9706.7776	0.0000	1.0000
HKY	15	-4801.6585	9710.5776	0.0000	1.0000
K2P	12	-4851.3578	9788.5241	0.0000	1.0000
SYM	16	-4844.5073	9803.4258	0.0000	1.0000
F81G	15	-4877.6448	9862.5501	0.0000	1.0000
F81I	15	-4879.0531	9865.3666	0.0000	1.0000
F81IG	16	-4877.4336	9869.2784	0.0000	1.0000
JC69G	12	-4911.8774	9909.5632	0.0000	1.0000
JC69I	12	-4913.1714	9912.1512	0.0000	1.0000
JC69IG	13	-4911.6152	9916.1896	0.0000	1.0000
F81	14	-4922.4417	9944.9933	0.0000	1.0000
JC69	11	-4958.5383	9995.7343	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 591-tomodel-phylip/chr5_14720.phylip = 1 - 1275;
 Partition Dummy = 1:591-tomodel-phylip/chr5_14720.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 591-tomodel-phylip/chr5_14720.phylip = 1 - 1275;
 Partition Dummy = 1:591-tomodel-phylip/chr5_14720.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 (GTRG)]

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

