# create noisy cylcic series
set.seed(2020)
x <- rnorm(1000)
x <- runif(1000) * sin(x) + rnorm(1000)
w <- x
# introduce 10 percent missingness
x[sample(1000, 100)] <- NA
# make predictions with lag 3 using k=4 nearest neighbors
towerTS(x, lag=10, k=4)
# -0.1685019
Run the code above in your browser using DataLab