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

Input data from file "261-aggressive-phylip/6297_hg18.phylip" (29 taxa, 998 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	-6545.1283	13229.4625	0.9490	0.9490
GTRG	64	-6549.1999	13235.3172	0.0508	0.9999
HKYIG	61	-6558.9839	13248.0490	0.0001	0.9999
HKYG	60	-6560.6965	13249.2052	0.0000	1.0000
SYMIG	62	-6561.1249	13254.6049	0.0000	1.0000
SYMG	61	-6563.0753	13256.2319	0.0000	1.0000
K2PIG	58	-6577.3178	13277.9242	0.0000	1.0000
K2PG	57	-6580.3770	13281.7880	0.0000	1.0000
GTRI	64	-6659.4367	13455.7908	0.0000	1.0000
HKYI	60	-6672.8511	13473.5144	0.0000	1.0000
SYMI	61	-6679.7420	13489.5652	0.0000	1.0000
K2PI	57	-6695.5326	13512.0993	0.0000	1.0000
F81G	59	-6835.9391	13797.4261	0.0000	1.0000
F81IG	60	-6836.1049	13800.0219	0.0000	1.0000
JC69IG	57	-6873.9902	13869.0145	0.0000	1.0000
JC69G	56	-6875.1938	13869.1719	0.0000	1.0000
F81I	59	-6925.6924	13976.9328	0.0000	1.0000
JC69I	56	-6976.3709	14071.5262	0.0000	1.0000
GTR	63	-7046.1779	14226.9897	0.0000	1.0000
HKY	59	-7073.5823	14272.7126	0.0000	1.0000
SYM	60	-7075.1650	14278.1421	0.0000	1.0000
K2P	56	-7106.1249	14331.0341	0.0000	1.0000
F81	58	-7324.7702	14772.8290	0.0000	1.0000
JC69	55	-7358.0131	14832.5655	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	65	-6545.1283	13220.2565	0.9556	0.9556
GTRG	64	-6549.1999	13226.3998	0.0443	0.9999
HKYIG	61	-6558.9839	13239.9678	0.0001	1.0000
HKYG	60	-6560.6965	13241.3930	0.0000	1.0000
SYMIG	62	-6561.1249	13246.2498	0.0000	1.0000
SYMG	61	-6563.0753	13248.1507	0.0000	1.0000
K2PIG	58	-6577.3178	13270.6356	0.0000	1.0000
K2PG	57	-6580.3770	13274.7540	0.0000	1.0000
GTRI	64	-6659.4367	13446.8734	0.0000	1.0000
HKYI	60	-6672.8511	13465.7022	0.0000	1.0000
SYMI	61	-6679.7420	13481.4840	0.0000	1.0000
K2PI	57	-6695.5326	13505.0652	0.0000	1.0000
F81G	59	-6835.9391	13789.8782	0.0000	1.0000
F81IG	60	-6836.1049	13792.2097	0.0000	1.0000
JC69IG	57	-6873.9902	13861.9805	0.0000	1.0000
JC69G	56	-6875.1938	13862.3876	0.0000	1.0000
F81I	59	-6925.6924	13969.3848	0.0000	1.0000
JC69I	56	-6976.3709	14064.7419	0.0000	1.0000
GTR	63	-7046.1779	14218.3559	0.0000	1.0000
HKY	59	-7073.5823	14265.1646	0.0000	1.0000
SYM	60	-7075.1650	14270.3299	0.0000	1.0000
K2P	56	-7106.1249	14324.2498	0.0000	1.0000
F81	58	-7324.7702	14765.5404	0.0000	1.0000
JC69	55	-7358.0131	14826.0262	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-6560.6965	13535.7382	0.6843	0.6843
GTRIG	65	-6545.1283	13539.1305	0.1255	0.8098
HKYIG	61	-6558.9839	13539.2187	0.1201	0.9299
GTRG	64	-6549.1999	13540.3680	0.0676	0.9974
SYMG	61	-6563.0753	13547.4016	0.0020	0.9995
SYMIG	62	-6561.1249	13550.4065	0.0004	0.9999
K2PG	57	-6580.3770	13554.3819	0.0001	1.0000
K2PIG	58	-6577.3178	13555.1693	0.0000	1.0000
HKYI	60	-6672.8511	13760.0474	0.0000	1.0000
GTRI	64	-6659.4367	13760.8416	0.0000	1.0000
SYMI	61	-6679.7420	13780.7350	0.0000	1.0000
K2PI	57	-6695.5326	13784.6932	0.0000	1.0000
F81G	59	-6835.9391	14079.3176	0.0000	1.0000
F81IG	60	-6836.1049	14086.5549	0.0000	1.0000
JC69G	56	-6875.1938	14137.1098	0.0000	1.0000
JC69IG	57	-6873.9902	14141.6084	0.0000	1.0000
F81I	59	-6925.6924	14258.8242	0.0000	1.0000
JC69I	56	-6976.3709	14339.4641	0.0000	1.0000
GTR	63	-7046.1779	14527.4183	0.0000	1.0000
HKY	59	-7073.5823	14554.6040	0.0000	1.0000
SYM	60	-7075.1650	14564.6751	0.0000	1.0000
K2P	56	-7106.1249	14598.9720	0.0000	1.0000
F81	58	-7324.7702	15050.0741	0.0000	1.0000
JC69	55	-7358.0131	15095.8426	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/6297_hg18.phylip = 1 - 998;
 Partition Dummy = 1:261-aggressive-phylip/6297_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/6297_hg18.phylip = 1 - 998;
 Partition Dummy = 1:261-aggressive-phylip/6297_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/6297_hg18.phylip = 1 - 998;
 Partition Dummy = 1:261-aggressive-phylip/6297_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

