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

Input data from file "591-tomodel-phylip/chr1_32388.phylip" (7 taxa, 655 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	-1924.7472	3890.8194	0.5529	0.5529
GTRI	20	-1925.7479	3892.8208	0.2032	0.7561
GTRIG	21	-1924.7264	3892.9125	0.1941	0.9503
HKYG	16	-1931.9748	3896.8022	0.0278	0.9780
HKYI	16	-1932.7999	3898.4524	0.0122	0.9902
HKYIG	17	-1931.9613	3898.8834	0.0098	1.0000
GTR	19	-1945.8204	3930.8376	0.0000	1.0000
HKY	15	-1953.1391	3937.0294	0.0000	1.0000
SYMG	17	-1961.8920	3958.7448	0.0000	1.0000
SYMI	17	-1962.3221	3959.6050	0.0000	1.0000
SYMIG	18	-1961.8967	3960.8689	0.0000	1.0000
K2PG	13	-1967.5437	3961.6552	0.0000	1.0000
K2PI	13	-1967.9628	3962.4934	0.0000	1.0000
K2PIG	14	-1967.5496	3963.7555	0.0000	1.0000
F81G	15	-1982.3484	3995.4480	0.0000	1.0000
SYM	16	-1981.4646	3995.7819	0.0000	1.0000
F81I	15	-1982.8348	3996.4209	0.0000	1.0000
F81IG	16	-1982.3673	3997.5873	0.0000	1.0000
K2P	12	-1987.6358	3999.7576	0.0000	1.0000
F81	14	-2000.8021	4030.2605	0.0000	1.0000
JC69G	12	-2010.4900	4045.4660	0.0000	1.0000
JC69I	12	-2010.7967	4046.0794	0.0000	1.0000
JC69IG	13	-2010.5012	4047.5703	0.0000	1.0000
JC69	11	-2028.6764	4079.7633	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-1924.7472	3889.4944	0.5509	0.5509
GTRIG	21	-1924.7264	3891.4527	0.2070	0.7579
GTRI	20	-1925.7479	3891.4959	0.2025	0.9604
HKYG	16	-1931.9748	3895.9495	0.0218	0.9823
HKYI	16	-1932.7999	3897.5997	0.0096	0.9919
HKYIG	17	-1931.9613	3897.9226	0.0081	1.0000
GTR	19	-1945.8204	3929.6408	0.0000	1.0000
HKY	15	-1953.1391	3936.2782	0.0000	1.0000
SYMG	17	-1961.8920	3957.7841	0.0000	1.0000
SYMI	17	-1962.3221	3958.6443	0.0000	1.0000
SYMIG	18	-1961.8967	3959.7935	0.0000	1.0000
K2PG	13	-1967.5437	3961.0873	0.0000	1.0000
K2PI	13	-1967.9628	3961.9255	0.0000	1.0000
K2PIG	14	-1967.5496	3963.0992	0.0000	1.0000
F81G	15	-1982.3484	3994.6968	0.0000	1.0000
SYM	16	-1981.4646	3994.9292	0.0000	1.0000
F81I	15	-1982.8348	3995.6697	0.0000	1.0000
F81IG	16	-1982.3673	3996.7346	0.0000	1.0000
K2P	12	-1987.6358	3999.2716	0.0000	1.0000
F81	14	-2000.8021	4029.6043	0.0000	1.0000
JC69G	12	-2010.4900	4044.9800	0.0000	1.0000
JC69I	12	-2010.7967	4045.5934	0.0000	1.0000
JC69IG	13	-2010.5012	4047.0025	0.0000	1.0000
JC69	11	-2028.6764	4079.3527	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1931.9748	3967.7037	0.6746	0.6746
HKYI	16	-1932.7999	3969.3539	0.2956	0.9702
HKYIG	17	-1931.9613	3974.1614	0.0267	0.9970
GTRG	20	-1924.7472	3979.1871	0.0022	0.9991
GTRI	20	-1925.7479	3981.1886	0.0008	0.9999
GTRIG	21	-1924.7264	3985.6301	0.0001	1.0000
HKY	15	-1953.1391	4003.5478	0.0000	1.0000
GTR	19	-1945.8204	4014.8488	0.0000	1.0000
K2PG	13	-1967.5437	4019.3876	0.0000	1.0000
K2PI	13	-1967.9628	4020.2258	0.0000	1.0000
K2PIG	14	-1967.5496	4025.8841	0.0000	1.0000
SYMG	17	-1961.8920	4034.0229	0.0000	1.0000
SYMI	17	-1962.3221	4034.8831	0.0000	1.0000
SYMIG	18	-1961.8967	4040.5169	0.0000	1.0000
K2P	12	-1987.6358	4053.0873	0.0000	1.0000
F81G	15	-1982.3484	4061.9664	0.0000	1.0000
F81I	15	-1982.8348	4062.9392	0.0000	1.0000
SYM	16	-1981.4646	4066.6834	0.0000	1.0000
F81IG	16	-1982.3673	4068.4888	0.0000	1.0000
F81	14	-2000.8021	4092.3892	0.0000	1.0000
JC69G	12	-2010.4900	4098.7956	0.0000	1.0000
JC69I	12	-2010.7967	4099.4091	0.0000	1.0000
JC69IG	13	-2010.5012	4105.3027	0.0000	1.0000
JC69	11	-2028.6764	4128.6837	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/chr1_32388.phylip = 1 - 655;
 Partition Dummy = 1:591-tomodel-phylip/chr1_32388.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/chr1_32388.phylip = 1 - 655;
 Partition Dummy = 1:591-tomodel-phylip/chr1_32388.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/chr1_32388.phylip = 1 - 655;
 Partition Dummy = 1:591-tomodel-phylip/chr1_32388.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

