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

Input data from file "261-aggressive-phylip/7325_hg18.phylip" (29 taxa, 286 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	-1291.3637	2748.3744	0.5961	0.5961
HKYG	60	-1298.4009	2749.3351	0.3687	0.9648
GTRIG	65	-1293.1547	2755.3095	0.0186	0.9834
HKYIG	61	-1299.9520	2755.6718	0.0155	0.9989
HKYI	60	-1304.8811	2762.2956	0.0006	0.9994
GTRI	64	-1298.3817	2762.4104	0.0005	1.0000
SYMG	61	-1306.5803	2768.9284	0.0000	1.0000
SYMIG	62	-1306.2842	2771.5999	0.0000	1.0000
SYMI	61	-1310.5364	2776.8406	0.0000	1.0000
K2PG	57	-1316.9770	2776.9540	0.0000	1.0000
K2PIG	58	-1316.6586	2779.4669	0.0000	1.0000
F81G	59	-1315.1999	2779.7273	0.0000	1.0000
K2PI	57	-1320.9789	2784.9579	0.0000	1.0000
F81IG	60	-1316.7687	2786.0707	0.0000	1.0000
F81I	59	-1321.4573	2792.2420	0.0000	1.0000
JC69G	56	-1333.6965	2807.2708	0.0000	1.0000
JC69IG	57	-1333.4185	2809.8370	0.0000	1.0000
JC69I	56	-1337.6470	2815.1716	0.0000	1.0000
HKY	59	-1338.1248	2825.5771	0.0000	1.0000
GTR	63	-1332.2465	2826.8172	0.0000	1.0000
SYM	60	-1343.9760	2840.4854	0.0000	1.0000
K2P	56	-1353.9608	2847.7993	0.0000	1.0000
F81	58	-1354.9430	2856.0358	0.0000	1.0000
JC69	55	-1370.5113	2877.8052	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1291.3637	2710.7273	0.8976	0.8976
GTRIG	65	-1293.1547	2716.3095	0.0551	0.9527
HKYG	60	-1298.4009	2716.8018	0.0431	0.9958
HKYIG	61	-1299.9520	2721.9039	0.0034	0.9991
GTRI	64	-1298.3817	2724.7634	0.0008	0.9999
HKYI	60	-1304.8811	2729.7622	0.0001	1.0000
SYMG	61	-1306.5803	2735.1605	0.0000	1.0000
SYMIG	62	-1306.2842	2736.5685	0.0000	1.0000
SYMI	61	-1310.5364	2743.0728	0.0000	1.0000
K2PG	57	-1316.9770	2747.9540	0.0000	1.0000
F81G	59	-1315.1999	2748.3999	0.0000	1.0000
K2PIG	58	-1316.6586	2749.3172	0.0000	1.0000
F81IG	60	-1316.7687	2753.5373	0.0000	1.0000
K2PI	57	-1320.9789	2755.9579	0.0000	1.0000
F81I	59	-1321.4573	2760.9146	0.0000	1.0000
JC69G	56	-1333.6965	2779.3931	0.0000	1.0000
JC69IG	57	-1333.4185	2780.8370	0.0000	1.0000
JC69I	56	-1337.6470	2787.2939	0.0000	1.0000
GTR	63	-1332.2465	2790.4929	0.0000	1.0000
HKY	59	-1338.1248	2794.2497	0.0000	1.0000
SYM	60	-1343.9760	2807.9521	0.0000	1.0000
K2P	56	-1353.9608	2819.9216	0.0000	1.0000
F81	58	-1354.9430	2825.8860	0.0000	1.0000
JC69	55	-1370.5113	2851.0226	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1298.4009	2936.1613	0.9726	0.9726
GTRG	64	-1291.3637	2944.7108	0.0135	0.9861
HKYIG	61	-1299.9520	2944.9194	0.0122	0.9983
HKYI	60	-1304.8811	2949.1217	0.0015	0.9998
GTRIG	65	-1293.1547	2953.9489	0.0001	0.9999
K2PG	57	-1316.9770	2956.3456	0.0000	1.0000
SYMG	61	-1306.5803	2958.1760	0.0000	1.0000
GTRI	64	-1298.3817	2958.7469	0.0000	1.0000
K2PIG	58	-1316.6586	2961.3647	0.0000	1.0000
SYMIG	62	-1306.2842	2963.2400	0.0000	1.0000
F81G	59	-1315.1999	2964.1034	0.0000	1.0000
K2PI	57	-1320.9789	2964.3494	0.0000	1.0000
SYMI	61	-1310.5364	2966.0883	0.0000	1.0000
F81IG	60	-1316.7687	2972.8968	0.0000	1.0000
F81I	59	-1321.4573	2976.6181	0.0000	1.0000
JC69G	56	-1333.6965	2984.1286	0.0000	1.0000
JC69IG	57	-1333.4185	2989.2285	0.0000	1.0000
JC69I	56	-1337.6470	2992.0294	0.0000	1.0000
HKY	59	-1338.1248	3009.9532	0.0000	1.0000
GTR	63	-1332.2465	3020.8204	0.0000	1.0000
K2P	56	-1353.9608	3024.6571	0.0000	1.0000
SYM	60	-1343.9760	3027.3116	0.0000	1.0000
F81	58	-1354.9430	3037.9335	0.0000	1.0000
JC69	55	-1370.5113	3052.1022	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/7325_hg18.phylip = 1 - 286;
 Partition Dummy = 1:261-aggressive-phylip/7325_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/7325_hg18.phylip = 1 - 286;
 Partition Dummy = 1:261-aggressive-phylip/7325_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/7325_hg18.phylip = 1 - 286;
 Partition Dummy = 1:261-aggressive-phylip/7325_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

