matlab 误差修正模型,求问大神用R做VEC（误差修正）模型的程序包或者代码是什么啊？...

#向量误差修正模型案例分析

> ############################

> #1.生成数据

> set.seed(12345)

> u1

> u2

> u3

> y1

> y2

> y3

> #调用urca包中的ca.jo()对时间序列y1 y2 y3进行Jonhansen协整检验

> #2.Jonhansen协整检验

> library(urca)

> data

## ca.jo(x, type = c("eigen", "trace"), ecdet = c("none", "const", "trend"), K = 2,spec=c("longrun", "transitory"), season = NULL, dumvar = NULL) 注意这里只是用默认设置。 ##

> model.vecm

> head(model.vecm@x) #ca.jo使用S4方法，故用@提取变量

y1          y2         y3

[1,]  0.5855288 -0.31135095 -1.0377854

[2,]  1.2949948  0.59430322 -0.5116634

[3,]  1.1856915  1.28751444 -0.1316301

[4,]  0.7321943  1.64792194  0.7132483

[5,]  1.3380818  0.09367809  1.3288