The derivation of species response curves with Gaussian logistic regression is sensitive to sampling intensity and curve characteristics

We investigated quantitatively the sensitivity of plant species response curves to sampling characteristics (number of plots, occurrence and frequency of species), along a simulated pH gradient. We defined 54 theoretical unimodal response curves, issued from combinations of six values for optimum (opt = 3, 4, …, 8), three values for tolerance (tol = 0.5, 1.0, and 1.5, sensu ter Braak and Looman [ter Braak, C.J.F., Looman, C.W.N., 1986. Weighted averaging, logistic regression and the Gaussian response model. Vegetatio 65, 3–11]), and three values for maximum probability of presence (pmax = 0.05, 0.20, and 0.50). For each of these 54 theoretical response curves, we built artificial binary data sets (presence/absence) to test the influence of species occurrence, frequency, or number of available plots. With real data extracted from EcoPlant, a phytoecological database for French forests [Gégout, J.-C., Coudun, Ch., Bailly, G., Jabiol, B., 2005. EcoPlant: a forest sites database linking floristic data with soil characteristics and climatic conditions. J. Veg. Sci. 16, 257–260], we compared the ecological response of 50 plant species to soil pH, based first on a small data set (100 randomly sampled plots), and then based on the whole data set available (3810 plots). s on simulated data showed that the curve optimum, amplitude, or maximum probability of presence, cannot be assessed reliably with logistic regression, when a species is too rare, or when its theoretical optimum lies near an extreme of the gradient. Those theoretical results were illustrated by real data extracted from EcoPlant. We suggest a general minimum value of 50 occurrences for species to derive acceptable ecological response curves with logistic regression.