R: Прогнозируемый набор данных тестирования вручную вместо использования функции прогноза

Question

У меня есть следующая модель, полученная из функции auto.arima:

BV = -3851.172 * TMAX + 0.9912 * TOURIST + eta

, где

eta = error - 0.6899 * error[t-1]

, где error ~ NID(0, 7.66e+10), где sigma^2 = 7.66e+10.

Я хочу сделать прогноз вручную в наборе данных тестирования вместо использования функции forecast. Вот мой исходный код:

f_test <- eta <- error <- c(); 
eta[1] <- error[1] <- rnorm(1, 0, sqrt(7.66e+10))
for (i in 2:nrow(testing)) {
    error[i] <- rnorm(1, 0, sqrt(7.66e+10))
    eta[i] <- error[i] - 0.6899 * error[i - 1]
    f_test[i] <- -3851.172 * testing[i, "TMAX"] + 0.9912 * testing[i, "Tourist"] + eta[i] - error[i]
}

Это правильно? Я не подхожу близко (учитывая, что ошибка в этом случае генерируется случайным образом) с прогнозом forecast.