I would like to do inverse distance interpolation with Trend removal. When doing: gstat(formula=z~Co,data=values,newdata=grid,nmax=10) will gstat first remove the trend and then interpolate using inverse distance? Here is the output I got: [ordinary or weighted least squares prediction]