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

Input data from file "591-tomodel-phylip/chr1_28846.phylip" (7 taxa, 1109 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	-2589.4368	5210.8736	0.4543	0.4543
HKYI	16	-2589.8787	5211.7574	0.2920	0.7463
HKYIG	17	-2589.4368	5212.8736	0.1671	0.9134
GTRG	20	-2587.9881	5215.9762	0.0354	0.9488
GTRI	20	-2588.4115	5216.8231	0.0232	0.9720
HKY	15	-2593.9235	5217.8470	0.0139	0.9859
GTRIG	21	-2588.0031	5218.0062	0.0128	0.9988
GTR	19	-2592.3740	5222.7479	0.0012	1.0000
K2PG	13	-2602.9503	5231.9006	0.0000	1.0000
K2PI	13	-2603.3838	5232.7676	0.0000	1.0000
K2PIG	14	-2602.9642	5233.9284	0.0000	1.0000
SYMG	17	-2601.8295	5237.6589	0.0000	1.0000
SYMI	17	-2602.2380	5238.4761	0.0000	1.0000
K2P	12	-2607.6489	5239.2978	0.0000	1.0000
SYMIG	18	-2601.8383	5239.6765	0.0000	1.0000
SYM	16	-2606.4344	5244.8687	0.0000	1.0000
F81G	15	-2623.9335	5277.8671	0.0000	1.0000
F81I	15	-2624.3525	5278.7051	0.0000	1.0000
F81IG	16	-2623.9621	5279.9242	0.0000	1.0000
F81	14	-2628.2865	5284.5730	0.0000	1.0000
JC69G	12	-2636.3790	5296.7581	0.0000	1.0000
JC69I	12	-2636.7786	5297.5572	0.0000	1.0000
JC69IG	13	-2636.4097	5298.8195	0.0000	1.0000
JC69	11	-2640.8071	5303.6142	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2589.4368	5291.0530	0.4658	0.4658
HKYI	16	-2589.8787	5291.9369	0.2994	0.7653
HKY	15	-2593.9235	5293.0152	0.1746	0.9399
K2PG	13	-2602.9503	5297.0464	0.0233	0.9632
K2PI	13	-2603.3838	5297.9134	0.0151	0.9783
HKYIG	17	-2589.4368	5298.0643	0.0140	0.9922
K2P	12	-2607.6489	5299.4323	0.0071	0.9993
K2PIG	14	-2602.9642	5304.0854	0.0007	1.0000
GTRG	20	-2587.9881	5316.2005	0.0000	1.0000
GTRI	20	-2588.4115	5317.0474	0.0000	1.0000
GTR	19	-2592.3740	5317.9610	0.0000	1.0000
SYMG	17	-2601.8295	5322.8496	0.0000	1.0000
GTRIG	21	-2588.0031	5323.2417	0.0000	1.0000
SYMI	17	-2602.2380	5323.6667	0.0000	1.0000
SYM	16	-2606.4344	5325.0482	0.0000	1.0000
SYMIG	18	-2601.8383	5329.8784	0.0000	1.0000
F81G	15	-2623.9335	5353.0353	0.0000	1.0000
F81I	15	-2624.3525	5353.8733	0.0000	1.0000
F81	14	-2628.2865	5354.7300	0.0000	1.0000
JC69G	12	-2636.3790	5356.8926	0.0000	1.0000
JC69I	12	-2636.7786	5357.6917	0.0000	1.0000
JC69	11	-2640.8071	5358.7376	0.0000	1.0000
F81IG	16	-2623.9621	5360.1036	0.0000	1.0000
JC69IG	13	-2636.4097	5363.9653	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/chr1_28846.phylip = 1 - 1109;
 Partition Dummy = 1:591-tomodel-phylip/chr1_28846.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/chr1_28846.phylip = 1 - 1109;
 Partition Dummy = 1:591-tomodel-phylip/chr1_28846.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/chr1_28846.phylip = 1 - 1109;
 Partition Dummy = 1:591-tomodel-phylip/chr1_28846.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

