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

Input data from file "261-aggressive-phylip/12887_hg18.phylip" (29 taxa, 859 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	-4223.3658	8587.5514	0.8928	0.8928
HKYIG	61	-4230.4287	8592.3479	0.0811	0.9740
GTRG	64	-4228.1335	8594.7455	0.0245	0.9984
HKYG	60	-4235.5321	8600.2372	0.0016	1.0000
SYMIG	62	-4243.6722	8621.1585	0.0000	1.0000
K2PIG	58	-4250.9227	8626.4004	0.0000	1.0000
SYMG	61	-4248.9755	8629.4416	0.0000	1.0000
K2PG	57	-4256.4786	8635.2118	0.0000	1.0000
GTRI	64	-4255.4134	8649.3054	0.0000	1.0000
HKYI	60	-4262.6324	8654.4377	0.0000	1.0000
SYMI	61	-4274.2146	8679.9197	0.0000	1.0000
K2PI	57	-4282.4358	8687.1262	0.0000	1.0000
F81IG	60	-4353.8640	8836.9010	0.0000	1.0000
F81G	59	-4358.8313	8844.5236	0.0000	1.0000
JC69IG	57	-4370.3211	8862.8970	0.0000	1.0000
JC69G	56	-4375.2872	8870.5344	0.0000	1.0000
F81I	59	-4384.4821	8895.8253	0.0000	1.0000
JC69I	56	-4399.7920	8919.5441	0.0000	1.0000
GTR	63	-4404.8655	8945.8744	0.0000	1.0000
HKY	59	-4411.5113	8949.8836	0.0000	1.0000
SYM	60	-4422.6128	8974.3985	0.0000	1.0000
K2P	56	-4431.0601	8982.0803	0.0000	1.0000
F81	58	-4530.5493	9185.6537	0.0000	1.0000
JC69	55	-4545.2328	9208.1368	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	65	-4223.3658	8576.7317	0.9340	0.9340
HKYIG	61	-4230.4287	8582.8573	0.0437	0.9777
GTRG	64	-4228.1335	8584.2670	0.0216	0.9993
HKYG	60	-4235.5321	8591.0642	0.0007	1.0000
SYMIG	62	-4243.6722	8611.3444	0.0000	1.0000
K2PIG	58	-4250.9227	8617.8454	0.0000	1.0000
SYMG	61	-4248.9755	8619.9510	0.0000	1.0000
K2PG	57	-4256.4786	8626.9572	0.0000	1.0000
GTRI	64	-4255.4134	8638.8268	0.0000	1.0000
HKYI	60	-4262.6324	8645.2648	0.0000	1.0000
SYMI	61	-4274.2146	8670.4291	0.0000	1.0000
K2PI	57	-4282.4358	8678.8715	0.0000	1.0000
F81IG	60	-4353.8640	8827.7280	0.0000	1.0000
F81G	59	-4358.8313	8835.6625	0.0000	1.0000
JC69IG	57	-4370.3211	8854.6423	0.0000	1.0000
JC69G	56	-4375.2872	8862.5743	0.0000	1.0000
F81I	59	-4384.4821	8886.9643	0.0000	1.0000
JC69I	56	-4399.7920	8911.5840	0.0000	1.0000
GTR	63	-4404.8655	8935.7310	0.0000	1.0000
HKY	59	-4411.5113	8941.0225	0.0000	1.0000
SYM	60	-4422.6128	8965.2255	0.0000	1.0000
K2P	56	-4431.0601	8974.1202	0.0000	1.0000
F81	58	-4530.5493	9177.0987	0.0000	1.0000
JC69	55	-4545.2328	9200.4656	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYIG	61	-4230.4287	8872.9592	0.8474	0.8474
HKYG	60	-4235.5321	8876.4104	0.1509	0.9983
GTRIG	65	-4223.3658	8885.8567	0.0013	0.9996
GTRG	64	-4228.1335	8888.6362	0.0003	1.0000
K2PIG	58	-4250.9227	8893.6800	0.0000	1.0000
K2PG	57	-4256.4786	8898.0360	0.0000	1.0000
SYMIG	62	-4243.6722	8906.2021	0.0000	1.0000
SYMG	61	-4248.9755	8910.0529	0.0000	1.0000
HKYI	60	-4262.6324	8930.6109	0.0000	1.0000
GTRI	64	-4255.4134	8943.1960	0.0000	1.0000
K2PI	57	-4282.4358	8949.9503	0.0000	1.0000
SYMI	61	-4274.2146	8960.5310	0.0000	1.0000
F81IG	60	-4353.8640	9113.0742	0.0000	1.0000
F81G	59	-4358.8313	9116.2529	0.0000	1.0000
JC69IG	57	-4370.3211	9125.7211	0.0000	1.0000
JC69G	56	-4375.2872	9128.8974	0.0000	1.0000
F81I	59	-4384.4821	9167.5546	0.0000	1.0000
JC69I	56	-4399.7920	9177.9071	0.0000	1.0000
HKY	59	-4411.5113	9221.6129	0.0000	1.0000
GTR	63	-4404.8655	9235.3444	0.0000	1.0000
K2P	56	-4431.0601	9240.4433	0.0000	1.0000
SYM	60	-4422.6128	9250.5717	0.0000	1.0000
F81	58	-4530.5493	9452.9333	0.0000	1.0000
JC69	55	-4545.2328	9462.0329	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/12887_hg18.phylip = 1 - 859;
 Partition Dummy = 1:261-aggressive-phylip/12887_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/12887_hg18.phylip = 1 - 859;
 Partition Dummy = 1:261-aggressive-phylip/12887_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/12887_hg18.phylip = 1 - 859;
 Partition Dummy = 1:261-aggressive-phylip/12887_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

