library(vars) vecm.level <- vec2var(vecm, r = 2) arch.test(vecm.level) normality.test(vecm.level) serial.test(vecm.level) predict(vecm.level) irf(vecm.level, boot = FALSE) fevd(vecm.level) class(vecm.level) methods(class = "vec2var")