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

Input data from file "261-aggressive-phylip/9935_hg18.phylip" (29 taxa, 686 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	-2360.8552	4863.1081	0.9275	0.9275
GTRIG	65	-2362.2380	4868.3147	0.0687	0.9961
HKYG	60	-2371.2890	4874.2900	0.0035	0.9996
HKYIG	61	-2372.6360	4879.3938	0.0003	0.9999
GTRI	64	-2370.0057	4881.4092	0.0001	1.0000
HKYI	60	-2376.2613	4884.2346	0.0000	1.0000
SYMG	61	-2379.0845	4892.2907	0.0000	1.0000
K2PG	57	-2387.1661	4898.8609	0.0000	1.0000
SYMIG	62	-2384.2845	4905.1083	0.0000	1.0000
K2PIG	58	-2390.5083	4907.9320	0.0000	1.0000
SYMI	61	-2387.9172	4909.9563	0.0000	1.0000
K2PI	57	-2393.6552	4911.8392	0.0000	1.0000
GTR	63	-2397.0492	4933.0630	0.0000	1.0000
F81G	59	-2405.2220	4939.7539	0.0000	1.0000
F81IG	60	-2405.0342	4941.7804	0.0000	1.0000
SYM	60	-2405.6432	4942.9984	0.0000	1.0000
HKY	59	-2407.5041	4944.3182	0.0000	1.0000
F81I	59	-2408.2277	4945.7654	0.0000	1.0000
JC69G	56	-2414.8018	4951.7531	0.0000	1.0000
K2P	56	-2417.4028	4956.9550	0.0000	1.0000
JC69IG	57	-2421.1163	4966.7612	0.0000	1.0000
JC69I	56	-2424.1057	4970.3608	0.0000	1.0000
F81	58	-2425.3134	4977.5423	0.0000	1.0000
JC69	55	-2443.0662	5005.9102	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-2360.8552	4849.7104	0.9139	0.9139
GTRIG	65	-2362.2380	4854.4760	0.0843	0.9983
HKYG	60	-2371.2890	4862.5780	0.0015	0.9998
HKYIG	61	-2372.6360	4867.2720	0.0001	0.9999
GTRI	64	-2370.0057	4868.0115	0.0001	1.0000
HKYI	60	-2376.2613	4872.5226	0.0000	1.0000
SYMG	61	-2379.0845	4880.1689	0.0000	1.0000
K2PG	57	-2387.1661	4888.3322	0.0000	1.0000
SYMIG	62	-2384.2845	4892.5690	0.0000	1.0000
K2PIG	58	-2390.5083	4897.0165	0.0000	1.0000
SYMI	61	-2387.9172	4897.8345	0.0000	1.0000
K2PI	57	-2393.6552	4901.3105	0.0000	1.0000
GTR	63	-2397.0492	4920.0984	0.0000	1.0000
F81G	59	-2405.2220	4928.4440	0.0000	1.0000
F81IG	60	-2405.0342	4930.0684	0.0000	1.0000
SYM	60	-2405.6432	4931.2864	0.0000	1.0000
HKY	59	-2407.5041	4933.0083	0.0000	1.0000
F81I	59	-2408.2277	4934.4555	0.0000	1.0000
JC69G	56	-2414.8018	4941.6037	0.0000	1.0000
K2P	56	-2417.4028	4946.8056	0.0000	1.0000
JC69IG	57	-2421.1163	4956.2325	0.0000	1.0000
JC69I	56	-2424.1057	4960.2113	0.0000	1.0000
F81	58	-2425.3134	4966.6268	0.0000	1.0000
JC69	55	-2443.0662	4996.1325	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-2371.2890	5134.4307	0.9157	0.9157
GTRG	64	-2360.8552	5139.6865	0.0661	0.9818
HKYIG	61	-2372.6360	5143.6556	0.0091	0.9909
HKYI	60	-2376.2613	5144.3753	0.0063	0.9972
K2PG	57	-2387.1661	5146.5922	0.0021	0.9993
GTRIG	65	-2362.2380	5148.9830	0.0006	1.0000
SYMG	61	-2379.0845	5156.5525	0.0000	1.0000
GTRI	64	-2370.0057	5157.9877	0.0000	1.0000
K2PI	57	-2393.6552	5159.5705	0.0000	1.0000
K2PIG	58	-2390.5083	5159.8074	0.0000	1.0000
SYMIG	62	-2384.2845	5173.4834	0.0000	1.0000
SYMI	61	-2387.9172	5174.2180	0.0000	1.0000
JC69G	56	-2414.8018	5195.3328	0.0000	1.0000
F81G	59	-2405.2220	5195.7658	0.0000	1.0000
HKY	59	-2407.5041	5200.3301	0.0000	1.0000
K2P	56	-2417.4028	5200.5347	0.0000	1.0000
F81I	59	-2408.2277	5201.7773	0.0000	1.0000
F81IG	60	-2405.0342	5201.9211	0.0000	1.0000
SYM	60	-2405.6432	5203.1391	0.0000	1.0000
GTR	63	-2397.0492	5205.5437	0.0000	1.0000
JC69I	56	-2424.1057	5213.9405	0.0000	1.0000
JC69IG	57	-2421.1163	5214.4925	0.0000	1.0000
F81	58	-2425.3134	5229.4177	0.0000	1.0000
JC69	55	-2443.0662	5245.3307	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/9935_hg18.phylip = 1 - 686;
 Partition Dummy = 1:261-aggressive-phylip/9935_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/9935_hg18.phylip = 1 - 686;
 Partition Dummy = 1:261-aggressive-phylip/9935_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/9935_hg18.phylip = 1 - 686;
 Partition Dummy = 1:261-aggressive-phylip/9935_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

