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, 2015 at 0:09

1 Answer 1


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 and acknowledge that you have read and understand our privacy policy and code of conduct.

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