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

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

Minimum AIC  model: GTRIG
Minimum AICc model: GTRIG
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
GTRIG	21	-1914.3650	3872.0731	0.4255	0.4255
GTRG	20	-1916.0644	3873.3479	0.2250	0.6505
HKYG	16	-1920.3104	3873.4058	0.2185	0.8690
HKYIG	17	-1920.6712	3876.2268	0.0533	0.9223
GTRI	20	-1917.5386	3876.2963	0.0515	0.9739
HKYI	16	-1922.4353	3877.6556	0.0261	1.0000
K2PG	13	-1932.7846	3892.0923	0.0000	1.0000
K2PIG	14	-1932.1515	3892.9073	0.0000	1.0000
K2PI	13	-1934.6847	3895.8923	0.0000	1.0000
SYMG	17	-1930.6110	3896.1063	0.0000	1.0000
SYMIG	18	-1930.0238	3897.0375	0.0000	1.0000
SYMI	17	-1931.6684	3898.2211	0.0000	1.0000
GTR	19	-1946.2553	3931.6120	0.0000	1.0000
F81G	15	-1952.1367	3934.9650	0.0000	1.0000
F81IG	16	-1951.6550	3936.0951	0.0000	1.0000
F81I	15	-1953.1848	3937.0613	0.0000	1.0000
HKY	15	-1953.4736	3937.6389	0.0000	1.0000
JC69G	12	-1963.1627	3950.7731	0.0000	1.0000
JC69IG	13	-1962.6343	3951.7916	0.0000	1.0000
JC69I	12	-1964.1383	3952.7242	0.0000	1.0000
SYM	16	-1961.4506	3955.6862	0.0000	1.0000
K2P	12	-1966.6533	3957.7543	0.0000	1.0000
F81	14	-1982.2929	3993.1900	0.0000	1.0000
JC69	11	-1993.7759	4009.9300	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	21	-1914.3650	3870.7301	0.4661	0.4661
GTRG	20	-1916.0644	3872.1288	0.2316	0.6978
HKYG	16	-1920.3104	3872.6208	0.1811	0.8789
GTRI	20	-1917.5386	3875.0771	0.0530	0.9319
HKYIG	17	-1920.6712	3875.3424	0.0464	0.9783
HKYI	16	-1922.4353	3876.8706	0.0216	1.0000
K2PG	13	-1932.7846	3891.5693	0.0000	1.0000
K2PIG	14	-1932.1515	3892.3030	0.0000	1.0000
SYMG	17	-1930.6110	3895.2219	0.0000	1.0000
K2PI	13	-1934.6847	3895.3693	0.0000	1.0000
SYMIG	18	-1930.0238	3896.0477	0.0000	1.0000
SYMI	17	-1931.6684	3897.3367	0.0000	1.0000
GTR	19	-1946.2553	3930.5106	0.0000	1.0000
F81G	15	-1952.1367	3934.2734	0.0000	1.0000
F81IG	16	-1951.6550	3935.3101	0.0000	1.0000
F81I	15	-1953.1848	3936.3697	0.0000	1.0000
HKY	15	-1953.4736	3936.9472	0.0000	1.0000
JC69G	12	-1963.1627	3950.3255	0.0000	1.0000
JC69IG	13	-1962.6343	3951.2687	0.0000	1.0000
JC69I	12	-1964.1383	3952.2766	0.0000	1.0000
SYM	16	-1961.4506	3954.9012	0.0000	1.0000
K2P	12	-1966.6533	3957.3067	0.0000	1.0000
F81	14	-1982.2929	3992.5857	0.0000	1.0000
JC69	11	-1993.7759	4009.5518	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1920.3104	3945.6651	0.8103	0.8103
HKYI	16	-1922.4353	3949.9148	0.0968	0.9071
K2PG	13	-1932.7846	3950.9177	0.0586	0.9657
HKYIG	17	-1920.6712	3952.9519	0.0212	0.9869
K2PI	13	-1934.6847	3954.7178	0.0088	0.9957
K2PIG	14	-1932.1515	3956.2167	0.0041	0.9998
GTRG	20	-1916.0644	3963.4341	0.0001	0.9999
GTRI	20	-1917.5386	3966.3824	0.0000	1.0000
GTRIG	21	-1914.3650	3966.6006	0.0000	1.0000
SYMG	17	-1930.6110	3972.8314	0.0000	1.0000
SYMI	17	-1931.6684	3974.9462	0.0000	1.0000
SYMIG	18	-1930.0238	3978.2224	0.0000	1.0000
F81G	15	-1952.1367	4002.7524	0.0000	1.0000
F81I	15	-1953.1848	4004.8486	0.0000	1.0000
JC69G	12	-1963.1627	4005.1087	0.0000	1.0000
HKY	15	-1953.4736	4005.4262	0.0000	1.0000
JC69I	12	-1964.1383	4007.0598	0.0000	1.0000
F81IG	16	-1951.6550	4008.3543	0.0000	1.0000
JC69IG	13	-1962.6343	4010.6171	0.0000	1.0000
K2P	12	-1966.6533	4012.0899	0.0000	1.0000
GTR	19	-1946.2553	4017.2506	0.0000	1.0000
SYM	16	-1961.4506	4027.9455	0.0000	1.0000
F81	14	-1982.2929	4056.4994	0.0000	1.0000
JC69	11	-1993.7759	4059.7697	0.0000	1.0000



To specify default values for best AIC, AICc, or BIC models in MrBayes:

[Mrbayes block for the best AIC model (GTRIG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr8_4241.phylip = 1 - 710;
 Partition Dummy = 1:591-tomodel-phylip/chr8_4241.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


[Mrbayes block for the best AICc model (GTRIG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr8_4241.phylip = 1 - 710;
 Partition Dummy = 1:591-tomodel-phylip/chr8_4241.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


[Mrbayes block for the best BIC model (HKYG)]

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

