
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
16:11 04/04/2011

Input data from file "683-tomodel-phylip/chr2_22024.phylip" (7 taxa, 754 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	-1934.7155	3910.5770	0.4828	0.4828
GTRI	20	-1935.2142	3911.5744	0.2932	0.7760
GTRIG	21	-1934.4341	3912.1304	0.2221	0.9981
HKYG	16	-1945.1990	3923.1362	0.0009	0.9990
HKYI	16	-1945.6539	3924.0460	0.0006	0.9996
HKYIG	17	-1944.9067	3924.6450	0.0004	1.0000
SYMG	17	-1955.0828	3944.9970	0.0000	1.0000
SYMI	17	-1955.4485	3945.7284	0.0000	1.0000
SYMIG	18	-1954.8537	3946.6380	0.0000	1.0000
K2PG	13	-1961.7201	3949.9321	0.0000	1.0000
K2PI	13	-1962.2929	3951.0778	0.0000	1.0000
K2PIG	14	-1961.4199	3951.4082	0.0000	1.0000
GTR	19	-1958.8897	3956.8147	0.0000	1.0000
HKY	15	-1967.1109	3964.8722	0.0000	1.0000
F81G	15	-1971.6009	3973.8521	0.0000	1.0000
F81I	15	-1972.1618	3974.9739	0.0000	1.0000
F81IG	16	-1971.3889	3975.5159	0.0000	1.0000
SYM	16	-1976.9068	3986.5516	0.0000	1.0000
K2P	12	-1982.6449	3989.7108	0.0000	1.0000
JC69G	12	-1985.3583	3995.1377	0.0000	1.0000
JC69I	12	-1985.9515	3996.3240	0.0000	1.0000
JC69IG	13	-1985.1362	3996.7643	0.0000	1.0000
F81	14	-1990.9094	4010.3872	0.0000	1.0000
JC69	11	-2004.4241	4031.2039	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-1934.7155	3909.4310	0.4766	0.4766
GTRI	20	-1935.2142	3910.4284	0.2895	0.7661
GTRIG	21	-1934.4341	3910.8681	0.2323	0.9985
HKYG	16	-1945.1990	3922.3980	0.0007	0.9992
HKYI	16	-1945.6539	3923.3078	0.0005	0.9996
HKYIG	17	-1944.9067	3923.8135	0.0004	1.0000
SYMG	17	-1955.0828	3944.1655	0.0000	1.0000
SYMI	17	-1955.4485	3944.8969	0.0000	1.0000
SYMIG	18	-1954.8537	3945.7074	0.0000	1.0000
K2PG	13	-1961.7201	3949.4402	0.0000	1.0000
K2PI	13	-1962.2929	3950.5859	0.0000	1.0000
K2PIG	14	-1961.4199	3950.8399	0.0000	1.0000
GTR	19	-1958.8897	3955.7793	0.0000	1.0000
HKY	15	-1967.1109	3964.2218	0.0000	1.0000
F81G	15	-1971.6009	3973.2017	0.0000	1.0000
F81I	15	-1972.1618	3974.3235	0.0000	1.0000
F81IG	16	-1971.3889	3974.7778	0.0000	1.0000
SYM	16	-1976.9068	3985.8135	0.0000	1.0000
K2P	12	-1982.6449	3989.2897	0.0000	1.0000
JC69G	12	-1985.3583	3994.7167	0.0000	1.0000
JC69I	12	-1985.9515	3995.9030	0.0000	1.0000
JC69IG	13	-1985.1362	3996.2724	0.0000	1.0000
F81	14	-1990.9094	4009.8189	0.0000	1.0000
JC69	11	-2004.4241	4030.8481	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1945.1990	3996.4043	0.5588	0.5588
HKYI	16	-1945.6539	3997.3141	0.3546	0.9134
GTRG	20	-1934.7155	4001.9388	0.0351	0.9485
HKYIG	17	-1944.9067	4002.4452	0.0273	0.9757
GTRI	20	-1935.2142	4002.9362	0.0213	0.9971
GTRIG	21	-1934.4341	4008.0014	0.0017	0.9988
K2PG	13	-1961.7201	4009.5703	0.0008	0.9995
K2PI	13	-1962.2929	4010.7160	0.0004	1.0000
K2PIG	14	-1961.4199	4015.5954	0.0000	1.0000
SYMG	17	-1955.0828	4022.7972	0.0000	1.0000
SYMI	17	-1955.4485	4023.5286	0.0000	1.0000
SYMIG	18	-1954.8537	4028.9645	0.0000	1.0000
HKY	15	-1967.1109	4033.6027	0.0000	1.0000
F81G	15	-1971.6009	4042.5826	0.0000	1.0000
GTR	19	-1958.8897	4043.6618	0.0000	1.0000
F81I	15	-1972.1618	4043.7044	0.0000	1.0000
K2P	12	-1982.6449	4044.7944	0.0000	1.0000
F81IG	16	-1971.3889	4048.7841	0.0000	1.0000
JC69G	12	-1985.3583	4050.2214	0.0000	1.0000
JC69I	12	-1985.9515	4051.4077	0.0000	1.0000
JC69IG	13	-1985.1362	4056.4025	0.0000	1.0000
SYM	16	-1976.9068	4059.8198	0.0000	1.0000
F81	14	-1990.9094	4074.5744	0.0000	1.0000
JC69	11	-2004.4241	4081.7274	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 683-tomodel-phylip/chr2_22024.phylip = 1 - 754;
 Partition Dummy = 1:683-tomodel-phylip/chr2_22024.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 683-tomodel-phylip/chr2_22024.phylip = 1 - 754;
 Partition Dummy = 1:683-tomodel-phylip/chr2_22024.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 683-tomodel-phylip/chr2_22024.phylip = 1 - 754;
 Partition Dummy = 1:683-tomodel-phylip/chr2_22024.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

