
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
11:50 08/01/2011

Input data from file "261-aggressive-phylip/6839_hg18.phylip" (29 taxa, 589 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	64	-1884.2073	3912.2924	0.7397	0.7397
GTRIG	65	-1884.2067	3914.8188	0.2091	0.9488
GTRI	64	-1886.9232	3917.7242	0.0489	0.9977
HKYG	60	-1895.3353	3924.5343	0.0016	0.9993
HKYIG	61	-1895.3348	3927.0226	0.0005	0.9998
HKYI	60	-1897.4626	3928.7889	0.0002	1.0000
GTR	63	-1912.5138	3966.3876	0.0000	1.0000
HKY	59	-1920.6853	3972.7544	0.0000	1.0000
K2PG	57	-1923.5307	3973.5133	0.0000	1.0000
SYMIG	62	-1917.5441	3973.9399	0.0000	1.0000
SYMG	61	-1919.3734	3975.0997	0.0000	1.0000
K2PIG	58	-1923.5302	3975.9736	0.0000	1.0000
K2PI	57	-1925.3382	3977.1284	0.0000	1.0000
SYMI	61	-1920.4594	3977.2717	0.0000	1.0000
F81G	59	-1936.8707	4005.1252	0.0000	1.0000
F81I	59	-1937.9759	4007.3355	0.0000	1.0000
F81IG	60	-1936.8700	4007.6036	0.0000	1.0000
K2P	56	-1947.0842	4018.1683	0.0000	1.0000
SYM	60	-1942.2324	4018.3284	0.0000	1.0000
JC69G	56	-1958.1627	4040.3254	0.0000	1.0000
JC69I	56	-1959.0116	4042.0232	0.0000	1.0000
JC69IG	57	-1958.1620	4042.7760	0.0000	1.0000
F81	58	-1957.4088	4043.7308	0.0000	1.0000
JC69	55	-1978.2912	4078.1397	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1884.2073	3896.4146	0.6967	0.6967
GTRIG	65	-1884.2067	3898.4135	0.2564	0.9531
GTRI	64	-1886.9232	3901.8463	0.0461	0.9992
HKYG	60	-1895.3353	3910.6707	0.0006	0.9997
HKYIG	61	-1895.3348	3912.6697	0.0002	0.9999
HKYI	60	-1897.4626	3914.9253	0.0001	1.0000
GTR	63	-1912.5138	3951.0276	0.0000	1.0000
SYMIG	62	-1917.5441	3959.0882	0.0000	1.0000
HKY	59	-1920.6853	3959.3707	0.0000	1.0000
SYMG	61	-1919.3734	3960.7467	0.0000	1.0000
K2PG	57	-1923.5307	3961.0613	0.0000	1.0000
SYMI	61	-1920.4594	3962.9187	0.0000	1.0000
K2PIG	58	-1923.5302	3963.0604	0.0000	1.0000
K2PI	57	-1925.3382	3964.6764	0.0000	1.0000
F81G	59	-1936.8707	3991.7415	0.0000	1.0000
F81IG	60	-1936.8700	3993.7400	0.0000	1.0000
F81I	59	-1937.9759	3993.9518	0.0000	1.0000
SYM	60	-1942.2324	4004.4648	0.0000	1.0000
K2P	56	-1947.0842	4006.1683	0.0000	1.0000
JC69G	56	-1958.1627	4028.3254	0.0000	1.0000
JC69I	56	-1959.0116	4030.0232	0.0000	1.0000
JC69IG	57	-1958.1620	4030.3240	0.0000	1.0000
F81	58	-1957.4088	4030.8176	0.0000	1.0000
JC69	55	-1978.2912	4066.5825	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1895.3353	4173.3762	0.7259	0.7259
GTRG	64	-1884.2073	4176.6339	0.1424	0.8683
HKYI	60	-1897.4626	4177.6309	0.0865	0.9548
HKYIG	61	-1895.3348	4179.7537	0.0299	0.9847
GTRI	64	-1886.9232	4182.0656	0.0094	0.9941
GTRIG	65	-1884.2067	4183.0112	0.0059	1.0000
K2PG	57	-1923.5307	4210.6316	0.0000	1.0000
K2PI	57	-1925.3382	4214.2467	0.0000	1.0000
K2PIG	58	-1923.5302	4217.0091	0.0000	1.0000
HKY	59	-1920.6853	4217.6978	0.0000	1.0000
GTR	63	-1912.5138	4226.8684	0.0000	1.0000
SYMG	61	-1919.3734	4227.8307	0.0000	1.0000
SYMI	61	-1920.4594	4230.0027	0.0000	1.0000
SYMIG	62	-1917.5441	4230.5506	0.0000	1.0000
F81G	59	-1936.8707	4250.0686	0.0000	1.0000
K2P	56	-1947.0842	4251.3602	0.0000	1.0000
F81I	59	-1937.9759	4252.2789	0.0000	1.0000
F81IG	60	-1936.8700	4256.4455	0.0000	1.0000
SYM	60	-1942.2324	4267.1703	0.0000	1.0000
JC69G	56	-1958.1627	4273.5172	0.0000	1.0000
JC69I	56	-1959.0116	4275.2150	0.0000	1.0000
JC69IG	57	-1958.1620	4279.8943	0.0000	1.0000
F81	58	-1957.4088	4284.7663	0.0000	1.0000
JC69	55	-1978.2912	4307.3959	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 261-aggressive-phylip/6839_hg18.phylip = 1 - 589;
 Partition Dummy = 1:261-aggressive-phylip/6839_hg18.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 261-aggressive-phylip/6839_hg18.phylip = 1 - 589;
 Partition Dummy = 1:261-aggressive-phylip/6839_hg18.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 261-aggressive-phylip/6839_hg18.phylip = 1 - 589;
 Partition Dummy = 1:261-aggressive-phylip/6839_hg18.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

