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

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

Minimum AIC  model: GTRI
Minimum AICc model: GTRI
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
GTRI	20	-988.5288	2019.0529	0.3589	0.3589
GTRG	20	-989.0183	2020.0318	0.2200	0.5789
HKYI	16	-993.8318	2020.9436	0.1394	0.7183
GTRIG	21	-988.5289	2021.2577	0.1192	0.8375
HKYG	16	-994.1411	2021.5621	0.1023	0.9398
HKYIG	17	-993.8319	2023.1072	0.0473	0.9871
K2PI	13	-1000.4338	2027.7181	0.0047	0.9918
K2PG	13	-1000.6174	2028.0853	0.0039	0.9957
K2PIG	14	-1000.4338	2029.8511	0.0016	0.9973
SYMI	17	-997.4589	2030.3612	0.0013	0.9986
SYMG	17	-997.7062	2030.8558	0.0010	0.9996
SYMIG	18	-997.4589	2032.5348	0.0004	1.0000
F81I	15	-1006.9471	2045.0209	0.0000	1.0000
F81G	15	-1007.1053	2045.3374	0.0000	1.0000
F81IG	16	-1006.9471	2047.1743	0.0000	1.0000
JC69I	12	-1012.3618	2049.4508	0.0000	1.0000
JC69G	12	-1012.4799	2049.6870	0.0000	1.0000
HKY	15	-1009.9827	2051.0922	0.0000	1.0000
JC69IG	13	-1012.3619	2051.5742	0.0000	1.0000
GTR	19	-1006.0512	2051.9033	0.0000	1.0000
K2P	12	-1015.7479	2056.2230	0.0000	1.0000
SYM	16	-1013.1930	2059.6660	0.0000	1.0000
F81	14	-1022.0697	2073.1230	0.0000	1.0000
JC69	11	-1027.1435	2076.9010	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-988.5288	2017.0577	0.3885	0.3885
GTRG	20	-989.0183	2018.0366	0.2381	0.6266
GTRIG	21	-988.5289	2019.0577	0.1429	0.7696
HKYI	16	-993.8318	2019.6636	0.1056	0.8751
HKYG	16	-994.1411	2020.2821	0.0775	0.9526
HKYIG	17	-993.8319	2021.6638	0.0388	0.9914
K2PI	13	-1000.4338	2026.8677	0.0029	0.9943
K2PG	13	-1000.6174	2027.2349	0.0024	0.9967
K2PIG	14	-1000.4338	2028.8675	0.0011	0.9978
SYMI	17	-997.4589	2028.9178	0.0010	0.9988
SYMG	17	-997.7062	2029.4124	0.0008	0.9996
SYMIG	18	-997.4589	2030.9177	0.0004	1.0000
F81I	15	-1006.9471	2043.8942	0.0000	1.0000
F81G	15	-1007.1053	2044.2106	0.0000	1.0000
F81IG	16	-1006.9471	2045.8943	0.0000	1.0000
JC69I	12	-1012.3618	2048.7236	0.0000	1.0000
JC69G	12	-1012.4799	2048.9597	0.0000	1.0000
HKY	15	-1009.9827	2049.9655	0.0000	1.0000
GTR	19	-1006.0512	2050.1023	0.0000	1.0000
JC69IG	13	-1012.3619	2050.7237	0.0000	1.0000
K2P	12	-1015.7479	2055.4958	0.0000	1.0000
SYM	16	-1013.1930	2058.3860	0.0000	1.0000
F81	14	-1022.0697	2072.1394	0.0000	1.0000
JC69	11	-1027.1435	2076.2871	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PI	13	-1000.4338	2080.0547	0.4946	0.4946
K2PG	13	-1000.6174	2080.4219	0.4117	0.9063
HKYI	16	-993.8318	2085.1245	0.0392	0.9455
HKYG	16	-994.1411	2085.7431	0.0288	0.9743
K2PIG	14	-1000.4338	2086.1459	0.0235	0.9978
HKYIG	17	-993.8319	2091.2160	0.0019	0.9997
JC69I	12	-1012.3618	2097.8193	0.0001	0.9998
JC69G	12	-1012.4799	2098.0554	0.0001	0.9998
SYMI	17	-997.4589	2098.4701	0.0000	0.9999
GTRI	20	-988.5288	2098.8839	0.0000	0.9999
SYMG	17	-997.7062	2098.9646	0.0000	1.0000
GTRG	20	-989.0183	2099.8628	0.0000	1.0000
JC69IG	13	-1012.3619	2103.9107	0.0000	1.0000
SYMIG	18	-997.4589	2104.5613	0.0000	1.0000
K2P	12	-1015.7479	2104.5915	0.0000	1.0000
GTRIG	21	-988.5289	2104.9752	0.0000	1.0000
F81I	15	-1006.9471	2105.2638	0.0000	1.0000
F81G	15	-1007.1053	2105.5802	0.0000	1.0000
HKY	15	-1009.9827	2111.3351	0.0000	1.0000
F81IG	16	-1006.9471	2111.3552	0.0000	1.0000
JC69	11	-1027.1435	2121.2915	0.0000	1.0000
SYM	16	-1013.1930	2123.8470	0.0000	1.0000
GTR	19	-1006.0512	2127.8372	0.0000	1.0000
F81	14	-1022.0697	2129.4178	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr2_12994.phylip = 1 - 442;
 Partition Dummy = 1:591-tomodel-phylip/chr2_12994.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr2_12994.phylip = 1 - 442;
 Partition Dummy = 1:591-tomodel-phylip/chr2_12994.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

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

