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

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

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


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-2815.9759	5671.9519	0.2585	0.2585
HKYG	16	-2820.0035	5672.0071	0.2515	0.5100
HKYI	16	-2820.4912	5672.9825	0.1544	0.6644
GTRI	20	-2816.5292	5673.0584	0.1487	0.8130
GTRIG	21	-2815.9843	5673.9686	0.0943	0.9073
HKYIG	17	-2820.0105	5674.0210	0.0919	0.9992
HKY	15	-2827.3470	5684.6941	0.0004	0.9996
GTR	19	-2823.5059	5685.0118	0.0004	1.0000
K2PG	13	-2877.0328	5780.0657	0.0000	1.0000
K2PI	13	-2877.3241	5780.6483	0.0000	1.0000
K2PIG	14	-2877.0242	5782.0484	0.0000	1.0000
SYMG	17	-2876.5365	5787.0729	0.0000	1.0000
SYMI	17	-2876.8411	5787.6821	0.0000	1.0000
SYMIG	18	-2876.5279	5789.0557	0.0000	1.0000
K2P	12	-2885.2842	5794.5685	0.0000	1.0000
SYM	16	-2884.8658	5801.7317	0.0000	1.0000
F81G	15	-2899.8829	5829.7657	0.0000	1.0000
F81I	15	-2900.0693	5830.1385	0.0000	1.0000
F81IG	16	-2899.8683	5831.7365	0.0000	1.0000
F81	14	-2906.1454	5840.2909	0.0000	1.0000
JC69G	12	-2945.5405	5915.0810	0.0000	1.0000
JC69I	12	-2945.6923	5915.3846	0.0000	1.0000
JC69IG	13	-2945.5187	5917.0373	0.0000	1.0000
JC69	11	-2952.5186	5927.0373	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2820.0035	5749.1099	0.6000	0.6000
HKYI	16	-2820.4912	5750.0852	0.3685	0.9685
HKYIG	17	-2820.0105	5755.9427	0.0197	0.9882
HKY	15	-2827.3470	5756.9779	0.0117	0.9999
GTRG	20	-2815.9759	5768.3303	0.0000	1.0000
GTRI	20	-2816.5292	5769.4369	0.0000	1.0000
GTRIG	21	-2815.9843	5775.1660	0.0000	1.0000
GTR	19	-2823.5059	5776.5714	0.0000	1.0000
K2PG	13	-2877.0328	5842.7117	0.0000	1.0000
K2PI	13	-2877.3241	5843.2943	0.0000	1.0000
K2PIG	14	-2877.0242	5849.5133	0.0000	1.0000
K2P	12	-2885.2842	5852.3955	0.0000	1.0000
SYMG	17	-2876.5365	5868.9946	0.0000	1.0000
SYMI	17	-2876.8411	5869.6038	0.0000	1.0000
SYMIG	18	-2876.5279	5875.7964	0.0000	1.0000
SYM	16	-2884.8658	5878.8345	0.0000	1.0000
F81G	15	-2899.8829	5902.0496	0.0000	1.0000
F81I	15	-2900.0693	5902.4224	0.0000	1.0000
F81	14	-2906.1454	5907.7558	0.0000	1.0000
F81IG	16	-2899.8683	5908.8393	0.0000	1.0000
JC69G	12	-2945.5405	5972.9080	0.0000	1.0000
JC69I	12	-2945.6923	5973.2117	0.0000	1.0000
JC69IG	13	-2945.5187	5979.6834	0.0000	1.0000
JC69	11	-2952.5186	5980.0454	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/chr13_2771.phylip = 1 - 915;
 Partition Dummy = 1:683-tomodel-phylip/chr13_2771.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 (HKYG)]

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

