Проблема в том, что colnames
из data.frame
, который вы генерируете в первом ряду, это macv.ABN
, macv.DAN
и macv.CRL
. Поэтому в model1
нет переменных ABN
и DAN
. Таким образом, рабочий код должен быть:
library(vars)
model1 = VECM(data.frame(ABN=macv$ABN, DAN=macv$DAN,CRL=macv$CRL), lag=10, estim ="ML")
summary(model1)
irf(model1,impulse = "ABN", response = "DAN", boot=FALSE, n.ahead = 10)