认真看了,他用的 是auto.arima模型,不是我说的TSA里面的arimax包。。。
然后用的forecast函数进行预测,这里没问题的。不过auto.arima有一个问题就是外生变量x只能有一个
fit_vars_0 <- auto.arima(df$cons, xreg = vars_matrix[, 1:2]) #xreg只有一列向量
expected_temp_income <- matrix(c(fcast_temp, 91, 91, 93, 96, 96, 96), ncol = 2, nrow = 6)
fcast_cons_temp_income <- forecast(fit_vars_0, xreg = expected_temp_income,h = 6)