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

Input data from file "261-aggressive-phylip/5754_hg18.phylip" (29 taxa, 1263 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	-4247.3426	8629.6302	0.6223	0.6223
HKYG	60	-4252.3073	8630.7044	0.3637	0.9860
GTRIG	65	-4250.3895	8637.9469	0.0097	0.9957
HKYIG	61	-4255.6421	8639.5822	0.0043	1.0000
SYMG	61	-4266.6843	8661.6666	0.0000	1.0000
K2PG	57	-4273.0166	8665.5204	0.0000	1.0000
SYMIG	62	-4269.8509	8670.2119	0.0000	1.0000
K2PIG	58	-4276.8593	8675.4030	0.0000	1.0000
GTR	63	-4351.0574	8834.8403	0.0000	1.0000
HKY	59	-4355.8471	8835.5796	0.0000	1.0000
GTRI	64	-4351.0604	8837.0658	0.0000	1.0000
HKYI	60	-4355.8502	8837.7902	0.0000	1.0000
SYM	60	-4367.1639	8860.4177	0.0000	1.0000
SYMI	61	-4367.1673	8862.6327	0.0000	1.0000
K2P	56	-4374.3868	8866.0671	0.0000	1.0000
K2PI	57	-4374.3901	8868.2674	0.0000	1.0000
F81G	59	-4376.3044	8876.4940	0.0000	1.0000
F81IG	60	-4377.4744	8881.0386	0.0000	1.0000
JC69G	56	-4391.7453	8900.7841	0.0000	1.0000
JC69IG	57	-4396.3934	8912.2738	0.0000	1.0000
F81I	59	-4472.5281	9068.9415	0.0000	1.0000
F81	58	-4475.8805	9073.4454	0.0000	1.0000
JC69	55	-4490.1853	9095.4742	0.0000	1.0000
JC69I	56	-4490.1890	9097.6716	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-4247.3426	8622.6853	0.7125	0.7125
HKYG	60	-4252.3073	8624.6145	0.2715	0.9840
GTRIG	65	-4250.3895	8630.7790	0.0125	0.9964
HKYIG	61	-4255.6421	8633.2842	0.0036	1.0000
SYMG	61	-4266.6843	8655.3685	0.0000	1.0000
K2PG	57	-4273.0166	8660.0332	0.0000	1.0000
SYMIG	62	-4269.8509	8663.7019	0.0000	1.0000
K2PIG	58	-4276.8593	8669.7186	0.0000	1.0000
GTR	63	-4351.0574	8828.1147	0.0000	1.0000
HKY	59	-4355.8471	8829.6943	0.0000	1.0000
GTRI	64	-4351.0604	8830.1209	0.0000	1.0000
HKYI	60	-4355.8502	8831.7004	0.0000	1.0000
SYM	60	-4367.1639	8854.3279	0.0000	1.0000
SYMI	61	-4367.1673	8856.3346	0.0000	1.0000
K2P	56	-4374.3868	8860.7735	0.0000	1.0000
K2PI	57	-4374.3901	8862.7802	0.0000	1.0000
F81G	59	-4376.3044	8870.6087	0.0000	1.0000
F81IG	60	-4377.4744	8874.9487	0.0000	1.0000
JC69G	56	-4391.7453	8895.4906	0.0000	1.0000
JC69IG	57	-4396.3934	8906.7867	0.0000	1.0000
F81I	59	-4472.5281	9063.0562	0.0000	1.0000
F81	58	-4475.8805	9067.7610	0.0000	1.0000
JC69	55	-4490.1853	9090.3706	0.0000	1.0000
JC69I	56	-4490.1890	9092.3780	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-4252.3073	8933.0892	0.9989	0.9989
HKYIG	61	-4255.6421	8946.9001	0.0010	0.9999
GTRG	64	-4247.3426	8951.7250	0.0001	1.0000
K2PG	57	-4273.0166	8953.0842	0.0000	1.0000
GTRIG	65	-4250.3895	8964.9600	0.0000	1.0000
K2PIG	58	-4276.8593	8967.9108	0.0000	1.0000
SYMG	61	-4266.6843	8968.9845	0.0000	1.0000
SYMIG	62	-4269.8509	8982.4591	0.0000	1.0000
HKY	59	-4355.8471	9133.0278	0.0000	1.0000
HKYI	60	-4355.8502	9140.1751	0.0000	1.0000
K2P	56	-4374.3868	9148.6833	0.0000	1.0000
GTR	63	-4351.0574	9152.0132	0.0000	1.0000
K2PI	57	-4374.3901	9155.8312	0.0000	1.0000
GTRI	64	-4351.0604	9159.1605	0.0000	1.0000
SYM	60	-4367.1639	9162.8026	0.0000	1.0000
SYMI	61	-4367.1673	9169.9505	0.0000	1.0000
F81G	59	-4376.3044	9173.9422	0.0000	1.0000
JC69G	56	-4391.7453	9183.4003	0.0000	1.0000
F81IG	60	-4377.4744	9183.4234	0.0000	1.0000
JC69IG	57	-4396.3934	9199.8377	0.0000	1.0000
F81	58	-4475.8805	9365.9533	0.0000	1.0000
F81I	59	-4472.5281	9366.3896	0.0000	1.0000
JC69	55	-4490.1853	9373.1391	0.0000	1.0000
JC69I	56	-4490.1890	9380.2878	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/5754_hg18.phylip = 1 - 1263;
 Partition Dummy = 1:261-aggressive-phylip/5754_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/5754_hg18.phylip = 1 - 1263;
 Partition Dummy = 1:261-aggressive-phylip/5754_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/5754_hg18.phylip = 1 - 1263;
 Partition Dummy = 1:261-aggressive-phylip/5754_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

