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

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

Minimum AIC  model: HKYI
Minimum AICc model: HKYI
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
HKYI	16	-1043.0930	2119.3895	0.2680	0.2680
HKYG	16	-1043.3118	2119.8272	0.2153	0.4834
GTRI	20	-1039.1663	2120.2077	0.1780	0.6614
GTRG	20	-1039.4270	2120.7290	0.1372	0.7986
HKYIG	17	-1043.3119	2121.9807	0.0734	0.8719
GTRIG	21	-1039.4270	2122.9210	0.0458	0.9178
K2PI	13	-1048.7003	2124.2005	0.0242	0.9420
K2PG	13	-1048.8934	2124.5868	0.0199	0.9619
SYMI	17	-1045.2157	2125.7884	0.0109	0.9728
SYMG	17	-1045.4648	2126.2865	0.0085	0.9813
K2PIG	14	-1048.8935	2126.7120	0.0069	0.9882
HKY	15	-1048.0187	2127.0971	0.0057	0.9939
SYMIG	18	-1045.4648	2128.4496	0.0029	0.9968
GTR	19	-1044.7185	2129.1296	0.0021	0.9989
K2P	12	-1053.3868	2131.4579	0.0006	0.9995
SYM	16	-1050.2857	2133.7750	0.0002	0.9997
F81I	15	-1051.8708	2134.8011	0.0001	0.9998
F81G	15	-1052.0562	2135.1721	0.0001	0.9999
F81IG	16	-1052.0563	2137.3162	0.0000	1.0000
JC69I	12	-1057.0345	2138.7533	0.0000	1.0000
JC69G	12	-1057.2100	2139.1041	0.0000	1.0000
JC69IG	13	-1057.2100	2141.2201	0.0000	1.0000
F81	14	-1056.5994	2142.1239	0.0000	1.0000
JC69	11	-1061.6398	2145.8573	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	16	-1043.0930	2118.1860	0.2332	0.2332
GTRI	20	-1039.1663	2118.3327	0.2167	0.4498
HKYG	16	-1043.3118	2118.6237	0.1873	0.6372
GTRG	20	-1039.4270	2118.8540	0.1670	0.8041
HKYIG	17	-1043.3119	2120.6237	0.0689	0.8731
GTRIG	21	-1039.4270	2120.8539	0.0614	0.9345
K2PI	13	-1048.7003	2123.4005	0.0172	0.9517
K2PG	13	-1048.8934	2123.7868	0.0142	0.9658
SYMI	17	-1045.2157	2124.4314	0.0103	0.9761
SYMG	17	-1045.4648	2124.9295	0.0080	0.9841
K2PIG	14	-1048.8935	2125.7869	0.0052	0.9893
HKY	15	-1048.0187	2126.0375	0.0046	0.9939
SYMIG	18	-1045.4648	2126.9296	0.0029	0.9969
GTR	19	-1044.7185	2127.4369	0.0023	0.9992
K2P	12	-1053.3868	2130.7736	0.0004	0.9996
SYM	16	-1050.2857	2132.5714	0.0002	0.9998
F81I	15	-1051.8708	2133.7415	0.0001	0.9999
F81G	15	-1052.0562	2134.1125	0.0001	0.9999
F81IG	16	-1052.0563	2136.1127	0.0000	1.0000
JC69I	12	-1057.0345	2138.0691	0.0000	1.0000
JC69G	12	-1057.2100	2138.4199	0.0000	1.0000
JC69IG	13	-1057.2100	2140.4201	0.0000	1.0000
F81	14	-1056.5994	2141.1988	0.0000	1.0000
JC69	11	-1061.6398	2145.2796	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PI	13	-1048.7003	2177.3584	0.4701	0.4701
K2PG	13	-1048.8934	2177.7447	0.3876	0.8577
K2P	12	-1053.3868	2180.5809	0.0939	0.9516
K2PIG	14	-1048.8935	2183.8953	0.0179	0.9695
HKYI	16	-1043.0930	2184.5956	0.0126	0.9821
HKYG	16	-1043.3118	2185.0333	0.0101	0.9922
JC69I	12	-1057.0345	2187.8763	0.0024	0.9946
JC69G	12	-1057.2100	2188.2272	0.0021	0.9967
HKY	15	-1048.0187	2188.2965	0.0020	0.9987
JC69	11	-1061.6398	2190.9363	0.0005	0.9992
HKYIG	17	-1043.3119	2191.1839	0.0005	0.9997
JC69IG	13	-1057.2100	2194.3779	0.0001	0.9998
SYMI	17	-1045.2157	2194.9916	0.0001	0.9998
SYMG	17	-1045.4648	2195.4898	0.0001	0.9999
F81I	15	-1051.8708	2196.0006	0.0000	0.9999
F81G	15	-1052.0562	2196.3715	0.0000	1.0000
SYM	16	-1050.2857	2198.9811	0.0000	1.0000
F81	14	-1056.5994	2199.3072	0.0000	1.0000
GTRI	20	-1039.1663	2201.3448	0.0000	1.0000
SYMIG	18	-1045.4648	2201.6404	0.0000	1.0000
GTRG	20	-1039.4270	2201.8661	0.0000	1.0000
F81IG	16	-1052.0563	2202.5223	0.0000	1.0000
GTR	19	-1044.7185	2206.2984	0.0000	1.0000
GTRIG	21	-1039.4270	2208.0166	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr2_22128.phylip = 1 - 469;
 Partition Dummy = 1:591-tomodel-phylip/chr2_22128.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr2_22128.phylip = 1 - 469;
 Partition Dummy = 1:591-tomodel-phylip/chr2_22128.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

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

