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

Input data from file "683-tomodel-phylip/chr4_15363.phylip" (7 taxa, 1146 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	-2901.8257	5835.6514	0.2690	0.2690
HKYIG	17	-2900.8972	5835.7944	0.2505	0.5195
GTRG	20	-2898.4638	5836.9276	0.1421	0.6616
HKYI	16	-2902.4817	5836.9633	0.1396	0.8012
GTRIG	21	-2897.6297	5837.2594	0.1204	0.9216
GTRI	20	-2899.0588	5838.1177	0.0784	1.0000
F81G	15	-2936.1502	5902.3004	0.0000	1.0000
F81IG	16	-2935.4864	5902.9728	0.0000	1.0000
F81I	15	-2936.6963	5903.3926	0.0000	1.0000
GTR	19	-2947.8887	5933.7774	0.0000	1.0000
HKY	15	-2952.0162	5934.0325	0.0000	1.0000
SYMG	17	-2954.9670	5943.9340	0.0000	1.0000
SYMIG	18	-2954.1313	5944.2625	0.0000	1.0000
SYMI	17	-2955.4541	5944.9081	0.0000	1.0000
K2PG	13	-2964.9000	5955.8001	0.0000	1.0000
K2PIG	14	-2964.1489	5956.2977	0.0000	1.0000
K2PI	13	-2965.3848	5956.7695	0.0000	1.0000
F81	14	-2984.4686	5996.9372	0.0000	1.0000
JC69G	12	-2993.4384	6010.8768	0.0000	1.0000
JC69IG	13	-2992.8068	6011.6135	0.0000	1.0000
JC69I	12	-2993.9247	6011.8494	0.0000	1.0000
SYM	16	-3005.7831	6043.5662	0.0000	1.0000
K2P	12	-3014.9103	6053.8206	0.0000	1.0000
JC69	11	-3042.0820	6106.1640	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2901.8257	5916.3559	0.6275	0.6275
HKYI	16	-2902.4817	5917.6678	0.3256	0.9531
HKYIG	17	-2900.8972	5921.5429	0.0469	1.0000
GTRG	20	-2898.4638	5937.8083	0.0000	1.0000
GTRI	20	-2899.0588	5938.9983	0.0000	1.0000
GTRIG	21	-2897.6297	5943.1841	0.0000	1.0000
F81G	15	-2936.1502	5977.9609	0.0000	1.0000
F81I	15	-2936.6963	5979.0531	0.0000	1.0000
F81IG	16	-2935.4864	5983.6774	0.0000	1.0000
HKY	15	-2952.0162	6009.6930	0.0000	1.0000
K2PG	13	-2964.9000	6021.3725	0.0000	1.0000
K2PI	13	-2965.3848	6022.3420	0.0000	1.0000
K2PIG	14	-2964.1489	6026.9142	0.0000	1.0000
GTR	19	-2947.8887	6029.6140	0.0000	1.0000
SYMG	17	-2954.9670	6029.6825	0.0000	1.0000
SYMI	17	-2955.4541	6030.6567	0.0000	1.0000
SYMIG	18	-2954.1313	6035.0551	0.0000	1.0000
F81	14	-2984.4686	6067.5536	0.0000	1.0000
JC69G	12	-2993.4384	6071.4052	0.0000	1.0000
JC69I	12	-2993.9247	6072.3778	0.0000	1.0000
JC69IG	13	-2992.8068	6077.1859	0.0000	1.0000
K2P	12	-3014.9103	6114.3490	0.0000	1.0000
SYM	16	-3005.7831	6124.2707	0.0000	1.0000
JC69	11	-3042.0820	6161.6484	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 683-tomodel-phylip/chr4_15363.phylip = 1 - 1146;
 Partition Dummy = 1:683-tomodel-phylip/chr4_15363.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 683-tomodel-phylip/chr4_15363.phylip = 1 - 1146;
 Partition Dummy = 1:683-tomodel-phylip/chr4_15363.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/chr4_15363.phylip = 1 - 1146;
 Partition Dummy = 1:683-tomodel-phylip/chr4_15363.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

