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

Input data from file "683-tomodel-phylip/chr5_10191.phylip" (7 taxa, 781 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	-2016.5157	4065.7435	0.4331	0.4331
HKYI	16	-2016.6899	4066.0919	0.3639	0.7970
HKYIG	17	-2016.5070	4067.8161	0.1537	0.9507
GTRG	20	-2015.2756	4071.6564	0.0225	0.9732
GTRI	20	-2015.4522	4072.0097	0.0189	0.9921
GTRIG	21	-2015.2688	4073.7549	0.0079	1.0000
K2PG	13	-2030.3751	4087.2247	0.0000	1.0000
K2PI	13	-2030.4942	4087.4630	0.0000	1.0000
K2PIG	14	-2030.3743	4089.2969	0.0000	1.0000
HKY	15	-2031.5769	4093.7812	0.0000	1.0000
SYMG	17	-2030.1381	4095.0784	0.0000	1.0000
SYMI	17	-2030.2793	4095.3607	0.0000	1.0000
SYMIG	18	-2030.1315	4097.1607	0.0000	1.0000
GTR	19	-2030.4057	4099.8102	0.0000	1.0000
K2P	12	-2045.2568	4114.9198	0.0000	1.0000
F81G	15	-2042.6467	4115.9209	0.0000	1.0000
F81I	15	-2042.6969	4116.0213	0.0000	1.0000
F81IG	16	-2042.6534	4118.0189	0.0000	1.0000
SYM	16	-2044.7310	4122.1741	0.0000	1.0000
JC69G	12	-2054.8713	4134.1489	0.0000	1.0000
JC69I	12	-2054.9071	4134.2205	0.0000	1.0000
JC69IG	13	-2054.8782	4136.2309	0.0000	1.0000
F81	14	-2056.8391	4142.2265	0.0000	1.0000
JC69	11	-2069.0348	4160.4128	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-2016.5157	4065.0314	0.4253	0.4253
HKYI	16	-2016.6899	4065.3799	0.3573	0.7827
HKYIG	17	-2016.5070	4067.0140	0.1578	0.9405
GTRG	20	-2015.2756	4070.5511	0.0269	0.9674
GTRI	20	-2015.4522	4070.9045	0.0226	0.9900
GTRIG	21	-2015.2688	4072.5375	0.0100	1.0000
K2PG	13	-2030.3751	4086.7501	0.0000	1.0000
K2PI	13	-2030.4942	4086.9884	0.0000	1.0000
K2PIG	14	-2030.3743	4088.7486	0.0000	1.0000
HKY	15	-2031.5769	4093.1537	0.0000	1.0000
SYMG	17	-2030.1381	4094.2763	0.0000	1.0000
SYMI	17	-2030.2793	4094.5586	0.0000	1.0000
SYMIG	18	-2030.1315	4096.2631	0.0000	1.0000
GTR	19	-2030.4057	4098.8115	0.0000	1.0000
K2P	12	-2045.2568	4114.5135	0.0000	1.0000
F81G	15	-2042.6467	4115.2935	0.0000	1.0000
F81I	15	-2042.6969	4115.3939	0.0000	1.0000
F81IG	16	-2042.6534	4117.3069	0.0000	1.0000
SYM	16	-2044.7310	4121.4621	0.0000	1.0000
JC69G	12	-2054.8713	4133.7427	0.0000	1.0000
JC69I	12	-2054.9071	4133.8143	0.0000	1.0000
JC69IG	13	-2054.8782	4135.7563	0.0000	1.0000
F81	14	-2056.8391	4141.6782	0.0000	1.0000
JC69	11	-2069.0348	4160.0695	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2016.5157	4139.6006	0.5218	0.5218
HKYI	16	-2016.6899	4139.9491	0.4384	0.9602
HKYIG	17	-2016.5070	4146.2437	0.0188	0.9790
K2PG	13	-2030.3751	4147.3376	0.0109	0.9899
K2PI	13	-2030.4942	4147.5759	0.0097	0.9996
K2PIG	14	-2030.3743	4153.9967	0.0004	1.0000
HKY	15	-2031.5769	4163.0624	0.0000	1.0000
GTRG	20	-2015.2756	4163.7626	0.0000	1.0000
GTRI	20	-2015.4522	4164.1160	0.0000	1.0000
GTRIG	21	-2015.2688	4170.4096	0.0000	1.0000
K2P	12	-2045.2568	4170.4404	0.0000	1.0000
SYMG	17	-2030.1381	4173.5060	0.0000	1.0000
SYMI	17	-2030.2793	4173.7884	0.0000	1.0000
SYMIG	18	-2030.1315	4180.1534	0.0000	1.0000
F81G	15	-2042.6467	4185.2021	0.0000	1.0000
F81I	15	-2042.6969	4185.3025	0.0000	1.0000
GTR	19	-2030.4057	4187.3624	0.0000	1.0000
JC69G	12	-2054.8713	4189.6696	0.0000	1.0000
JC69I	12	-2054.9071	4189.7412	0.0000	1.0000
F81IG	16	-2042.6534	4191.8761	0.0000	1.0000
SYM	16	-2044.7310	4196.0313	0.0000	1.0000
JC69IG	13	-2054.8782	4196.3438	0.0000	1.0000
F81	14	-2056.8391	4206.9263	0.0000	1.0000
JC69	11	-2069.0348	4211.3358	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/chr5_10191.phylip = 1 - 781;
 Partition Dummy = 1:683-tomodel-phylip/chr5_10191.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/chr5_10191.phylip = 1 - 781;
 Partition Dummy = 1:683-tomodel-phylip/chr5_10191.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/chr5_10191.phylip = 1 - 781;
 Partition Dummy = 1:683-tomodel-phylip/chr5_10191.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

