Analysed non-linear data with GAM regression, but reviewer has suggested fitting exponential or logarithmic curves instead. Which to use? I am trying to publish my Masters thesis which looked at elephant impacts on vegetation, with a focus on the effect of artificial waterholes. I have made my initial submission and am now performing major revision, which involves some questioning of my statistical approaches.
Some of the elephant impact measures I looked at showed a non-linear relationship with distance to water which led me to use GAMs. However, a reviewer has said that they feel this is statistical overkill and achieves little more than "fitting complicated nonlinear functions that demonstrate that distance matters". They suggest I should rather describe non-linear relationships between elephant impacts and distance to water by fitting simple functions which others can use (e.g. exponential decay curves). I think the point being made is that these simpler functions would allow me to say something more generalisable about the elephant impacts (e.g. they resemble a logarithmic curve), while the GAMs are only useful for describing the impacts I found in my particular study area and is thus of less general usefulness.
I'm still a novice at statistics so I would appreciate some thoughts on which approach seems more sensible. I have included some examples from my data below.
The top graph shows canopy volume/ha in one vegetation type at different distances from water(the circles represent individual sampling plots). The line shows the predicted relationship from the GAM regression. This sort of relationship appears to resemble a logarithmic curve so the reviewer would perhaps recommend I fit that instead in a case like this.


The second graph shows canopy volume/ha in relation to distance to water, but in a second vegetation type. This vegetation type was more heavily impacted by elephants. Again the line on the graph is derived from the GAM prediction, but this could be one where I could look to fit an exponential curve?
I'd appreciate any thoughts on the merits of my use of GAMs versus the reviewer's recommendation to use simpler functions. As an aside, I also have not really explored fitting logarithmic or exponential curves to data in R, so I would also much appreciate any pointers in that regard.

That GAMs are somehow “statistical overkill” and simpler functions are “more generalizable” is a contentious claims. A priori, if you had no idea about the functional relationship between your inputs and outputs, a GAM is a sensible approach which subsumes exponential and logarithmic shapes as a special case. To fit a model and decide post facto that “hey this looks like this other function, I’m going to fit that instead” is not a proper inferential procedure anyway.  Besides, the first plot is clearly non-monotonic so I don’t see why a log would be suggested.
Assuming you fit the models appropriately (looks like these are count or ratio data e.g. something per unit something meaning you should use an offset of some sort) I would probably politely respond to the reviewer to say that GAMs would be more generalizable because they freely estimate the effect as compared to specifying the functional form to be exponential or logarithmic. Specifying the functional form is a higher bias model than estimating the functional form with penalized basis functions. This means however that GAMs are more variable, which may or may not be an issue. You can at least be honest with the uncertainty by bootstrapping.
