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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
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
GTRG	20	-2209.8380	4460.8917	0.6066	0.6066
GTRIG	21	-2209.7996	4462.9383	0.2180	0.8247
GTRI	20	-2211.0822	4463.3800	0.1748	0.9995
HKYG	16	-2221.6890	4476.1606	0.0003	0.9998
HKYI	16	-2222.5953	4477.9734	0.0001	0.9999
HKYIG	17	-2221.6843	4478.2505	0.0001	1.0000
GTR	19	-2223.8970	4486.8923	0.0000	1.0000
SYMG	17	-2227.1835	4489.2489	0.0000	1.0000
SYMI	17	-2228.0925	4491.0667	0.0000	1.0000
SYMIG	18	-2227.1714	4491.3298	0.0000	1.0000
K2PG	13	-2234.5163	4495.5541	0.0000	1.0000
K2PI	13	-2235.3054	4497.1322	0.0000	1.0000
K2PIG	14	-2234.5146	4497.6319	0.0000	1.0000
HKY	15	-2233.5666	4497.8228	0.0000	1.0000
SYM	16	-2239.5283	4511.8394	0.0000	1.0000
K2P	12	-2245.9605	4516.3674	0.0000	1.0000
F81G	15	-2245.4393	4521.5682	0.0000	1.0000
F81I	15	-2246.2852	4523.2601	0.0000	1.0000
F81IG	16	-2245.4404	4523.6636	0.0000	1.0000
JC69G	12	-2256.1463	4536.7389	0.0000	1.0000
JC69I	12	-2256.9628	4538.3719	0.0000	1.0000
JC69IG	13	-2256.1473	4538.8160	0.0000	1.0000
F81	14	-2256.4853	4541.5733	0.0000	1.0000
JC69	11	-2267.0007	4556.3786	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-2209.8380	4459.6760	0.5984	0.5984
GTRIG	21	-2209.7996	4461.5992	0.2288	0.8271
GTRI	20	-2211.0822	4462.1644	0.1724	0.9996
HKYG	16	-2221.6890	4475.3779	0.0002	0.9998
HKYI	16	-2222.5953	4477.1906	0.0001	0.9999
HKYIG	17	-2221.6843	4477.3687	0.0001	1.0000
GTR	19	-2223.8970	4485.7941	0.0000	1.0000
SYMG	17	-2227.1835	4488.3670	0.0000	1.0000
SYMI	17	-2228.0925	4490.1849	0.0000	1.0000
SYMIG	18	-2227.1714	4490.3428	0.0000	1.0000
K2PG	13	-2234.5163	4495.0326	0.0000	1.0000
K2PI	13	-2235.3054	4496.6107	0.0000	1.0000
K2PIG	14	-2234.5146	4497.0293	0.0000	1.0000
HKY	15	-2233.5666	4497.1331	0.0000	1.0000
SYM	16	-2239.5283	4511.0567	0.0000	1.0000
K2P	12	-2245.9605	4515.9210	0.0000	1.0000
F81G	15	-2245.4393	4520.8786	0.0000	1.0000
F81I	15	-2246.2852	4522.5705	0.0000	1.0000
F81IG	16	-2245.4404	4522.8808	0.0000	1.0000
JC69G	12	-2256.1463	4536.2926	0.0000	1.0000
JC69I	12	-2256.9628	4537.9256	0.0000	1.0000
JC69IG	13	-2256.1473	4538.2945	0.0000	1.0000
F81	14	-2256.4853	4540.9707	0.0000	1.0000
JC69	11	-2267.0007	4556.0015	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2221.6890	4548.4671	0.5304	0.5304
HKYI	16	-2222.5953	4550.2799	0.2143	0.7447
GTRG	20	-2209.8380	4551.0376	0.1467	0.8914
GTRI	20	-2211.0822	4553.5259	0.0423	0.9337
K2PG	13	-2234.5163	4554.4176	0.0271	0.9608
HKYIG	17	-2221.6843	4555.0260	0.0200	0.9807
K2PI	13	-2235.3054	4555.9957	0.0123	0.9930
GTRIG	21	-2209.7996	4557.5288	0.0057	0.9988
K2PIG	14	-2234.5146	4560.9824	0.0010	0.9998
HKY	15	-2233.5666	4565.6543	0.0001	0.9999
SYMG	17	-2227.1835	4566.0244	0.0001	1.0000
SYMI	17	-2228.0925	4567.8422	0.0000	1.0000
K2P	12	-2245.9605	4570.7379	0.0000	1.0000
SYMIG	18	-2227.1714	4572.5682	0.0000	1.0000
GTR	19	-2223.8970	4572.5875	0.0000	1.0000
SYM	16	-2239.5283	4584.1459	0.0000	1.0000
F81G	15	-2245.4393	4589.3997	0.0000	1.0000
F81I	15	-2246.2852	4591.0916	0.0000	1.0000
JC69G	12	-2256.1463	4591.1095	0.0000	1.0000
JC69I	12	-2256.9628	4592.7425	0.0000	1.0000
F81IG	16	-2245.4404	4595.9701	0.0000	1.0000
JC69IG	13	-2256.1473	4597.6795	0.0000	1.0000
F81	14	-2256.4853	4604.9238	0.0000	1.0000
JC69	11	-2267.0007	4606.2503	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 591-tomodel-phylip/chr2_14247.phylip = 1 - 712;
 Partition Dummy = 1:591-tomodel-phylip/chr2_14247.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 591-tomodel-phylip/chr2_14247.phylip = 1 - 712;
 Partition Dummy = 1:591-tomodel-phylip/chr2_14247.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 (HKYG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr2_14247.phylip = 1 - 712;
 Partition Dummy = 1:591-tomodel-phylip/chr2_14247.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

