I want to perform regression kriging in R using gstat. I have a SpatialPointsDataFrame which is of this form:

enter image description here

After I have run my glm model, in order to see which variables are statistically significant, I end up trying to apply regression kriging with:

k = krige(Jan~prop+Lat+dem, st, grd, fit_jan3)

krige can not see my auxiliary variables, giving me this error message:

Error in eval(expr, envir, enclos) : object 'prop' not found.

Any help please?

  • please replace this picture by a the output of str(st) – agstudy Jul 11 '15 at 0:09

For this to work, prop, Lat and dem need to be present in grd as well. Also, from your question we cannot infer whether dem is present in st.

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.