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

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

Minimum AIC  model: GTRG
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	-755.3773	1686.6324	0.7582	0.7582
HKYIG	61	-755.3401	1690.8649	0.0913	0.8495
HKYI	60	-757.6759	1691.2297	0.0761	0.9256
K2PG	57	-764.2975	1691.9382	0.0534	0.9790
GTRG	64	-751.0731	1695.6581	0.0083	0.9873
K2PIG	58	-764.2128	1695.8842	0.0074	0.9948
K2PI	57	-767.2299	1697.8030	0.0028	0.9976
GTRI	64	-752.8627	1699.2372	0.0014	0.9990
GTRIG	65	-751.0713	1700.2378	0.0008	0.9998
SYMG	61	-762.2481	1704.6807	0.0001	0.9999
HKY	59	-767.3370	1706.3103	0.0000	1.0000
SYMIG	62	-762.2168	1708.9918	0.0000	1.0000
SYMI	61	-764.7023	1709.5892	0.0000	1.0000
K2P	56	-776.3930	1712.0749	0.0000	1.0000
GTR	63	-763.1408	1715.2816	0.0000	1.0000
SYM	60	-773.5950	1723.0678	0.0000	1.0000
F81G	59	-782.0857	1735.8078	0.0000	1.0000
JC69G	56	-788.4648	1736.2184	0.0000	1.0000
F81IG	60	-782.0346	1739.9471	0.0000	1.0000
JC69IG	57	-788.3985	1740.1402	0.0000	1.0000
F81I	59	-784.7196	1741.0755	0.0000	1.0000
JC69I	56	-791.3369	1741.9627	0.0000	1.0000
F81	58	-793.6654	1754.7894	0.0000	1.0000
JC69	55	-800.1986	1755.6913	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-751.0731	1630.1462	0.3799	0.3799
HKYG	60	-755.3773	1630.7546	0.2803	0.6601
GTRIG	65	-751.0713	1632.1426	0.1400	0.8001
HKYIG	61	-755.3401	1632.6802	0.1070	0.9072
GTRI	64	-752.8627	1633.7254	0.0635	0.9706
HKYI	60	-757.6759	1635.3519	0.0281	0.9987
K2PG	57	-764.2975	1642.5949	0.0008	0.9995
K2PIG	58	-764.2128	1644.4256	0.0003	0.9998
SYMG	61	-762.2481	1646.4961	0.0001	0.9999
SYMIG	62	-762.2168	1648.4336	0.0000	0.9999
K2PI	57	-767.2299	1648.4597	0.0000	1.0000
SYMI	61	-764.7023	1651.4045	0.0000	1.0000
GTR	63	-763.1408	1652.2816	0.0000	1.0000
HKY	59	-767.3370	1652.6739	0.0000	1.0000
K2P	56	-776.3930	1664.7860	0.0000	1.0000
SYM	60	-773.5950	1667.1900	0.0000	1.0000
F81G	59	-782.0857	1682.1714	0.0000	1.0000
F81IG	60	-782.0346	1684.0693	0.0000	1.0000
F81I	59	-784.7196	1687.4392	0.0000	1.0000
JC69G	56	-788.4648	1688.9295	0.0000	1.0000
JC69IG	57	-788.3985	1690.7969	0.0000	1.0000
JC69I	56	-791.3369	1694.6738	0.0000	1.0000
F81	58	-793.6654	1703.3307	0.0000	1.0000
JC69	55	-800.1986	1710.3972	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-755.3773	1826.2043	0.6327	0.6327
K2PG	57	-764.2975	1828.2722	0.2250	0.8577
HKYI	60	-757.6759	1830.8016	0.0635	0.9212
HKYIG	61	-755.3401	1831.3875	0.0474	0.9686
K2PIG	58	-764.2128	1833.3603	0.0177	0.9863
K2PI	57	-767.2299	1834.1370	0.0120	0.9983
GTRG	64	-751.0731	1838.6259	0.0013	0.9996
GTRI	64	-752.8627	1842.2051	0.0002	0.9998
GTRIG	65	-751.0713	1843.8798	0.0001	0.9999
HKY	59	-767.3370	1844.8661	0.0001	0.9999
SYMG	61	-762.2481	1845.2033	0.0000	1.0000
K2P	56	-776.3930	1847.2057	0.0000	1.0000
SYMI	61	-764.7023	1850.1118	0.0000	1.0000
SYMIG	62	-762.2168	1850.3984	0.0000	1.0000
GTR	63	-763.1408	1857.5038	0.0000	1.0000
SYM	60	-773.5950	1862.6397	0.0000	1.0000
JC69G	56	-788.4648	1871.3493	0.0000	1.0000
F81G	59	-782.0857	1874.3636	0.0000	1.0000
JC69IG	57	-788.3985	1876.4742	0.0000	1.0000
JC69I	56	-791.3369	1877.0936	0.0000	1.0000
F81IG	60	-782.0346	1879.5190	0.0000	1.0000
F81I	59	-784.7196	1879.6314	0.0000	1.0000
JC69	55	-800.1986	1889.5594	0.0000	1.0000
F81	58	-793.6654	1892.2655	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 183-phylip/chr2_8593.phylip = 1 - 192;
 Partition Dummy = 1:183-phylip/chr2_8593.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 Prset applyto=(1) revmatpr=Dirichlet(1.0,1.0,1.0,1.0,1.0,1.0) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

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

