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

Input data from file "683-tomodel-phylip/chr7_2690.phylip" (7 taxa, 902 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	-2906.8346	5853.6692	0.6070	0.6070
GTRIG	21	-2906.8356	5855.6712	0.2231	0.8300
GTRI	20	-2908.2607	5856.5214	0.1458	0.9759
HKYG	16	-2914.6347	5861.2694	0.0136	0.9894
HKYI	16	-2915.5276	5863.0552	0.0056	0.9950
HKYIG	17	-2914.6363	5863.2727	0.0050	1.0000
GTR	19	-2918.9454	5875.8909	0.0000	1.0000
HKY	15	-2926.4331	5882.8662	0.0000	1.0000
K2PG	13	-2948.1683	5922.3365	0.0000	1.0000
K2PI	13	-2948.7719	5923.5437	0.0000	1.0000
SYMG	17	-2944.8172	5923.6344	0.0000	1.0000
K2PIG	14	-2948.1691	5924.3382	0.0000	1.0000
SYMI	17	-2945.5384	5925.0768	0.0000	1.0000
SYMIG	18	-2944.8243	5925.6487	0.0000	1.0000
K2P	12	-2960.9580	5945.9159	0.0000	1.0000
SYM	16	-2957.5250	5947.0499	0.0000	1.0000
F81G	15	-3018.6890	6067.3779	0.0000	1.0000
F81I	15	-3019.3097	6068.6194	0.0000	1.0000
F81IG	16	-3018.6980	6069.3959	0.0000	1.0000
F81	14	-3028.9073	6085.8146	0.0000	1.0000
JC69G	12	-3046.2552	6116.5104	0.0000	1.0000
JC69I	12	-3046.7111	6117.4222	0.0000	1.0000
JC69IG	13	-3046.2508	6118.5016	0.0000	1.0000
JC69	11	-3056.9708	6135.9415	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2914.6347	5938.1433	0.6912	0.6912
HKYI	16	-2915.5276	5939.9290	0.2830	0.9742
HKYIG	17	-2914.6363	5944.9511	0.0230	0.9972
GTRG	20	-2906.8346	5949.7615	0.0021	0.9993
GTRI	20	-2908.2607	5952.6137	0.0005	0.9998
HKY	15	-2926.4331	5954.9355	0.0002	0.9999
GTRIG	21	-2906.8356	5956.5681	0.0001	1.0000
GTR	19	-2918.9454	5967.1785	0.0000	1.0000
K2PG	13	-2948.1683	5984.7965	0.0000	1.0000
K2PI	13	-2948.7719	5986.0037	0.0000	1.0000
K2PIG	14	-2948.1691	5991.6028	0.0000	1.0000
K2P	12	-2960.9580	6003.5713	0.0000	1.0000
SYMG	17	-2944.8172	6005.3128	0.0000	1.0000
SYMI	17	-2945.5384	6006.7552	0.0000	1.0000
SYMIG	18	-2944.8243	6012.1317	0.0000	1.0000
SYM	16	-2957.5250	6023.9238	0.0000	1.0000
F81G	15	-3018.6890	6139.4472	0.0000	1.0000
F81I	15	-3019.3097	6140.6886	0.0000	1.0000
F81IG	16	-3018.6980	6146.2698	0.0000	1.0000
F81	14	-3028.9073	6153.0792	0.0000	1.0000
JC69G	12	-3046.2552	6174.1658	0.0000	1.0000
JC69I	12	-3046.7111	6175.0776	0.0000	1.0000
JC69IG	13	-3046.2508	6180.9616	0.0000	1.0000
JC69	11	-3056.9708	6188.7923	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/chr7_2690.phylip = 1 - 902;
 Partition Dummy = 1:683-tomodel-phylip/chr7_2690.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/chr7_2690.phylip = 1 - 902;
 Partition Dummy = 1:683-tomodel-phylip/chr7_2690.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/chr7_2690.phylip = 1 - 902;
 Partition Dummy = 1:683-tomodel-phylip/chr7_2690.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

