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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: K2PI


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	-957.8546	2083.3283	0.5731	0.5731
HKYI	60	-965.4892	2085.1839	0.2266	0.7997
HKYG	60	-965.9889	2086.1834	0.1375	0.9372
HKYIG	61	-965.6247	2088.7611	0.0379	0.9751
K2PI	57	-973.4046	2091.2793	0.0108	0.9858
K2PG	57	-973.8182	2092.1065	0.0071	0.9930
GTRI	64	-962.9850	2093.5891	0.0034	0.9963
GTRIG	65	-961.9945	2095.0417	0.0016	0.9980
K2PIG	58	-973.8176	2095.3204	0.0014	0.9994
SYMI	61	-970.1469	2097.8056	0.0004	0.9998
SYMG	61	-971.5455	2100.6027	0.0001	0.9999
SYMIG	62	-970.3245	2101.4981	0.0001	1.0000
GTR	63	-970.3424	2104.9029	0.0000	1.0000
HKY	59	-979.6124	2110.1550	0.0000	1.0000
K2P	56	-988.9503	2119.1850	0.0000	1.0000
SYM	60	-985.7002	2125.6060	0.0000	1.0000
F81I	59	-989.6707	2130.2717	0.0000	1.0000
F81G	59	-990.4137	2131.7576	0.0000	1.0000
JC69I	56	-996.1079	2133.5003	0.0000	1.0000
F81IG	60	-989.8964	2133.9984	0.0000	1.0000
JC69G	56	-996.6978	2134.6800	0.0000	1.0000
JC69IG	57	-996.3552	2137.1804	0.0000	1.0000
F81	58	-1004.4370	2156.5593	0.0000	1.0000
JC69	55	-1011.4906	2161.1091	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-957.8546	2043.7093	0.9406	0.9406
HKYI	60	-965.4892	2050.9783	0.0248	0.9654
HKYG	60	-965.9889	2051.9778	0.0151	0.9805
HKYIG	61	-965.6247	2053.2493	0.0080	0.9884
GTRI	64	-962.9850	2053.9701	0.0056	0.9940
GTRIG	65	-961.9945	2053.9891	0.0055	0.9995
K2PI	57	-973.4046	2060.8093	0.0002	0.9997
K2PG	57	-973.8182	2061.6365	0.0001	0.9998
SYMI	61	-970.1469	2062.2938	0.0001	0.9999
K2PIG	58	-973.8176	2063.6352	0.0000	0.9999
SYMIG	62	-970.3245	2064.6491	0.0000	1.0000
SYMG	61	-971.5455	2065.0910	0.0000	1.0000
GTR	63	-970.3424	2066.6849	0.0000	1.0000
HKY	59	-979.6124	2077.2247	0.0000	1.0000
K2P	56	-988.9503	2089.9006	0.0000	1.0000
SYM	60	-985.7002	2091.4004	0.0000	1.0000
F81I	59	-989.6707	2097.3414	0.0000	1.0000
F81G	59	-990.4137	2098.8273	0.0000	1.0000
F81IG	60	-989.8964	2099.7928	0.0000	1.0000
JC69I	56	-996.1079	2104.2159	0.0000	1.0000
JC69G	56	-996.6978	2105.3956	0.0000	1.0000
JC69IG	57	-996.3552	2106.7104	0.0000	1.0000
F81	58	-1004.4370	2124.8741	0.0000	1.0000
JC69	55	-1011.4906	2132.9812	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PI	57	-973.4046	2266.9652	0.3683	0.3683
K2PG	57	-973.8182	2267.7924	0.2436	0.6119
HKYI	60	-965.4892	2267.9846	0.2213	0.8332
HKYG	60	-965.9889	2268.9841	0.1342	0.9674
K2PIG	58	-973.8176	2273.4079	0.0147	0.9821
HKYIG	61	-965.6247	2273.8724	0.0117	0.9937
GTRG	64	-957.8546	2275.1826	0.0061	0.9998
SYMI	61	-970.1469	2282.9169	0.0001	0.9999
GTRI	64	-962.9850	2285.4434	0.0000	1.0000
SYMG	61	-971.5455	2285.7140	0.0000	1.0000
SYMIG	62	-970.3245	2288.8889	0.0000	1.0000
GTRIG	65	-961.9945	2289.0792	0.0000	1.0000
HKY	59	-979.6124	2290.6142	0.0000	1.0000
K2P	56	-988.9503	2292.4397	0.0000	1.0000
GTR	63	-970.3424	2294.5414	0.0000	1.0000
JC69I	56	-996.1079	2306.7551	0.0000	1.0000
JC69G	56	-996.6978	2307.9348	0.0000	1.0000
SYM	60	-985.7002	2308.4066	0.0000	1.0000
F81I	59	-989.6707	2310.7309	0.0000	1.0000
F81G	59	-990.4137	2312.2168	0.0000	1.0000
JC69IG	57	-996.3552	2312.8663	0.0000	1.0000
F81IG	60	-989.8964	2316.7991	0.0000	1.0000
JC69	55	-1011.4906	2331.9036	0.0000	1.0000
F81	58	-1004.4370	2334.6468	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/10555_hg18.phylip = 1 - 275;
 Partition Dummy = 1:261-aggressive-phylip/10555_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/10555_hg18.phylip = 1 - 275;
 Partition Dummy = 1:261-aggressive-phylip/10555_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 (K2PI)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/10555_hg18.phylip = 1 - 275;
 Partition Dummy = 1:261-aggressive-phylip/10555_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

