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

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

Minimum AIC  model: GTRIG
Minimum AICc model: GTRIG
Minimum BIC  model: HKYIG


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	-4355.2508	8854.8255	0.9230	0.9230
GTRG	64	-4358.9696	8859.8058	0.0765	0.9995
HKYIG	61	-4367.7021	8869.9481	0.0005	1.0000
HKYG	60	-4373.9603	8880.0399	0.0000	1.0000
GTRI	64	-4374.5624	8890.9914	0.0000	1.0000
HKYI	60	-4387.1339	8906.3870	0.0000	1.0000
SYMG	61	-4392.5992	8919.7422	0.0000	1.0000
K2PIG	58	-4397.2778	8921.8494	0.0000	1.0000
SYMIG	62	-4397.2822	8931.5411	0.0000	1.0000
K2PG	57	-4403.5535	8931.9998	0.0000	1.0000
K2PI	57	-4411.7142	8948.3214	0.0000	1.0000
SYMI	61	-4409.9678	8954.4795	0.0000	1.0000
F81G	59	-4456.6922	9043.0869	0.0000	1.0000
F81IG	60	-4460.0564	9052.2320	0.0000	1.0000
GTR	63	-4467.1488	9073.7152	0.0000	1.0000
F81I	59	-4476.7940	9083.2905	0.0000	1.0000
HKY	59	-4479.2440	9088.1904	0.0000	1.0000
JC69G	56	-4485.6182	9093.7363	0.0000	1.0000
JC69IG	57	-4487.3359	9099.5648	0.0000	1.0000
K2P	56	-4507.9184	9138.3367	0.0000	1.0000
JC69I	56	-4508.2919	9139.0838	0.0000	1.0000
SYM	60	-4504.1661	9140.4513	0.0000	1.0000
F81	58	-4567.2033	9261.7003	0.0000	1.0000
JC69	55	-4595.9459	9312.0066	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	65	-4355.2508	8840.5017	0.9379	0.9379
GTRG	64	-4358.9696	8845.9392	0.0619	0.9998
HKYIG	61	-4367.7021	8857.4041	0.0002	1.0000
HKYG	60	-4373.9603	8867.9207	0.0000	1.0000
GTRI	64	-4374.5624	8877.1248	0.0000	1.0000
HKYI	60	-4387.1339	8894.2678	0.0000	1.0000
SYMG	61	-4392.5992	8907.1983	0.0000	1.0000
K2PIG	58	-4397.2778	8910.5557	0.0000	1.0000
SYMIG	62	-4397.2822	8918.5644	0.0000	1.0000
K2PG	57	-4403.5535	8921.1069	0.0000	1.0000
K2PI	57	-4411.7142	8937.4285	0.0000	1.0000
SYMI	61	-4409.9678	8941.9355	0.0000	1.0000
F81G	59	-4456.6922	9031.3844	0.0000	1.0000
F81IG	60	-4460.0564	9040.1128	0.0000	1.0000
GTR	63	-4467.1488	9060.2975	0.0000	1.0000
F81I	59	-4476.7940	9071.5880	0.0000	1.0000
HKY	59	-4479.2440	9076.4880	0.0000	1.0000
JC69G	56	-4485.6182	9083.2363	0.0000	1.0000
JC69IG	57	-4487.3359	9088.6718	0.0000	1.0000
K2P	56	-4507.9184	9127.8367	0.0000	1.0000
SYM	60	-4504.1661	9128.3321	0.0000	1.0000
JC69I	56	-4508.2919	9128.5838	0.0000	1.0000
F81	58	-4567.2033	9250.4066	0.0000	1.0000
JC69	55	-4595.9459	9301.8917	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYIG	61	-4367.7021	9131.8911	0.5028	0.5028
GTRIG	65	-4355.2508	9132.9878	0.2906	0.7934
GTRG	64	-4358.9696	9133.9256	0.1818	0.9752
HKYG	60	-4373.9603	9137.9079	0.0248	1.0000
HKYI	60	-4387.1339	9164.2550	0.0000	1.0000
GTRI	64	-4374.5624	9165.1112	0.0000	1.0000
K2PIG	58	-4397.2778	9171.5433	0.0000	1.0000
K2PG	57	-4403.5535	9177.5948	0.0000	1.0000
SYMG	61	-4392.5992	9181.6853	0.0000	1.0000
K2PI	57	-4411.7142	9193.9164	0.0000	1.0000
SYMIG	62	-4397.2822	9197.5512	0.0000	1.0000
SYMI	61	-4409.9678	9216.4225	0.0000	1.0000
F81G	59	-4456.6922	9296.8718	0.0000	1.0000
F81IG	60	-4460.0564	9310.1000	0.0000	1.0000
JC69G	56	-4485.6182	9335.2244	0.0000	1.0000
F81I	59	-4476.7940	9337.0754	0.0000	1.0000
HKY	59	-4479.2440	9341.9754	0.0000	1.0000
GTR	63	-4467.1488	9343.7841	0.0000	1.0000
JC69IG	57	-4487.3359	9345.1597	0.0000	1.0000
K2P	56	-4507.9184	9379.8248	0.0000	1.0000
JC69I	56	-4508.2919	9380.5718	0.0000	1.0000
SYM	60	-4504.1661	9398.3194	0.0000	1.0000
F81	58	-4567.2033	9511.3942	0.0000	1.0000
JC69	55	-4595.9459	9549.3800	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/3313_hg18.phylip = 1 - 665;
 Partition Dummy = 1:261-aggressive-phylip/3313_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/3313_hg18.phylip = 1 - 665;
 Partition Dummy = 1:261-aggressive-phylip/3313_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 (HKYIG)]

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


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

