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

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

Minimum AIC  model: SYMG
Minimum AICc model: SYMG
Minimum BIC  model: SYMG


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
SYMG	61	-2066.2128	4273.0561	0.6054	0.6054
SYMIG	62	-2065.8524	4274.9938	0.2298	0.8351
GTRG	64	-2063.5067	4275.6586	0.1648	0.9999
GTRIG	65	-2070.2241	4291.7914	0.0001	1.0000
K2PG	57	-2081.7658	4293.6583	0.0000	1.0000
HKYG	60	-2079.8392	4297.6637	0.0000	1.0000
K2PIG	58	-2092.0397	4316.8129	0.0000	1.0000
HKYIG	61	-2090.3994	4321.4294	0.0000	1.0000
JC69G	56	-2139.2726	4406.0781	0.0000	1.0000
JC69IG	57	-2139.2781	4408.6831	0.0000	1.0000
F81G	59	-2137.5002	4410.3534	0.0000	1.0000
F81IG	60	-2137.5057	4412.9967	0.0000	1.0000
SYMI	61	-2144.8522	4430.3350	0.0000	1.0000
GTRI	64	-2141.2208	4431.0868	0.0000	1.0000
K2PI	57	-2173.8019	4477.7307	0.0000	1.0000
HKYI	60	-2172.6490	4483.2832	0.0000	1.0000
JC69I	56	-2212.3048	4552.1424	0.0000	1.0000
F81I	59	-2210.3787	4556.1104	0.0000	1.0000
SYM	60	-2217.6765	4573.3382	0.0000	1.0000
GTR	63	-2214.8972	4575.7549	0.0000	1.0000
K2P	56	-2238.6896	4604.9120	0.0000	1.0000
HKY	59	-2237.6716	4610.6960	0.0000	1.0000
JC69	55	-2281.8416	4688.6346	0.0000	1.0000
F81	58	-2280.1764	4693.0862	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
SYMG	61	-2066.2128	4254.4255	0.4399	0.4399
GTRG	64	-2063.5067	4255.0134	0.3279	0.7678
SYMIG	62	-2065.8524	4255.7049	0.2320	0.9998
GTRIG	65	-2070.2241	4270.4481	0.0001	1.0000
K2PG	57	-2081.7658	4277.5315	0.0000	1.0000
HKYG	60	-2079.8392	4279.6784	0.0000	1.0000
K2PIG	58	-2092.0397	4300.0794	0.0000	1.0000
HKYIG	61	-2090.3994	4302.7988	0.0000	1.0000
JC69G	56	-2139.2726	4390.5452	0.0000	1.0000
JC69IG	57	-2139.2781	4392.5562	0.0000	1.0000
F81G	59	-2137.5002	4393.0004	0.0000	1.0000
F81IG	60	-2137.5057	4395.0115	0.0000	1.0000
GTRI	64	-2141.2208	4410.4417	0.0000	1.0000
SYMI	61	-2144.8522	4411.7044	0.0000	1.0000
K2PI	57	-2173.8019	4461.6039	0.0000	1.0000
HKYI	60	-2172.6490	4465.2979	0.0000	1.0000
JC69I	56	-2212.3048	4536.6096	0.0000	1.0000
F81I	59	-2210.3787	4538.7574	0.0000	1.0000
SYM	60	-2217.6765	4555.3530	0.0000	1.0000
GTR	63	-2214.8972	4555.7945	0.0000	1.0000
K2P	56	-2238.6896	4589.3791	0.0000	1.0000
HKY	59	-2237.6716	4593.3431	0.0000	1.0000
JC69	55	-2281.8416	4673.6832	0.0000	1.0000
F81	58	-2280.1764	4676.3527	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
SYMG	61	-2066.2128	4507.4821	0.9039	0.9039
SYMIG	62	-2065.8524	4512.9099	0.0599	0.9638
K2PG	57	-2081.7658	4513.9942	0.0348	0.9986
GTRG	64	-2063.5067	4520.5154	0.0013	1.0000
HKYG	60	-2079.8392	4528.5865	0.0000	1.0000
GTRIG	65	-2070.2241	4540.0986	0.0000	1.0000
K2PIG	58	-2092.0397	4540.6906	0.0000	1.0000
HKYIG	61	-2090.3994	4555.8554	0.0000	1.0000
JC69G	56	-2139.2726	4622.8595	0.0000	1.0000
JC69IG	57	-2139.2781	4629.0189	0.0000	1.0000
F81G	59	-2137.5002	4637.7601	0.0000	1.0000
F81IG	60	-2137.5057	4643.9196	0.0000	1.0000
SYMI	61	-2144.8522	4664.7610	0.0000	1.0000
GTRI	64	-2141.2208	4675.9436	0.0000	1.0000
K2PI	57	-2173.8019	4698.0666	0.0000	1.0000
HKYI	60	-2172.6490	4714.2060	0.0000	1.0000
JC69I	56	-2212.3048	4768.9238	0.0000	1.0000
F81I	59	-2210.3787	4783.5171	0.0000	1.0000
SYM	60	-2217.6765	4804.2611	0.0000	1.0000
GTR	63	-2214.8972	4817.1480	0.0000	1.0000
K2P	56	-2238.6896	4821.6933	0.0000	1.0000
HKY	59	-2237.6716	4838.1027	0.0000	1.0000
JC69	55	-2281.8416	4901.8489	0.0000	1.0000
F81	58	-2280.1764	4916.9639	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/11593_hg18.phylip = 1 - 468;
 Partition Dummy = 1:261-aggressive-phylip/11593_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=Fixed(Equal);
END;


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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/11593_hg18.phylip = 1 - 468;
 Partition Dummy = 1:261-aggressive-phylip/11593_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=Fixed(Equal);
END;


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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/11593_hg18.phylip = 1 - 468;
 Partition Dummy = 1:261-aggressive-phylip/11593_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=Fixed(Equal);
END;


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

