Nearest neighbours method

ˆZ(s0)=ki=1(Z(si)wi)

wi=d0i/ki=1d0i=1/k

nn <- gstat(
  formula = vble ~ 1, 
  data = d2, 
  locations = ~ x + y, 
  nmax = 5, 
  set = list(idp = 0)
)