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

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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
Minimum BIC  model: HKYG


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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-2585.8468	5203.6935	0.3930	0.3930
GTRG	20	-2582.7414	5205.4828	0.1607	0.5537
HKYIG	17	-2585.7734	5205.5468	0.1556	0.7093
HKYI	16	-2586.9443	5205.8886	0.1312	0.8404
GTRIG	21	-2582.0933	5206.1866	0.1130	0.9534
GTRI	20	-2583.9798	5207.9597	0.0466	1.0000
K2PG	13	-2635.2254	5296.4508	0.0000	1.0000
K2PI	13	-2635.5096	5297.0193	0.0000	1.0000
SYMG	17	-2631.9670	5297.9340	0.0000	1.0000
SYMI	17	-2632.1272	5298.2544	0.0000	1.0000
K2PIG	14	-2635.2332	5298.4663	0.0000	1.0000
SYMIG	18	-2631.9723	5299.9447	0.0000	1.0000
F81G	15	-2657.5246	5345.0492	0.0000	1.0000
F81I	15	-2657.8114	5345.6228	0.0000	1.0000
F81IG	16	-2657.5214	5347.0427	0.0000	1.0000
HKY	15	-2667.2513	5364.5026	0.0000	1.0000
GTR	19	-2664.1265	5366.2529	0.0000	1.0000
JC69I	12	-2694.3151	5412.6301	0.0000	1.0000
JC69G	12	-2694.8273	5413.6546	0.0000	1.0000
JC69IG	13	-2695.7784	5417.5569	0.0000	1.0000
K2P	12	-2713.5408	5451.0817	0.0000	1.0000
SYM	16	-2709.7702	5451.5405	0.0000	1.0000
F81	14	-2732.0620	5492.1239	0.0000	1.0000
JC69	11	-2768.8875	5559.7749	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2585.8468	5280.7262	0.7303	0.7303
HKYI	16	-2586.9443	5282.9213	0.2437	0.9739
HKYIG	17	-2585.7734	5287.3940	0.0260	1.0000
GTRG	20	-2582.7414	5301.7737	0.0000	1.0000
GTRI	20	-2583.9798	5304.2505	0.0000	1.0000
GTRIG	21	-2582.0933	5307.2920	0.0000	1.0000
K2PG	13	-2635.2254	5359.0399	0.0000	1.0000
K2PI	13	-2635.5096	5359.6083	0.0000	1.0000
K2PIG	14	-2635.2332	5365.8699	0.0000	1.0000
SYMG	17	-2631.9670	5379.7812	0.0000	1.0000
SYMI	17	-2632.1272	5380.1016	0.0000	1.0000
SYMIG	18	-2631.9723	5386.6064	0.0000	1.0000
F81G	15	-2657.5246	5417.2673	0.0000	1.0000
F81I	15	-2657.8114	5417.8410	0.0000	1.0000
F81IG	16	-2657.5214	5424.0754	0.0000	1.0000
HKY	15	-2667.2513	5436.7207	0.0000	1.0000
GTR	19	-2664.1265	5457.7292	0.0000	1.0000
JC69I	12	-2694.3151	5470.4046	0.0000	1.0000
JC69G	12	-2694.8273	5471.4292	0.0000	1.0000
JC69IG	13	-2695.7784	5480.1460	0.0000	1.0000
K2P	12	-2713.5408	5508.8562	0.0000	1.0000
SYM	16	-2709.7702	5528.5732	0.0000	1.0000
F81	14	-2732.0620	5559.5275	0.0000	1.0000
JC69	11	-2768.8875	5612.7349	0.0000	1.0000



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

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

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


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

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


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

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

