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

Input data from file "591-tomodel-phylip/chr1_28767.phylip" (7 taxa, 943 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	-2614.4699	5268.9398	0.5452	0.5452
GTRI	20	-2615.2477	5270.4955	0.2505	0.7957
GTRIG	21	-2614.4518	5270.9035	0.2043	1.0000
HKYG	16	-2628.6307	5289.2613	0.0000	1.0000
HKYI	16	-2628.9246	5289.8491	0.0000	1.0000
HKYIG	17	-2628.6209	5291.2417	0.0000	1.0000
GTR	19	-2628.0298	5294.0597	0.0000	1.0000
HKY	15	-2640.8056	5311.6111	0.0000	1.0000
SYMG	17	-2639.6533	5313.3066	0.0000	1.0000
SYMI	17	-2640.1886	5314.3772	0.0000	1.0000
SYMIG	18	-2639.6409	5315.2817	0.0000	1.0000
K2PG	13	-2648.0408	5322.0817	0.0000	1.0000
K2PI	13	-2648.3534	5322.7069	0.0000	1.0000
K2PIG	14	-2648.0290	5324.0580	0.0000	1.0000
F81G	15	-2652.2388	5334.4777	0.0000	1.0000
F81I	15	-2652.5385	5335.0770	0.0000	1.0000
SYM	16	-2652.0606	5336.1211	0.0000	1.0000
F81IG	16	-2652.2369	5336.4738	0.0000	1.0000
K2P	12	-2660.1548	5344.3096	0.0000	1.0000
F81	14	-2663.8035	5355.6070	0.0000	1.0000
JC69G	12	-2672.5975	5369.1950	0.0000	1.0000
JC69I	12	-2672.8825	5369.7650	0.0000	1.0000
JC69IG	13	-2672.5926	5371.1852	0.0000	1.0000
JC69	11	-2684.2323	5390.4646	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	20	-2614.4699	5365.9211	0.3828	0.3828
HKYG	16	-2628.6307	5366.8464	0.2410	0.6238
HKYI	16	-2628.9246	5367.4342	0.1796	0.8034
GTRI	20	-2615.2477	5367.4768	0.1759	0.9793
GTRIG	21	-2614.4518	5372.7339	0.0127	0.9920
HKYIG	17	-2628.6209	5373.6759	0.0079	0.9999
HKY	15	-2640.8056	5384.3471	0.0000	0.9999
K2PG	13	-2648.0408	5385.1195	0.0000	1.0000
K2PI	13	-2648.3534	5385.7447	0.0000	1.0000
GTR	19	-2628.0298	5386.1919	0.0000	1.0000
K2PIG	14	-2648.0290	5391.9449	0.0000	1.0000
SYMG	17	-2639.6533	5395.7407	0.0000	1.0000
SYMI	17	-2640.1886	5396.8113	0.0000	1.0000
K2P	12	-2660.1548	5402.4984	0.0000	1.0000
SYMIG	18	-2639.6409	5402.5649	0.0000	1.0000
F81G	15	-2652.2388	5407.2137	0.0000	1.0000
F81I	15	-2652.5385	5407.8130	0.0000	1.0000
SYM	16	-2652.0606	5413.7062	0.0000	1.0000
F81IG	16	-2652.2369	5414.0589	0.0000	1.0000
F81	14	-2663.8035	5423.4939	0.0000	1.0000
JC69G	12	-2672.5975	5427.3838	0.0000	1.0000
JC69I	12	-2672.8825	5427.9538	0.0000	1.0000
JC69IG	13	-2672.5926	5434.2231	0.0000	1.0000
JC69	11	-2684.2323	5443.8043	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/chr1_28767.phylip = 1 - 943;
 Partition Dummy = 1:591-tomodel-phylip/chr1_28767.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/chr1_28767.phylip = 1 - 943;
 Partition Dummy = 1:591-tomodel-phylip/chr1_28767.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/chr1_28767.phylip = 1 - 943;
 Partition Dummy = 1:591-tomodel-phylip/chr1_28767.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

