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

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

Minimum AIC  model: GTRIG
Minimum AICc model: GTRIG
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
GTRIG	65	-2638.7877	5423.0071	0.5399	0.5399
HKYIG	61	-2644.4557	5424.4186	0.2666	0.8065
GTRG	64	-2641.3911	5425.7193	0.1391	0.9456
HKYG	60	-2647.2741	5427.5963	0.0544	1.0000
K2PIG	58	-2658.3375	5444.8313	0.0000	1.0000
GTRI	64	-2652.4442	5447.8255	0.0000	1.0000
K2PG	57	-2661.4224	5448.5681	0.0000	1.0000
HKYI	60	-2657.9573	5448.9628	0.0000	1.0000
SYMIG	62	-2656.2824	5450.5398	0.0000	1.0000
SYMG	61	-2659.3033	5454.1138	0.0000	1.0000
K2PI	57	-2672.4010	5470.5253	0.0000	1.0000
SYMI	61	-2670.2209	5475.9489	0.0000	1.0000
F81IG	60	-2715.1751	5563.3982	0.0000	1.0000
F81G	59	-2717.6953	5565.9884	0.0000	1.0000
JC69IG	57	-2726.9346	5579.5926	0.0000	1.0000
JC69G	56	-2729.5547	5582.4086	0.0000	1.0000
F81I	59	-2727.5684	5585.7347	0.0000	1.0000
JC69I	56	-2739.8193	5602.9378	0.0000	1.0000
GTR	63	-2734.3647	5609.1810	0.0000	1.0000
HKY	59	-2742.5775	5615.7528	0.0000	1.0000
K2P	56	-2760.6153	5644.5298	0.0000	1.0000
SYM	60	-2756.2089	5645.4659	0.0000	1.0000
F81	58	-2812.0582	5752.2728	0.0000	1.0000
JC69	55	-2825.8017	5772.4867	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	65	-2638.7877	5407.5755	0.7040	0.7040
GTRG	64	-2641.3911	5410.7821	0.1417	0.8457
HKYIG	61	-2644.4557	5410.9115	0.1328	0.9784
HKYG	60	-2647.2741	5414.5482	0.0216	1.0000
K2PIG	58	-2658.3375	5432.6750	0.0000	1.0000
GTRI	64	-2652.4442	5432.8884	0.0000	1.0000
HKYI	60	-2657.9573	5435.9147	0.0000	1.0000
SYMIG	62	-2656.2824	5436.5649	0.0000	1.0000
K2PG	57	-2661.4224	5436.8447	0.0000	1.0000
SYMG	61	-2659.3033	5440.6067	0.0000	1.0000
K2PI	57	-2672.4010	5458.8019	0.0000	1.0000
SYMI	61	-2670.2209	5462.4417	0.0000	1.0000
F81IG	60	-2715.1751	5550.3501	0.0000	1.0000
F81G	59	-2717.6953	5553.3906	0.0000	1.0000
JC69IG	57	-2726.9346	5567.8692	0.0000	1.0000
JC69G	56	-2729.5547	5571.1095	0.0000	1.0000
F81I	59	-2727.5684	5573.1368	0.0000	1.0000
JC69I	56	-2739.8193	5591.6387	0.0000	1.0000
GTR	63	-2734.3647	5594.7294	0.0000	1.0000
HKY	59	-2742.5775	5603.1549	0.0000	1.0000
SYM	60	-2756.2089	5632.4178	0.0000	1.0000
K2P	56	-2760.6153	5633.2306	0.0000	1.0000
F81	58	-2812.0582	5740.1165	0.0000	1.0000
JC69	55	-2825.8017	5761.6033	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-2647.2741	5680.5246	0.5904	0.5904
HKYIG	61	-2644.4557	5681.3208	0.3965	0.9868
K2PG	57	-2661.4224	5689.5223	0.0066	0.9934
K2PIG	58	-2658.3375	5689.7855	0.0058	0.9991
GTRG	64	-2641.3911	5694.4903	0.0005	0.9997
GTRIG	65	-2638.7877	5695.7166	0.0003	1.0000
HKYI	60	-2657.9573	5701.8911	0.0000	1.0000
SYMG	61	-2659.3033	5711.0160	0.0000	1.0000
SYMIG	62	-2656.2824	5711.4072	0.0000	1.0000
K2PI	57	-2672.4010	5711.4795	0.0000	1.0000
GTRI	64	-2652.4442	5716.5965	0.0000	1.0000
SYMI	61	-2670.2209	5732.8511	0.0000	1.0000
F81G	59	-2717.6953	5814.9340	0.0000	1.0000
F81IG	60	-2715.1751	5816.3265	0.0000	1.0000
JC69G	56	-2729.5547	5819.3541	0.0000	1.0000
JC69IG	57	-2726.9346	5820.5468	0.0000	1.0000
F81I	59	-2727.5684	5834.6803	0.0000	1.0000
JC69I	56	-2739.8193	5839.8833	0.0000	1.0000
HKY	59	-2742.5775	5864.6984	0.0000	1.0000
GTR	63	-2734.3647	5874.0046	0.0000	1.0000
K2P	56	-2760.6153	5881.4753	0.0000	1.0000
SYM	60	-2756.2089	5898.3942	0.0000	1.0000
F81	58	-2812.0582	5997.2270	0.0000	1.0000
JC69	55	-2825.8017	6005.4150	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/689_hg18.phylip = 1 - 622;
 Partition Dummy = 1:261-aggressive-phylip/689_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/689_hg18.phylip = 1 - 622;
 Partition Dummy = 1:261-aggressive-phylip/689_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 (HKYG)]

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

