Inverse Distance Weighting method (IDW)

res <- grid |> interpolate(idw)

res

## class       : SpatRaster 
## dimensions  : 100, 100, 2  (nrow, ncol, nlyr)
## resolution  : 91.04224, 92.98719  (x, y)
## extent      : 472365.2, 481469.4, 4199888, 4209186  (xmin, xmax, ymin, ymax)
## coord. ref. : GGRS87 / Greek Grid (EPSG:2100) 
## source(s)   : memory
## names       : var1.pred, var1.var 
## min values  :  939.5008,      NaN 
## max values  : 3266.4193,      NaN