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

Input data from file "683-tomodel-phylip/chr8_4130.phylip" (7 taxa, 883 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	-2696.4216	5424.8432	0.3525	0.3525
GTRG	20	-2692.5272	5425.0545	0.3172	0.6697
HKYIG	17	-2696.3731	5426.7462	0.1361	0.8058
GTRIG	21	-2692.4777	5426.9555	0.1226	0.9284
HKYI	16	-2698.4024	5428.8049	0.0486	0.9770
GTRI	20	-2695.1532	5430.3065	0.0230	1.0000
HKY	15	-2720.9639	5471.9278	0.0000	1.0000
GTR	19	-2717.1483	5472.2966	0.0000	1.0000
K2PG	13	-2746.3172	5518.6345	0.0000	1.0000
SYMG	17	-2742.5059	5519.0118	0.0000	1.0000
K2PIG	14	-2746.2949	5520.5898	0.0000	1.0000
SYMIG	18	-2742.4761	5520.9523	0.0000	1.0000
K2PI	13	-2747.4942	5520.9885	0.0000	1.0000
SYMI	17	-2743.9997	5521.9994	0.0000	1.0000
F81G	15	-2767.0608	5564.1217	0.0000	1.0000
F81I	15	-2768.0447	5566.0895	0.0000	1.0000
F81IG	16	-2767.0494	5566.0988	0.0000	1.0000
SYM	16	-2767.8926	5567.7853	0.0000	1.0000
K2P	12	-2772.1322	5568.2644	0.0000	1.0000
F81	14	-2790.1376	5608.2753	0.0000	1.0000
JC69G	12	-2809.7858	5643.5715	0.0000	1.0000
JC69I	12	-2810.4540	5644.9080	0.0000	1.0000
JC69IG	13	-2809.7947	5645.5894	0.0000	1.0000
JC69	11	-2833.0263	5688.0525	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2696.4216	5501.3764	0.8523	0.8523
HKYI	16	-2698.4024	5505.3381	0.1176	0.9698
HKYIG	17	-2696.3731	5508.0627	0.0301	0.9999
GTRG	20	-2692.5272	5520.7210	0.0001	1.0000
GTRI	20	-2695.1532	5525.9730	0.0000	1.0000
GTRIG	21	-2692.4777	5527.4053	0.0000	1.0000
HKY	15	-2720.9639	5543.6777	0.0000	1.0000
GTR	19	-2717.1483	5563.1798	0.0000	1.0000
K2PG	13	-2746.3172	5580.8177	0.0000	1.0000
K2PI	13	-2747.4942	5583.1717	0.0000	1.0000
K2PIG	14	-2746.2949	5587.5563	0.0000	1.0000
SYMG	17	-2742.5059	5600.3283	0.0000	1.0000
SYMI	17	-2743.9997	5603.3160	0.0000	1.0000
SYMIG	18	-2742.4761	5607.0522	0.0000	1.0000
K2P	12	-2772.1322	5625.6643	0.0000	1.0000
F81G	15	-2767.0608	5635.8715	0.0000	1.0000
F81I	15	-2768.0447	5637.8394	0.0000	1.0000
F81IG	16	-2767.0494	5642.6320	0.0000	1.0000
SYM	16	-2767.8926	5644.3185	0.0000	1.0000
F81	14	-2790.1376	5675.2418	0.0000	1.0000
JC69G	12	-2809.7858	5700.9714	0.0000	1.0000
JC69I	12	-2810.4540	5702.3079	0.0000	1.0000
JC69IG	13	-2809.7947	5707.7726	0.0000	1.0000
JC69	11	-2833.0263	5740.6691	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/chr8_4130.phylip = 1 - 883;
 Partition Dummy = 1:683-tomodel-phylip/chr8_4130.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/chr8_4130.phylip = 1 - 883;
 Partition Dummy = 1:683-tomodel-phylip/chr8_4130.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/chr8_4130.phylip = 1 - 883;
 Partition Dummy = 1:683-tomodel-phylip/chr8_4130.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

