协整:使用VECM /ca.jo和p值的不同系数

时间:2019-03-28 20:41:25

标签: r p-value standard-error

我正在尝试获取我的协整向量的p值。我读到许多有关它的问题,大多数答案都依赖于urca软件包中的ca.jo函数(Bernhard Pfaff的书-第156页-是参考)。因此,我决定尝试一下(我短期内使用外生变量,因此tsDyn软件包是我最好的选择,因为用户可以轻松设置它)。

以下是我在VECM(tsDyn)中使用的功能。请注意,我排除了外生变量:

X = VECM(sub_target_1, lag=3, r=1, include="const", estim="ML", LRinclude="const")

这些是我的结果:

#############
###Model VECM 
#############
Full sample size: 431   End sample size: 427
Number of variables: 4  Number of estimated slope parameters 52
AIC -15665.99   BIC -15442.86   SSR 1.146335
Cointegrating vector (estimated by ML):
   Clp_Fx_Nom      Copper       Oil    Ratio     const
r1          1 -0.03725038 0.1276867 1.018937 -8.932703


                    ECT                Clp_Fx_Nom -1       Copper -1          
Equation Clp_Fx_Nom 0.0042(0.0173)     -0.0316(0.0568)     -0.0427(0.0258).   
Equation Copper     0.0249(0.0379)     -0.2829(0.1244)*    -0.1045(0.0565). 

这是我用来比较的ca.jo(urca)函数:

test=ca.jo(sub_target_1, type="trace",K=3,ecdet="const", spec="longrun")

结果如下:

###################### 
# Johansen-Procedure # 
###################### 

Test type: trace statistic , without linear trend and constant in cointegration 

Eigenvalues (lambda):
[1]  9.466481e-02  2.847656e-02  1.536945e-02  9.699295e-03 -2.795648e-18

Values of teststatistic and critical values of test:

          test 10pct  5pct  1pct
r <= 3 |  4.17  7.52  9.24 12.97
r <= 2 | 10.80 17.85 19.96 24.60
r <= 1 | 23.17 32.00 34.91 41.07
r = 0  | 65.73 49.65 53.12 60.16

Eigenvectors, normalised to first column:
(These are the cointegration relations)

              Clp_Fx_Nom.l3   Copper.l3     Oil.l3   Ratio.l3   constant
Clp_Fx_Nom.l3    1.00000000   1.0000000  1.0000000   1.000000   1.000000
Copper.l3       -0.03333813   0.5059200  0.3900929  12.301568   2.931141
Oil.l3           0.11958746  -0.2209898  0.2773398   2.166322  -1.252026
Ratio.l3         1.03753074   1.2162162 -0.2086412 -17.503155 -10.584605
constant        -8.96705858 -11.0296017 -9.3375875 -45.202912   6.269255

如您所见,结果略有不同。我想知道为什么会这样。设置其中一个时,我会丢失某些东西吗?

此外,如何仅使用VECM(tsDyn包)中的输入来查找p值?

0 个答案:

没有答案