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...