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

Input data from file "591-tomodel-phylip/chr9_5246.phylip" (7 taxa, 1024 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	-2495.0721	5030.1441	0.4646	0.4646
GTRI	20	-2495.5285	5031.0571	0.2943	0.7589
GTRIG	21	-2495.0031	5032.0063	0.1831	0.9420
HKYG	16	-2501.8648	5035.7297	0.0285	0.9705
HKYI	16	-2502.3089	5036.6177	0.0183	0.9887
HKYIG	17	-2501.7928	5037.5857	0.0113	1.0000
GTR	19	-2512.5779	5063.1558	0.0000	1.0000
HKY	15	-2518.9552	5067.9105	0.0000	1.0000
SYMG	17	-2517.4180	5068.8360	0.0000	1.0000
SYMI	17	-2517.9114	5069.8228	0.0000	1.0000
K2PG	13	-2521.9586	5069.9172	0.0000	1.0000
SYMIG	18	-2517.3306	5070.6613	0.0000	1.0000
K2PI	13	-2522.4169	5070.8339	0.0000	1.0000
K2PIG	14	-2521.8751	5071.7502	0.0000	1.0000
F81G	15	-2532.7933	5095.5867	0.0000	1.0000
F81I	15	-2533.1988	5096.3977	0.0000	1.0000
F81IG	16	-2532.7245	5097.4490	0.0000	1.0000
SYM	16	-2534.7339	5101.4677	0.0000	1.0000
K2P	12	-2539.1981	5102.3962	0.0000	1.0000
JC69G	12	-2550.7419	5125.4837	0.0000	1.0000
JC69I	12	-2551.1692	5126.3383	0.0000	1.0000
JC69IG	13	-2550.6617	5127.3233	0.0000	1.0000
F81	14	-2549.8049	5127.6099	0.0000	1.0000
JC69	11	-2567.9487	5157.8974	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2501.8648	5114.6332	0.5965	0.5965
HKYI	16	-2502.3089	5115.5213	0.3826	0.9791
HKYIG	17	-2501.7928	5121.4207	0.0200	0.9991
GTRG	20	-2495.0721	5128.7735	0.0005	0.9996
GTRI	20	-2495.5285	5129.6865	0.0003	0.9999
K2PG	13	-2521.9586	5134.0263	0.0000	1.0000
K2PI	13	-2522.4169	5134.9430	0.0000	1.0000
GTRIG	21	-2495.0031	5135.5672	0.0000	1.0000
K2PIG	14	-2521.8751	5140.7908	0.0000	1.0000
HKY	15	-2518.9552	5141.8826	0.0000	1.0000
SYMG	17	-2517.4180	5152.6710	0.0000	1.0000
SYMI	17	-2517.9114	5153.6579	0.0000	1.0000
GTR	19	-2512.5779	5156.8538	0.0000	1.0000
SYMIG	18	-2517.3306	5159.4278	0.0000	1.0000
K2P	12	-2539.1981	5161.5739	0.0000	1.0000
F81G	15	-2532.7933	5169.5588	0.0000	1.0000
F81I	15	-2533.1988	5170.3698	0.0000	1.0000
F81IG	16	-2532.7245	5176.3526	0.0000	1.0000
SYM	16	-2534.7339	5180.3713	0.0000	1.0000
JC69G	12	-2550.7419	5184.6614	0.0000	1.0000
JC69I	12	-2551.1692	5185.5160	0.0000	1.0000
JC69IG	13	-2550.6617	5191.4325	0.0000	1.0000
F81	14	-2549.8049	5196.6505	0.0000	1.0000
JC69	11	-2567.9487	5212.1436	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/chr9_5246.phylip = 1 - 1024;
 Partition Dummy = 1:591-tomodel-phylip/chr9_5246.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/chr9_5246.phylip = 1 - 1024;
 Partition Dummy = 1:591-tomodel-phylip/chr9_5246.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/chr9_5246.phylip = 1 - 1024;
 Partition Dummy = 1:591-tomodel-phylip/chr9_5246.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

