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

Input data from file "683-tomodel-phylip/chr20_1391.phylip" (7 taxa, 833 characters)

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


Note: number of parameters compared to sample size
is low (Nchar/Nparams < 40). Consider the use of AICc.


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

Model	df	lnL		AICc		wAICc	cwAICc
HKYG	16	-3177.2395	6387.1458	0.5635	0.5635
HKYIG	17	-3177.1651	6389.0811	0.2141	0.7776
HKYI	16	-3178.4284	6389.5236	0.1716	0.9492
GTRG	20	-3175.9495	6392.9335	0.0312	0.9804
GTRIG	21	-3175.8883	6394.9160	0.0116	0.9920
GTRI	20	-3177.3122	6395.6589	0.0080	1.0000
K2PG	13	-3193.0918	6412.6281	0.0000	1.0000
K2PIG	14	-3192.9998	6414.5130	0.0000	1.0000
K2PI	13	-3194.0426	6414.5296	0.0000	1.0000
SYMG	17	-3191.6218	6417.9945	0.0000	1.0000
SYMI	17	-3192.3082	6419.3672	0.0000	1.0000
SYMIG	18	-3191.4489	6419.7381	0.0000	1.0000
HKY	15	-3204.7076	6440.0026	0.0000	1.0000
GTR	19	-3203.7411	6446.4171	0.0000	1.0000
K2P	12	-3220.4769	6465.3343	0.0000	1.0000
SYM	16	-3218.8636	6470.3938	0.0000	1.0000
F81G	15	-3288.4914	6607.5704	0.0000	1.0000
F81I	15	-3289.3452	6609.2780	0.0000	1.0000
F81IG	16	-3288.4402	6609.5470	0.0000	1.0000
JC69G	12	-3302.4330	6629.2465	0.0000	1.0000
JC69I	12	-3303.0532	6630.4868	0.0000	1.0000
JC69IG	13	-3302.3330	6631.1104	0.0000	1.0000
F81	14	-3313.0009	6654.5152	0.0000	1.0000
JC69	11	-3326.7485	6675.8185	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-3177.2395	6386.4791	0.5523	0.5523
HKYIG	17	-3177.1651	6388.3302	0.2189	0.7712
HKYI	16	-3178.4284	6388.8569	0.1682	0.9395
GTRG	20	-3175.9495	6391.8990	0.0368	0.9762
GTRIG	21	-3175.8883	6393.7767	0.0144	0.9906
GTRI	20	-3177.3122	6394.6244	0.0094	1.0000
K2PG	13	-3193.0918	6412.1837	0.0000	1.0000
K2PIG	14	-3192.9998	6413.9995	0.0000	1.0000
K2PI	13	-3194.0426	6414.0851	0.0000	1.0000
SYMG	17	-3191.6218	6417.2435	0.0000	1.0000
SYMI	17	-3192.3082	6418.6163	0.0000	1.0000
SYMIG	18	-3191.4489	6418.8978	0.0000	1.0000
HKY	15	-3204.7076	6439.4151	0.0000	1.0000
GTR	19	-3203.7411	6445.4823	0.0000	1.0000
K2P	12	-3220.4769	6464.9538	0.0000	1.0000
SYM	16	-3218.8636	6469.7271	0.0000	1.0000
F81G	15	-3288.4914	6606.9829	0.0000	1.0000
F81I	15	-3289.3452	6608.6905	0.0000	1.0000
F81IG	16	-3288.4402	6608.8803	0.0000	1.0000
JC69G	12	-3302.4330	6628.8660	0.0000	1.0000
JC69I	12	-3303.0532	6630.1063	0.0000	1.0000
JC69IG	13	-3302.3330	6630.6660	0.0000	1.0000
F81	14	-3313.0009	6654.0018	0.0000	1.0000
JC69	11	-3326.7485	6675.4970	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-3177.2395	6462.0796	0.7427	0.7427
HKYI	16	-3178.4284	6464.4574	0.2262	0.9690
HKYIG	17	-3177.1651	6468.6558	0.0277	0.9967
K2PG	13	-3193.0918	6473.6091	0.0023	0.9990
K2PI	13	-3194.0426	6475.5106	0.0009	0.9999
K2PIG	14	-3192.9998	6480.1500	0.0001	1.0000
GTRG	20	-3175.9495	6486.3997	0.0000	1.0000
GTRI	20	-3177.3122	6489.1251	0.0000	1.0000
GTRIG	21	-3175.8883	6493.0024	0.0000	1.0000
SYMG	17	-3191.6218	6497.5691	0.0000	1.0000
SYMI	17	-3192.3082	6498.9419	0.0000	1.0000
SYMIG	18	-3191.4489	6503.9484	0.0000	1.0000
HKY	15	-3204.7076	6510.2906	0.0000	1.0000
K2P	12	-3220.4769	6521.6542	0.0000	1.0000
GTR	19	-3203.7411	6535.2579	0.0000	1.0000
SYM	16	-3218.8636	6545.3277	0.0000	1.0000
F81G	15	-3288.4914	6677.8584	0.0000	1.0000
F81I	15	-3289.3452	6679.5660	0.0000	1.0000
F81IG	16	-3288.4402	6684.4809	0.0000	1.0000
JC69G	12	-3302.4330	6685.5664	0.0000	1.0000
JC69I	12	-3303.0532	6686.8067	0.0000	1.0000
JC69IG	13	-3302.3330	6692.0914	0.0000	1.0000
F81	14	-3313.0009	6720.1522	0.0000	1.0000
JC69	11	-3326.7485	6727.4724	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/chr20_1391.phylip = 1 - 833;
 Partition Dummy = 1:683-tomodel-phylip/chr20_1391.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/chr20_1391.phylip = 1 - 833;
 Partition Dummy = 1:683-tomodel-phylip/chr20_1391.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/chr20_1391.phylip = 1 - 833;
 Partition Dummy = 1:683-tomodel-phylip/chr20_1391.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

