I am looking to correlate crop area with climate variables and then predict for future if suitability of crop areas will change/remain same under different climate scenarios (in line with species distribution modelling where presence/absence of species is correlated with the environmental variables and hence gives an idea of what climatic conditions do species prefer). My response variable is area under crop (a particular crop) for 600 administrative units and I have some environmental variables as explanatory variables.

I was wondering if someone has experience using biomod2 package in R. biomod2 package allows user to run an ensemble of models (in contrast with running a single model) for species distribution modelling. It takes in species presence (coded as 1) and absence (coded as 0) and correlate it with the environmental variables giving an indication of what environmental variables species like based on where it is found. Since the biomod2 package needs species presence/absence data as response variable, in my case, what do you suggest I should do considering my response variable is a continuous data: Should I develop a threshold crop area value below which it will be considered 0 (i.e. absence of crop area) and above which it will be given the value of 1 (presence of crop area). Or is there any way I can use my crop areas as it is (I haven't tried it yet)

  • 1
    $\begingroup$ There seems to be a statistical question here but it also depends on people understanding exactly what biomod2 does. Also, a name and date reference for almost everybody either means nothing or obliges people to search for a paper and then read it. You may hit the target you want of someone who knows this package and what you are asking, but in general revising the question so it is self-contained is strongly advised. Asking for code-specific advice without an underlying statistical question would be off-topic here. $\endgroup$ – Nick Cox Mar 3 '15 at 13:44
  • $\begingroup$ I have tried editing the question to make it more clear. Hope this explains more clearly I or maybe I will try posting it in stackoverflow. $\endgroup$ – user53020 Mar 3 '15 at 13:55
  • $\begingroup$ I don't think this is more suited to Stack Overflow. $\endgroup$ – Nick Cox Mar 3 '15 at 17:42
  • $\begingroup$ Ok. I have posted it there $\endgroup$ – user53020 Mar 4 '15 at 20:18
  • $\begingroup$ I don't think that.... $\endgroup$ – Nick Cox Mar 4 '15 at 20:22

biomod2 is not handling continuous response variables.

You will have to switch to abundance/count based models (not supported in biomod2) or to convert your continuous data into binary variable if you want to work within biomod2 framework.

Both approaches are interesting and complementary. Comparing the outputs should give you a better idea of your species abundance and distribution across space.

Another interesting idea should be to weight presences by the abundances (Yweight argument) but this is still a under development part of biomod2 and should be used with cautious.


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.