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

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

Minimum AIC  model: HKYI
Minimum AICc model: HKYI
Minimum BIC  model: HKYI


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
HKYI	16	-808.3342	1650.1038	0.4378	0.4378
HKYG	16	-808.4879	1650.4112	0.3754	0.8132
HKYIG	17	-808.3341	1652.2873	0.1469	0.9601
GTRI	20	-807.4535	1657.1470	0.0129	0.9730
GTRG	20	-807.6196	1657.4792	0.0110	0.9840
F81I	15	-813.8751	1659.0134	0.0051	0.9891
F81G	15	-813.9728	1659.2088	0.0046	0.9937
GTRIG	21	-807.4533	1659.3771	0.0042	0.9980
F81IG	16	-813.8751	1661.1856	0.0017	0.9997
HKY	15	-816.6520	1664.5671	0.0003	1.0000
GTR	19	-815.8165	1671.6542	0.0000	1.0000
F81	14	-821.7852	1672.6727	0.0000	1.0000
K2PI	13	-831.0800	1689.1129	0.0000	1.0000
K2PG	13	-831.3522	1689.6573	0.0000	1.0000
K2PIG	14	-831.0798	1691.2620	0.0000	1.0000
SYMI	17	-829.9001	1695.4193	0.0000	1.0000
SYMG	17	-830.1756	1695.9702	0.0000	1.0000
JC69I	12	-835.7242	1696.2631	0.0000	1.0000
JC69G	12	-835.8660	1696.5467	0.0000	1.0000
SYMIG	18	-829.9002	1697.6147	0.0000	1.0000
JC69IG	13	-835.7242	1698.4012	0.0000	1.0000
K2P	12	-839.8384	1704.4915	0.0000	1.0000
JC69	11	-843.8763	1710.4401	0.0000	1.0000
SYM	16	-838.6260	1710.6874	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	16	-808.3342	1648.6685	0.4258	0.4258
HKYG	16	-808.4879	1648.9758	0.3652	0.7910
HKYIG	17	-808.3341	1650.6683	0.1567	0.9477
GTRI	20	-807.4535	1654.9070	0.0188	0.9665
GTRG	20	-807.6196	1655.2392	0.0159	0.9824
GTRIG	21	-807.4533	1656.9065	0.0069	0.9894
F81I	15	-813.8751	1657.7502	0.0045	0.9939
F81G	15	-813.9728	1657.9456	0.0041	0.9980
F81IG	16	-813.8751	1659.7502	0.0017	0.9997
HKY	15	-816.6520	1663.3040	0.0003	1.0000
GTR	19	-815.8165	1669.6329	0.0000	1.0000
F81	14	-821.7852	1671.5704	0.0000	1.0000
K2PI	13	-831.0800	1688.1600	0.0000	1.0000
K2PG	13	-831.3522	1688.7044	0.0000	1.0000
K2PIG	14	-831.0798	1690.1597	0.0000	1.0000
SYMI	17	-829.9001	1693.8002	0.0000	1.0000
SYMG	17	-830.1756	1694.3511	0.0000	1.0000
JC69I	12	-835.7242	1695.4485	0.0000	1.0000
JC69G	12	-835.8660	1695.7320	0.0000	1.0000
SYMIG	18	-829.9002	1695.8004	0.0000	1.0000
JC69IG	13	-835.7242	1697.4484	0.0000	1.0000
K2P	12	-839.8384	1703.6769	0.0000	1.0000
SYM	16	-838.6260	1709.2520	0.0000	1.0000
JC69	11	-843.8763	1709.7526	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-808.3342	1712.3711	0.4840	0.4840
HKYG	16	-808.4879	1712.6784	0.4151	0.8991
F81I	15	-813.8751	1717.4714	0.0378	0.9369
F81G	15	-813.9728	1717.6669	0.0343	0.9711
HKYIG	17	-808.3341	1718.3523	0.0243	0.9955
HKY	15	-816.6520	1723.0252	0.0024	0.9978
F81IG	16	-813.8751	1723.4528	0.0019	0.9997
F81	14	-821.7852	1727.3102	0.0003	1.0000
GTRI	20	-807.4535	1734.5352	0.0000	1.0000
GTRG	20	-807.6196	1734.8675	0.0000	1.0000
K2PI	13	-831.0800	1739.9184	0.0000	1.0000
K2PG	13	-831.3522	1740.4628	0.0000	1.0000
GTRIG	21	-807.4533	1740.5162	0.0000	1.0000
JC69I	12	-835.7242	1743.2255	0.0000	1.0000
JC69G	12	-835.8660	1743.5090	0.0000	1.0000
GTR	19	-815.8165	1745.2798	0.0000	1.0000
K2PIG	14	-831.0798	1745.8995	0.0000	1.0000
JC69IG	13	-835.7242	1749.2067	0.0000	1.0000
K2P	12	-839.8384	1751.4538	0.0000	1.0000
JC69	11	-843.8763	1753.5482	0.0000	1.0000
SYMI	17	-829.9001	1761.4843	0.0000	1.0000
SYMG	17	-830.1756	1762.0352	0.0000	1.0000
SYMIG	18	-829.9002	1767.4658	0.0000	1.0000
SYM	16	-838.6260	1772.9546	0.0000	1.0000



To specify default values for best AIC, AICc, or BIC models in MrBayes:

[Mrbayes block for the best AIC model (HKYI)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr3_6248.phylip = 1 - 396;
 Partition Dummy = 1:591-tomodel-phylip/chr3_6248.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


[Mrbayes block for the best AICc model (HKYI)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr3_6248.phylip = 1 - 396;
 Partition Dummy = 1:591-tomodel-phylip/chr3_6248.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


[Mrbayes block for the best BIC model (HKYI)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr3_6248.phylip = 1 - 396;
 Partition Dummy = 1:591-tomodel-phylip/chr3_6248.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

