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

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

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


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	-4094.6283	8331.9865	0.9999	0.9999
HKYIG	61	-4109.4890	8352.1344	0.0000	1.0000
GTRG	64	-4106.5146	8353.3551	0.0000	1.0000
HKYG	60	-4115.6089	8361.9983	0.0000	1.0000
K2PIG	58	-4137.9924	8402.0347	0.0000	1.0000
SYMIG	62	-4137.8984	8411.3359	0.0000	1.0000
K2PG	57	-4145.7564	8415.2079	0.0000	1.0000
SYMG	61	-4145.5091	8424.1745	0.0000	1.0000
GTRI	64	-4166.0526	8472.4311	0.0000	1.0000
HKYI	60	-4171.2818	8473.3441	0.0000	1.0000
K2PI	57	-4208.1025	8539.9000	0.0000	1.0000
SYMI	61	-4209.0908	8551.3380	0.0000	1.0000
F81IG	60	-4229.7639	8590.3084	0.0000	1.0000
F81G	59	-4234.2051	8596.8220	0.0000	1.0000
JC69IG	57	-4248.9583	8621.6117	0.0000	1.0000
JC69G	56	-4252.1998	8625.7466	0.0000	1.0000
F81I	59	-4289.7963	8708.0043	0.0000	1.0000
JC69I	56	-4305.4617	8732.2703	0.0000	1.0000
GTR	63	-4477.0213	9091.9717	0.0000	1.0000
HKY	59	-4484.9633	9098.3383	0.0000	1.0000
K2P	56	-4513.4504	9148.2477	0.0000	1.0000
SYM	60	-4514.9503	9160.6811	0.0000	1.0000
F81	58	-4590.6113	9307.2726	0.0000	1.0000
JC69	55	-4608.9622	9336.9303	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	65	-4094.6283	8319.2566	1.0000	1.0000
HKYIG	61	-4109.4890	8340.9780	0.0000	1.0000
GTRG	64	-4106.5146	8341.0292	0.0000	1.0000
HKYG	60	-4115.6089	8351.2177	0.0000	1.0000
K2PIG	58	-4137.9924	8391.9848	0.0000	1.0000
SYMIG	62	-4137.8984	8399.7968	0.0000	1.0000
K2PG	57	-4145.7564	8405.5128	0.0000	1.0000
SYMG	61	-4145.5091	8413.0182	0.0000	1.0000
GTRI	64	-4166.0526	8460.1051	0.0000	1.0000
HKYI	60	-4171.2818	8462.5635	0.0000	1.0000
K2PI	57	-4208.1025	8530.2050	0.0000	1.0000
SYMI	61	-4209.0908	8540.1816	0.0000	1.0000
F81IG	60	-4229.7639	8579.5279	0.0000	1.0000
F81G	59	-4234.2051	8586.4102	0.0000	1.0000
JC69IG	57	-4248.9583	8611.9167	0.0000	1.0000
JC69G	56	-4252.1998	8616.3996	0.0000	1.0000
F81I	59	-4289.7963	8697.5925	0.0000	1.0000
JC69I	56	-4305.4617	8722.9233	0.0000	1.0000
GTR	63	-4477.0213	9080.0427	0.0000	1.0000
HKY	59	-4484.9633	9087.9266	0.0000	1.0000
K2P	56	-4513.4504	9138.9007	0.0000	1.0000
SYM	60	-4514.9503	9149.9006	0.0000	1.0000
F81	58	-4590.6113	9297.2227	0.0000	1.0000
JC69	55	-4608.9622	9327.9245	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRIG	65	-4094.6283	8618.6888	0.8304	0.8304
HKYIG	61	-4109.4890	8621.9837	0.1599	0.9903
HKYG	60	-4115.6089	8627.6168	0.0096	0.9998
GTRG	64	-4106.5146	8635.8548	0.0002	1.0000
K2PIG	58	-4137.9924	8659.1705	0.0000	1.0000
K2PG	57	-4145.7564	8668.0919	0.0000	1.0000
SYMIG	62	-4137.8984	8685.4091	0.0000	1.0000
SYMG	61	-4145.5091	8694.0239	0.0000	1.0000
HKYI	60	-4171.2818	8738.9626	0.0000	1.0000
GTRI	64	-4166.0526	8754.9308	0.0000	1.0000
K2PI	57	-4208.1025	8792.7841	0.0000	1.0000
SYMI	61	-4209.0908	8821.1873	0.0000	1.0000
F81IG	60	-4229.7639	8855.9269	0.0000	1.0000
F81G	59	-4234.2051	8858.2026	0.0000	1.0000
JC69G	56	-4252.1998	8874.3720	0.0000	1.0000
JC69IG	57	-4248.9583	8874.4958	0.0000	1.0000
F81I	59	-4289.7963	8969.3849	0.0000	1.0000
JC69I	56	-4305.4617	8980.8958	0.0000	1.0000
HKY	59	-4484.9633	9359.7189	0.0000	1.0000
GTR	63	-4477.0213	9370.2616	0.0000	1.0000
K2P	56	-4513.4504	9396.8732	0.0000	1.0000
SYM	60	-4514.9503	9426.2996	0.0000	1.0000
F81	58	-4590.6113	9564.4084	0.0000	1.0000
JC69	55	-4608.9622	9581.2902	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/8721_hg18.phylip = 1 - 740;
 Partition Dummy = 1:261-aggressive-phylip/8721_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/8721_hg18.phylip = 1 - 740;
 Partition Dummy = 1:261-aggressive-phylip/8721_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 (GTRIG)]

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


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

