I apply Thiessen Polygons interpolation using the krige command in the gstat R package. Using krige
without a variogram model sets automatically krige function to perform IDW, and when nmax parameter is set to 1, it is actually performing thiessen polygons interpolation.
What I have is this:
dec_vor = krige(Dec ~ 1, st, grd, nmax = 1)
Then, when I want to perform cross validation, i do this:
dec_vor_cv = krige.cv(Dec ~ 1, st, nfold = 52, dec_vor)
and I get this error message:
"Error in load.variogram.model(object$model[[name]], c(i - 1, i - 1)) :
model should be of mode variogramModel (use function vgm)"