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

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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
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
HKYG	60	-4041.4011	8209.4933	0.6892	0.6892
GTRG	64	-4038.2513	8212.1356	0.1839	0.8731
HKYIG	61	-4042.5975	8214.1155	0.0683	0.9414
GTRIG	65	-4038.2726	8214.4241	0.0586	1.0000
SYMG	61	-4102.4599	8333.8402	0.0000	1.0000
F81G	59	-4104.8685	8334.2027	0.0000	1.0000
F81IG	60	-4104.8899	8336.4708	0.0000	1.0000
HKY	59	-4106.5435	8337.5528	0.0000	1.0000
SYMIG	62	-4103.2094	8337.5727	0.0000	1.0000
K2PG	57	-4109.2969	8338.6212	0.0000	1.0000
HKYI	60	-4106.5481	8339.7871	0.0000	1.0000
K2PIG	58	-4109.3184	8340.8814	0.0000	1.0000
GTR	63	-4103.7540	8340.8994	0.0000	1.0000
GTRI	64	-4103.7582	8343.1493	0.0000	1.0000
JC69G	56	-4159.4449	8436.7041	0.0000	1.0000
JC69IG	57	-4159.4664	8438.9602	0.0000	1.0000
SYM	60	-4163.1966	8453.0843	0.0000	1.0000
K2PI	57	-4167.4387	8454.9047	0.0000	1.0000
F81	58	-4166.4228	8455.0902	0.0000	1.0000
SYMI	61	-4163.2019	8455.3242	0.0000	1.0000
K2P	56	-4168.9082	8455.6306	0.0000	1.0000
F81I	59	-4166.4280	8457.3217	0.0000	1.0000
JC69	55	-4218.7319	8553.0688	0.0000	1.0000
JC69I	56	-4218.7373	8555.2888	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-4041.4011	8202.8022	0.5909	0.5909
GTRG	64	-4038.2513	8204.5026	0.2525	0.8434
GTRIG	65	-4038.2726	8206.5453	0.0909	0.9343
HKYIG	61	-4042.5975	8207.1951	0.0657	1.0000
SYMG	61	-4102.4599	8326.9198	0.0000	1.0000
F81G	59	-4104.8685	8327.7370	0.0000	1.0000
F81IG	60	-4104.8899	8329.7798	0.0000	1.0000
SYMIG	62	-4103.2094	8330.4189	0.0000	1.0000
HKY	59	-4106.5435	8331.0871	0.0000	1.0000
K2PG	57	-4109.2969	8332.5939	0.0000	1.0000
HKYI	60	-4106.5481	8333.0961	0.0000	1.0000
GTR	63	-4103.7540	8333.5080	0.0000	1.0000
K2PIG	58	-4109.3184	8334.6368	0.0000	1.0000
GTRI	64	-4103.7582	8335.5163	0.0000	1.0000
JC69G	56	-4159.4449	8430.8899	0.0000	1.0000
JC69IG	57	-4159.4664	8432.9329	0.0000	1.0000
SYM	60	-4163.1966	8446.3933	0.0000	1.0000
SYMI	61	-4163.2019	8448.4038	0.0000	1.0000
F81	58	-4166.4228	8448.8457	0.0000	1.0000
K2PI	57	-4167.4387	8448.8774	0.0000	1.0000
K2P	56	-4168.9082	8449.8164	0.0000	1.0000
F81I	59	-4166.4280	8450.8559	0.0000	1.0000
JC69	55	-4218.7319	8547.4637	0.0000	1.0000
JC69I	56	-4218.7373	8549.4746	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-4041.4011	8505.9136	0.9912	0.9912
HKYIG	61	-4042.5975	8515.3583	0.0088	1.0000
GTRG	64	-4038.2513	8527.8213	0.0000	1.0000
GTRIG	65	-4038.2726	8534.9159	0.0000	1.0000
K2PG	57	-4109.2969	8620.5497	0.0000	1.0000
F81G	59	-4104.8685	8625.7964	0.0000	1.0000
K2PIG	58	-4109.3184	8627.6445	0.0000	1.0000
HKY	59	-4106.5435	8629.1465	0.0000	1.0000
F81IG	60	-4104.8899	8632.8911	0.0000	1.0000
SYMG	61	-4102.4599	8635.0830	0.0000	1.0000
HKYI	60	-4106.5481	8636.2074	0.0000	1.0000
SYMIG	62	-4103.2094	8643.6339	0.0000	1.0000
GTR	63	-4103.7540	8651.7749	0.0000	1.0000
GTRI	64	-4103.7582	8658.8351	0.0000	1.0000
JC69G	56	-4159.4449	8713.7938	0.0000	1.0000
JC69IG	57	-4159.4664	8720.8886	0.0000	1.0000
K2P	56	-4168.9082	8732.7204	0.0000	1.0000
K2PI	57	-4167.4387	8736.8332	0.0000	1.0000
F81	58	-4166.4228	8741.8533	0.0000	1.0000
F81I	59	-4166.4280	8748.9154	0.0000	1.0000
SYM	60	-4163.1966	8749.5046	0.0000	1.0000
SYMI	61	-4163.2019	8756.5670	0.0000	1.0000
JC69	55	-4218.7319	8825.3158	0.0000	1.0000
JC69I	56	-4218.7373	8832.3785	0.0000	1.0000



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

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

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


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

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


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

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

