Model	Parameter as a function of covariates			Number of parameters	Comparison to model	Log Likelihood	−2 (logLR)	p-value
Model	f	b	g	Number of parameters	Comparison to model	Log Likelihood	−2 (logLR)	p-value
1	int	int	int	5		−7087.924
2	int + pH	int	int	6	1	−7051.742	72.364	<0.0001
3	int	int + pH	int	6	1	−7088.574	0	0.9999
4	int	int	int + pH	6	1	−7087.469	0.91	0.340
5	int + S	int	int	6	1	−7083.645	8.558	0.003
6	int	int + S	int	6	1	−7088.068	0	0.9999
7	int	int	int + S	6	1	−7087.771	0.306	0.580
8	int + TP	int	int	6	1	−7100.969	0	0.9999
9	int	int + TP	int	6	1	−7041.934	91.98	<0.0001
10	int	int	int + TP	6	1	−7084.677	6.494	0.011
11	int + N	int	int	6	1	−7020.683	134.482	<0.0001
12	int	int + N	int	6	1	−7112.818	0	0.9999
13	int	int	int + N	6	1	−7113.944	0	0.9999
14	int + pH + S	int	int	7	2	−7039.702	24.08	<0.0001
15	int + pH + N	int	int	7	2	−6958.700	186.084	<0.0001
16	int + S + N	int	int	7	2	−7014.577	74.33	<0.0001
17	int + pH + N + S	int	int	8	15	−6952.366	12.668	<0.001
18	int + pH + N	int + TP	int	8	15	−6881.009	155.382	<0.0001
19	int + pH + S	int + TP	int	8	15	−6967.294	0	0.9999
20	int + pH + N + S	int + TP	int	9	18	−6875.410	11.198	<0.001