I want to perform regression kriging in R using gstat. I have a SpatialPointsDataFrame which is of this form:
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?