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

Input data from file "261-aggressive-phylip/11075_hg18.phylip" (29 taxa, 659 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	-2444.7741	5031.5550	0.9874	0.9874
GTRIG	65	-2447.9805	5040.4298	0.0117	0.9991
HKYG	60	-2456.6709	5045.5827	0.0009	1.0000
HKYIG	61	-2460.3983	5055.4666	0.0000	1.0000
SYMG	61	-2471.7430	5078.1561	0.0000	1.0000
SYMIG	62	-2471.7514	5080.6102	0.0000	1.0000
K2PG	57	-2481.2194	5087.4404	0.0000	1.0000
K2PIG	58	-2484.8330	5097.0727	0.0000	1.0000
F81G	59	-2530.1285	5190.0767	0.0000	1.0000
F81IG	60	-2531.1559	5194.5526	0.0000	1.0000
GTRI	64	-2526.9601	5195.9270	0.0000	1.0000
HKYI	60	-2540.7249	5213.6905	0.0000	1.0000
JC69G	56	-2550.5681	5223.7408	0.0000	1.0000
GTR	63	-2543.4171	5226.3871	0.0000	1.0000
JC69IG	57	-2551.5878	5228.1774	0.0000	1.0000
HKY	59	-2556.2254	5242.2706	0.0000	1.0000
SYMI	61	-2554.1269	5242.9237	0.0000	1.0000
K2PI	57	-2566.1536	5257.3088	0.0000	1.0000
SYM	60	-2569.6535	5271.5479	0.0000	1.0000
K2P	56	-2582.5090	5287.6227	0.0000	1.0000
F81I	59	-2611.7503	5353.3203	0.0000	1.0000
JC69I	56	-2631.6226	5385.8499	0.0000	1.0000
F81	58	-2629.4979	5386.4024	0.0000	1.0000
JC69	55	-2649.5903	5419.3963	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-2444.7741	5017.5483	0.9850	0.9850
GTRIG	65	-2447.9805	5025.9610	0.0147	0.9996
HKYG	60	-2456.6709	5033.3419	0.0004	1.0000
HKYIG	61	-2460.3983	5042.7966	0.0000	1.0000
SYMG	61	-2471.7430	5065.4861	0.0000	1.0000
SYMIG	62	-2471.7514	5067.5028	0.0000	1.0000
K2PG	57	-2481.2194	5076.4387	0.0000	1.0000
K2PIG	58	-2484.8330	5085.6660	0.0000	1.0000
F81G	59	-2530.1285	5178.2570	0.0000	1.0000
GTRI	64	-2526.9601	5181.9202	0.0000	1.0000
F81IG	60	-2531.1559	5182.3118	0.0000	1.0000
HKYI	60	-2540.7249	5201.4497	0.0000	1.0000
GTR	63	-2543.4171	5212.8342	0.0000	1.0000
JC69G	56	-2550.5681	5213.1362	0.0000	1.0000
JC69IG	57	-2551.5878	5217.1757	0.0000	1.0000
SYMI	61	-2554.1269	5230.2537	0.0000	1.0000
HKY	59	-2556.2254	5230.4509	0.0000	1.0000
K2PI	57	-2566.1536	5246.3072	0.0000	1.0000
SYM	60	-2569.6535	5259.3071	0.0000	1.0000
K2P	56	-2582.5090	5277.0180	0.0000	1.0000
F81I	59	-2611.7503	5341.5006	0.0000	1.0000
F81	58	-2629.4979	5374.9957	0.0000	1.0000
JC69I	56	-2631.6226	5375.2453	0.0000	1.0000
JC69	55	-2649.5903	5409.1807	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-2456.6709	5302.7853	0.7465	0.7465
GTRG	64	-2444.7741	5304.9546	0.2524	0.9989
HKYIG	61	-2460.3983	5316.7307	0.0007	0.9996
GTRIG	65	-2447.9805	5317.8580	0.0004	1.0000
K2PG	57	-2481.2194	5332.4100	0.0000	1.0000
SYMG	61	-2471.7430	5339.4202	0.0000	1.0000
SYMIG	62	-2471.7514	5345.9276	0.0000	1.0000
K2PIG	58	-2484.8330	5346.1280	0.0000	1.0000
F81G	59	-2530.1285	5443.2097	0.0000	1.0000
F81IG	60	-2531.1559	5451.7552	0.0000	1.0000
JC69G	56	-2550.5681	5464.6167	0.0000	1.0000
GTRI	64	-2526.9601	5469.3265	0.0000	1.0000
HKYI	60	-2540.7249	5470.8931	0.0000	1.0000
JC69IG	57	-2551.5878	5473.1469	0.0000	1.0000
HKY	59	-2556.2254	5495.4036	0.0000	1.0000
GTR	63	-2543.4171	5495.7497	0.0000	1.0000
K2PI	57	-2566.1536	5502.2784	0.0000	1.0000
SYMI	61	-2554.1269	5504.1879	0.0000	1.0000
K2P	56	-2582.5090	5528.4985	0.0000	1.0000
SYM	60	-2569.6535	5528.7505	0.0000	1.0000
F81I	59	-2611.7503	5606.4533	0.0000	1.0000
JC69I	56	-2631.6226	5626.7258	0.0000	1.0000
F81	58	-2629.4979	5635.4577	0.0000	1.0000
JC69	55	-2649.5903	5656.1705	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/11075_hg18.phylip = 1 - 659;
 Partition Dummy = 1:261-aggressive-phylip/11075_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/11075_hg18.phylip = 1 - 659;
 Partition Dummy = 1:261-aggressive-phylip/11075_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/11075_hg18.phylip = 1 - 659;
 Partition Dummy = 1:261-aggressive-phylip/11075_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

