Я пытаюсь подобрать модель VARX, используя пакет mts в R.
Однако при подборе фактической модели я получаю следующее сообщение об ошибке:
«Ошибка в бета-версии [, (jdx + 1) :( jdx + kx)]: индекс за пределами».
В связи с этим я не могу сохранить выходные данные VARX в объекте, которые мне понадобятся для получения IRF экзогенной переменной (xt) по эндогенной переменной (zt).
Ваша помощь высоко ценится!
datar = ts(datar,frequency = 4, start = c(1973, 1))
zt=datar[,1]; xt=datar[,-1]
#Lag Order selection
VARXorder(zt,xt,maxp=4,maxm=8) ## order selection
#model estimation varx
m1=VARX(zt,1,xt,8)
#getting the IRFs
VARXirf(m1, lag=12, orth=FALSE)