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

Input data from file "261-aggressive-phylip/8720_hg18.phylip" (29 taxa, 783 characters)

Minimum AIC  model: GTRIG
Minimum AICc model: GTRIG
Minimum BIC  model: GTRG


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
GTRIG	65	-2704.2525	5550.4715	0.6598	0.6598
GTRG	64	-2706.1166	5551.8210	0.3360	0.9958
GTRI	64	-2710.5082	5560.6042	0.0042	1.0000
HKYIG	61	-2721.6089	5575.7087	0.0000	1.0000
HKYG	60	-2722.7921	5575.7226	0.0000	1.0000
SYMIG	62	-2723.6957	5582.2414	0.0000	1.0000
HKYI	60	-2726.1544	5582.4473	0.0000	1.0000
SYMG	61	-2725.2790	5583.0490	0.0000	1.0000
SYMI	61	-2729.2195	5590.9299	0.0000	1.0000
K2PG	57	-2740.6312	5604.3823	0.0000	1.0000
K2PIG	58	-2739.6070	5604.6670	0.0000	1.0000
K2PI	57	-2743.7286	5610.5773	0.0000	1.0000
F81G	59	-2779.5601	5686.9128	0.0000	1.0000
F81IG	60	-2778.5098	5687.1582	0.0000	1.0000
F81I	59	-2782.6765	5693.1455	0.0000	1.0000
GTR	63	-2790.3788	5717.9732	0.0000	1.0000
JC69G	56	-2798.6719	5718.1371	0.0000	1.0000
JC69IG	57	-2797.7983	5718.7166	0.0000	1.0000
JC69I	56	-2801.6714	5724.1361	0.0000	1.0000
HKY	59	-2803.2563	5734.3051	0.0000	1.0000
SYM	60	-2807.0817	5744.3018	0.0000	1.0000
K2P	56	-2819.6461	5760.0855	0.0000	1.0000
F81	58	-2856.4944	5838.4418	0.0000	1.0000
JC69	55	-2874.8241	5868.1213	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	65	-2704.2525	5538.5050	0.7010	0.7010
GTRG	64	-2706.1166	5540.2333	0.2954	0.9963
GTRI	64	-2710.5082	5549.0164	0.0037	1.0000
HKYIG	61	-2721.6089	5565.2177	0.0000	1.0000
HKYG	60	-2722.7921	5565.5841	0.0000	1.0000
SYMIG	62	-2723.6957	5571.3914	0.0000	1.0000
HKYI	60	-2726.1544	5572.3088	0.0000	1.0000
SYMG	61	-2725.2790	5572.5580	0.0000	1.0000
SYMI	61	-2729.2195	5580.4389	0.0000	1.0000
K2PIG	58	-2739.6070	5595.2139	0.0000	1.0000
K2PG	57	-2740.6312	5595.2623	0.0000	1.0000
K2PI	57	-2743.7286	5601.4573	0.0000	1.0000
F81IG	60	-2778.5098	5677.0197	0.0000	1.0000
F81G	59	-2779.5601	5677.1202	0.0000	1.0000
F81I	59	-2782.6765	5683.3529	0.0000	1.0000
GTR	63	-2790.3788	5706.7576	0.0000	1.0000
JC69G	56	-2798.6719	5709.3437	0.0000	1.0000
JC69IG	57	-2797.7983	5709.5966	0.0000	1.0000
JC69I	56	-2801.6714	5715.3428	0.0000	1.0000
HKY	59	-2803.2563	5724.5126	0.0000	1.0000
SYM	60	-2807.0817	5734.1633	0.0000	1.0000
K2P	56	-2819.6461	5751.2921	0.0000	1.0000
F81	58	-2856.4944	5828.9887	0.0000	1.0000
JC69	55	-2874.8241	5859.6481	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	64	-2706.1166	5838.6738	0.7791	0.7791
GTRIG	65	-2704.2525	5841.6086	0.1796	0.9587
HKYG	60	-2722.7921	5845.3721	0.0274	0.9861
GTRI	64	-2710.5082	5847.4569	0.0096	0.9958
HKYIG	61	-2721.6089	5849.6688	0.0032	0.9989
HKYI	60	-2726.1544	5852.0968	0.0009	0.9999
SYMG	61	-2725.2790	5857.0091	0.0001	1.0000
SYMIG	62	-2723.6957	5860.5056	0.0000	1.0000
K2PG	57	-2740.6312	5861.0609	0.0000	1.0000
SYMI	61	-2729.2195	5864.8900	0.0000	1.0000
K2PIG	58	-2739.6070	5865.6756	0.0000	1.0000
K2PI	57	-2743.7286	5867.2558	0.0000	1.0000
F81G	59	-2779.5601	5952.2450	0.0000	1.0000
F81IG	60	-2778.5098	5956.8076	0.0000	1.0000
F81I	59	-2782.6765	5958.4777	0.0000	1.0000
JC69G	56	-2798.6719	5970.4791	0.0000	1.0000
JC69IG	57	-2797.7983	5975.3952	0.0000	1.0000
JC69I	56	-2801.6714	5976.4782	0.0000	1.0000
HKY	59	-2803.2563	5999.6374	0.0000	1.0000
GTR	63	-2790.3788	6000.5350	0.0000	1.0000
K2P	56	-2819.6461	6012.4276	0.0000	1.0000
SYM	60	-2807.0817	6013.9513	0.0000	1.0000
F81	58	-2856.4944	6099.4504	0.0000	1.0000
JC69	55	-2874.8241	6116.1204	0.0000	1.0000



To specify default values for best AIC, AICc, or BIC models in MrBayes:

[Mrbayes block for the best AIC model (GTRIG)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/8720_hg18.phylip = 1 - 783;
 Partition Dummy = 1:261-aggressive-phylip/8720_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


[Mrbayes block for the best AICc model (GTRIG)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/8720_hg18.phylip = 1 - 783;
 Partition Dummy = 1:261-aggressive-phylip/8720_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


[Mrbayes block for the best BIC model (GTRG)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/8720_hg18.phylip = 1 - 783;
 Partition Dummy = 1:261-aggressive-phylip/8720_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;


-------------------------------------------------------------
End of Output

