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

Input data from file "683-tomodel-phylip/chr9_6300.phylip" (7 taxa, 1067 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	-2750.5551	5533.1101	0.4345	0.4345
HKYIG	17	-2750.1431	5534.2862	0.2413	0.6758
HKYI	16	-2751.1711	5534.3422	0.2347	0.9105
GTRG	20	-2748.8656	5537.7312	0.0431	0.9536
GTRIG	21	-2748.4751	5538.9502	0.0234	0.9770
GTRI	20	-2749.4937	5538.9873	0.0230	1.0000
F81G	15	-2790.9598	5611.9197	0.0000	1.0000
F81I	15	-2791.3225	5612.6450	0.0000	1.0000
F81IG	16	-2790.8762	5613.7523	0.0000	1.0000
HKY	15	-2792.9802	5615.9604	0.0000	1.0000
GTR	19	-2790.8276	5619.6553	0.0000	1.0000
K2PG	13	-2813.4074	5652.8148	0.0000	1.0000
K2PIG	14	-2813.0156	5654.0312	0.0000	1.0000
K2PI	13	-2814.1215	5654.2431	0.0000	1.0000
SYMG	17	-2810.8856	5655.7713	0.0000	1.0000
SYMIG	18	-2810.4172	5656.8343	0.0000	1.0000
SYMI	17	-2811.8421	5657.6842	0.0000	1.0000
F81	14	-2830.6594	5689.3189	0.0000	1.0000
JC69G	12	-2849.6377	5723.2754	0.0000	1.0000
JC69I	12	-2850.1651	5724.3302	0.0000	1.0000
JC69IG	13	-2849.4817	5724.9634	0.0000	1.0000
K2P	12	-2853.3992	5730.7984	0.0000	1.0000
SYM	16	-2850.4063	5732.8126	0.0000	1.0000
JC69	11	-2887.7578	5797.5156	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2750.5551	5612.6718	0.6304	0.6304
HKYI	16	-2751.1711	5613.9039	0.3405	0.9709
HKYIG	17	-2750.1431	5618.8205	0.0291	1.0000
GTRG	20	-2748.8656	5637.1834	0.0000	1.0000
GTRI	20	-2749.4937	5638.4394	0.0000	1.0000
GTRIG	21	-2748.4751	5643.3750	0.0000	1.0000
F81G	15	-2790.9598	5686.5088	0.0000	1.0000
F81I	15	-2791.3225	5687.2341	0.0000	1.0000
HKY	15	-2792.9802	5690.5495	0.0000	1.0000
F81IG	16	-2790.8762	5693.3140	0.0000	1.0000
GTR	19	-2790.8276	5714.1348	0.0000	1.0000
K2PG	13	-2813.4074	5717.4587	0.0000	1.0000
K2PI	13	-2814.1215	5718.8869	0.0000	1.0000
K2PIG	14	-2813.0156	5723.6477	0.0000	1.0000
SYMG	17	-2810.8856	5740.3056	0.0000	1.0000
SYMI	17	-2811.8421	5742.2185	0.0000	1.0000
SYMIG	18	-2810.4172	5746.3412	0.0000	1.0000
F81	14	-2830.6594	5758.9353	0.0000	1.0000
JC69G	12	-2849.6377	5782.9467	0.0000	1.0000
JC69I	12	-2850.1651	5784.0015	0.0000	1.0000
JC69IG	13	-2849.4817	5789.6073	0.0000	1.0000
K2P	12	-2853.3992	5790.4696	0.0000	1.0000
SYM	16	-2850.4063	5812.3743	0.0000	1.0000
JC69	11	-2887.7578	5852.2142	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 683-tomodel-phylip/chr9_6300.phylip = 1 - 1067;
 Partition Dummy = 1:683-tomodel-phylip/chr9_6300.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 683-tomodel-phylip/chr9_6300.phylip = 1 - 1067;
 Partition Dummy = 1:683-tomodel-phylip/chr9_6300.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 683-tomodel-phylip/chr9_6300.phylip = 1 - 1067;
 Partition Dummy = 1:683-tomodel-phylip/chr9_6300.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

