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

Input data from file "591-tomodel-phylip/chr20_5160.phylip" (7 taxa, 638 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	-1890.2369	3821.8352	0.5875	0.5875
GTRIG	21	-1889.9076	3823.3152	0.2803	0.8678
GTRI	20	-1891.7971	3824.9556	0.1234	0.9913
HKYG	16	-1899.2624	3831.4007	0.0049	0.9962
HKYIG	17	-1898.7709	3832.5289	0.0028	0.9990
HKYI	16	-1900.8534	3834.5829	0.0010	1.0000
SYMG	17	-1908.8205	3852.6282	0.0000	1.0000
SYMIG	18	-1908.3528	3853.8106	0.0000	1.0000
K2PG	13	-1914.4082	3855.3998	0.0000	1.0000
SYMI	17	-1910.3428	3855.6727	0.0000	1.0000
K2PIG	14	-1913.8896	3856.4534	0.0000	1.0000
K2PI	13	-1915.9719	3858.5272	0.0000	1.0000
GTR	19	-1925.1461	3889.5219	0.0000	1.0000
HKY	15	-1934.9217	3900.6151	0.0000	1.0000
SYM	16	-1942.9068	3918.6896	0.0000	1.0000
K2P	12	-1949.4536	3923.4065	0.0000	1.0000
F81G	15	-1950.3235	3931.4186	0.0000	1.0000
F81IG	16	-1950.1057	3933.0873	0.0000	1.0000
F81I	15	-1951.3611	3933.4940	0.0000	1.0000
JC69G	12	-1962.5621	3949.6234	0.0000	1.0000
JC69IG	13	-1962.3617	3951.3067	0.0000	1.0000
JC69I	12	-1963.5752	3951.6496	0.0000	1.0000
F81	14	-1978.9104	3986.4950	0.0000	1.0000
JC69	11	-1991.4009	4005.2234	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-1890.2369	3820.4738	0.5769	0.5769
GTRIG	21	-1889.9076	3821.8152	0.2950	0.8719
GTRI	20	-1891.7971	3823.5942	0.1212	0.9932
HKYG	16	-1899.2624	3830.5247	0.0038	0.9969
HKYIG	17	-1898.7709	3831.5418	0.0023	0.9992
HKYI	16	-1900.8534	3833.7069	0.0008	1.0000
SYMG	17	-1908.8205	3851.6411	0.0000	1.0000
SYMIG	18	-1908.3528	3852.7056	0.0000	1.0000
SYMI	17	-1910.3428	3854.6856	0.0000	1.0000
K2PG	13	-1914.4082	3854.8165	0.0000	1.0000
K2PIG	14	-1913.8896	3855.7792	0.0000	1.0000
K2PI	13	-1915.9719	3857.9438	0.0000	1.0000
GTR	19	-1925.1461	3888.2921	0.0000	1.0000
HKY	15	-1934.9217	3899.8434	0.0000	1.0000
SYM	16	-1942.9068	3917.8136	0.0000	1.0000
K2P	12	-1949.4536	3922.9073	0.0000	1.0000
F81G	15	-1950.3235	3930.6469	0.0000	1.0000
F81IG	16	-1950.1057	3932.2113	0.0000	1.0000
F81I	15	-1951.3611	3932.7223	0.0000	1.0000
JC69G	12	-1962.5621	3949.1242	0.0000	1.0000
JC69IG	13	-1962.3617	3950.7234	0.0000	1.0000
JC69I	12	-1963.5752	3951.1504	0.0000	1.0000
F81	14	-1978.9104	3985.8208	0.0000	1.0000
JC69	11	-1991.4009	4004.8017	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1899.2624	3901.8581	0.7694	0.7694
HKYI	16	-1900.8534	3905.0403	0.1567	0.9261
HKYIG	17	-1898.7709	3907.3335	0.0498	0.9759
GTRG	20	-1890.2369	3909.6405	0.0157	0.9916
GTRI	20	-1891.7971	3912.7609	0.0033	0.9950
K2PG	13	-1914.4082	3912.7749	0.0033	0.9982
GTRIG	21	-1889.9076	3915.4403	0.0009	0.9991
K2PI	13	-1915.9719	3915.9022	0.0007	0.9998
K2PIG	14	-1913.8896	3918.1960	0.0002	1.0000
SYMG	17	-1908.8205	3927.4329	0.0000	1.0000
SYMI	17	-1910.3428	3930.4774	0.0000	1.0000
SYMIG	18	-1908.3528	3932.9557	0.0000	1.0000
HKY	15	-1934.9217	3966.7185	0.0000	1.0000
GTR	19	-1925.1461	3973.0005	0.0000	1.0000
K2P	12	-1949.4536	3976.4073	0.0000	1.0000
SYM	16	-1942.9068	3989.1470	0.0000	1.0000
F81G	15	-1950.3235	3997.5220	0.0000	1.0000
F81I	15	-1951.3611	3999.5973	0.0000	1.0000
JC69G	12	-1962.5621	4002.6242	0.0000	1.0000
F81IG	16	-1950.1057	4003.5447	0.0000	1.0000
JC69I	12	-1963.5752	4004.6504	0.0000	1.0000
JC69IG	13	-1962.3617	4008.6818	0.0000	1.0000
F81	14	-1978.9104	4048.2376	0.0000	1.0000
JC69	11	-1991.4009	4053.8434	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/chr20_5160.phylip = 1 - 638;
 Partition Dummy = 1:591-tomodel-phylip/chr20_5160.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/chr20_5160.phylip = 1 - 638;
 Partition Dummy = 1:591-tomodel-phylip/chr20_5160.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/chr20_5160.phylip = 1 - 638;
 Partition Dummy = 1:591-tomodel-phylip/chr20_5160.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

