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

Input data from file "591-tomodel-phylip/chr7_2725.phylip" (7 taxa, 520 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	-1036.0278	2105.1372	0.3147	0.3147
HKYG	16	-1036.3648	2105.8111	0.2247	0.5394
GTRI	20	-1032.4961	2106.6756	0.1458	0.6853
HKYIG	17	-1035.7913	2106.8018	0.1369	0.8222
GTRG	20	-1032.8794	2107.4421	0.0994	0.9216
GTRIG	21	-1032.2352	2108.3258	0.0639	0.9855
K2PI	13	-1043.7632	2114.2458	0.0033	0.9888
F81I	15	-1041.7696	2114.4916	0.0029	0.9918
K2PG	13	-1044.0794	2114.8782	0.0024	0.9942
F81G	15	-1041.9749	2114.9022	0.0024	0.9966
K2PIG	14	-1043.5696	2115.9709	0.0014	0.9980
F81IG	16	-1041.5949	2116.2713	0.0012	0.9992
SYMI	17	-1042.0136	2119.2463	0.0003	0.9994
SYMG	17	-1042.3560	2119.9311	0.0002	0.9996
HKY	15	-1044.8433	2120.6390	0.0001	0.9998
SYMIG	18	-1041.8255	2121.0163	0.0001	0.9999
JC69I	12	-1049.2321	2123.0796	0.0000	0.9999
GTR	19	-1041.8540	2123.2281	0.0000	0.9999
JC69G	12	-1049.4250	2123.4653	0.0000	1.0000
JC69IG	13	-1049.0921	2124.9035	0.0000	1.0000
F81	14	-1050.0180	2128.8677	0.0000	1.0000
K2P	12	-1052.4623	2129.5400	0.0000	1.0000
SYM	16	-1050.8899	2134.8613	0.0000	1.0000
JC69	11	-1057.4148	2137.3493	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	16	-1036.0278	2104.0556	0.2799	0.2799
HKYG	16	-1036.3648	2104.7296	0.1998	0.4797
GTRI	20	-1032.4961	2104.9923	0.1752	0.6549
HKYIG	17	-1035.7913	2105.5827	0.1304	0.7853
GTRG	20	-1032.8794	2105.7587	0.1194	0.9047
GTRIG	21	-1032.2352	2106.4703	0.0837	0.9884
K2PI	13	-1043.7632	2113.5265	0.0025	0.9908
F81I	15	-1041.7696	2113.5392	0.0024	0.9933
F81G	15	-1041.9749	2113.9498	0.0020	0.9953
K2PG	13	-1044.0794	2114.1588	0.0018	0.9971
K2PIG	14	-1043.5696	2115.1392	0.0011	0.9982
F81IG	16	-1041.5949	2115.1898	0.0011	0.9992
SYMI	17	-1042.0136	2118.0272	0.0003	0.9995
SYMG	17	-1042.3560	2118.7120	0.0002	0.9997
SYMIG	18	-1041.8255	2119.6510	0.0001	0.9998
HKY	15	-1044.8433	2119.6867	0.0001	0.9999
GTR	19	-1041.8540	2121.7081	0.0000	0.9999
JC69I	12	-1049.2321	2122.4642	0.0000	1.0000
JC69G	12	-1049.4250	2122.8499	0.0000	1.0000
JC69IG	13	-1049.0921	2124.1841	0.0000	1.0000
F81	14	-1050.0180	2128.0361	0.0000	1.0000
K2P	12	-1052.4623	2128.9246	0.0000	1.0000
SYM	16	-1050.8899	2133.7798	0.0000	1.0000
JC69	11	-1057.4148	2136.8296	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PI	13	-1043.7632	2168.8262	0.4282	0.4282
K2PG	13	-1044.0794	2169.4586	0.3121	0.7403
HKYI	16	-1036.0278	2172.1169	0.0826	0.8229
HKYG	16	-1036.3648	2172.7908	0.0590	0.8819
JC69I	12	-1049.2321	2173.5102	0.0412	0.9231
JC69G	12	-1049.4250	2173.8958	0.0339	0.9570
K2PIG	14	-1043.5696	2174.6928	0.0228	0.9798
F81I	15	-1041.7696	2177.3466	0.0060	0.9859
F81G	15	-1041.9749	2177.7573	0.0049	0.9908
HKYIG	17	-1035.7913	2177.8977	0.0046	0.9954
JC69IG	13	-1049.0921	2179.4839	0.0021	0.9974
K2P	12	-1052.4623	2179.9706	0.0016	0.9991
F81IG	16	-1041.5949	2183.2510	0.0003	0.9994
HKY	15	-1044.8433	2183.4941	0.0003	0.9997
JC69	11	-1057.4148	2183.6217	0.0003	0.9999
F81	14	-1050.0180	2187.5897	0.0000	1.0000
GTRI	20	-1032.4961	2190.0689	0.0000	1.0000
SYMI	17	-1042.0136	2190.3423	0.0000	1.0000
GTRG	20	-1032.8794	2190.8353	0.0000	1.0000
SYMG	17	-1042.3560	2191.0271	0.0000	1.0000
GTRIG	21	-1032.2352	2195.8007	0.0000	1.0000
SYMIG	18	-1041.8255	2196.2199	0.0000	1.0000
SYM	16	-1050.8899	2201.8410	0.0000	1.0000
GTR	19	-1041.8540	2202.5308	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/chr7_2725.phylip = 1 - 520;
 Partition Dummy = 1:591-tomodel-phylip/chr7_2725.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/chr7_2725.phylip = 1 - 520;
 Partition Dummy = 1:591-tomodel-phylip/chr7_2725.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/chr7_2725.phylip = 1 - 520;
 Partition Dummy = 1:591-tomodel-phylip/chr7_2725.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

