1
$\begingroup$

Let’s say I have the hypothesis that soil fertility affects the relationship between weed biomass and crop biomass. One way to go about testing that hypothesis might be to model the relationship between weed and crop biomass at various levels of soil fertility and then deciding whether to group or separate these data (for example, through an F-test). Example 1

This seems less than ideal though because it treats soil fertility, which is really continuous and numerical, like it's categorical variable. Also, what if the data looks like this? Example 2

That is, instead of having data at three distinct levels of soil fertility, I have observations at a wide range of values. What can I do to ask the question, “does soil fertility affect the relationship between weed and crop biomass?” Or specifically, since this relationship is modelled by a rectangular hyperbola of the equation...

Yield Loss = I * Weed Biomass / (1 + (I/A)* Weed Biomass )

How do I ask "does soil fertility affect the slope (I) or asymptote parameter (A) of this relationship?"

I feel that I am overlooking something simple. If this were a linear regression, it would be straightforward to see if soil fertility affects the slope or intercept of the relationship, but I am having trouble understanding how I would do something similar in a nonlinear regression situation.

Many thanks!

Angela

$\endgroup$
2
  • $\begingroup$ Is there some reason why you need to model this in terms of percent crop yield loss instead of as crop yield itself? $\endgroup$
    – EdM
    Jun 17, 2015 at 18:10
  • $\begingroup$ I should say that this is not actually my specific research question. It was just the simplest example that I could think of (one that was fairly intuitive and wouldn’t require too much background information). That being said, I chose to model in terms of percent yield loss for a few reasons of convenience, but no reasons of necessity. For example, percent yield loss is nice because it lets me combine/compare data across different sites and years (where maximum crop yield might be very different). $\endgroup$
    – Angela
    Jun 17, 2015 at 18:37

1 Answer 1

1
$\begingroup$

The general approach to evaluating the addition of a variable to a model isn't that much different for linear and non-linear models. In both cases you examine whether the reduction of residual error when you add the variable more than balances off the associated loss of a degree of freedom. This page, from the guide to a commercial software package, describes pretty simply how to use an F-test or the Akaike Information Criterion to make this assessment in non-linear regression.

Your bigger problem may be how to incorporate a variable like soil fertility into the model to begin with. For the example you show, you would have to specify a particular functional relation between that variable and the intercepts and asymptotes, and the validity of your results would depend on the choice of that relation. (Of course, linear regression also makes such a choice.)

Your choice of example raises two additional concerns, related to the use of percent crop yield loss as the dependent variable. I understand that this example is not your specific research question, but your choice of this example may indicate some additional issues you need to consider as you proceed with your own research.

First, in linear regression, a dependent variable with a necessarily restricted range like this might be better examined with a generalized linear model that transforms the dependent variable to cover the range of ($-\infty, \infty)$ before the regression. I am not sure how this is handled in non-linear regression. If you are studying agronomy there presumably is competent statistical help at your institution that can point you in the correct direction.

Second, and maybe more important, your choice of this dependent variable throws away all the information you have about maximum crop yield, and does not explicitly include the inevitable errors in measuring maximum crop yield. Throwing away data before analysis is almost never a good idea. In my opinion this example would be much better analyzed with actual crop yield as the dependent variable, and maximum yield included as a predictor variable along with weed biomass and fertility; it might even then reduce to a linear model. There presumably are strong relations among soil fertility, maximum crop yield, and weed biomass that would be more readily explored in this way.

$\endgroup$
1
  • $\begingroup$ Great—thank you! Your explanation helped me to identify the crux of my question as well as answer it. I was (1) struggling with where even to put soil fertility into the model and (2) getting hung up on the idea that I would determine the usefulness of soil fertility as a variable based solely on the significance of the parameter estimate coupled with it. As for the concerns you brought up about the choice in dependent variable, I’ll have to think on that more in depth with regard to my particular situation. Thanks for bringing them up! $\endgroup$
    – Angela
    Jun 18, 2015 at 16:27

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.