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

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

Minimum AIC  model: HKY
Minimum AICc model: HKY
Minimum BIC  model: K2P


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
HKY	59	-788.0161	1727.9078	0.5729	0.5729
HKYI	60	-787.4976	1730.1875	0.1832	0.7561
HKYG	60	-787.7335	1730.6593	0.1447	0.9009
HKYIG	61	-787.5069	1733.5549	0.0340	0.9349
K2P	56	-796.0253	1734.1638	0.0251	0.9600
GTR	63	-784.9215	1735.1796	0.0151	0.9751
K2PI	57	-795.7492	1736.8349	0.0066	0.9817
K2PG	57	-795.9150	1737.1665	0.0056	0.9873
GTRI	64	-784.4259	1737.6362	0.0044	0.9917
GTRG	64	-784.6802	1738.1447	0.0034	0.9951
SYM	60	-792.0667	1739.3258	0.0019	0.9970
K2PIG	58	-795.7584	1740.1072	0.0013	0.9983
GTRIG	65	-784.4360	1741.1380	0.0008	0.9991
SYMI	61	-791.8447	1742.2304	0.0004	0.9995
SYMG	61	-791.9979	1742.5369	0.0004	0.9999
SYMIG	62	-791.8538	1745.6299	0.0001	1.0000
F81	58	-804.7597	1758.1099	0.0000	1.0000
F81I	59	-804.3579	1760.5914	0.0000	1.0000
F81G	59	-804.5717	1761.0190	0.0000	1.0000
JC69	55	-811.9162	1762.7526	0.0000	1.0000
F81IG	60	-804.3681	1763.9285	0.0000	1.0000
JC69I	56	-811.6496	1765.4124	0.0000	1.0000
JC69G	56	-811.8195	1765.7521	0.0000	1.0000
JC69IG	57	-811.6590	1768.6544	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKY	59	-788.0161	1694.0322	0.3044	0.3044
HKYI	60	-787.4976	1694.9952	0.1881	0.4925
HKYG	60	-787.7335	1695.4670	0.1486	0.6411
GTR	63	-784.9215	1695.8430	0.1231	0.7642
GTRI	64	-784.4259	1696.8519	0.0743	0.8386
HKYIG	61	-787.5069	1697.0139	0.0686	0.9071
GTRG	64	-784.6802	1697.3604	0.0576	0.9648
GTRIG	65	-784.4360	1698.8720	0.0271	0.9918
K2P	56	-796.0253	1704.0506	0.0020	0.9939
SYM	60	-792.0667	1704.1335	0.0019	0.9958
K2PI	57	-795.7492	1705.4984	0.0010	0.9968
SYMI	61	-791.8447	1705.6893	0.0009	0.9977
K2PG	57	-795.9150	1705.8300	0.0008	0.9985
SYMG	61	-791.9979	1705.9958	0.0008	0.9993
K2PIG	58	-795.7584	1707.5168	0.0004	0.9997
SYMIG	62	-791.8538	1707.7076	0.0003	1.0000
F81	58	-804.7597	1725.5195	0.0000	1.0000
F81I	59	-804.3579	1726.7158	0.0000	1.0000
F81G	59	-804.5717	1727.1434	0.0000	1.0000
F81IG	60	-804.3681	1728.7362	0.0000	1.0000
JC69	55	-811.9162	1733.8324	0.0000	1.0000
JC69I	56	-811.6496	1735.2992	0.0000	1.0000
JC69G	56	-811.8195	1735.6389	0.0000	1.0000
JC69IG	57	-811.6590	1737.3179	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2P	56	-796.0253	1905.3544	0.5088	0.5088
HKY	59	-788.0161	1906.1201	0.3469	0.8557
K2PI	57	-795.7492	1910.3970	0.0409	0.8966
HKYI	60	-787.4976	1910.6778	0.0355	0.9321
K2PG	57	-795.9150	1910.7285	0.0346	0.9668
HKYG	60	-787.7335	1911.1496	0.0281	0.9948
K2PIG	58	-795.7584	1916.0100	0.0025	0.9973
HKYIG	61	-787.5069	1916.2913	0.0021	0.9995
SYM	60	-792.0667	1919.8162	0.0004	0.9998
GTR	63	-784.9215	1922.3098	0.0001	0.9999
SYMI	61	-791.8447	1924.9667	0.0000	1.0000
SYMG	61	-791.9979	1925.2732	0.0000	1.0000
GTRI	64	-784.4259	1926.9134	0.0000	1.0000
GTRG	64	-784.6802	1927.4219	0.0000	1.0000
SYMIG	62	-791.8538	1930.5797	0.0000	1.0000
JC69	55	-811.9162	1931.5415	0.0000	1.0000
GTRIG	65	-784.4360	1932.5283	0.0000	1.0000
F81	58	-804.7597	1934.0127	0.0000	1.0000
JC69I	56	-811.6496	1936.6031	0.0000	1.0000
JC69G	56	-811.8195	1936.9428	0.0000	1.0000
F81I	59	-804.3579	1938.8038	0.0000	1.0000
F81G	59	-804.5717	1939.2314	0.0000	1.0000
JC69IG	57	-811.6590	1942.2164	0.0000	1.0000
F81IG	60	-804.3681	1944.4188	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 183-phylip/chr5_11231.phylip = 1 - 269;
 Partition Dummy = 1:183-phylip/chr5_11231.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 183-phylip/chr5_11231.phylip = 1 - 269;
 Partition Dummy = 1:183-phylip/chr5_11231.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 183-phylip/chr5_11231.phylip = 1 - 269;
 Partition Dummy = 1:183-phylip/chr5_11231.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0);
END;


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

