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

Input data from file "683-tomodel-phylip/chr11_6038.phylip" (7 taxa, 668 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	-2266.3389	4565.5135	0.4751	0.4751
HKYIG	17	-2266.3237	4567.5889	0.1683	0.6435
K2PG	13	-2270.6005	4567.7576	0.1547	0.7982
HKYI	16	-2268.2202	4569.2760	0.0724	0.8706
K2PIG	14	-2270.5875	4569.8182	0.0552	0.9258
GTRG	20	-2265.2559	4571.8102	0.0204	0.9462
K2PI	13	-2272.6379	4571.8324	0.0202	0.9664
SYMG	17	-2268.6611	4572.2636	0.0163	0.9826
GTRIG	21	-2265.2429	4573.9161	0.0071	0.9897
SYMIG	18	-2268.6463	4574.3466	0.0057	0.9955
GTRI	20	-2267.2580	4575.8143	0.0028	0.9982
SYMI	17	-2270.8700	4576.6815	0.0018	1.0000
HKY	15	-2289.6549	4610.0460	0.0000	1.0000
K2P	12	-2293.5562	4611.5887	0.0000	1.0000
SYM	16	-2291.6295	4616.0946	0.0000	1.0000
GTR	19	-2288.4904	4616.1535	0.0000	1.0000
F81G	15	-2313.3371	4657.4104	0.0000	1.0000
JC69G	12	-2317.0988	4658.6739	0.0000	1.0000
F81IG	16	-2313.3288	4659.4933	0.0000	1.0000
F81I	15	-2314.8633	4660.4629	0.0000	1.0000
JC69IG	13	-2317.0930	4660.7427	0.0000	1.0000
JC69I	12	-2318.6404	4661.7570	0.0000	1.0000
F81	14	-2334.7002	4698.0435	0.0000	1.0000
JC69	11	-2338.1290	4698.6604	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-2266.3389	4564.6779	0.4792	0.4792
HKYIG	17	-2266.3237	4566.6473	0.1790	0.6582
K2PG	13	-2270.6005	4567.2010	0.1357	0.7940
HKYI	16	-2268.2202	4568.4404	0.0730	0.8670
K2PIG	14	-2270.5875	4569.1751	0.0506	0.9176
GTRG	20	-2265.2559	4570.5119	0.0259	0.9435
K2PI	13	-2272.6379	4571.2758	0.0177	0.9612
SYMG	17	-2268.6611	4571.3221	0.0173	0.9785
GTRIG	21	-2265.2429	4572.4858	0.0097	0.9881
SYMIG	18	-2268.6463	4573.2926	0.0065	0.9946
GTRI	20	-2267.2580	4574.5160	0.0035	0.9981
SYMI	17	-2270.8700	4575.7399	0.0019	1.0000
HKY	15	-2289.6549	4609.3098	0.0000	1.0000
K2P	12	-2293.5562	4611.1124	0.0000	1.0000
GTR	19	-2288.4904	4614.9807	0.0000	1.0000
SYM	16	-2291.6295	4615.2589	0.0000	1.0000
F81G	15	-2313.3371	4656.6742	0.0000	1.0000
JC69G	12	-2317.0988	4658.1975	0.0000	1.0000
F81IG	16	-2313.3288	4658.6576	0.0000	1.0000
F81I	15	-2314.8633	4659.7267	0.0000	1.0000
JC69IG	13	-2317.0930	4660.1861	0.0000	1.0000
JC69I	12	-2318.6404	4661.2807	0.0000	1.0000
F81	14	-2334.7002	4697.4004	0.0000	1.0000
JC69	11	-2338.1290	4698.2580	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	13	-2270.6005	4625.7567	0.8514	0.8514
K2PI	13	-2272.6379	4629.8315	0.1110	0.9624
K2PIG	14	-2270.5875	4632.2351	0.0334	0.9958
HKYG	16	-2266.3389	4636.7465	0.0035	0.9993
HKYI	16	-2268.2202	4640.5090	0.0005	0.9998
HKYIG	17	-2266.3237	4643.2202	0.0001	1.0000
SYMG	17	-2268.6611	4647.8950	0.0000	1.0000
SYMI	17	-2270.8700	4652.3128	0.0000	1.0000
SYMIG	18	-2268.6463	4654.3698	0.0000	1.0000
GTRG	20	-2265.2559	4660.5976	0.0000	1.0000
GTRI	20	-2267.2580	4664.6017	0.0000	1.0000
K2P	12	-2293.5562	4665.1638	0.0000	1.0000
GTRIG	21	-2265.2429	4667.0758	0.0000	1.0000
HKY	15	-2289.6549	4676.8742	0.0000	1.0000
SYM	16	-2291.6295	4687.3276	0.0000	1.0000
GTR	19	-2288.4904	4700.5622	0.0000	1.0000
JC69G	12	-2317.0988	4712.2490	0.0000	1.0000
JC69I	12	-2318.6404	4715.3322	0.0000	1.0000
JC69IG	13	-2317.0930	4718.7418	0.0000	1.0000
F81G	15	-2313.3371	4724.2385	0.0000	1.0000
F81I	15	-2314.8633	4727.2910	0.0000	1.0000
F81IG	16	-2313.3288	4730.7263	0.0000	1.0000
JC69	11	-2338.1290	4747.8052	0.0000	1.0000
F81	14	-2334.7002	4760.4604	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/chr11_6038.phylip = 1 - 668;
 Partition Dummy = 1:683-tomodel-phylip/chr11_6038.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/chr11_6038.phylip = 1 - 668;
 Partition Dummy = 1:683-tomodel-phylip/chr11_6038.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 683-tomodel-phylip/chr11_6038.phylip = 1 - 668;
 Partition Dummy = 1:683-tomodel-phylip/chr11_6038.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

