Dear List, dear Edzer, is it correct if I use the term "regression kriging" when I run this function?: kriged_value <- gstat::krige(z ~ x + y, [...])@data$var1.pred
Or should I call it "kriging with external drift" (or "universal kriging" if x and y are coordinates), and use the term "regression kriging" only in the case of running this?: linear_model <- lm(z ~ x + y, [...]) residuals <- linear_model$residuals kriged_residuals <- gstat::krige(residuals ~ 1, [...])@data$var1.pred kriged_value <- linear_model$fitted.values + kriged_residuals Thank you in advance, Ákos Bede-Fazekas Budapest, Hungary _______________________________________________ R-sig-Geo mailing list [email protected] https://stat.ethz.ch/mailman/listinfo/r-sig-geo
