), function(x) order(x, decreasing = FALSE))
znn <- data.frame(z)
z <- lapply(data.frame(D), function

. It now forces interpolation with NN.
if (nnas > 0)
nn <- knnidw(1, rmax = .Machine$double.xmax)

<- zu[1:2]
znn <- zu[(l-1):l]
} else {

= .Machine$double.xmax)
sub_grid <- data.frame(X = grid$X[isna],  Y = grid$Y[isna])
znn <- nn(ground, sub_grid)

= "ZNN", damped = FALSE)) %>%
map_dfc(~forecast::forecast(.,h=h)$mean) %>%
mutate(ts = "ETSNTNS", idx = row_number

<- ets(training, model= "ZNN", damped = FALSE)
forecast(fit_ets,h, level=c(95))
}, error=function(e

<- ets(training, model= "ZNN", damped = FALSE)
forecast(fit_ets,h, level=c(95))
}, error=function(e

(training, model= "ZNN", damped = FALSE)
fcast <- forecast(fit_ets,h, level=c(95))
} else if (predictions[i

(training, model= "ZNN", damped = FALSE)
fcast <- forecast(fit_ets,h, level=c(95))
} else if (predictions[i

,"MNN");
modelDo <- "select";
model <- "ZNN

");
modelDo <- "select";
model <- "ZNN";