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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
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
GTRG	64	-1931.0110	4008.3078	0.4243	0.4243
HKYG	60	-1936.3419	4008.6316	0.3609	0.7853
GTRIG	65	-1931.0135	4010.9256	0.1146	0.8999
HKYIG	61	-1936.3444	4011.2042	0.0997	0.9996
K2PG	57	-1947.4665	4023.2447	0.0002	0.9999
K2PIG	58	-1947.4689	4025.7838	0.0001	0.9999
SYMG	61	-1943.6924	4025.9000	0.0001	1.0000
SYMIG	62	-1943.6948	4028.4838	0.0000	1.0000
GTRI	64	-1953.4708	4053.2273	0.0000	1.0000
HKYI	60	-1959.9799	4055.9075	0.0000	1.0000
K2PI	57	-1971.0678	4070.4473	0.0000	1.0000
SYMI	61	-1967.2968	4073.1088	0.0000	1.0000
GTR	63	-1969.5122	4082.7087	0.0000	1.0000
HKY	59	-1976.0463	4085.4839	0.0000	1.0000
F81G	59	-1976.3002	4085.9916	0.0000	1.0000
F81IG	60	-1976.3157	4088.5792	0.0000	1.0000
JC69G	56	-1984.4927	4094.7737	0.0000	1.0000
JC69IG	57	-1984.5074	4097.3264	0.0000	1.0000
K2P	56	-1987.4231	4100.6346	0.0000	1.0000
SYM	60	-1983.4349	4102.8174	0.0000	1.0000
F81I	59	-1998.0713	4129.5340	0.0000	1.0000
JC69I	56	-2006.3441	4138.4765	0.0000	1.0000
F81	58	-2013.2005	4157.2469	0.0000	1.0000
JC69	55	-2021.7961	4166.8680	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1931.0110	3990.0221	0.5786	0.5786
GTRIG	65	-1931.0135	3992.0270	0.2123	0.7909
HKYG	60	-1936.3419	3992.6839	0.1529	0.9438
HKYIG	61	-1936.3444	3994.6889	0.0561	0.9999
K2PG	57	-1947.4665	4008.9331	0.0000	0.9999
SYMG	61	-1943.6924	4009.3847	0.0000	1.0000
K2PIG	58	-1947.4689	4010.9378	0.0000	1.0000
SYMIG	62	-1943.6948	4011.3897	0.0000	1.0000
GTRI	64	-1953.4708	4034.9416	0.0000	1.0000
HKYI	60	-1959.9799	4039.9598	0.0000	1.0000
K2PI	57	-1971.0678	4056.1356	0.0000	1.0000
SYMI	61	-1967.2968	4056.5936	0.0000	1.0000
GTR	63	-1969.5122	4065.0245	0.0000	1.0000
HKY	59	-1976.0463	4070.0926	0.0000	1.0000
F81G	59	-1976.3002	4070.6003	0.0000	1.0000
F81IG	60	-1976.3157	4072.6314	0.0000	1.0000
JC69G	56	-1984.4927	4080.9854	0.0000	1.0000
JC69IG	57	-1984.5074	4083.0147	0.0000	1.0000
K2P	56	-1987.4231	4086.8463	0.0000	1.0000
SYM	60	-1983.4349	4086.8697	0.0000	1.0000
F81I	59	-1998.0713	4114.1427	0.0000	1.0000
JC69I	56	-2006.3441	4124.6882	0.0000	1.0000
F81	58	-2013.2005	4142.4010	0.0000	1.0000
JC69	55	-2021.7961	4153.5922	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1936.3419	4247.9136	0.8150	0.8150
K2PG	57	-1947.4665	4251.4013	0.1425	0.9574
HKYIG	61	-1936.3444	4254.1724	0.0356	0.9931
K2PIG	58	-1947.4689	4257.6599	0.0062	0.9993
GTRG	64	-1931.0110	4262.2671	0.0006	0.9999
GTRIG	65	-1931.0135	4268.5258	0.0000	1.0000
SYMG	61	-1943.6924	4268.8683	0.0000	1.0000
SYMIG	62	-1943.6948	4275.1270	0.0000	1.0000
HKYI	60	-1959.9799	4295.1895	0.0000	1.0000
K2PI	57	-1971.0678	4298.6038	0.0000	1.0000
GTRI	64	-1953.4708	4307.1866	0.0000	1.0000
SYMI	61	-1967.2968	4316.0771	0.0000	1.0000
JC69G	56	-1984.4927	4319.1998	0.0000	1.0000
HKY	59	-1976.0463	4321.0685	0.0000	1.0000
F81G	59	-1976.3002	4321.5762	0.0000	1.0000
K2P	56	-1987.4231	4325.0607	0.0000	1.0000
JC69IG	57	-1984.5074	4325.4830	0.0000	1.0000
F81IG	60	-1976.3157	4327.8612	0.0000	1.0000
GTR	63	-1969.5122	4333.0157	0.0000	1.0000
SYM	60	-1983.4349	4342.0994	0.0000	1.0000
JC69I	56	-2006.3441	4362.9026	0.0000	1.0000
F81I	59	-1998.0713	4365.1186	0.0000	1.0000
JC69	55	-2021.7961	4387.5528	0.0000	1.0000
F81	58	-2013.2005	4389.1230	0.0000	1.0000



To specify default values for best AIC, AICc, or BIC models in MrBayes:

[Mrbayes block for the best AIC model (GTRG)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/11371_hg18.phylip = 1 - 520;
 Partition Dummy = 1:261-aggressive-phylip/11371_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 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);
END;


[Mrbayes block for the best AICc model (GTRG)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/11371_hg18.phylip = 1 - 520;
 Partition Dummy = 1:261-aggressive-phylip/11371_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 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);
END;


[Mrbayes block for the best BIC model (HKYG)]

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

