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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: K2PG


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	-988.1360	2018.3158	0.3701	0.3701
GTRI	20	-988.3243	2018.6923	0.3066	0.6767
GTRIG	21	-987.8351	2019.9239	0.1656	0.8423
SYMG	17	-993.5952	2022.6686	0.0420	0.8843
SYMI	17	-993.8138	2023.1058	0.0337	0.9180
HKYG	16	-995.2404	2023.7915	0.0239	0.9420
HKYI	16	-995.3988	2024.1084	0.0204	0.9624
SYMIG	18	-993.3604	2024.3770	0.0179	0.9803
HKYIG	17	-994.9481	2025.3745	0.0109	0.9911
K2PG	13	-1000.2810	2027.4329	0.0039	0.9950
K2PI	13	-1000.4515	2027.7738	0.0033	0.9983
K2PIG	14	-1000.0501	2029.1074	0.0017	1.0000
GTR	19	-998.7313	2037.3073	0.0000	1.0000
SYM	16	-1003.3680	2040.0469	0.0000	1.0000
HKY	15	-1005.9931	2043.1399	0.0000	1.0000
K2P	12	-1010.4721	2045.6888	0.0000	1.0000
F81G	15	-1017.5311	2066.2161	0.0000	1.0000
F81I	15	-1017.6845	2066.5228	0.0000	1.0000
F81IG	16	-1017.3718	2068.0543	0.0000	1.0000
JC69G	12	-1021.9455	2068.6357	0.0000	1.0000
JC69I	12	-1022.0958	2068.9363	0.0000	1.0000
JC69IG	13	-1021.8047	2070.4801	0.0000	1.0000
F81	14	-1027.0192	2083.0456	0.0000	1.0000
JC69	11	-1031.2139	2085.0564	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-988.1360	2016.2720	0.3791	0.3791
GTRI	20	-988.3243	2016.6485	0.3140	0.6931
GTRIG	21	-987.8351	2017.6702	0.1884	0.8814
SYMG	17	-993.5952	2021.1904	0.0324	0.9139
SYMI	17	-993.8138	2021.6276	0.0260	0.9399
HKYG	16	-995.2404	2022.4807	0.0170	0.9569
SYMIG	18	-993.3604	2022.7208	0.0151	0.9720
HKYI	16	-995.3988	2022.7975	0.0145	0.9865
HKYIG	17	-994.9481	2023.8962	0.0084	0.9949
K2PG	13	-1000.2810	2026.5620	0.0022	0.9971
K2PI	13	-1000.4515	2026.9030	0.0019	0.9989
K2PIG	14	-1000.0501	2028.1002	0.0010	1.0000
GTR	19	-998.7313	2035.4627	0.0000	1.0000
SYM	16	-1003.3680	2038.7361	0.0000	1.0000
HKY	15	-1005.9931	2041.9861	0.0000	1.0000
K2P	12	-1010.4721	2044.9442	0.0000	1.0000
F81G	15	-1017.5311	2065.0622	0.0000	1.0000
F81I	15	-1017.6845	2065.3690	0.0000	1.0000
F81IG	16	-1017.3718	2066.7435	0.0000	1.0000
JC69G	12	-1021.9455	2067.8911	0.0000	1.0000
JC69I	12	-1022.0958	2068.1917	0.0000	1.0000
JC69IG	13	-1021.8047	2069.6093	0.0000	1.0000
F81	14	-1027.0192	2082.0384	0.0000	1.0000
JC69	11	-1031.2139	2084.4279	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	13	-1000.2810	2079.4516	0.5139	0.5139
K2PI	13	-1000.4515	2079.7925	0.4334	0.9473
K2PIG	14	-1000.0501	2085.0582	0.0311	0.9784
HKYG	16	-995.2404	2087.5755	0.0088	0.9873
HKYI	16	-995.3988	2087.8923	0.0076	0.9948
SYMG	17	-993.5952	2090.3536	0.0022	0.9970
SYMI	17	-993.8138	2090.7908	0.0018	0.9988
HKYIG	17	-994.9481	2093.0594	0.0006	0.9993
K2P	12	-1010.4721	2093.7653	0.0004	0.9997
SYMIG	18	-993.3604	2095.9525	0.0001	0.9999
GTRG	20	-988.1360	2097.6405	0.0001	0.9999
GTRI	20	-988.3243	2098.0170	0.0000	1.0000
HKY	15	-1005.9931	2103.0125	0.0000	1.0000
GTRIG	21	-987.8351	2103.1072	0.0000	1.0000
SYM	16	-1003.3680	2103.8309	0.0000	1.0000
GTR	19	-998.7313	2112.7628	0.0000	1.0000
JC69G	12	-1021.9455	2116.7122	0.0000	1.0000
JC69I	12	-1022.0958	2117.0128	0.0000	1.0000
JC69IG	13	-1021.8047	2122.4988	0.0000	1.0000
F81G	15	-1017.5311	2126.0886	0.0000	1.0000
F81I	15	-1017.6845	2126.3954	0.0000	1.0000
JC69	11	-1031.2139	2129.1805	0.0000	1.0000
F81IG	16	-1017.3718	2131.8383	0.0000	1.0000
F81	14	-1027.0192	2138.9963	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 591-tomodel-phylip/chr7_4994.phylip = 1 - 432;
 Partition Dummy = 1:591-tomodel-phylip/chr7_4994.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 591-tomodel-phylip/chr7_4994.phylip = 1 - 432;
 Partition Dummy = 1:591-tomodel-phylip/chr7_4994.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 (K2PG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr7_4994.phylip = 1 - 432;
 Partition Dummy = 1:591-tomodel-phylip/chr7_4994.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

