
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
00:52 02/17/2011

Input data from file "183-phylip/chr5_10180.phylip" (29 taxa, 209 characters)

Minimum AIC  model: K2PG
Minimum AICc model: K2PG
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
K2PG	57	-794.9437	1747.6754	0.5097	0.5097
K2PI	57	-795.1494	1748.0868	0.4150	0.9247
K2PIG	58	-794.9434	1751.5134	0.0748	0.9995
HKYG	60	-797.1877	1763.8349	0.0002	0.9996
K2P	56	-805.0711	1764.1422	0.0001	0.9998
HKYI	60	-797.5917	1764.6429	0.0001	0.9999
SYMG	61	-796.5122	1766.4801	0.0000	0.9999
SYMI	61	-796.9955	1767.4467	0.0000	1.0000
HKYIG	61	-797.1878	1767.8313	0.0000	1.0000
GTRG	64	-791.7500	1769.2778	0.0000	1.0000
GTRI	64	-792.2011	1770.1799	0.0000	1.0000
SYMIG	62	-796.5115	1770.5299	0.0000	1.0000
GTRIG	65	-791.7494	1773.4988	0.0000	1.0000
JC69G	56	-812.7982	1779.5965	0.0000	1.0000
JC69I	56	-813.1612	1780.3223	0.0000	1.0000
HKY	59	-808.2928	1782.1024	0.0000	1.0000
JC69IG	57	-812.7981	1783.3844	0.0000	1.0000
F81G	59	-809.0847	1783.6861	0.0000	1.0000
F81I	59	-809.3057	1784.1282	0.0000	1.0000
SYM	60	-807.6886	1784.8366	0.0000	1.0000
GTR	63	-802.3139	1786.2416	0.0000	1.0000
F81IG	60	-809.0845	1787.6284	0.0000	1.0000
JC69	55	-823.6747	1797.6108	0.0000	1.0000
F81	58	-819.0858	1799.7983	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
K2PG	57	-794.9437	1703.8873	0.4454	0.4454
K2PI	57	-795.1494	1704.2987	0.3626	0.8080
K2PIG	58	-794.9434	1705.8868	0.1639	0.9719
GTRG	64	-791.7500	1711.5001	0.0099	0.9818
GTRI	64	-792.2011	1712.4021	0.0063	0.9881
GTRIG	65	-791.7494	1713.4988	0.0036	0.9918
HKYG	60	-797.1877	1714.3754	0.0024	0.9941
SYMG	61	-796.5122	1715.0243	0.0017	0.9958
HKYI	60	-797.5917	1715.1835	0.0016	0.9974
SYMI	61	-796.9955	1715.9909	0.0010	0.9985
HKYIG	61	-797.1878	1716.3755	0.0009	0.9993
SYMIG	62	-796.5115	1717.0230	0.0006	1.0000
K2P	56	-805.0711	1722.1422	0.0000	1.0000
GTR	63	-802.3139	1730.6278	0.0000	1.0000
HKY	59	-808.2928	1734.5856	0.0000	1.0000
SYM	60	-807.6886	1735.3771	0.0000	1.0000
F81G	59	-809.0847	1736.1693	0.0000	1.0000
F81I	59	-809.3057	1736.6114	0.0000	1.0000
JC69G	56	-812.7982	1737.5965	0.0000	1.0000
F81IG	60	-809.0845	1738.1690	0.0000	1.0000
JC69I	56	-813.1612	1738.3223	0.0000	1.0000
JC69IG	57	-812.7981	1739.5963	0.0000	1.0000
F81	58	-819.0858	1754.1716	0.0000	1.0000
JC69	55	-823.6747	1757.3494	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-794.9437	1894.4004	0.5308	0.5308
K2PI	57	-795.1494	1894.8118	0.4321	0.9629
K2PIG	58	-794.9434	1899.7422	0.0367	0.9997
K2P	56	-805.0711	1909.3129	0.0003	1.0000
HKYG	60	-797.1877	1914.9155	0.0000	1.0000
HKYI	60	-797.5917	1915.7235	0.0000	1.0000
SYMG	61	-796.5122	1918.9067	0.0000	1.0000
SYMI	61	-796.9955	1919.8733	0.0000	1.0000
HKYIG	61	-797.1878	1920.2579	0.0000	1.0000
SYMIG	62	-796.5115	1924.2477	0.0000	1.0000
JC69G	56	-812.7982	1924.7672	0.0000	1.0000
GTRG	64	-791.7500	1925.4095	0.0000	1.0000
JC69I	56	-813.1612	1925.4930	0.0000	1.0000
GTRI	64	-792.2011	1926.3115	0.0000	1.0000
JC69IG	57	-812.7981	1930.1093	0.0000	1.0000
GTRIG	65	-791.7494	1930.7506	0.0000	1.0000
HKY	59	-808.2928	1931.7833	0.0000	1.0000
F81G	59	-809.0847	1933.3670	0.0000	1.0000
F81I	59	-809.3057	1933.8091	0.0000	1.0000
SYM	60	-807.6886	1935.9172	0.0000	1.0000
F81IG	60	-809.0845	1938.7090	0.0000	1.0000
JC69	55	-823.6747	1941.1778	0.0000	1.0000
GTR	63	-802.3139	1941.1949	0.0000	1.0000
F81	58	-819.0858	1948.0270	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 183-phylip/chr5_10180.phylip = 1 - 209;
 Partition Dummy = 1:183-phylip/chr5_10180.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;


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

BEGIN MRBAYES;
 Charset 183-phylip/chr5_10180.phylip = 1 - 209;
 Partition Dummy = 1:183-phylip/chr5_10180.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;


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

BEGIN MRBAYES;
 Charset 183-phylip/chr5_10180.phylip = 1 - 209;
 Partition Dummy = 1:183-phylip/chr5_10180.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

