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

Input data from file "591-tomodel-phylip/chr3_14142.phylip" (7 taxa, 624 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	-1910.6608	3862.7147	0.5385	0.5385
GTRI	20	-1911.4347	3864.2625	0.2484	0.7869
GTRIG	21	-1910.5324	3864.5997	0.2098	0.9967
HKYG	16	-1920.5647	3874.0256	0.0019	0.9986
HKYI	16	-1921.2544	3875.4050	0.0009	0.9996
HKYIG	17	-1920.9945	3876.9990	0.0004	1.0000
SYMG	17	-1925.5201	3886.0500	0.0000	1.0000
SYMI	17	-1926.2542	3887.5184	0.0000	1.0000
SYMIG	18	-1926.0764	3889.2833	0.0000	1.0000
K2PG	13	-1934.3857	3895.3681	0.0000	1.0000
K2PI	13	-1934.9913	3896.5793	0.0000	1.0000
K2PIG	14	-1934.7555	3898.2006	0.0000	1.0000
GTR	19	-1941.9670	3923.1924	0.0000	1.0000
HKY	15	-1951.7520	3934.2935	0.0000	1.0000
SYM	16	-1952.0802	3937.0567	0.0000	1.0000
F81G	15	-1959.3529	3949.4954	0.0000	1.0000
K2P	12	-1962.6140	3949.7387	0.0000	1.0000
F81I	15	-1959.8033	3950.3961	0.0000	1.0000
F81IG	16	-1959.2832	3951.4625	0.0000	1.0000
JC69G	12	-1971.3206	3967.1519	0.0000	1.0000
JC69I	12	-1971.5675	3967.6457	0.0000	1.0000
JC69IG	13	-1973.3477	3973.2921	0.0000	1.0000
F81	14	-1985.6448	3999.9792	0.0000	1.0000
JC69	11	-1995.9427	4014.3167	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-1910.6608	3861.3216	0.5307	0.5307
GTRI	20	-1911.4347	3862.8695	0.2448	0.7755
GTRIG	21	-1910.5324	3863.0648	0.2220	0.9975
HKYG	16	-1920.5647	3873.1293	0.0014	0.9989
HKYI	16	-1921.2544	3874.5088	0.0007	0.9996
HKYIG	17	-1920.9945	3875.9891	0.0003	1.0000
SYMG	17	-1925.5201	3885.0401	0.0000	1.0000
SYMI	17	-1926.2542	3886.5085	0.0000	1.0000
SYMIG	18	-1926.0764	3888.1527	0.0000	1.0000
K2PG	13	-1934.3857	3894.7713	0.0000	1.0000
K2PI	13	-1934.9913	3895.9826	0.0000	1.0000
K2PIG	14	-1934.7555	3897.5109	0.0000	1.0000
GTR	19	-1941.9670	3921.9341	0.0000	1.0000
HKY	15	-1951.7520	3933.5040	0.0000	1.0000
SYM	16	-1952.0802	3936.1605	0.0000	1.0000
F81G	15	-1959.3529	3948.7059	0.0000	1.0000
K2P	12	-1962.6140	3949.2281	0.0000	1.0000
F81I	15	-1959.8033	3949.6066	0.0000	1.0000
F81IG	16	-1959.2832	3950.5663	0.0000	1.0000
JC69G	12	-1971.3206	3966.6413	0.0000	1.0000
JC69I	12	-1971.5675	3967.1351	0.0000	1.0000
JC69IG	13	-1973.3477	3972.6954	0.0000	1.0000
F81	14	-1985.6448	3999.2896	0.0000	1.0000
JC69	11	-1995.9427	4013.8853	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1920.5647	3944.1077	0.6135	0.6135
HKYI	16	-1921.2544	3945.4872	0.3078	0.9213
GTRG	20	-1910.6608	3950.0446	0.0315	0.9528
HKYIG	17	-1920.9945	3951.4036	0.0160	0.9688
GTRI	20	-1911.4347	3951.5925	0.0145	0.9833
K2PG	13	-1934.3857	3952.4413	0.0095	0.9929
K2PI	13	-1934.9913	3953.6525	0.0052	0.9980
GTRIG	21	-1910.5324	3956.2240	0.0014	0.9995
K2PIG	14	-1934.7555	3959.6170	0.0003	0.9997
SYMG	17	-1925.5201	3960.4547	0.0002	0.9999
SYMI	17	-1926.2542	3961.9230	0.0001	1.0000
SYMIG	18	-1926.0764	3968.0034	0.0000	1.0000
HKY	15	-1951.7520	4000.0463	0.0000	1.0000
K2P	12	-1962.6140	4002.4619	0.0000	1.0000
GTR	19	-1941.9670	4006.2210	0.0000	1.0000
SYM	16	-1952.0802	4007.1389	0.0000	1.0000
F81G	15	-1959.3529	4015.2481	0.0000	1.0000
F81I	15	-1959.8033	4016.1489	0.0000	1.0000
JC69G	12	-1971.3206	4019.8751	0.0000	1.0000
JC69I	12	-1971.5675	4020.3689	0.0000	1.0000
F81IG	16	-1959.2832	4021.5447	0.0000	1.0000
JC69IG	13	-1973.3477	4030.3654	0.0000	1.0000
F81	14	-1985.6448	4061.3957	0.0000	1.0000
JC69	11	-1995.9427	4062.6830	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_14142.phylip = 1 - 624;
 Partition Dummy = 1:591-tomodel-phylip/chr3_14142.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_14142.phylip = 1 - 624;
 Partition Dummy = 1:591-tomodel-phylip/chr3_14142.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/chr3_14142.phylip = 1 - 624;
 Partition Dummy = 1:591-tomodel-phylip/chr3_14142.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

