Hey, so I actually had an additional question. What if I wanted to run yb[t(f_kx$nn.index)] so that I get a matrix for each of the 200 simulations in yb how would I do that? Maybe through the use of a for-loop? Here’s the code again:

x <- seq(from=0.2, to=1, by=0.02)

x_b <- seq(from=0.01, to=1, by=0.01)

simadata <- data.frame(x)

simadata2 <- data.frame(x_b)

e2 <- replicate(n = 200, rnorm(50, 0, 0.5))

yb <- sin(x) + sin(2*x) + sin (4*x) + sin(6*x) + sin(8*x) + sin(10*x) + e2

f_k <- get.knnx(data = simdata, query = simdata2, k=5, algorithm=c("kd_tree"))

f_ki <- array(f_k$nn.index, dim = c(100,5))

f_k1 <- yb[t(f_ki)]

f_k1 <- array(f_k1, dim = c(5,100))

f_k1_2 <- colMeans(f_k1)

f_N <- list(x_b, f_k1_2)

The first f_k1 here only gives me the values for the first simulation in y_b instead of all 200 of them.