数据检验后Y,X1,X2,X3,X4都是一阶协整

长期传导线性回归结果：
Y=-2.037657+0.546322X1+0.470159X2-0.405363X3-0.537549X4               
    (0.005325)(0.001907)    (0.002680)  (0.003261)   (0.001730)

如何用R语言作出短期传导机制的修正模型？[类似问题也有别人问过，但是我没找到准确的答案，好像例子都是基于只有X1的，而且解释不清楚，理解上有困难= =]
麻烦顺便把修正模型的形式写一下，我不太确定自己有没有写错...

万分感谢~
[不要告诉我用EVIEWS很容易之类的，同一个数据其他处理都是用R语言，想保持前后一致，就这样~]