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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: HKYI


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-3750.2093	7540.4186	0.4840	0.4840
GTRI	20	-3750.5848	7541.1696	0.3325	0.8166
GTRIG	21	-3750.2093	7542.4186	0.1781	0.9946
HKYI	16	-3759.4763	7550.9526	0.0025	0.9971
HKYG	16	-3759.6516	7551.3033	0.0021	0.9992
HKYIG	17	-3759.6518	7553.3035	0.0008	1.0000
SYMI	17	-3783.4073	7600.8146	0.0000	1.0000
SYMG	17	-3784.4891	7602.9782	0.0000	1.0000
K2PI	13	-3788.9214	7603.8429	0.0000	1.0000
SYMIG	18	-3784.4893	7604.9786	0.0000	1.0000
K2PG	13	-3789.7023	7605.4045	0.0000	1.0000
K2PIG	14	-3789.7024	7607.4049	0.0000	1.0000
F81I	15	-3824.0601	7678.1201	0.0000	1.0000
F81G	15	-3824.9601	7679.9201	0.0000	1.0000
F81IG	16	-3824.9602	7681.9204	0.0000	1.0000
GTR	19	-3823.6943	7685.3887	0.0000	1.0000
HKY	15	-3832.1586	7694.3173	0.0000	1.0000
JC69I	12	-3853.9364	7731.8729	0.0000	1.0000
JC69G	12	-3855.2904	7734.5808	0.0000	1.0000
JC69IG	13	-3855.2906	7736.5813	0.0000	1.0000
SYM	16	-3856.8545	7745.7089	0.0000	1.0000
K2P	12	-3861.3941	7746.7883	0.0000	1.0000
F81	14	-3891.2157	7810.4314	0.0000	1.0000
JC69	11	-3921.7079	7865.4157	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-3759.4763	7632.6058	0.5332	0.5332
HKYG	16	-3759.6516	7632.9564	0.4474	0.9806
HKYIG	17	-3759.6518	7640.0600	0.0128	0.9935
GTRG	20	-3750.2093	7642.4851	0.0038	0.9973
GTRI	20	-3750.5848	7643.2360	0.0026	0.9999
GTRIG	21	-3750.2093	7649.5884	0.0001	1.0000
K2PI	13	-3788.9214	7670.1860	0.0000	1.0000
K2PG	13	-3789.7023	7671.7477	0.0000	1.0000
K2PIG	14	-3789.7024	7678.8514	0.0000	1.0000
SYMI	17	-3783.4073	7687.5710	0.0000	1.0000
SYMG	17	-3784.4891	7689.7347	0.0000	1.0000
SYMIG	18	-3784.4893	7696.8384	0.0000	1.0000
F81I	15	-3824.0601	7754.6700	0.0000	1.0000
F81G	15	-3824.9601	7756.4700	0.0000	1.0000
F81IG	16	-3824.9602	7763.5736	0.0000	1.0000
HKY	15	-3832.1586	7770.8671	0.0000	1.0000
GTR	19	-3823.6943	7782.3518	0.0000	1.0000
JC69I	12	-3853.9364	7793.1127	0.0000	1.0000
JC69G	12	-3855.2904	7795.8207	0.0000	1.0000
JC69IG	13	-3855.2906	7802.9245	0.0000	1.0000
K2P	12	-3861.3941	7808.0281	0.0000	1.0000
SYM	16	-3856.8545	7827.3621	0.0000	1.0000
F81	14	-3891.2157	7881.8779	0.0000	1.0000
JC69	11	-3921.7079	7921.5523	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 683-tomodel-phylip/chr6_7254.phylip = 1 - 1216;
 Partition Dummy = 1:683-tomodel-phylip/chr6_7254.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 683-tomodel-phylip/chr6_7254.phylip = 1 - 1216;
 Partition Dummy = 1:683-tomodel-phylip/chr6_7254.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 (HKYI)]

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr6_7254.phylip = 1 - 1216;
 Partition Dummy = 1:683-tomodel-phylip/chr6_7254.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


-------------------------------------------------------------
End of Output

