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

Input data from file "683-tomodel-phylip/chr2_8757.phylip" (7 taxa, 294 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	-814.1754	1662.3147	0.4777	0.4777
HKYI	16	-814.5916	1663.1472	0.3150	0.7927
HKYIG	17	-814.1335	1664.4844	0.1614	0.9541
GTRG	20	-812.5960	1668.2689	0.0243	0.9785
GTRI	20	-813.2202	1669.5173	0.0130	0.9915
GTRIG	21	-812.5105	1670.4180	0.0083	0.9998
HKY	15	-823.5772	1678.8811	0.0001	0.9999
F81G	15	-824.9250	1681.5766	0.0000	1.0000
F81I	15	-825.1587	1682.0440	0.0000	1.0000
F81IG	16	-824.9105	1683.7849	0.0000	1.0000
GTR	19	-821.9953	1684.7643	0.0000	1.0000
K2PG	13	-833.5430	1694.3860	0.0000	1.0000
K2PI	13	-833.7513	1694.8026	0.0000	1.0000
K2PIG	14	-833.5361	1696.5775	0.0000	1.0000
F81	14	-834.3700	1698.2454	0.0000	1.0000
SYMG	17	-832.3605	1700.9385	0.0000	1.0000
SYMI	17	-832.7139	1701.6453	0.0000	1.0000
SYMIG	18	-832.3288	1703.1448	0.0000	1.0000
K2P	12	-843.5819	1712.2742	0.0000	1.0000
JC69G	12	-844.2401	1713.5905	0.0000	1.0000
JC69I	12	-844.3868	1713.8840	0.0000	1.0000
JC69IG	13	-844.2448	1715.7895	0.0000	1.0000
SYM	16	-842.1721	1718.3080	0.0000	1.0000
JC69	11	-853.7207	1730.3776	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-814.1754	1660.3508	0.4513	0.4513
HKYI	16	-814.5916	1661.1833	0.2977	0.7490
HKYIG	17	-814.1335	1662.2670	0.1731	0.9222
GTRG	20	-812.5960	1665.1920	0.0401	0.9623
GTRI	20	-813.2202	1666.4404	0.0215	0.9838
GTRIG	21	-812.5105	1667.0210	0.0161	0.9998
HKY	15	-823.5772	1677.1545	0.0001	0.9999
F81G	15	-824.9250	1679.8500	0.0000	1.0000
F81I	15	-825.1587	1680.3173	0.0000	1.0000
F81IG	16	-824.9105	1681.8210	0.0000	1.0000
GTR	19	-821.9953	1681.9905	0.0000	1.0000
K2PG	13	-833.5430	1693.0860	0.0000	1.0000
K2PI	13	-833.7513	1693.5026	0.0000	1.0000
K2PIG	14	-833.5361	1695.0722	0.0000	1.0000
F81	14	-834.3700	1696.7400	0.0000	1.0000
SYMG	17	-832.3605	1698.7211	0.0000	1.0000
SYMI	17	-832.7139	1699.4279	0.0000	1.0000
SYMIG	18	-832.3288	1700.6576	0.0000	1.0000
K2P	12	-843.5819	1711.1639	0.0000	1.0000
JC69G	12	-844.2401	1712.4802	0.0000	1.0000
JC69I	12	-844.3868	1712.7737	0.0000	1.0000
JC69IG	13	-844.2448	1714.4895	0.0000	1.0000
SYM	16	-842.1721	1716.3441	0.0000	1.0000
JC69	11	-853.7207	1729.4414	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-814.1754	1719.2881	0.5805	0.5805
HKYI	16	-814.5916	1720.1206	0.3829	0.9634
HKYIG	17	-814.1335	1724.8878	0.0353	0.9987
HKY	15	-823.5772	1732.4082	0.0008	0.9995
F81G	15	-824.9250	1735.1037	0.0002	0.9997
F81I	15	-825.1587	1735.5710	0.0002	0.9999
GTRG	20	-812.5960	1738.8636	0.0000	0.9999
GTRI	20	-813.2202	1740.1120	0.0000	1.0000
F81IG	16	-824.9105	1740.7583	0.0000	1.0000
K2PG	13	-833.5430	1740.9725	0.0000	1.0000
K2PI	13	-833.7513	1741.3891	0.0000	1.0000
GTRIG	21	-812.5105	1744.3762	0.0000	1.0000
K2PIG	14	-833.5361	1746.6423	0.0000	1.0000
F81	14	-834.3700	1748.3101	0.0000	1.0000
GTR	19	-821.9953	1751.9786	0.0000	1.0000
K2P	12	-843.5819	1755.3668	0.0000	1.0000
JC69G	12	-844.2401	1756.6832	0.0000	1.0000
JC69I	12	-844.3868	1756.9766	0.0000	1.0000
SYMG	17	-832.3605	1761.3419	0.0000	1.0000
SYMI	17	-832.7139	1762.0487	0.0000	1.0000
JC69IG	13	-844.2448	1762.3761	0.0000	1.0000
SYMIG	18	-832.3288	1766.9620	0.0000	1.0000
JC69	11	-853.7207	1769.9608	0.0000	1.0000
SYM	16	-842.1721	1775.2814	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/chr2_8757.phylip = 1 - 294;
 Partition Dummy = 1:683-tomodel-phylip/chr2_8757.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/chr2_8757.phylip = 1 - 294;
 Partition Dummy = 1:683-tomodel-phylip/chr2_8757.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/chr2_8757.phylip = 1 - 294;
 Partition Dummy = 1:683-tomodel-phylip/chr2_8757.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

