
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
21:54 04/04/2011

Input data from file "591-tomodel-phylip/chrUn_random_7197.phylip" (7 taxa, 618 characters)

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
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
HKYG	16	-1072.8884	2178.6820	0.3304	0.3304
HKYI	16	-1072.9402	2178.7856	0.3137	0.6441
HKYIG	17	-1072.7480	2180.5159	0.1321	0.7762
HKY	15	-1074.8648	2180.5269	0.1314	0.9076
GTRG	20	-1071.0691	2183.5451	0.0290	0.9366
GTRI	20	-1071.1340	2183.6750	0.0272	0.9638
GTR	19	-1072.9692	2185.2093	0.0126	0.9765
GTRIG	21	-1070.9398	2185.4299	0.0113	0.9878
F81G	15	-1078.2549	2187.3072	0.0044	0.9922
F81I	15	-1078.3054	2187.4081	0.0042	0.9964
F81IG	16	-1078.1081	2189.1213	0.0018	0.9982
F81	14	-1080.2480	2189.1926	0.0017	0.9999
K2PG	13	-1085.6433	2197.8892	0.0000	1.0000
K2PI	13	-1085.6810	2197.9647	0.0000	1.0000
K2PIG	14	-1085.4906	2199.6778	0.0000	1.0000
K2P	12	-1087.6940	2199.9038	0.0000	1.0000
SYMG	17	-1083.9080	2202.8361	0.0000	1.0000
SYMI	17	-1083.9546	2202.9292	0.0000	1.0000
SYMIG	18	-1083.7603	2204.6624	0.0000	1.0000
SYM	16	-1085.9181	2204.7414	0.0000	1.0000
JC69G	12	-1090.8391	2206.1938	0.0000	1.0000
JC69I	12	-1090.8791	2206.2740	0.0000	1.0000
JC69IG	13	-1090.6816	2207.9659	0.0000	1.0000
JC69	11	-1092.8875	2208.2107	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-1072.8884	2177.7769	0.3229	0.3229
HKYI	16	-1072.9402	2177.8804	0.3066	0.6296
HKYIG	17	-1072.7480	2179.4959	0.1367	0.7663
HKY	15	-1074.8648	2179.7295	0.1216	0.8879
GTRG	20	-1071.0691	2182.1381	0.0365	0.9244
GTRI	20	-1071.1340	2182.2679	0.0342	0.9586
GTRIG	21	-1070.9398	2183.8795	0.0153	0.9738
GTR	19	-1072.9692	2183.9384	0.0148	0.9887
F81G	15	-1078.2549	2186.5099	0.0041	0.9928
F81I	15	-1078.3054	2186.6107	0.0039	0.9967
F81IG	16	-1078.1081	2188.2161	0.0017	0.9984
F81	14	-1080.2480	2188.4961	0.0015	0.9999
K2PG	13	-1085.6433	2197.2865	0.0000	1.0000
K2PI	13	-1085.6810	2197.3621	0.0000	1.0000
K2PIG	14	-1085.4906	2198.9813	0.0000	1.0000
K2P	12	-1087.6940	2199.3881	0.0000	1.0000
SYMG	17	-1083.9080	2201.8161	0.0000	1.0000
SYMI	17	-1083.9546	2201.9092	0.0000	1.0000
SYMIG	18	-1083.7603	2203.5205	0.0000	1.0000
SYM	16	-1085.9181	2203.8362	0.0000	1.0000
JC69G	12	-1090.8391	2205.6781	0.0000	1.0000
JC69I	12	-1090.8791	2205.7583	0.0000	1.0000
JC69IG	13	-1090.6816	2207.3633	0.0000	1.0000
JC69	11	-1092.8875	2207.7751	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKY	15	-1074.8648	2246.1269	0.5443	0.5443
HKYG	16	-1072.8884	2248.6007	0.1580	0.7023
HKYI	16	-1072.9402	2248.7042	0.1500	0.8523
F81	14	-1080.2480	2250.4669	0.0621	0.9145
K2P	12	-1087.6940	2252.5059	0.0224	0.9369
F81G	15	-1078.2549	2252.9072	0.0183	0.9552
F81I	15	-1078.3054	2253.0080	0.0174	0.9727
HKYIG	17	-1072.7480	2254.7462	0.0073	0.9800
K2PG	13	-1085.6433	2254.8309	0.0070	0.9870
K2PI	13	-1085.6810	2254.9064	0.0068	0.9938
JC69	11	-1092.8875	2256.4665	0.0031	0.9969
JC69G	12	-1090.8391	2258.7960	0.0010	0.9978
JC69I	12	-1090.8791	2258.8762	0.0009	0.9988
F81IG	16	-1078.1081	2259.0399	0.0009	0.9996
K2PIG	14	-1085.4906	2260.9521	0.0003	0.9999
JC69IG	13	-1090.6816	2264.9076	0.0000	1.0000
GTR	19	-1072.9692	2268.0417	0.0000	1.0000
GTRG	20	-1071.0691	2270.6679	0.0000	1.0000
GTRI	20	-1071.1340	2270.7977	0.0000	1.0000
SYM	16	-1085.9181	2274.6601	0.0000	1.0000
GTRIG	21	-1070.9398	2276.8358	0.0000	1.0000
SYMG	17	-1083.9080	2277.0664	0.0000	1.0000
SYMI	17	-1083.9546	2277.1595	0.0000	1.0000
SYMIG	18	-1083.7603	2283.1973	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chrUn_random_7197.phylip = 1 - 618;
 Partition Dummy = 1:591-tomodel-phylip/chrUn_random_7197.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 AICc model (HKYG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chrUn_random_7197.phylip = 1 - 618;
 Partition Dummy = 1:591-tomodel-phylip/chrUn_random_7197.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 (HKY)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chrUn_random_7197.phylip = 1 - 618;
 Partition Dummy = 1:591-tomodel-phylip/chrUn_random_7197.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

