Перекодировать для цикла (рекурсии), используя Rcpp

Question

Я должен вычислить цикл for несколько раз по 100 000 с помощью R и хочу оптимизировать его с помощью Rcpp (или альтернатив):

for (t in 1:num) {
      Xt<- c(t(yt[.N:(.N-h)]))      
      yt<- rbindlist(list(yt, as.list(Xt %*% coeff)))
 }

Xt %*% coeff

- это умножение матриц.

После долгих попыток с Rcpp ... Есть ли подсказки, как перекодировать это с помощью Rcpp? Это легко возможно, потому что вычисление Xt не является для меня тривиальным.

Я также реализовал функцию рекурсии R, которая на самом деле не быстрее, чем этот цикл.

Спасибо