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

Input data from file "183-phylip/chr1_5365.phylip" (29 taxa, 281 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	60	-761.4659	1676.2046	0.4707	0.4707
HKYI	60	-761.6302	1676.5332	0.3994	0.8700
HKYIG	61	-761.1974	1678.9336	0.1203	0.9903
GTRG	64	-759.4493	1685.4172	0.0047	0.9950
GTRI	64	-759.7036	1685.9258	0.0036	0.9986
GTRIG	65	-759.0958	1688.0985	0.0012	0.9999
K2PG	57	-775.2289	1694.1079	0.0001	0.9999
K2PI	57	-775.3016	1694.2533	0.0001	1.0000
K2PIG	58	-774.9643	1696.7575	0.0000	1.0000
HKY	59	-776.4948	1703.0257	0.0000	1.0000
SYMG	61	-774.2249	1704.9886	0.0000	1.0000
SYMI	61	-774.3230	1705.1848	0.0000	1.0000
SYMIG	62	-773.9421	1707.7191	0.0000	1.0000
F81G	59	-778.8991	1707.8343	0.0000	1.0000
F81I	59	-779.0195	1708.0752	0.0000	1.0000
F81IG	60	-778.6430	1710.5588	0.0000	1.0000
GTR	63	-774.0380	1711.2373	0.0000	1.0000
JC69G	56	-791.1237	1722.7474	0.0000	1.0000
JC69I	56	-791.2071	1722.9143	0.0000	1.0000
K2P	56	-791.2302	1722.9604	0.0000	1.0000
JC69IG	57	-790.8804	1725.4110	0.0000	1.0000
SYM	60	-789.8453	1732.9632	0.0000	1.0000
F81	58	-794.1114	1735.0516	0.0000	1.0000
JC69	55	-806.8572	1751.0921	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-761.4659	1642.9319	0.3779	0.3779
HKYI	60	-761.6302	1643.2605	0.3207	0.6986
HKYIG	61	-761.1974	1644.3947	0.1819	0.8804
GTRG	64	-759.4493	1646.8987	0.0520	0.9324
GTRI	64	-759.7036	1647.4073	0.0403	0.9727
GTRIG	65	-759.0958	1648.1916	0.0272	1.0000
K2PG	57	-775.2289	1664.4577	0.0000	1.0000
K2PI	57	-775.3016	1664.6031	0.0000	1.0000
K2PIG	58	-774.9643	1665.9286	0.0000	1.0000
SYMG	61	-774.2249	1670.4497	0.0000	1.0000
SYMI	61	-774.3230	1670.6460	0.0000	1.0000
HKY	59	-776.4948	1670.9895	0.0000	1.0000
SYMIG	62	-773.9421	1671.8843	0.0000	1.0000
GTR	63	-774.0380	1674.0760	0.0000	1.0000
F81G	59	-778.8991	1675.7981	0.0000	1.0000
F81I	59	-779.0195	1676.0390	0.0000	1.0000
F81IG	60	-778.6430	1677.2861	0.0000	1.0000
JC69G	56	-791.1237	1694.2474	0.0000	1.0000
JC69I	56	-791.2071	1694.4143	0.0000	1.0000
K2P	56	-791.2302	1694.4604	0.0000	1.0000
JC69IG	57	-790.8804	1695.7608	0.0000	1.0000
SYM	60	-789.8453	1699.6905	0.0000	1.0000
F81	58	-794.1114	1704.2228	0.0000	1.0000
JC69	55	-806.8572	1723.7144	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-761.4659	1861.2332	0.5163	0.5163
HKYI	60	-761.6302	1861.5617	0.4381	0.9545
HKYIG	61	-761.1974	1866.3344	0.0403	0.9948
K2PG	57	-775.2289	1871.8439	0.0026	0.9973
K2PI	57	-775.3016	1871.9893	0.0024	0.9997
K2PIG	58	-774.9643	1876.9532	0.0002	0.9999
GTRG	64	-759.4493	1879.7534	0.0000	1.0000
GTRI	64	-759.7036	1880.2620	0.0000	1.0000
GTRIG	65	-759.0958	1884.6846	0.0000	1.0000
HKY	59	-776.4948	1885.6524	0.0000	1.0000
F81G	59	-778.8991	1890.4611	0.0000	1.0000
F81I	59	-779.0195	1890.7019	0.0000	1.0000
SYMG	61	-774.2249	1892.3894	0.0000	1.0000
SYMI	61	-774.3230	1892.5856	0.0000	1.0000
F81IG	60	-778.6430	1895.5874	0.0000	1.0000
SYMIG	62	-773.9421	1897.4622	0.0000	1.0000
JC69G	56	-791.1237	1897.9953	0.0000	1.0000
JC69I	56	-791.2071	1898.1621	0.0000	1.0000
K2P	56	-791.2302	1898.2082	0.0000	1.0000
JC69IG	57	-790.8804	1903.1470	0.0000	1.0000
GTR	63	-774.0380	1903.2923	0.0000	1.0000
F81	58	-794.1114	1915.2474	0.0000	1.0000
SYM	60	-789.8453	1917.9918	0.0000	1.0000
JC69	55	-806.8572	1923.8239	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 183-phylip/chr1_5365.phylip = 1 - 281;
 Partition Dummy = 1:183-phylip/chr1_5365.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 183-phylip/chr1_5365.phylip = 1 - 281;
 Partition Dummy = 1:183-phylip/chr1_5365.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 183-phylip/chr1_5365.phylip = 1 - 281;
 Partition Dummy = 1:183-phylip/chr1_5365.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

