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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
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
GTRG	20	-1879.5164	3800.1029	0.4656	0.4656
GTRI	20	-1879.9420	3800.9541	0.3042	0.7698
GTRIG	21	-1879.1679	3801.5144	0.2299	0.9997
HKYG	16	-1891.7729	3816.2353	0.0001	0.9998
HKYI	16	-1892.3032	3817.2959	0.0001	0.9999
HKYIG	17	-1891.3142	3817.4051	0.0001	1.0000
GTR	19	-1896.6002	3832.1674	0.0000	1.0000
HKY	15	-1908.6211	3847.8499	0.0000	1.0000
SYMG	17	-1913.7058	3862.1883	0.0000	1.0000
SYMI	17	-1914.1287	3863.0341	0.0000	1.0000
SYMIG	18	-1913.4889	3863.8469	0.0000	1.0000
K2PG	13	-1921.6006	3869.6608	0.0000	1.0000
K2PI	13	-1922.0876	3870.6348	0.0000	1.0000
K2PIG	14	-1921.2786	3871.0881	0.0000	1.0000
F81G	15	-1924.7314	3880.0704	0.0000	1.0000
F81I	15	-1925.1882	3880.9839	0.0000	1.0000
F81IG	16	-1924.4659	3881.6212	0.0000	1.0000
SYM	16	-1928.1225	3888.9346	0.0000	1.0000
K2P	12	-1936.8965	3898.1865	0.0000	1.0000
F81	14	-1939.7245	3907.9800	0.0000	1.0000
JC69G	12	-1950.3754	3925.1443	0.0000	1.0000
JC69I	12	-1950.8236	3926.0406	0.0000	1.0000
JC69IG	13	-1950.1418	3926.7431	0.0000	1.0000
JC69	11	-1964.8107	3951.9540	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-1879.5164	3799.0328	0.4597	0.4597
GTRI	20	-1879.9420	3799.8841	0.3004	0.7601
GTRIG	21	-1879.1679	3800.3358	0.2396	0.9997
HKYG	16	-1891.7729	3815.5458	0.0001	0.9999
HKYI	16	-1892.3032	3816.6064	0.0001	0.9999
HKYIG	17	-1891.3142	3816.6285	0.0001	1.0000
GTR	19	-1896.6002	3831.2005	0.0000	1.0000
HKY	15	-1908.6211	3847.2423	0.0000	1.0000
SYMG	17	-1913.7058	3861.4117	0.0000	1.0000
SYMI	17	-1914.1287	3862.2574	0.0000	1.0000
SYMIG	18	-1913.4889	3862.9778	0.0000	1.0000
K2PG	13	-1921.6006	3869.2012	0.0000	1.0000
K2PI	13	-1922.0876	3870.1752	0.0000	1.0000
K2PIG	14	-1921.2786	3870.5572	0.0000	1.0000
F81G	15	-1924.7314	3879.4628	0.0000	1.0000
F81I	15	-1925.1882	3880.3763	0.0000	1.0000
F81IG	16	-1924.4659	3880.9317	0.0000	1.0000
SYM	16	-1928.1225	3888.2451	0.0000	1.0000
K2P	12	-1936.8965	3897.7931	0.0000	1.0000
F81	14	-1939.7245	3907.4490	0.0000	1.0000
JC69G	12	-1950.3754	3924.7508	0.0000	1.0000
JC69I	12	-1950.8236	3925.6471	0.0000	1.0000
JC69IG	13	-1950.1418	3926.2835	0.0000	1.0000
JC69	11	-1964.8107	3951.6215	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1891.7729	3890.6191	0.4555	0.4555
HKYI	16	-1892.3032	3891.6798	0.2680	0.7234
GTRG	20	-1879.5164	3892.8745	0.1475	0.8709
GTRI	20	-1879.9420	3893.7258	0.0963	0.9673
HKYIG	17	-1891.3142	3896.3939	0.0254	0.9926
GTRIG	21	-1879.1679	3898.8696	0.0074	1.0000
HKY	15	-1908.6211	3917.6236	0.0000	1.0000
GTR	19	-1896.6002	3920.3501	0.0000	1.0000
K2PG	13	-1921.6006	3930.1983	0.0000	1.0000
K2PI	13	-1922.0876	3931.1723	0.0000	1.0000
K2PIG	14	-1921.2786	3936.2463	0.0000	1.0000
SYMG	17	-1913.7058	3941.1771	0.0000	1.0000
SYMI	17	-1914.1287	3942.0228	0.0000	1.0000
SYMIG	18	-1913.4889	3947.4353	0.0000	1.0000
F81G	15	-1924.7314	3949.8441	0.0000	1.0000
F81I	15	-1925.1882	3950.7576	0.0000	1.0000
K2P	12	-1936.8965	3954.0981	0.0000	1.0000
F81IG	16	-1924.4659	3956.0051	0.0000	1.0000
SYM	16	-1928.1225	3963.3184	0.0000	1.0000
F81	14	-1939.7245	3973.1382	0.0000	1.0000
JC69G	12	-1950.3754	3981.0558	0.0000	1.0000
JC69I	12	-1950.8236	3981.9521	0.0000	1.0000
JC69IG	13	-1950.1418	3987.2806	0.0000	1.0000
JC69	11	-1964.8107	4003.2344	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 591-tomodel-phylip/chr9_6343.phylip = 1 - 806;
 Partition Dummy = 1:591-tomodel-phylip/chr9_6343.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 (GTRG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr9_6343.phylip = 1 - 806;
 Partition Dummy = 1:591-tomodel-phylip/chr9_6343.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 BIC model (HKYG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr9_6343.phylip = 1 - 806;
 Partition Dummy = 1:591-tomodel-phylip/chr9_6343.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

