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

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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
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
HKYG	16	-974.8581	1983.3696	0.5608	0.5608
HKYI	16	-975.7783	1985.2101	0.2234	0.7842
HKYIG	17	-974.8221	1985.5100	0.1923	0.9765
GTRG	20	-974.1939	1990.9724	0.0125	0.9890
GTRI	20	-974.8461	1992.2769	0.0065	0.9956
GTRIG	21	-974.1813	1993.2146	0.0041	0.9996
HKY	15	-983.3431	1998.1407	0.0003	1.0000
GTR	19	-982.2434	2004.8181	0.0000	1.0000
F81G	15	-997.3294	2026.1133	0.0000	1.0000
F81I	15	-997.7929	2027.0404	0.0000	1.0000
F81IG	16	-997.3214	2028.2962	0.0000	1.0000
K2PG	13	-1003.9584	2035.0132	0.0000	1.0000
K2PI	13	-1004.7940	2036.6844	0.0000	1.0000
K2PIG	14	-1003.8948	2037.0586	0.0000	1.0000
F81	14	-1005.2078	2039.6845	0.0000	1.0000
SYMI	17	-1004.1722	2044.2102	0.0000	1.0000
SYMG	17	-1004.6710	2045.2079	0.0000	1.0000
SYMIG	18	-1004.6351	2047.3619	0.0000	1.0000
K2P	12	-1013.2185	2051.3739	0.0000	1.0000
SYM	16	-1012.9778	2059.6091	0.0000	1.0000
JC69G	12	-1022.5694	2070.0757	0.0000	1.0000
JC69I	12	-1022.9979	2070.9328	0.0000	1.0000
JC69IG	13	-1022.5543	2072.2049	0.0000	1.0000
JC69	11	-1030.9637	2084.7179	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-974.8581	1981.7161	0.5412	0.5412
HKYI	16	-975.7783	1983.5566	0.2156	0.7568
HKYIG	17	-974.8221	1983.6442	0.2064	0.9632
GTRG	20	-974.1939	1988.3877	0.0193	0.9825
GTRI	20	-974.8461	1989.6923	0.0100	0.9925
GTRIG	21	-974.1813	1990.3627	0.0072	0.9997
HKY	15	-983.3431	1996.6862	0.0003	1.0000
GTR	19	-982.2434	2002.4868	0.0000	1.0000
F81G	15	-997.3294	2024.6588	0.0000	1.0000
F81I	15	-997.7929	2025.5859	0.0000	1.0000
F81IG	16	-997.3214	2026.6427	0.0000	1.0000
K2PG	13	-1003.9584	2033.9169	0.0000	1.0000
K2PI	13	-1004.7940	2035.5880	0.0000	1.0000
K2PIG	14	-1003.8948	2035.7897	0.0000	1.0000
F81	14	-1005.2078	2038.4156	0.0000	1.0000
SYMI	17	-1004.1722	2042.3444	0.0000	1.0000
SYMG	17	-1004.6710	2043.3420	0.0000	1.0000
SYMIG	18	-1004.6351	2045.2702	0.0000	1.0000
K2P	12	-1013.2185	2050.4369	0.0000	1.0000
SYM	16	-1012.9778	2057.9556	0.0000	1.0000
JC69G	12	-1022.5694	2069.1388	0.0000	1.0000
JC69I	12	-1022.9979	2069.9958	0.0000	1.0000
JC69IG	13	-1022.5543	2071.1086	0.0000	1.0000
JC69	11	-1030.9637	2083.9274	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-974.8581	2043.2592	0.6859	0.6859
HKYI	16	-975.7783	2045.0997	0.2733	0.9591
HKYIG	17	-974.8221	2049.0336	0.0382	0.9973
HKY	15	-983.3431	2054.3828	0.0026	1.0000
GTRG	20	-974.1939	2065.3165	0.0000	1.0000
GTRI	20	-974.8461	2066.6210	0.0000	1.0000
GTRIG	21	-974.1813	2071.1379	0.0000	1.0000
GTR	19	-982.2434	2075.5692	0.0000	1.0000
F81G	15	-997.3294	2082.3554	0.0000	1.0000
F81I	15	-997.7929	2083.2825	0.0000	1.0000
K2PG	13	-1003.9584	2083.9206	0.0000	1.0000
K2PI	13	-1004.7940	2085.5917	0.0000	1.0000
F81IG	16	-997.3214	2088.1858	0.0000	1.0000
K2PIG	14	-1003.8948	2089.6398	0.0000	1.0000
F81	14	-1005.2078	2092.2658	0.0000	1.0000
K2P	12	-1013.2185	2096.5942	0.0000	1.0000
SYMI	17	-1004.1722	2107.7338	0.0000	1.0000
SYMG	17	-1004.6710	2108.7315	0.0000	1.0000
SYMIG	18	-1004.6351	2114.5061	0.0000	1.0000
JC69G	12	-1022.5694	2115.2960	0.0000	1.0000
JC69I	12	-1022.9979	2116.1531	0.0000	1.0000
SYM	16	-1012.9778	2119.4986	0.0000	1.0000
JC69IG	13	-1022.5543	2121.1123	0.0000	1.0000
JC69	11	-1030.9637	2126.2383	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/chr1_32461.phylip = 1 - 346;
 Partition Dummy = 1:591-tomodel-phylip/chr1_32461.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/chr1_32461.phylip = 1 - 346;
 Partition Dummy = 1:591-tomodel-phylip/chr1_32461.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 (HKYG)]

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

