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

Input data from file "591-tomodel-phylip/chr5_11231.phylip" (7 taxa, 844 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	-2718.7234	5469.4468	0.5583	0.5583
HKYIG	17	-2718.7407	5471.4814	0.2019	0.7602
GTRG	20	-2716.2095	5472.4190	0.1263	0.8865
HKYI	16	-2721.0053	5474.0106	0.0570	0.9435
GTRIG	21	-2716.2190	5474.4379	0.0460	0.9895
GTRI	20	-2718.6976	5477.3953	0.0105	1.0000
HKY	15	-2748.6834	5527.3668	0.0000	1.0000
GTR	19	-2746.9239	5531.8479	0.0000	1.0000
F81G	15	-2763.7840	5557.5680	0.0000	1.0000
F81IG	16	-2763.7896	5559.5792	0.0000	1.0000
F81I	15	-2765.5170	5561.0340	0.0000	1.0000
K2PG	13	-2769.6944	5565.3888	0.0000	1.0000
K2PIG	14	-2769.6919	5567.3838	0.0000	1.0000
K2PI	13	-2772.0065	5570.0130	0.0000	1.0000
SYMG	17	-2768.6733	5571.3466	0.0000	1.0000
SYMIG	18	-2768.6702	5573.3404	0.0000	1.0000
SYMI	17	-2771.0739	5576.1478	0.0000	1.0000
F81	14	-2791.0411	5610.0821	0.0000	1.0000
K2P	12	-2798.9049	5621.8098	0.0000	1.0000
SYM	16	-2798.1274	5628.2548	0.0000	1.0000
JC69G	12	-2806.1388	5636.2776	0.0000	1.0000
JC69IG	13	-2806.1446	5638.2892	0.0000	1.0000
JC69I	12	-2807.8896	5639.7792	0.0000	1.0000
JC69	11	-2833.2719	5688.5438	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2718.7234	5545.2572	0.8803	0.8803
HKYI	16	-2721.0053	5549.8210	0.0899	0.9702
HKYIG	17	-2718.7407	5552.0300	0.0298	1.0000
GTRG	20	-2716.2095	5567.1820	0.0000	1.0000
GTRI	20	-2718.6976	5572.1583	0.0000	1.0000
GTRIG	21	-2716.2190	5573.9391	0.0000	1.0000
HKY	15	-2748.6834	5598.4391	0.0000	1.0000
GTR	19	-2746.9239	5621.8728	0.0000	1.0000
K2PG	13	-2769.6944	5626.9847	0.0000	1.0000
F81G	15	-2763.7840	5628.6403	0.0000	1.0000
K2PI	13	-2772.0065	5631.6090	0.0000	1.0000
F81I	15	-2765.5170	5632.1063	0.0000	1.0000
K2PIG	14	-2769.6919	5633.7179	0.0000	1.0000
F81IG	16	-2763.7896	5635.3896	0.0000	1.0000
SYMG	17	-2768.6733	5651.8952	0.0000	1.0000
SYMI	17	-2771.0739	5656.6964	0.0000	1.0000
SYMIG	18	-2768.6702	5658.6272	0.0000	1.0000
F81	14	-2791.0411	5676.4163	0.0000	1.0000
K2P	12	-2798.9049	5678.6676	0.0000	1.0000
JC69G	12	-2806.1388	5693.1354	0.0000	1.0000
JC69I	12	-2807.8896	5696.6370	0.0000	1.0000
JC69IG	13	-2806.1446	5699.8852	0.0000	1.0000
SYM	16	-2798.1274	5704.0652	0.0000	1.0000
JC69	11	-2833.2719	5740.6635	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_11231.phylip = 1 - 844;
 Partition Dummy = 1:591-tomodel-phylip/chr5_11231.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_11231.phylip = 1 - 844;
 Partition Dummy = 1:591-tomodel-phylip/chr5_11231.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_11231.phylip = 1 - 844;
 Partition Dummy = 1:591-tomodel-phylip/chr5_11231.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

