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

Input data from file "591-tomodel-phylip/chrZ_10013.phylip" (7 taxa, 851 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	-3233.0681	6506.1363	0.5433	0.5433
GTRIG	21	-3233.0682	6508.1364	0.1999	0.7432
HKYG	16	-3238.1509	6508.3018	0.1840	0.9272
HKYIG	17	-3238.1510	6510.3020	0.0677	0.9949
GTRI	20	-3238.1738	6516.3476	0.0033	0.9982
HKYI	16	-3242.7543	6517.5086	0.0018	1.0000
GTR	19	-3260.9892	6559.9784	0.0000	1.0000
HKY	15	-3266.0109	6562.0218	0.0000	1.0000
SYMG	17	-3269.4005	6572.8009	0.0000	1.0000
SYMIG	18	-3269.4005	6574.8010	0.0000	1.0000
K2PG	13	-3274.5204	6575.0409	0.0000	1.0000
K2PIG	14	-3274.5205	6577.0410	0.0000	1.0000
SYMI	17	-3273.7115	6581.4230	0.0000	1.0000
K2PI	13	-3278.6248	6583.2495	0.0000	1.0000
SYM	16	-3297.4925	6626.9850	0.0000	1.0000
K2P	12	-3302.1307	6628.2614	0.0000	1.0000
F81G	15	-3305.2910	6640.5820	0.0000	1.0000
F81IG	16	-3305.2910	6642.5820	0.0000	1.0000
F81I	15	-3309.1091	6648.2181	0.0000	1.0000
F81	14	-3331.1221	6690.2443	0.0000	1.0000
JC69G	12	-3333.6871	6691.3742	0.0000	1.0000
JC69IG	13	-3333.6871	6693.3742	0.0000	1.0000
JC69I	12	-3337.2943	6698.5886	0.0000	1.0000
JC69	11	-3359.3496	6740.6993	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-3238.1509	6584.2444	0.9574	0.9574
HKYIG	17	-3238.1510	6590.9910	0.0328	0.9902
HKYI	16	-3242.7543	6593.4512	0.0096	0.9998
GTRG	20	-3233.0681	6601.0645	0.0002	1.0000
GTRIG	21	-3233.0682	6607.8111	0.0000	1.0000
GTRI	20	-3238.1738	6611.2758	0.0000	1.0000
HKY	15	-3266.0109	6633.2180	0.0000	1.0000
K2PG	13	-3274.5204	6636.7443	0.0000	1.0000
K2PIG	14	-3274.5205	6643.4907	0.0000	1.0000
K2PI	13	-3278.6248	6644.9529	0.0000	1.0000
GTR	19	-3260.9892	6650.1602	0.0000	1.0000
SYMG	17	-3269.4005	6653.4899	0.0000	1.0000
SYMIG	18	-3269.4005	6660.2364	0.0000	1.0000
SYMI	17	-3273.7115	6662.1120	0.0000	1.0000
K2P	12	-3302.1307	6685.2184	0.0000	1.0000
SYM	16	-3297.4925	6702.9276	0.0000	1.0000
F81G	15	-3305.2910	6711.7782	0.0000	1.0000
F81IG	16	-3305.2910	6718.5246	0.0000	1.0000
F81I	15	-3309.1091	6719.4143	0.0000	1.0000
JC69G	12	-3333.6871	6748.3311	0.0000	1.0000
JC69IG	13	-3333.6871	6755.0776	0.0000	1.0000
JC69I	12	-3337.2943	6755.5456	0.0000	1.0000
F81	14	-3331.1221	6756.6940	0.0000	1.0000
JC69	11	-3359.3496	6792.9098	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 591-tomodel-phylip/chrZ_10013.phylip = 1 - 851;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_10013.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 591-tomodel-phylip/chrZ_10013.phylip = 1 - 851;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_10013.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 591-tomodel-phylip/chrZ_10013.phylip = 1 - 851;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_10013.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

