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

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

Minimum AIC  model: GTR
Minimum AICc model: GTR
Minimum BIC  model: HKY


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
GTR	63	-954.8954	2068.7050	0.3152	0.3152
HKY	59	-961.3998	2069.2333	0.2421	0.5573
GTRG	64	-954.2398	2070.5780	0.1236	0.6809
GTRI	64	-954.4688	2071.0359	0.0983	0.7791
HKYG	60	-960.8054	2071.1269	0.0939	0.8731
HKYI	60	-961.0064	2071.5290	0.0768	0.9499
GTRIG	65	-954.2391	2073.7868	0.0248	0.9747
HKYIG	61	-960.8045	2074.2325	0.0199	0.9946
K2P	56	-970.7420	2078.8173	0.0020	0.9966
K2PG	57	-970.0839	2080.5104	0.0009	0.9975
SYM	60	-965.6077	2080.7315	0.0008	0.9982
K2PI	57	-970.2844	2080.9114	0.0007	0.9989
SYMG	61	-964.8853	2082.3942	0.0003	0.9993
SYMI	61	-965.0788	2082.7811	0.0003	0.9995
K2PIG	58	-970.0828	2083.5416	0.0002	0.9997
F81	58	-970.6140	2084.6040	0.0001	0.9998
SYMIG	62	-964.8848	2085.5256	0.0001	0.9999
F81G	59	-970.0547	2086.5430	0.0000	1.0000
F81I	59	-970.2640	2086.9617	0.0000	1.0000
F81IG	60	-970.0538	2089.6238	0.0000	1.0000
JC69	55	-980.0414	2094.4306	0.0000	1.0000
JC69G	56	-979.3986	2096.1305	0.0000	1.0000
JC69I	56	-979.6070	2096.5474	0.0000	1.0000
JC69IG	57	-979.3976	2099.1379	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTR	63	-954.8954	2035.7908	0.3656	0.3656
GTRG	64	-954.2398	2036.4796	0.2590	0.6246
GTRI	64	-954.4688	2036.9375	0.2060	0.8306
GTRIG	65	-954.2391	2038.4782	0.0954	0.9260
HKY	59	-961.3998	2040.7995	0.0299	0.9559
HKYG	60	-960.8054	2041.6107	0.0199	0.9758
HKYI	60	-961.0064	2042.0129	0.0163	0.9921
HKYIG	61	-960.8045	2043.6090	0.0073	0.9994
SYM	60	-965.6077	2051.2154	0.0002	0.9996
SYMG	61	-964.8853	2051.7707	0.0001	0.9997
SYMI	61	-965.0788	2052.1576	0.0001	0.9998
K2P	56	-970.7420	2053.4840	0.0001	0.9999
SYMIG	62	-964.8848	2053.7695	0.0000	0.9999
K2PG	57	-970.0839	2054.1678	0.0000	0.9999
K2PI	57	-970.2844	2054.5688	0.0000	1.0000
K2PIG	58	-970.0828	2056.1656	0.0000	1.0000
F81	58	-970.6140	2057.2280	0.0000	1.0000
F81G	59	-970.0547	2058.1093	0.0000	1.0000
F81I	59	-970.2640	2058.5279	0.0000	1.0000
F81IG	60	-970.0538	2060.1077	0.0000	1.0000
JC69	55	-980.0414	2070.0828	0.0000	1.0000
JC69G	56	-979.3986	2070.7972	0.0000	1.0000
JC69I	56	-979.6070	2071.2141	0.0000	1.0000
JC69IG	57	-979.3976	2072.7953	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKY	59	-961.3998	2261.0667	0.5619	0.5619
K2P	56	-970.7420	2262.5511	0.2675	0.8293
HKYG	60	-960.8054	2265.6112	0.0579	0.8872
HKYI	60	-961.0064	2266.0134	0.0474	0.9346
K2PG	57	-970.0839	2266.9683	0.0294	0.9640
K2PI	57	-970.2844	2267.3693	0.0240	0.9880
GTR	63	-954.8954	2270.9913	0.0039	0.9920
HKYIG	61	-960.8045	2271.3428	0.0033	0.9953
K2PIG	58	-970.0828	2272.6994	0.0017	0.9969
F81	58	-970.6140	2273.7618	0.0010	0.9979
SYM	60	-965.6077	2275.2158	0.0005	0.9984
GTRG	64	-954.2398	2275.4135	0.0004	0.9988
JC69	55	-980.0414	2275.4166	0.0004	0.9993
GTRI	64	-954.4688	2275.8713	0.0003	0.9996
F81G	59	-970.0547	2278.3764	0.0001	0.9997
F81I	59	-970.2640	2278.7951	0.0001	0.9998
SYMG	61	-964.8853	2279.5045	0.0001	0.9998
JC69G	56	-979.3986	2279.8643	0.0000	0.9999
SYMI	61	-965.0788	2279.8914	0.0000	0.9999
JC69I	56	-979.6070	2280.2812	0.0000	1.0000
GTRIG	65	-954.2391	2281.1454	0.0000	1.0000
F81IG	60	-970.0538	2284.1082	0.0000	1.0000
SYMIG	62	-964.8848	2285.2367	0.0000	1.0000
JC69IG	57	-979.3976	2285.5958	0.0000	1.0000



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

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

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


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

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


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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/5864_hg18.phylip = 1 - 309;
 Partition Dummy = 1:261-aggressive-phylip/5864_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0);
END;


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

