tz <- zoo(cbind(Y =TR_full, x = DateMth))
tz[652:815, "Y"]<- NA

rr <- tz[1:652]
rr<-cbind(`lag(Y, -1)` = lag(rr$Y, -1),rr)


fit4 <- dyn$randomForest(Y ~ lag(Y,-1)+x, tz, na.action=na.omit)

fcast4 <- predict(fit4, newdata=rr)