Прогнозирование в Dynami c Регрессия

Question

Я пытаюсь понять, могу ли я использовать количество случаев COVID в качестве ковариаты для прогнозирования количества госпитализаций COVID.

Я не понимаю, как писать термин xreg в части прогноза код и выдает ошибку.

Мой код ниже.

library(tidyverse)
library(fpp2)

#obtaining data from the COVID tracking project
US <- read_csv(url("https://covidtracking.com/api/v1/states/daily.csv"))
# filtering out Ohio specific data
OH <- filter(US, state =="OH")
# selecting columns of interest
OH <- dplyr::select(OH, date, state, cases = positive, hosp = hospitalizedCumulative, icu = inIcuCumulative, death)
# coverting to XTS
OH <- xts(OH[,c(-1,-2)], order.by = strptime(OH$date, format = "%Y%m%d"))
# selecting observations with complete values
OH <- first(OH, "-20 days")

# forecasting for the next 14 days (and where things are not looking right)
OH_model <- auto.arima(OH$hosp, xreg = OH$cases) %>% forecast(xreg = rep(mean(OH$cases), 14)) %>% autoplot()

Спасибо