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

Input data from file "683-tomodel-phylip/chr19_354.phylip" (7 taxa, 1060 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	-2883.1267	5798.2534	0.3070	0.3070
GTRG	20	-2879.3293	5798.6585	0.2507	0.5577
HKYIG	17	-2882.9439	5799.8877	0.1356	0.6933
GTRIG	21	-2879.0900	5800.1800	0.1172	0.8104
HKYI	16	-2884.1311	5800.2622	0.1124	0.9229
GTRI	20	-2880.5078	5801.0155	0.0771	1.0000
GTR	19	-2910.7610	5859.5220	0.0000	1.0000
HKY	15	-2914.8433	5859.6866	0.0000	1.0000
K2PG	13	-2922.0109	5870.0217	0.0000	1.0000
K2PIG	14	-2921.8187	5871.6373	0.0000	1.0000
K2PI	13	-2922.9482	5871.8963	0.0000	1.0000
SYMG	17	-2920.3580	5874.7159	0.0000	1.0000
SYMIG	18	-2920.1430	5876.2861	0.0000	1.0000
SYMI	17	-2921.3836	5876.7671	0.0000	1.0000
K2P	12	-2955.4698	5934.9395	0.0000	1.0000
F81G	15	-2954.9980	5939.9960	0.0000	1.0000
SYM	16	-2954.0329	5940.0658	0.0000	1.0000
F81IG	16	-2954.8697	5941.7395	0.0000	1.0000
F81I	15	-2955.8962	5941.7924	0.0000	1.0000
F81	14	-2985.3269	5998.6537	0.0000	1.0000
JC69G	12	-2988.3382	6000.6763	0.0000	1.0000
JC69IG	13	-2988.2081	6002.4161	0.0000	1.0000
JC69I	12	-2989.2312	6002.4624	0.0000	1.0000
JC69	11	-3019.5952	6061.1905	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2883.1267	5877.7098	0.7127	0.7127
HKYI	16	-2884.1311	5879.7186	0.2610	0.9737
HKYIG	17	-2882.9439	5884.3102	0.0263	1.0000
GTRG	20	-2879.3293	5897.9790	0.0000	1.0000
GTRI	20	-2880.5078	5900.3360	0.0000	1.0000
GTRIG	21	-2879.0900	5904.4665	0.0000	1.0000
HKY	15	-2914.8433	5934.1770	0.0000	1.0000
K2PG	13	-2922.0109	5934.5800	0.0000	1.0000
K2PI	13	-2922.9482	5936.4546	0.0000	1.0000
K2PIG	14	-2921.8187	5941.1616	0.0000	1.0000
GTR	19	-2910.7610	5953.8765	0.0000	1.0000
SYMG	17	-2920.3580	5959.1384	0.0000	1.0000
SYMI	17	-2921.3836	5961.1896	0.0000	1.0000
SYMIG	18	-2920.1430	5965.6745	0.0000	1.0000
K2P	12	-2955.4698	5994.5318	0.0000	1.0000
F81G	15	-2954.9980	6014.4864	0.0000	1.0000
F81I	15	-2955.8962	6016.2828	0.0000	1.0000
SYM	16	-2954.0329	6019.5222	0.0000	1.0000
F81IG	16	-2954.8697	6021.1958	0.0000	1.0000
JC69G	12	-2988.3382	6060.2686	0.0000	1.0000
JC69I	12	-2989.2312	6062.0547	0.0000	1.0000
JC69IG	13	-2988.2081	6066.9744	0.0000	1.0000
F81	14	-2985.3269	6068.1781	0.0000	1.0000
JC69	11	-3019.5952	6115.8168	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/chr19_354.phylip = 1 - 1060;
 Partition Dummy = 1:683-tomodel-phylip/chr19_354.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/chr19_354.phylip = 1 - 1060;
 Partition Dummy = 1:683-tomodel-phylip/chr19_354.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/chr19_354.phylip = 1 - 1060;
 Partition Dummy = 1:683-tomodel-phylip/chr19_354.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

