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

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

Minimum AIC  model: GTRIG
Minimum AICc model: GTRIG
Minimum BIC  model: SYMG


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	21	-2856.9118	5755.8236	0.5144	0.5144
GTRG	20	-2857.9881	5755.9761	0.4767	0.9911
SYMG	17	-2865.8765	5765.7531	0.0036	0.9947
SYMIG	18	-2864.9411	5765.8822	0.0034	0.9981
GTRI	20	-2863.5036	5767.0073	0.0019	1.0000
SYMI	17	-2870.9213	5775.8425	0.0000	1.0000
HKYG	16	-2876.6737	5785.3474	0.0000	1.0000
HKYIG	17	-2875.9529	5785.9059	0.0000	1.0000
K2PG	13	-2882.0312	5790.0624	0.0000	1.0000
K2PIG	14	-2881.3240	5790.6480	0.0000	1.0000
HKYI	16	-2881.1965	5794.3930	0.0000	1.0000
K2PI	13	-2886.5450	5799.0899	0.0000	1.0000
GTR	19	-2916.3024	5870.6048	0.0000	1.0000
F81G	15	-2921.0153	5872.0306	0.0000	1.0000
F81IG	16	-2920.5093	5873.0186	0.0000	1.0000
JC69G	12	-2925.9345	5875.8690	0.0000	1.0000
SYM	16	-2922.0279	5876.0558	0.0000	1.0000
JC69IG	13	-2925.4498	5876.8995	0.0000	1.0000
F81I	15	-2924.9515	5879.9029	0.0000	1.0000
JC69I	12	-2929.8135	5883.6270	0.0000	1.0000
HKY	15	-2928.6825	5887.3650	0.0000	1.0000
K2P	12	-2933.7994	5891.5988	0.0000	1.0000
F81	14	-2970.3435	5968.6870	0.0000	1.0000
JC69	11	-2975.3276	5972.6553	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
SYMG	17	-2865.8765	5847.4504	0.7782	0.7782
GTRG	20	-2857.9881	5852.0906	0.0765	0.8547
SYMIG	18	-2864.9411	5852.3852	0.0660	0.9207
K2PG	13	-2882.0312	5852.5368	0.0612	0.9819
GTRIG	21	-2856.9118	5856.7438	0.0075	0.9894
SYMI	17	-2870.9213	5857.5398	0.0050	0.9944
K2PIG	14	-2881.3240	5857.9281	0.0041	0.9985
K2PI	13	-2886.5450	5861.5643	0.0007	0.9992
HKYG	16	-2876.6737	5862.2389	0.0005	0.9997
GTRI	20	-2863.5036	5863.1217	0.0003	1.0000
HKYIG	17	-2875.9529	5867.6032	0.0000	1.0000
HKYI	16	-2881.1965	5871.2845	0.0000	1.0000
JC69G	12	-2925.9345	5933.5377	0.0000	1.0000
JC69IG	13	-2925.4498	5939.3739	0.0000	1.0000
JC69I	12	-2929.8135	5941.2956	0.0000	1.0000
F81G	15	-2921.0153	5944.1165	0.0000	1.0000
K2P	12	-2933.7994	5949.2675	0.0000	1.0000
F81IG	16	-2920.5093	5949.9102	0.0000	1.0000
F81I	15	-2924.9515	5951.9887	0.0000	1.0000
SYM	16	-2922.0279	5952.9474	0.0000	1.0000
HKY	15	-2928.6825	5959.4508	0.0000	1.0000
GTR	19	-2916.3024	5961.9135	0.0000	1.0000
JC69	11	-2975.3276	6025.5182	0.0000	1.0000
F81	14	-2970.3435	6035.9671	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chrZ_6385.phylip = 1 - 903;
 Partition Dummy = 1:683-tomodel-phylip/chrZ_6385.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chrZ_6385.phylip = 1 - 903;
 Partition Dummy = 1:683-tomodel-phylip/chrZ_6385.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chrZ_6385.phylip = 1 - 903;
 Partition Dummy = 1:683-tomodel-phylip/chrZ_6385.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;


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

