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

Input data from file "261-aggressive-phylip/11775_hg18.phylip" (29 taxa, 656 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	-2894.9528	5931.9834	0.6649	0.6649
HKYG	60	-2900.5600	5933.4224	0.3238	0.9887
GTRIG	65	-2898.2474	5941.0371	0.0072	0.9958
HKYIG	61	-2903.8226	5942.3791	0.0037	0.9995
GTRI	64	-2902.6867	5947.4511	0.0003	0.9998
HKYI	60	-2908.0051	5948.3127	0.0002	1.0000
SYMG	61	-2934.9688	6004.6717	0.0000	1.0000
K2PG	57	-2940.4410	6005.9390	0.0000	1.0000
SYMIG	62	-2934.8603	6006.8943	0.0000	1.0000
K2PIG	58	-2941.7811	6011.0261	0.0000	1.0000
SYMI	61	-2941.8843	6018.5026	0.0000	1.0000
K2PI	57	-2947.2040	6019.4650	0.0000	1.0000
HKY	59	-2975.3126	6080.5045	0.0000	1.0000
GTR	63	-2972.3527	6084.3271	0.0000	1.0000
F81G	59	-2992.3989	6114.6769	0.0000	1.0000
F81IG	60	-2996.5429	6125.3884	0.0000	1.0000
F81I	59	-2999.1473	6128.1738	0.0000	1.0000
K2P	56	-3020.7412	6164.1402	0.0000	1.0000
SYM	60	-3016.2846	6164.8718	0.0000	1.0000
JC69G	56	-3028.0060	6178.6698	0.0000	1.0000
JC69IG	57	-3027.7276	6180.5121	0.0000	1.0000
JC69I	56	-3033.7229	6190.1036	0.0000	1.0000
F81	58	-3063.6265	6254.7170	0.0000	1.0000
JC69	55	-3102.0965	6324.4597	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-2894.9528	5917.9055	0.8214	0.8214
HKYG	60	-2900.5600	5921.1199	0.1646	0.9860
GTRIG	65	-2898.2474	5926.4948	0.0112	0.9972
HKYIG	61	-2903.8226	5929.6451	0.0023	0.9995
GTRI	64	-2902.6867	5933.3733	0.0004	0.9999
HKYI	60	-2908.0051	5936.0102	0.0001	1.0000
SYMG	61	-2934.9688	5991.9377	0.0000	1.0000
SYMIG	62	-2934.8603	5993.7206	0.0000	1.0000
K2PG	57	-2940.4410	5994.8821	0.0000	1.0000
K2PIG	58	-2941.7811	5999.5622	0.0000	1.0000
SYMI	61	-2941.8843	6005.7686	0.0000	1.0000
K2PI	57	-2947.2040	6008.4081	0.0000	1.0000
HKY	59	-2975.3126	6068.6253	0.0000	1.0000
GTR	63	-2972.3527	6070.7055	0.0000	1.0000
F81G	59	-2992.3989	6102.7977	0.0000	1.0000
F81IG	60	-2996.5429	6113.0859	0.0000	1.0000
F81I	59	-2999.1473	6116.2946	0.0000	1.0000
SYM	60	-3016.2846	6152.5692	0.0000	1.0000
K2P	56	-3020.7412	6153.4825	0.0000	1.0000
JC69G	56	-3028.0060	6168.0120	0.0000	1.0000
JC69IG	57	-3027.7276	6169.4552	0.0000	1.0000
JC69I	56	-3033.7229	6179.4458	0.0000	1.0000
F81	58	-3063.6265	6243.2531	0.0000	1.0000
JC69	55	-3102.0965	6314.1931	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-2900.5600	6190.2896	0.9973	0.9973
HKYIG	61	-2903.8226	6203.3009	0.0015	0.9988
GTRG	64	-2894.9528	6205.0198	0.0006	0.9994
HKYI	60	-2908.0051	6205.1799	0.0006	1.0000
GTRIG	65	-2898.2474	6218.0952	0.0000	1.0000
GTRI	64	-2902.6867	6220.4876	0.0000	1.0000
K2PG	57	-2940.4410	6250.5933	0.0000	1.0000
K2PIG	58	-2941.7811	6259.7595	0.0000	1.0000
K2PI	57	-2947.2040	6264.1193	0.0000	1.0000
SYMG	61	-2934.9688	6265.5935	0.0000	1.0000
SYMIG	62	-2934.8603	6271.8626	0.0000	1.0000
SYMI	61	-2941.8843	6279.4244	0.0000	1.0000
HKY	59	-2975.3126	6333.3088	0.0000	1.0000
GTR	63	-2972.3527	6353.3336	0.0000	1.0000
F81G	59	-2992.3989	6367.4812	0.0000	1.0000
F81I	59	-2999.1473	6380.9780	0.0000	1.0000
F81IG	60	-2996.5429	6382.2555	0.0000	1.0000
K2P	56	-3020.7412	6404.7075	0.0000	1.0000
JC69G	56	-3028.0060	6419.2370	0.0000	1.0000
SYM	60	-3016.2846	6421.7389	0.0000	1.0000
JC69IG	57	-3027.7276	6425.1664	0.0000	1.0000
JC69I	56	-3033.7229	6430.6708	0.0000	1.0000
F81	58	-3063.6265	6503.4504	0.0000	1.0000
JC69	55	-3102.0965	6560.9319	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/11775_hg18.phylip = 1 - 656;
 Partition Dummy = 1:261-aggressive-phylip/11775_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/11775_hg18.phylip = 1 - 656;
 Partition Dummy = 1:261-aggressive-phylip/11775_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/11775_hg18.phylip = 1 - 656;
 Partition Dummy = 1:261-aggressive-phylip/11775_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

