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

Input data from file "591-tomodel-phylip/chr2_10510.phylip" (7 taxa, 732 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	-1935.7609	3912.7033	0.5666	0.5666
GTRI	20	-1936.6894	3914.5602	0.2239	0.7905
GTRIG	21	-1935.6968	3914.6950	0.2093	0.9998
HKYG	16	-1948.6011	3929.9631	0.0001	0.9999
HKYI	16	-1949.2794	3931.3197	0.0001	1.0000
HKYIG	17	-1948.5194	3931.8959	0.0000	1.0000
GTR	19	-1951.5871	3942.2416	0.0000	1.0000
SYMG	17	-1959.6717	3954.2006	0.0000	1.0000
SYMI	17	-1960.2762	3955.4096	0.0000	1.0000
SYMIG	18	-1959.5906	3956.1404	0.0000	1.0000
HKY	15	-1964.2749	3959.2202	0.0000	1.0000
K2PG	13	-1967.2396	3960.9862	0.0000	1.0000
K2PI	13	-1967.8470	3962.2010	0.0000	1.0000
K2PIG	14	-1967.1370	3962.8597	0.0000	1.0000
SYM	16	-1975.5304	3983.8216	0.0000	1.0000
F81G	15	-1978.1115	3986.8933	0.0000	1.0000
F81I	15	-1978.6744	3988.0192	0.0000	1.0000
F81IG	16	-1978.0399	3988.8406	0.0000	1.0000
K2P	12	-1983.5401	3991.5141	0.0000	1.0000
JC69G	12	-1994.9009	4014.2357	0.0000	1.0000
JC69I	12	-1995.4191	4015.2721	0.0000	1.0000
F81	14	-1993.4305	4015.4468	0.0000	1.0000
JC69IG	13	-1994.8184	4016.1437	0.0000	1.0000
JC69	11	-2010.5989	4043.5644	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-1935.7609	3911.5218	0.5594	0.5594
GTRI	20	-1936.6894	3913.3788	0.2210	0.7804
GTRIG	21	-1935.6968	3913.3936	0.2194	0.9998
HKYG	16	-1948.6011	3929.2023	0.0001	0.9999
HKYI	16	-1949.2794	3930.5589	0.0000	1.0000
HKYIG	17	-1948.5194	3931.0388	0.0000	1.0000
GTR	19	-1951.5871	3941.1742	0.0000	1.0000
SYMG	17	-1959.6717	3953.3435	0.0000	1.0000
SYMI	17	-1960.2762	3954.5524	0.0000	1.0000
SYMIG	18	-1959.5906	3955.1811	0.0000	1.0000
HKY	15	-1964.2749	3958.5498	0.0000	1.0000
K2PG	13	-1967.2396	3960.4792	0.0000	1.0000
K2PI	13	-1967.8470	3961.6941	0.0000	1.0000
K2PIG	14	-1967.1370	3962.2740	0.0000	1.0000
SYM	16	-1975.5304	3983.0608	0.0000	1.0000
F81G	15	-1978.1115	3986.2229	0.0000	1.0000
F81I	15	-1978.6744	3987.3488	0.0000	1.0000
F81IG	16	-1978.0399	3988.0797	0.0000	1.0000
K2P	12	-1983.5401	3991.0801	0.0000	1.0000
JC69G	12	-1994.9009	4013.8018	0.0000	1.0000
JC69I	12	-1995.4191	4014.8381	0.0000	1.0000
F81	14	-1993.4305	4014.8610	0.0000	1.0000
JC69IG	13	-1994.8184	4015.6367	0.0000	1.0000
JC69	11	-2010.5989	4043.1977	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1948.6011	4002.7348	0.3910	0.3910
GTRG	20	-1935.7609	4003.4375	0.2752	0.6662
HKYI	16	-1949.2794	4004.0914	0.1984	0.8646
GTRI	20	-1936.6894	4005.2944	0.1087	0.9734
HKYIG	17	-1948.5194	4009.1670	0.0157	0.9891
GTRIG	21	-1935.6968	4009.9050	0.0108	0.9999
K2PG	13	-1967.2396	4020.2243	0.0001	1.0000
K2PI	13	-1967.8470	4021.4392	0.0000	1.0000
K2PIG	14	-1967.1370	4026.6149	0.0000	1.0000
HKY	15	-1964.2749	4027.4865	0.0000	1.0000
GTR	19	-1951.5871	4028.4940	0.0000	1.0000
SYMG	17	-1959.6717	4031.4718	0.0000	1.0000
SYMI	17	-1960.2762	4032.6807	0.0000	1.0000
SYMIG	18	-1959.5906	4037.9052	0.0000	1.0000
K2P	12	-1983.5401	4046.2295	0.0000	1.0000
F81G	15	-1978.1115	4055.1596	0.0000	1.0000
F81I	15	-1978.6744	4056.2855	0.0000	1.0000
SYM	16	-1975.5304	4056.5932	0.0000	1.0000
F81IG	16	-1978.0399	4061.6122	0.0000	1.0000
JC69G	12	-1994.9009	4068.9512	0.0000	1.0000
JC69I	12	-1995.4191	4069.9875	0.0000	1.0000
JC69IG	13	-1994.8184	4075.3819	0.0000	1.0000
F81	14	-1993.4305	4079.2019	0.0000	1.0000
JC69	11	-2010.5989	4093.7513	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/chr2_10510.phylip = 1 - 732;
 Partition Dummy = 1:591-tomodel-phylip/chr2_10510.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/chr2_10510.phylip = 1 - 732;
 Partition Dummy = 1:591-tomodel-phylip/chr2_10510.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/chr2_10510.phylip = 1 - 732;
 Partition Dummy = 1:591-tomodel-phylip/chr2_10510.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

