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

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

Minimum AIC  model: GTRG
Minimum AICc model: HKYG
Minimum BIC  model: K2PG


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	-988.1087	2125.9735	0.6920	0.6920
HKYIG	61	-988.1126	2129.0986	0.1451	0.8371
GTRG	64	-983.6456	2129.6714	0.1089	0.9460
K2PG	57	-996.0192	2132.5926	0.0253	0.9713
GTRIG	65	-983.6492	2132.9000	0.0217	0.9930
K2PIG	58	-996.0221	2135.6409	0.0055	0.9985
SYMG	61	-992.8779	2138.6293	0.0012	0.9997
SYMIG	62	-992.8812	2141.7788	0.0003	1.0000
HKY	59	-1019.6076	2185.8792	0.0000	1.0000
HKYI	60	-1019.6072	2188.9705	0.0000	1.0000
GTRI	64	-1014.8582	2192.0965	0.0000	1.0000
GTR	63	-1016.7293	2192.6437	0.0000	1.0000
K2P	56	-1027.6720	2192.8800	0.0000	1.0000
K2PI	57	-1027.6707	2195.8957	0.0000	1.0000
SYM	60	-1024.9845	2199.7252	0.0000	1.0000
SYMI	61	-1024.9834	2202.8403	0.0000	1.0000
F81G	59	-1037.2115	2221.0870	0.0000	1.0000
JC69G	56	-1041.8127	2221.1613	0.0000	1.0000
F81IG	60	-1036.8694	2223.4949	0.0000	1.0000
JC69IG	57	-1041.8164	2224.1870	0.0000	1.0000
F81	58	-1067.8459	2279.2887	0.0000	1.0000
JC69	55	-1072.6981	2279.9381	0.0000	1.0000
JC69I	56	-1072.6977	2282.9313	0.0000	1.0000
F81I	59	-1068.4485	2283.5610	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-983.6456	2095.2913	0.4473	0.4473
HKYG	60	-988.1087	2096.2174	0.2815	0.7288
GTRIG	65	-983.6492	2097.2984	0.1640	0.8928
HKYIG	61	-988.1126	2098.2251	0.1032	0.9960
K2PG	57	-996.0192	2106.0384	0.0021	0.9980
SYMG	61	-992.8779	2107.7558	0.0009	0.9989
K2PIG	58	-996.0221	2108.0441	0.0008	0.9997
SYMIG	62	-992.8812	2109.7624	0.0003	1.0000
HKY	59	-1019.6076	2157.2152	0.0000	1.0000
GTRI	64	-1014.8582	2157.7163	0.0000	1.0000
HKYI	60	-1019.6072	2159.2144	0.0000	1.0000
GTR	63	-1016.7293	2159.4585	0.0000	1.0000
K2P	56	-1027.6720	2167.3440	0.0000	1.0000
K2PI	57	-1027.6707	2169.3415	0.0000	1.0000
SYM	60	-1024.9845	2169.9691	0.0000	1.0000
SYMI	61	-1024.9834	2171.9669	0.0000	1.0000
F81G	59	-1037.2115	2192.4230	0.0000	1.0000
F81IG	60	-1036.8694	2193.7388	0.0000	1.0000
JC69G	56	-1041.8127	2195.6253	0.0000	1.0000
JC69IG	57	-1041.8164	2197.6327	0.0000	1.0000
F81	58	-1067.8459	2251.6919	0.0000	1.0000
F81I	59	-1068.4485	2254.8970	0.0000	1.0000
JC69	55	-1072.6981	2255.3962	0.0000	1.0000
JC69I	56	-1072.6977	2257.3953	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-996.0192	2318.4687	0.6277	0.6277
HKYG	60	-988.1087	2319.8283	0.3181	0.9457
K2PIG	58	-996.0221	2324.2013	0.0357	0.9814
HKYIG	61	-988.1126	2325.5628	0.0181	0.9995
GTRG	64	-983.6456	2333.8095	0.0003	0.9998
SYMG	61	-992.8779	2335.0936	0.0002	1.0000
GTRIG	65	-983.6492	2339.5435	0.0000	1.0000
SYMIG	62	-992.8812	2340.8270	0.0000	1.0000
K2P	56	-1027.6720	2376.0475	0.0000	1.0000
HKY	59	-1019.6076	2377.0992	0.0000	1.0000
K2PI	57	-1027.6707	2381.7718	0.0000	1.0000
HKYI	60	-1019.6072	2382.8253	0.0000	1.0000
SYM	60	-1024.9845	2393.5799	0.0000	1.0000
GTR	63	-1016.7293	2394.2499	0.0000	1.0000
GTRI	64	-1014.8582	2396.2346	0.0000	1.0000
SYMI	61	-1024.9834	2399.3046	0.0000	1.0000
JC69G	56	-1041.8127	2404.3288	0.0000	1.0000
JC69IG	57	-1041.8164	2410.0631	0.0000	1.0000
F81G	59	-1037.2115	2412.3071	0.0000	1.0000
F81IG	60	-1036.8694	2417.3497	0.0000	1.0000
JC69	55	-1072.6981	2460.3728	0.0000	1.0000
JC69I	56	-1072.6977	2466.0988	0.0000	1.0000
F81	58	-1067.8459	2467.8491	0.0000	1.0000
F81I	59	-1068.4485	2474.7811	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/11498_hg18.phylip = 1 - 307;
 Partition Dummy = 1:261-aggressive-phylip/11498_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 (HKYG)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/11498_hg18.phylip = 1 - 307;
 Partition Dummy = 1:261-aggressive-phylip/11498_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 (K2PG)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/11498_hg18.phylip = 1 - 307;
 Partition Dummy = 1:261-aggressive-phylip/11498_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

