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

Input data from file "683-tomodel-phylip/chrZ_4766.phylip" (7 taxa, 1187 characters)

Minimum AIC  model: SYMG
Minimum AICc model: K2PG
Minimum BIC  model: K2PG


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

Model	df	lnL		AIC		wAIC	cwAIC
SYMG	17	-3596.7638	7227.5277	0.2732	0.2732
K2PG	13	-3600.8177	7227.6355	0.2589	0.5321
SYMIG	18	-3596.7108	7229.4216	0.1060	0.6380
K2PIG	14	-3600.8038	7229.6076	0.0966	0.7346
K2PI	13	-3601.8351	7229.6701	0.0936	0.8282
SYMI	17	-3597.8694	7229.7388	0.0904	0.9186
GTRG	20	-3596.1318	7232.2637	0.0256	0.9442
HKYG	16	-3600.3208	7232.6417	0.0212	0.9654
GTRIG	21	-3596.0720	7234.1439	0.0100	0.9754
GTRI	20	-3597.1723	7234.3447	0.0090	0.9844
HKYIG	17	-3600.3054	7234.6108	0.0079	0.9923
HKYI	16	-3601.3381	7234.6762	0.0077	1.0000
K2P	12	-3633.1290	7290.2580	0.0000	1.0000
SYM	16	-3629.4753	7290.9507	0.0000	1.0000
HKY	15	-3632.7340	7295.4679	0.0000	1.0000
GTR	19	-3629.1550	7296.3101	0.0000	1.0000
JC69G	12	-3659.6896	7343.3792	0.0000	1.0000
JC69I	12	-3660.6650	7345.3300	0.0000	1.0000
JC69IG	13	-3659.7036	7345.4073	0.0000	1.0000
F81G	15	-3659.4158	7348.8315	0.0000	1.0000
F81I	15	-3660.3731	7350.7463	0.0000	1.0000
F81IG	16	-3659.4303	7350.8607	0.0000	1.0000
JC69	11	-3689.8816	7401.7632	0.0000	1.0000
F81	14	-3689.6838	7407.3677	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	13	-3600.8177	7293.6649	0.7189	0.7189
K2PI	13	-3601.8351	7295.6995	0.2599	0.9788
K2PIG	14	-3600.8038	7300.7162	0.0212	0.9999
SYMG	17	-3596.7638	7313.8738	0.0000	0.9999
HKYG	16	-3600.3208	7313.9086	0.0000	1.0000
HKYI	16	-3601.3381	7315.9432	0.0000	1.0000
SYMI	17	-3597.8694	7316.0850	0.0000	1.0000
SYMIG	18	-3596.7108	7320.8469	0.0000	1.0000
HKYIG	17	-3600.3054	7320.9569	0.0000	1.0000
GTRG	20	-3596.1318	7333.8474	0.0000	1.0000
GTRI	20	-3597.1723	7335.9284	0.0000	1.0000
GTRIG	21	-3596.0720	7340.8068	0.0000	1.0000
K2P	12	-3633.1290	7351.2082	0.0000	1.0000
HKY	15	-3632.7340	7371.6557	0.0000	1.0000
SYM	16	-3629.4753	7372.2176	0.0000	1.0000
GTR	19	-3629.1550	7392.8146	0.0000	1.0000
JC69G	12	-3659.6896	7404.3294	0.0000	1.0000
JC69I	12	-3660.6650	7406.2802	0.0000	1.0000
JC69IG	13	-3659.7036	7411.4367	0.0000	1.0000
F81G	15	-3659.4158	7425.0193	0.0000	1.0000
F81I	15	-3660.3731	7426.9341	0.0000	1.0000
F81IG	16	-3659.4303	7432.1276	0.0000	1.0000
JC69	11	-3689.8816	7457.6342	0.0000	1.0000
F81	14	-3689.6838	7478.4763	0.0000	1.0000



To specify default values for best AIC, AICc, or BIC models in MrBayes:

[Mrbayes block for the best AIC model (SYMG)]

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chrZ_4766.phylip = 1 - 1187;
 Partition Dummy = 1:683-tomodel-phylip/chrZ_4766.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=Fixed(Equal);
END;


[Mrbayes block for the best AICc model (K2PG)]

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chrZ_4766.phylip = 1 - 1187;
 Partition Dummy = 1:683-tomodel-phylip/chrZ_4766.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


[Mrbayes block for the best BIC model (K2PG)]

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chrZ_4766.phylip = 1 - 1187;
 Partition Dummy = 1:683-tomodel-phylip/chrZ_4766.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


-------------------------------------------------------------
End of Output

