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

Input data from file "261-aggressive-phylip/11141_hg18.phylip" (29 taxa, 605 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	-2607.3648	5358.1370	0.5837	0.5837
GTRIG	65	-2606.7734	5359.4651	0.3005	0.8842
HKYG	60	-2614.4153	5362.2864	0.0733	0.9575
HKYIG	61	-2613.9259	5363.7818	0.0347	0.9922
GTRI	64	-2611.9669	5367.3412	0.0059	0.9980
HKYI	60	-2618.0366	5369.5291	0.0020	1.0000
K2PG	57	-2634.0196	5394.1269	0.0000	1.0000
K2PIG	58	-2633.6983	5395.9315	0.0000	1.0000
SYMG	61	-2632.1377	5400.2054	0.0000	1.0000
K2PI	57	-2637.3926	5400.8730	0.0000	1.0000
SYMIG	62	-2631.7780	5401.9694	0.0000	1.0000
SYMI	61	-2635.6539	5407.2377	0.0000	1.0000
GTR	63	-2658.1725	5457.2507	0.0000	1.0000
HKY	59	-2665.8727	5462.7362	0.0000	1.0000
K2P	56	-2685.8130	5495.2757	0.0000	1.0000
SYM	60	-2682.4017	5498.2593	0.0000	1.0000
F81G	59	-2699.6136	5530.2181	0.0000	1.0000
F81IG	60	-2699.2969	5532.0498	0.0000	1.0000
F81I	59	-2702.8056	5536.6019	0.0000	1.0000
JC69G	56	-2712.4151	5548.4798	0.0000	1.0000
JC69IG	57	-2712.1584	5550.4045	0.0000	1.0000
JC69I	56	-2715.3760	5554.4017	0.0000	1.0000
F81	58	-2748.4761	5625.4869	0.0000	1.0000
JC69	55	-2761.5253	5644.2710	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-2607.3648	5342.7296	0.5709	0.5709
GTRIG	65	-2606.7734	5343.5468	0.3794	0.9503
HKYG	60	-2614.4153	5348.8305	0.0270	0.9773
HKYIG	61	-2613.9259	5349.8518	0.0162	0.9936
GTRI	64	-2611.9669	5351.9338	0.0057	0.9993
HKYI	60	-2618.0366	5356.0732	0.0007	1.0000
K2PG	57	-2634.0196	5382.0392	0.0000	1.0000
K2PIG	58	-2633.6983	5383.3967	0.0000	1.0000
SYMG	61	-2632.1377	5386.2754	0.0000	1.0000
SYMIG	62	-2631.7780	5387.5561	0.0000	1.0000
K2PI	57	-2637.3926	5388.7852	0.0000	1.0000
SYMI	61	-2635.6539	5393.3077	0.0000	1.0000
GTR	63	-2658.1725	5442.3450	0.0000	1.0000
HKY	59	-2665.8727	5449.7453	0.0000	1.0000
K2P	56	-2685.8130	5483.6261	0.0000	1.0000
SYM	60	-2682.4017	5484.8034	0.0000	1.0000
F81G	59	-2699.6136	5517.2273	0.0000	1.0000
F81IG	60	-2699.2969	5518.5939	0.0000	1.0000
F81I	59	-2702.8056	5523.6111	0.0000	1.0000
JC69G	56	-2712.4151	5536.8301	0.0000	1.0000
JC69IG	57	-2712.1584	5538.3168	0.0000	1.0000
JC69I	56	-2715.3760	5542.7521	0.0000	1.0000
F81	58	-2748.4761	5612.9521	0.0000	1.0000
JC69	55	-2761.5253	5633.0506	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-2614.4153	5613.1442	0.9120	0.9120
HKYIG	61	-2613.9259	5618.5707	0.0605	0.9724
HKYI	60	-2618.0366	5620.3869	0.0244	0.9968
GTRG	64	-2607.3648	5624.6642	0.0029	0.9997
GTRIG	65	-2606.7734	5629.8866	0.0002	0.9999
K2PG	57	-2634.0196	5633.1372	0.0000	1.0000
GTRI	64	-2611.9669	5633.8684	0.0000	1.0000
K2PIG	58	-2633.6983	5638.8999	0.0000	1.0000
K2PI	57	-2637.3926	5639.8832	0.0000	1.0000
SYMG	61	-2632.1377	5654.9943	0.0000	1.0000
SYMIG	62	-2631.7780	5660.6802	0.0000	1.0000
SYMI	61	-2635.6539	5662.0266	0.0000	1.0000
HKY	59	-2665.8727	5709.6538	0.0000	1.0000
GTR	63	-2658.1725	5719.8744	0.0000	1.0000
K2P	56	-2685.8130	5730.3189	0.0000	1.0000
SYM	60	-2682.4017	5749.1171	0.0000	1.0000
F81G	59	-2699.6136	5777.1358	0.0000	1.0000
F81IG	60	-2699.2969	5782.9076	0.0000	1.0000
F81I	59	-2702.8056	5783.5196	0.0000	1.0000
JC69G	56	-2712.4151	5783.5229	0.0000	1.0000
JC69IG	57	-2712.1584	5789.4148	0.0000	1.0000
JC69I	56	-2715.3760	5789.4449	0.0000	1.0000
F81	58	-2748.4761	5868.4554	0.0000	1.0000
JC69	55	-2761.5253	5875.3381	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/11141_hg18.phylip = 1 - 605;
 Partition Dummy = 1:261-aggressive-phylip/11141_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/11141_hg18.phylip = 1 - 605;
 Partition Dummy = 1:261-aggressive-phylip/11141_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/11141_hg18.phylip = 1 - 605;
 Partition Dummy = 1:261-aggressive-phylip/11141_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

