
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
16:18 04/04/2011

Input data from file "683-tomodel-phylip/chr8_9141.phylip" (7 taxa, 437 characters)

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


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	20	-1022.7997	2087.6186	0.2109	0.2109
HKYG	16	-1027.1958	2087.6869	0.2038	0.4148
HKYI	16	-1027.3853	2088.0657	0.1687	0.5834
GTRI	20	-1023.1509	2088.3210	0.1485	0.7319
GTRIG	21	-1022.1439	2088.5143	0.1348	0.8666
HKYIG	17	-1026.5375	2088.5356	0.1333	1.0000
K2PG	13	-1041.8206	2110.5017	0.0000	1.0000
K2PI	13	-1041.9481	2110.7567	0.0000	1.0000
K2PIG	14	-1041.1152	2111.2257	0.0000	1.0000
GTR	19	-1037.2576	2114.3378	0.0000	1.0000
HKY	15	-1041.7389	2114.6179	0.0000	1.0000
SYMG	17	-1040.1164	2115.6934	0.0000	1.0000
SYMI	17	-1040.2756	2116.0119	0.0000	1.0000
SYMIG	18	-1039.6049	2116.8461	0.0000	1.0000
F81G	15	-1044.6013	2120.3428	0.0000	1.0000
F81I	15	-1044.8084	2120.7570	0.0000	1.0000
F81IG	16	-1044.0213	2121.3378	0.0000	1.0000
K2P	12	-1056.8739	2138.4837	0.0000	1.0000
JC69G	12	-1057.8754	2140.4865	0.0000	1.0000
JC69I	12	-1058.0466	2140.8290	0.0000	1.0000
JC69IG	13	-1057.3050	2141.4705	0.0000	1.0000
SYM	16	-1054.5025	2142.3002	0.0000	1.0000
F81	14	-1058.1023	2145.1999	0.0000	1.0000
JC69	11	-1071.6706	2165.9624	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-1022.7997	2085.5994	0.2427	0.2427
GTRIG	21	-1022.1439	2086.2878	0.1720	0.4147
GTRI	20	-1023.1509	2086.3018	0.1708	0.5855
HKYG	16	-1027.1958	2086.3917	0.1633	0.7488
HKYI	16	-1027.3853	2086.7705	0.1351	0.8840
HKYIG	17	-1026.5375	2087.0750	0.1160	1.0000
K2PG	13	-1041.8206	2109.6412	0.0000	1.0000
K2PI	13	-1041.9481	2109.8962	0.0000	1.0000
K2PIG	14	-1041.1152	2110.2304	0.0000	1.0000
GTR	19	-1037.2576	2112.5153	0.0000	1.0000
HKY	15	-1041.7389	2113.4777	0.0000	1.0000
SYMG	17	-1040.1164	2114.2327	0.0000	1.0000
SYMI	17	-1040.2756	2114.5513	0.0000	1.0000
SYMIG	18	-1039.6049	2115.2098	0.0000	1.0000
F81G	15	-1044.6013	2119.2027	0.0000	1.0000
F81I	15	-1044.8084	2119.6169	0.0000	1.0000
F81IG	16	-1044.0213	2120.0426	0.0000	1.0000
K2P	12	-1056.8739	2137.7478	0.0000	1.0000
JC69G	12	-1057.8754	2139.7507	0.0000	1.0000
JC69I	12	-1058.0466	2140.0931	0.0000	1.0000
JC69IG	13	-1057.3050	2140.6100	0.0000	1.0000
SYM	16	-1054.5025	2141.0050	0.0000	1.0000
F81	14	-1058.1023	2144.2047	0.0000	1.0000
JC69	11	-1071.6706	2165.3412	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1027.1958	2151.6706	0.5185	0.5185
HKYI	16	-1027.3853	2152.0494	0.4290	0.9475
HKYIG	17	-1026.5375	2156.4339	0.0479	0.9954
K2PG	13	-1041.8206	2162.6804	0.0021	0.9975
K2PI	13	-1041.9481	2162.9353	0.0019	0.9994
GTRG	20	-1022.7997	2167.1981	0.0002	0.9996
K2PIG	14	-1041.1152	2167.3495	0.0002	0.9998
GTRI	20	-1023.1509	2167.9004	0.0002	1.0000
GTRIG	21	-1022.1439	2171.9664	0.0000	1.0000
HKY	15	-1041.7389	2174.6767	0.0000	1.0000
F81G	15	-1044.6013	2180.4017	0.0000	1.0000
F81I	15	-1044.8084	2180.8159	0.0000	1.0000
SYMG	17	-1040.1164	2183.5916	0.0000	1.0000
SYMI	17	-1040.2756	2183.9101	0.0000	1.0000
F81IG	16	-1044.0213	2185.3215	0.0000	1.0000
K2P	12	-1056.8739	2186.7070	0.0000	1.0000
SYMIG	18	-1039.6049	2188.6486	0.0000	1.0000
JC69G	12	-1057.8754	2188.7099	0.0000	1.0000
JC69I	12	-1058.0466	2189.0523	0.0000	1.0000
GTR	19	-1037.2576	2190.0340	0.0000	1.0000
JC69IG	13	-1057.3050	2193.6491	0.0000	1.0000
F81	14	-1058.1023	2201.3237	0.0000	1.0000
SYM	16	-1054.5025	2206.2839	0.0000	1.0000
JC69	11	-1071.6706	2210.2205	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 683-tomodel-phylip/chr8_9141.phylip = 1 - 437;
 Partition Dummy = 1:683-tomodel-phylip/chr8_9141.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 683-tomodel-phylip/chr8_9141.phylip = 1 - 437;
 Partition Dummy = 1:683-tomodel-phylip/chr8_9141.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 (HKYG)]

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr8_9141.phylip = 1 - 437;
 Partition Dummy = 1:683-tomodel-phylip/chr8_9141.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;


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

