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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: HKYI


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-2892.0165	5824.0331	0.4212	0.4212
GTRI	20	-2892.0209	5824.0418	0.4194	0.8407
GTRIG	21	-2892.0700	5826.1401	0.1469	0.9876
HKYI	16	-2900.3117	5832.6235	0.0057	0.9933
HKYG	16	-2900.4454	5832.8907	0.0050	0.9983
HKYIG	17	-2900.5424	5835.0849	0.0017	1.0000
SYMG	17	-2918.2353	5870.4705	0.0000	1.0000
SYMI	17	-2918.5793	5871.1586	0.0000	1.0000
SYMIG	18	-2918.2297	5872.4595	0.0000	1.0000
K2PI	13	-2928.2847	5882.5695	0.0000	1.0000
K2PG	13	-2928.4072	5882.8144	0.0000	1.0000
K2PIG	14	-2928.4747	5884.9493	0.0000	1.0000
GTR	19	-2924.7794	5887.5588	0.0000	1.0000
HKY	15	-2933.5293	5897.0585	0.0000	1.0000
SYM	16	-2950.0447	5932.0894	0.0000	1.0000
F81I	15	-2956.9632	5943.9264	0.0000	1.0000
K2P	12	-2960.3849	5944.7698	0.0000	1.0000
F81G	15	-2957.6783	5945.3565	0.0000	1.0000
F81IG	16	-2957.8117	5947.6233	0.0000	1.0000
JC69I	12	-2980.4904	5984.9809	0.0000	1.0000
JC69G	12	-2981.0533	5986.1066	0.0000	1.0000
JC69IG	13	-2980.4984	5986.9968	0.0000	1.0000
F81	14	-2987.6522	6003.3043	0.0000	1.0000
JC69	11	-3010.5325	6043.0649	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-2900.3117	5910.5941	0.5238	0.5238
HKYG	16	-2900.4454	5910.8614	0.4583	0.9821
HKYIG	17	-2900.5424	5917.9287	0.0134	0.9954
GTRG	20	-2892.0165	5921.4964	0.0022	0.9977
GTRI	20	-2892.0209	5921.5051	0.0022	0.9999
GTRIG	21	-2892.0700	5928.4765	0.0001	1.0000
K2PI	13	-2928.2847	5945.9206	0.0000	1.0000
K2PG	13	-2928.4072	5946.1655	0.0000	1.0000
K2PIG	14	-2928.4747	5953.1736	0.0000	1.0000
SYMG	17	-2918.2353	5953.3143	0.0000	1.0000
SYMI	17	-2918.5793	5954.0024	0.0000	1.0000
SYMIG	18	-2918.2297	5960.1764	0.0000	1.0000
HKY	15	-2933.5293	5970.1560	0.0000	1.0000
GTR	19	-2924.7794	5980.1489	0.0000	1.0000
K2P	12	-2960.3849	6003.2478	0.0000	1.0000
SYM	16	-2950.0447	6010.0600	0.0000	1.0000
F81I	15	-2956.9632	6017.0238	0.0000	1.0000
F81G	15	-2957.6783	6018.4540	0.0000	1.0000
F81IG	16	-2957.8117	6025.5939	0.0000	1.0000
JC69I	12	-2980.4904	6043.4588	0.0000	1.0000
JC69G	12	-2981.0533	6044.5846	0.0000	1.0000
JC69IG	13	-2980.4984	6050.3479	0.0000	1.0000
F81	14	-2987.6522	6071.5286	0.0000	1.0000
JC69	11	-3010.5325	6096.6697	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 683-tomodel-phylip/chr1_22335.phylip = 1 - 966;
 Partition Dummy = 1:683-tomodel-phylip/chr1_22335.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 683-tomodel-phylip/chr1_22335.phylip = 1 - 966;
 Partition Dummy = 1:683-tomodel-phylip/chr1_22335.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 (HKYI)]

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr1_22335.phylip = 1 - 966;
 Partition Dummy = 1:683-tomodel-phylip/chr1_22335.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

