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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: GTRG


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	-1926.5270	3894.3894	0.5447	0.5447
GTRIG	21	-1926.0927	3895.6568	0.2891	0.8338
GTRI	20	-1927.7146	3896.7646	0.1661	0.9999
HKYG	16	-1940.1359	3913.1313	0.0000	1.0000
HKYIG	17	-1939.8533	3914.6750	0.0000	1.0000
HKYI	16	-1941.0036	3914.8666	0.0000	1.0000
SYMG	17	-1943.5453	3922.0590	0.0000	1.0000
SYMI	17	-1944.2605	3923.4893	0.0000	1.0000
SYMIG	18	-1943.2617	3923.6074	0.0000	1.0000
K2PG	13	-1954.4035	3935.3793	0.0000	1.0000
K2PI	13	-1955.0443	3936.6609	0.0000	1.0000
K2PIG	14	-1954.2127	3937.0869	0.0000	1.0000
GTR	19	-1960.9246	3961.0555	0.0000	1.0000
F81G	15	-1969.7808	3970.3186	0.0000	1.0000
F81I	15	-1970.3693	3971.4957	0.0000	1.0000
F81IG	16	-1969.7034	3972.2662	0.0000	1.0000
HKY	15	-1973.9907	3978.7386	0.0000	1.0000
SYM	16	-1976.8694	3986.5982	0.0000	1.0000
JC69G	12	-1983.5992	3991.6881	0.0000	1.0000
JC69I	12	-1983.9910	3992.4717	0.0000	1.0000
JC69IG	13	-1983.5954	3993.7630	0.0000	1.0000
K2P	12	-1988.7121	4001.9141	0.0000	1.0000
F81	14	-2002.9913	4034.6440	0.0000	1.0000
JC69	11	-2016.4531	4055.3200	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-1926.5270	3893.0540	0.5339	0.5339
GTRIG	21	-1926.0927	3894.1854	0.3032	0.8371
GTRI	20	-1927.7146	3895.4291	0.1628	0.9999
HKYG	16	-1940.1359	3912.2719	0.0000	1.0000
HKYIG	17	-1939.8533	3913.7067	0.0000	1.0000
HKYI	16	-1941.0036	3914.0072	0.0000	1.0000
SYMG	17	-1943.5453	3921.0906	0.0000	1.0000
SYMI	17	-1944.2605	3922.5209	0.0000	1.0000
SYMIG	18	-1943.2617	3922.5234	0.0000	1.0000
K2PG	13	-1954.4035	3934.8070	0.0000	1.0000
K2PI	13	-1955.0443	3936.0886	0.0000	1.0000
K2PIG	14	-1954.2127	3936.4255	0.0000	1.0000
GTR	19	-1960.9246	3959.8491	0.0000	1.0000
F81G	15	-1969.7808	3969.5615	0.0000	1.0000
F81I	15	-1970.3693	3970.7386	0.0000	1.0000
F81IG	16	-1969.7034	3971.4068	0.0000	1.0000
HKY	15	-1973.9907	3977.9815	0.0000	1.0000
SYM	16	-1976.8694	3985.7388	0.0000	1.0000
JC69G	12	-1983.5992	3991.1983	0.0000	1.0000
JC69I	12	-1983.9910	3991.9819	0.0000	1.0000
JC69IG	13	-1983.5954	3993.1907	0.0000	1.0000
K2P	12	-1988.7121	4001.4243	0.0000	1.0000
F81	14	-2002.9913	4033.9826	0.0000	1.0000
JC69	11	-2016.4531	4054.9062	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	20	-1926.5270	3982.5934	0.4673	0.4673
HKYG	16	-1940.1359	3983.9034	0.2428	0.7101
GTRI	20	-1927.7146	3984.9686	0.1425	0.8526
HKYI	16	-1941.0036	3985.6388	0.1019	0.9546
GTRIG	21	-1926.0927	3988.2019	0.0283	0.9829
HKYIG	17	-1939.8533	3989.8152	0.0126	0.9955
K2PG	13	-1954.4035	3993.0076	0.0026	0.9980
K2PI	13	-1955.0443	3994.2892	0.0013	0.9994
SYMG	17	-1943.5453	3997.1991	0.0003	0.9997
SYMI	17	-1944.2605	3998.6295	0.0002	0.9999
K2PIG	14	-1954.2127	3999.1031	0.0001	1.0000
SYMIG	18	-1943.2617	4003.1089	0.0000	1.0000
F81G	15	-1969.7808	4036.7161	0.0000	1.0000
F81I	15	-1970.3693	4037.8932	0.0000	1.0000
F81IG	16	-1969.7034	4043.0383	0.0000	1.0000
GTR	19	-1960.9246	4044.9116	0.0000	1.0000
JC69G	12	-1983.5992	4044.9220	0.0000	1.0000
HKY	15	-1973.9907	4045.1361	0.0000	1.0000
JC69I	12	-1983.9910	4045.7056	0.0000	1.0000
JC69IG	13	-1983.5954	4051.3914	0.0000	1.0000
K2P	12	-1988.7121	4055.1479	0.0000	1.0000
SYM	16	-1976.8694	4057.3704	0.0000	1.0000
F81	14	-2002.9913	4096.6602	0.0000	1.0000
JC69	11	-2016.4531	4104.1529	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/chr3_3118.phylip = 1 - 650;
 Partition Dummy = 1:591-tomodel-phylip/chr3_3118.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/chr3_3118.phylip = 1 - 650;
 Partition Dummy = 1:591-tomodel-phylip/chr3_3118.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 (GTRG)]

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


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

