Возможно, я неправильно понял ваш вопрос, но ваш код уже включает соответствующий прогноз.
ARDL_Model <- dynamac::dynardl(Y ~ X1 + X2 + X3 + X4 , data = My_Data,
diffs = c("X1","X2", "X3","X4"),
lagdiffs = list("X1" = c(1:5), "X2" = 1, "X3" = c(1:2), "X4" = c(1:2)),
ec = TRUE, simulate = TRUE, shockvar = "X2", graph= TRUE)
Создает набор смоделированных результатов в ARDL_Model$simulate
.ARDL_Model$simulate$central
является ответом зависимой переменной.При shocktime
(поскольку здесь он не указан, по умолчанию t = 10), X2
будет шокировано стандартным отклонением его значения (поскольку shockval
не указано).Другие переменные будут храниться в своих средствах.Таким образом, вы сможете прогнозировать ответ зависимой переменной, но вам нужно будет запустить отдельное моделирование (т. Е. Для каждого shockvar
и shockval
, которое вы хотите спрогнозировать, ответ).
Виньетка в https://cran.r -project.org / web / packages / dynamac / vignettes / dynamac-vignette.html также может быть полезной.