
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
00:53 02/17/2011

Input data from file "183-phylip/chr2_25783.phylip" (29 taxa, 276 characters)

Minimum AIC  model: GTRG
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	-942.4551	2048.9480	0.4566	0.4566
GTRG	64	-940.9972	2049.4258	0.3596	0.8162
HKY	59	-951.0057	2052.7892	0.0669	0.8831
GTRI	64	-942.9940	2053.4193	0.0488	0.9319
GTRIG	65	-941.6515	2054.1601	0.0337	0.9656
HKYG	60	-950.4448	2054.9361	0.0229	0.9884
HKYI	60	-951.6067	2057.2600	0.0072	0.9956
HKYIG	61	-950.4444	2058.2345	0.0044	1.0000
F81	58	-976.3075	2100.1542	0.0000	1.0000
F81G	59	-975.5075	2101.7928	0.0000	1.0000
F81I	59	-976.6169	2104.0116	0.0000	1.0000
F81IG	60	-975.5073	2105.0611	0.0000	1.0000
K2PG	57	-982.8339	2109.9981	0.0000	1.0000
K2P	56	-984.9916	2111.1339	0.0000	1.0000
K2PI	57	-984.0895	2112.5093	0.0000	1.0000
K2PIG	58	-982.8332	2113.2057	0.0000	1.0000
SYMG	61	-980.7333	2118.8125	0.0000	1.0000
SYM	60	-982.5319	2119.1103	0.0000	1.0000
SYMI	61	-981.9553	2121.2564	0.0000	1.0000
SYMIG	62	-980.7331	2122.1422	0.0000	1.0000
JC69G	56	-1001.4327	2144.0161	0.0000	1.0000
JC69	55	-1003.4580	2144.9160	0.0000	1.0000
JC69I	56	-1002.6151	2146.3810	0.0000	1.0000
JC69IG	57	-1001.4323	2147.1949	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-940.9972	2009.9945	0.5067	0.5067
GTR	63	-942.4551	2010.9102	0.3205	0.8272
GTRIG	65	-941.6515	2013.3029	0.0969	0.9241
GTRI	64	-942.9940	2013.9880	0.0688	0.9929
HKY	59	-951.0057	2020.0114	0.0034	0.9963
HKYG	60	-950.4448	2020.8896	0.0022	0.9985
HKYIG	61	-950.4444	2022.8887	0.0008	0.9993
HKYI	60	-951.6067	2023.2134	0.0007	1.0000
F81	58	-976.3075	2068.6150	0.0000	1.0000
F81G	59	-975.5075	2069.0150	0.0000	1.0000
F81IG	60	-975.5073	2071.0146	0.0000	1.0000
F81I	59	-976.6169	2071.2339	0.0000	1.0000
K2PG	57	-982.8339	2079.6678	0.0000	1.0000
K2PIG	58	-982.8332	2081.6665	0.0000	1.0000
K2P	56	-984.9916	2081.9832	0.0000	1.0000
K2PI	57	-984.0895	2082.1790	0.0000	1.0000
SYMG	61	-980.7333	2083.4667	0.0000	1.0000
SYM	60	-982.5319	2085.0638	0.0000	1.0000
SYMIG	62	-980.7331	2085.4662	0.0000	1.0000
SYMI	61	-981.9553	2085.9106	0.0000	1.0000
JC69G	56	-1001.4327	2114.8654	0.0000	1.0000
JC69IG	57	-1001.4323	2116.8647	0.0000	1.0000
JC69	55	-1003.4580	2116.9160	0.0000	1.0000
JC69I	56	-1002.6151	2117.2303	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKY	59	-951.0057	2233.6151	0.8109	0.8109
HKYG	60	-950.4448	2238.1137	0.0855	0.8964
GTR	63	-942.4551	2238.9955	0.0550	0.9515
HKYI	60	-951.6067	2240.4375	0.0268	0.9782
GTRG	64	-940.9972	2241.7001	0.0142	0.9925
HKYIG	61	-950.4444	2243.7332	0.0052	0.9976
GTRI	64	-942.9940	2245.6937	0.0019	0.9996
GTRIG	65	-941.6515	2248.6290	0.0004	1.0000
F81	58	-976.3075	2278.5982	0.0000	1.0000
F81G	59	-975.5075	2282.6186	0.0000	1.0000
K2P	56	-984.9916	2284.7257	0.0000	1.0000
F81I	59	-976.6169	2284.8375	0.0000	1.0000
K2PG	57	-982.8339	2286.0306	0.0000	1.0000
F81IG	60	-975.5073	2288.2386	0.0000	1.0000
K2PI	57	-984.0895	2288.5419	0.0000	1.0000
K2PIG	58	-982.8332	2291.6497	0.0000	1.0000
SYM	60	-982.5319	2302.2878	0.0000	1.0000
SYMG	61	-980.7333	2304.3112	0.0000	1.0000
SYMI	61	-981.9553	2306.7550	0.0000	1.0000
SYMIG	62	-980.7331	2309.9310	0.0000	1.0000
JC69	55	-1003.4580	2316.0381	0.0000	1.0000
JC69G	56	-1001.4327	2317.6079	0.0000	1.0000
JC69I	56	-1002.6151	2319.9727	0.0000	1.0000
JC69IG	57	-1001.4323	2323.2275	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 183-phylip/chr2_25783.phylip = 1 - 276;
 Partition Dummy = 1:183-phylip/chr2_25783.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 (GTR)]

BEGIN MRBAYES;
 Charset 183-phylip/chr2_25783.phylip = 1 - 276;
 Partition Dummy = 1:183-phylip/chr2_25783.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 183-phylip/chr2_25783.phylip = 1 - 276;
 Partition Dummy = 1:183-phylip/chr2_25783.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

