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

Input data from file "261-aggressive-phylip/1996_hg18.phylip" (29 taxa, 536 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	-2070.4798	4289.2149	0.5602	0.5602
GTRG	64	-2072.0323	4289.7290	0.4332	0.9934
HKYIG	61	-2081.0488	4300.0555	0.0025	0.9958
HKYG	60	-2082.4815	4300.3736	0.0021	0.9979
GTRI	64	-2077.4063	4300.4772	0.0020	1.0000
HKYI	60	-2086.2989	4308.0082	0.0000	1.0000
SYMG	61	-2096.6177	4331.1931	0.0000	1.0000
SYMIG	62	-2096.0772	4332.6702	0.0000	1.0000
SYMI	61	-2099.7897	4337.5373	0.0000	1.0000
K2PG	57	-2105.9695	4339.7717	0.0000	1.0000
K2PIG	58	-2105.1727	4340.6933	0.0000	1.0000
K2PI	57	-2108.6114	4345.0555	0.0000	1.0000
F81G	59	-2159.8283	4452.5305	0.0000	1.0000
F81IG	60	-2158.8468	4453.1042	0.0000	1.0000
F81I	59	-2162.4237	4457.7213	0.0000	1.0000
GTR	63	-2157.9731	4459.0309	0.0000	1.0000
HKY	59	-2168.2347	4469.3434	0.0000	1.0000
JC69G	56	-2183.2546	4491.8370	0.0000	1.0000
JC69IG	57	-2182.4546	4492.7419	0.0000	1.0000
JC69I	56	-2185.8135	4496.9547	0.0000	1.0000
SYM	60	-2183.5257	4502.4618	0.0000	1.0000
K2P	56	-2193.6545	4512.6367	0.0000	1.0000
F81	58	-2243.1027	4616.5534	0.0000	1.0000
JC69	55	-2266.1024	4655.0380	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	65	-2070.4798	4270.9596	0.6327	0.6327
GTRG	64	-2072.0323	4272.0645	0.3641	0.9968
GTRI	64	-2077.4063	4282.8126	0.0017	0.9985
HKYIG	61	-2081.0488	4284.0977	0.0009	0.9994
HKYG	60	-2082.4815	4284.9630	0.0006	1.0000
HKYI	60	-2086.2989	4292.5977	0.0000	1.0000
SYMG	61	-2096.6177	4315.2353	0.0000	1.0000
SYMIG	62	-2096.0772	4316.1544	0.0000	1.0000
SYMI	61	-2099.7897	4321.5795	0.0000	1.0000
K2PG	57	-2105.9695	4325.9391	0.0000	1.0000
K2PIG	58	-2105.1727	4326.3453	0.0000	1.0000
K2PI	57	-2108.6114	4331.2228	0.0000	1.0000
F81G	59	-2159.8283	4437.6566	0.0000	1.0000
F81IG	60	-2158.8468	4437.6937	0.0000	1.0000
GTR	63	-2157.9731	4441.9462	0.0000	1.0000
F81I	59	-2162.4237	4442.8473	0.0000	1.0000
HKY	59	-2168.2347	4454.4695	0.0000	1.0000
JC69G	56	-2183.2546	4478.5093	0.0000	1.0000
JC69IG	57	-2182.4546	4478.9093	0.0000	1.0000
JC69I	56	-2185.8135	4483.6269	0.0000	1.0000
SYM	60	-2183.5257	4487.0513	0.0000	1.0000
K2P	56	-2193.6545	4499.3089	0.0000	1.0000
F81	58	-2243.1027	4602.2054	0.0000	1.0000
JC69	55	-2266.1024	4642.2047	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-2082.4815	4542.0111	0.7417	0.7417
HKYIG	61	-2081.0488	4545.4298	0.1342	0.8760
GTRG	64	-2072.0323	4546.2491	0.0891	0.9651
GTRIG	65	-2070.4798	4549.4283	0.0182	0.9833
HKYI	60	-2086.2989	4549.6458	0.0163	0.9996
GTRI	64	-2077.4063	4556.9972	0.0004	1.0000
K2PG	57	-2105.9695	4570.1347	0.0000	1.0000
K2PIG	58	-2105.1727	4574.8251	0.0000	1.0000
K2PI	57	-2108.6114	4575.4185	0.0000	1.0000
SYMG	61	-2096.6177	4576.5675	0.0000	1.0000
SYMIG	62	-2096.0772	4581.7707	0.0000	1.0000
SYMI	61	-2099.7897	4582.9117	0.0000	1.0000
F81G	59	-2159.8283	4690.4205	0.0000	1.0000
F81IG	60	-2158.8468	4694.7417	0.0000	1.0000
F81I	59	-2162.4237	4695.6113	0.0000	1.0000
HKY	59	-2168.2347	4707.2334	0.0000	1.0000
GTR	63	-2157.9731	4711.8466	0.0000	1.0000
JC69G	56	-2183.2546	4718.4208	0.0000	1.0000
JC69IG	57	-2182.4546	4723.1049	0.0000	1.0000
JC69I	56	-2185.8135	4723.5385	0.0000	1.0000
K2P	56	-2193.6545	4739.2205	0.0000	1.0000
SYM	60	-2183.5257	4744.0994	0.0000	1.0000
F81	58	-2243.1027	4850.6852	0.0000	1.0000
JC69	55	-2266.1024	4877.8321	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/1996_hg18.phylip = 1 - 536;
 Partition Dummy = 1:261-aggressive-phylip/1996_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/1996_hg18.phylip = 1 - 536;
 Partition Dummy = 1:261-aggressive-phylip/1996_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/1996_hg18.phylip = 1 - 536;
 Partition Dummy = 1:261-aggressive-phylip/1996_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

