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

Input data from file "591-tomodel-phylip/chr2_14216.phylip" (7 taxa, 1100 characters)

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


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-3826.9208	7693.8415	0.7218	0.7218
GTRIG	21	-3826.8994	7695.7987	0.2713	0.9931
GTRI	20	-3831.5727	7703.1455	0.0069	1.0000
SYMG	17	-3847.5313	7729.0627	0.0000	1.0000
SYMIG	18	-3847.5220	7731.0441	0.0000	1.0000
HKYG	16	-3851.1269	7734.2537	0.0000	1.0000
SYMI	17	-3851.0154	7736.0309	0.0000	1.0000
HKYIG	17	-3851.1579	7736.3158	0.0000	1.0000
HKYI	16	-3853.8124	7739.6248	0.0000	1.0000
K2PG	13	-3861.4880	7748.9760	0.0000	1.0000
K2PIG	14	-3861.5047	7751.0093	0.0000	1.0000
K2PI	13	-3864.1095	7754.2189	0.0000	1.0000
GTR	19	-3890.4835	7818.9670	0.0000	1.0000
SYM	16	-3911.1196	7854.2392	0.0000	1.0000
HKY	15	-3914.2184	7858.4367	0.0000	1.0000
K2P	12	-3926.3621	7876.7241	0.0000	1.0000
F81G	15	-3958.9904	7947.9809	0.0000	1.0000
F81IG	16	-3959.0616	7950.1231	0.0000	1.0000
F81I	15	-3960.4927	7950.9855	0.0000	1.0000
JC69G	12	-3966.3279	7956.6559	0.0000	1.0000
JC69IG	13	-3966.3817	7958.7633	0.0000	1.0000
JC69I	12	-3967.9017	7959.8035	0.0000	1.0000
F81	14	-4016.5723	8061.1446	0.0000	1.0000
JC69	11	-4024.9634	8071.9268	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	20	-3826.9208	7793.9028	0.9611	0.9611
GTRIG	21	-3826.8994	7800.8631	0.0296	0.9907
GTRI	20	-3831.5727	7803.2068	0.0092	0.9999
K2PG	13	-3861.4880	7814.0159	0.0000	0.9999
SYMG	17	-3847.5313	7814.1148	0.0000	1.0000
HKYG	16	-3851.1269	7814.3028	0.0000	1.0000
K2PI	13	-3864.1095	7819.2588	0.0000	1.0000
HKYI	16	-3853.8124	7819.6739	0.0000	1.0000
K2PIG	14	-3861.5047	7821.0522	0.0000	1.0000
SYMI	17	-3851.0154	7821.0830	0.0000	1.0000
SYMIG	18	-3847.5220	7821.0992	0.0000	1.0000
HKYIG	17	-3851.1579	7821.3679	0.0000	1.0000
GTR	19	-3890.4835	7914.0252	0.0000	1.0000
HKY	15	-3914.2184	7933.4827	0.0000	1.0000
SYM	16	-3911.1196	7934.2882	0.0000	1.0000
K2P	12	-3926.3621	7936.7609	0.0000	1.0000
JC69G	12	-3966.3279	8016.6927	0.0000	1.0000
JC69I	12	-3967.9017	8019.8403	0.0000	1.0000
F81G	15	-3958.9904	8023.0269	0.0000	1.0000
JC69IG	13	-3966.3817	8023.8032	0.0000	1.0000
F81I	15	-3960.4927	8026.0315	0.0000	1.0000
F81IG	16	-3959.0616	8030.1722	0.0000	1.0000
JC69	11	-4024.9634	8126.9605	0.0000	1.0000
F81	14	-4016.5723	8131.1875	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/chr2_14216.phylip = 1 - 1100;
 Partition Dummy = 1:591-tomodel-phylip/chr2_14216.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/chr2_14216.phylip = 1 - 1100;
 Partition Dummy = 1:591-tomodel-phylip/chr2_14216.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 (GTRG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr2_14216.phylip = 1 - 1100;
 Partition Dummy = 1:591-tomodel-phylip/chr2_14216.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;


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

