
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
21:53 04/04/2011

Input data from file "591-tomodel-phylip/chr8_1333.phylip" (7 taxa, 765 characters)

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


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	-2010.7431	4054.2134	0.4192	0.4192
HKYI	16	-2010.8106	4054.3485	0.3918	0.8109
HKYIG	17	-2010.7883	4056.3958	0.1408	0.9517
GTRI	20	-2009.7341	4060.5972	0.0172	0.9689
GTRG	20	-2009.7461	4060.6213	0.0170	0.9859
GTRIG	21	-2009.7883	4062.8203	0.0057	0.9916
K2PG	13	-2018.6552	4063.7951	0.0035	0.9951
K2PI	13	-2018.7903	4064.0653	0.0030	0.9981
K2PIG	14	-2018.6931	4065.9462	0.0012	0.9993
SYMI	17	-2016.8740	4068.5672	0.0003	0.9996
SYMG	17	-2016.9657	4068.7507	0.0003	0.9999
SYMIG	18	-2017.0093	4070.9355	0.0001	1.0000
HKY	15	-2027.8636	4086.3681	0.0000	1.0000
GTR	19	-2026.6297	4092.2795	0.0000	1.0000
K2P	12	-2034.8435	4094.1018	0.0000	1.0000
SYM	16	-2032.8723	4098.4718	0.0000	1.0000
F81G	15	-2066.9363	4164.5135	0.0000	1.0000
F81I	15	-2067.2836	4165.2080	0.0000	1.0000
F81IG	16	-2066.9534	4166.6341	0.0000	1.0000
JC69G	12	-2073.1465	4170.7080	0.0000	1.0000
JC69I	12	-2073.3989	4171.2126	0.0000	1.0000
JC69IG	13	-2073.1685	4172.8216	0.0000	1.0000
F81	14	-2081.5156	4191.5913	0.0000	1.0000
JC69	11	-2087.2136	4196.7778	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-2010.7431	4053.4861	0.4129	0.4129
HKYI	16	-2010.8106	4053.6213	0.3859	0.7988
HKYIG	17	-2010.7883	4055.5765	0.1452	0.9440
GTRI	20	-2009.7341	4059.4681	0.0207	0.9648
GTRG	20	-2009.7461	4059.4923	0.0205	0.9853
GTRIG	21	-2009.7883	4061.5767	0.0072	0.9925
K2PG	13	-2018.6552	4063.3104	0.0030	0.9955
K2PI	13	-2018.7903	4063.5806	0.0027	0.9982
K2PIG	14	-2018.6931	4065.3862	0.0011	0.9993
SYMI	17	-2016.8740	4067.7479	0.0003	0.9996
SYMG	17	-2016.9657	4067.9314	0.0003	0.9999
SYMIG	18	-2017.0093	4070.0186	0.0001	1.0000
HKY	15	-2027.8636	4085.7273	0.0000	1.0000
GTR	19	-2026.6297	4091.2594	0.0000	1.0000
K2P	12	-2034.8435	4093.6869	0.0000	1.0000
SYM	16	-2032.8723	4097.7445	0.0000	1.0000
F81G	15	-2066.9363	4163.8727	0.0000	1.0000
F81I	15	-2067.2836	4164.5671	0.0000	1.0000
F81IG	16	-2066.9534	4165.9068	0.0000	1.0000
JC69G	12	-2073.1465	4170.2931	0.0000	1.0000
JC69I	12	-2073.3989	4170.7977	0.0000	1.0000
JC69IG	13	-2073.1685	4172.3370	0.0000	1.0000
F81	14	-2081.5156	4191.0313	0.0000	1.0000
JC69	11	-2087.2136	4196.4272	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	13	-2018.6552	4123.6288	0.4624	0.4624
K2PI	13	-2018.7903	4123.8990	0.4040	0.8664
HKYG	16	-2010.7431	4127.7241	0.0597	0.9261
HKYI	16	-2010.8106	4127.8593	0.0558	0.9818
K2PIG	14	-2018.6931	4130.3445	0.0161	0.9979
HKYIG	17	-2010.7883	4134.4544	0.0021	1.0000
SYMI	17	-2016.8740	4146.6258	0.0000	1.0000
SYMG	17	-2016.9657	4146.8093	0.0000	1.0000
K2P	12	-2034.8435	4149.3655	0.0000	1.0000
GTRI	20	-2009.7341	4152.2657	0.0000	1.0000
GTRG	20	-2009.7461	4152.2898	0.0000	1.0000
SYMIG	18	-2017.0093	4153.5364	0.0000	1.0000
HKY	15	-2027.8636	4155.3254	0.0000	1.0000
GTRIG	21	-2009.7883	4159.0141	0.0000	1.0000
SYM	16	-2032.8723	4171.9826	0.0000	1.0000
GTR	19	-2026.6297	4179.4170	0.0000	1.0000
JC69G	12	-2073.1465	4225.9716	0.0000	1.0000
JC69I	12	-2073.3989	4226.4762	0.0000	1.0000
JC69IG	13	-2073.1685	4232.6553	0.0000	1.0000
F81G	15	-2066.9363	4233.4708	0.0000	1.0000
F81I	15	-2067.2836	4234.1653	0.0000	1.0000
F81IG	16	-2066.9534	4240.1448	0.0000	1.0000
JC69	11	-2087.2136	4247.4658	0.0000	1.0000
F81	14	-2081.5156	4255.9896	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 591-tomodel-phylip/chr8_1333.phylip = 1 - 765;
 Partition Dummy = 1:591-tomodel-phylip/chr8_1333.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 591-tomodel-phylip/chr8_1333.phylip = 1 - 765;
 Partition Dummy = 1:591-tomodel-phylip/chr8_1333.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 (K2PG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr8_1333.phylip = 1 - 765;
 Partition Dummy = 1:591-tomodel-phylip/chr8_1333.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

