この特定のRコマンドの作成者は、元のDickey-Fullerの式に精通していると想定しているため、値の解釈方法に関する関連ドキュメントを提供していません。私はエンダースが非常に役立つリソースであることを発見しました(Applied Econometric Time Series 3e、2010、p。206-209-他のエディションも問題ないと思います)。以下では、例としてデンマークの実質所得であるURCAパッケージのデータを使用します。
> income <- ts(denmark$LRY)
Dickey-Fullerがur.dfの「タイプ」オプションと一致するため、異なる仮説を取得するために使用される3つの異なる数式を最初に説明すると役立つ場合があります。エンダースは、これら3つのケースのすべてで、一貫した用語として、ガンマ、以前のyの値の係数、ラグ項を使用することを指定しています。gamma = 0の場合、単位根があります(ランダムウォーク、非定常)。帰無仮説がgamma = 0の場合、p <0.05の場合、帰無(95%レベル)を棄却し、単位根がないと仮定します。nullを拒否できなかった場合(p> 0.05)、単位根が存在すると推定します。ここから、タウとファイの解釈に進むことができます。
Δy(t)=γ∗y(t−1)+e(t)
e(t)γ=a−1y=a∗y(t−1)+e(t)y(t−1)
type = "none"の場合、tau(またはR出力のtau1)は、ガンマ= 0の帰無仮説です。デンマークの所得の例を使用すると、 "Value of test-statistic is 0.7944"となり、 "Test statisticsの重要な値は:tau1 -2.6 -1.95 -1.61。テスト統計がnullを拒否できない3つの領域すべて(1%、5%、10%)内にある場合、データはランダムウォークであると想定する必要があります。単位根が存在します。この場合、tau1はガンマ= 0の仮説を参照します。「z.lag1」はガンマ項、ラグ項の係数(y(t-1))です。これはp = 0.431は有意であるとして拒否できず、単にガンマがこのモデルに対して統計的に有意ではないことを意味します。Rからの出力は次のとおりです。
> summary(ur.df(y=income, type = "none",lags=1))
>
> ###############################################
> # Augmented Dickey-Fuller Test Unit Root Test #
> ###############################################
>
> Test regression none
>
>
> Call:
> lm(formula = z.diff ~ z.lag.1 - 1 + z.diff.lag)
>
> Residuals:
> Min 1Q Median 3Q Max
> -0.044067 -0.016747 -0.006596 0.010305 0.085688
>
> Coefficients:
> Estimate Std. Error t value Pr(>|t|)
> z.lag.1 0.0004636 0.0005836 0.794 0.431
> z.diff.lag 0.1724315 0.1362615 1.265 0.211
>
> Residual standard error: 0.0251 on 51 degrees of freedom
> Multiple R-squared: 0.04696, Adjusted R-squared: 0.009589
> F-statistic: 1.257 on 2 and 51 DF, p-value: 0.2933
>
>
> Value of test-statistic is: 0.7944
>
> Critical values for test statistics:
> 1pct 5pct 10pct
> tau1 -2.6 -1.95 -1.61
Δy(t)=a0+γ∗y(t−1)+e(t)
γ=0γ=0
phi1項は、2番目の仮説を指します。これは、a0 =ガンマ= 0の帰無仮説を組み合わせたものです。これは、値の両方が同時に0であることがテストされることを意味します。p <0.05の場合、nullを拒否し、これらの少なくとも1つが偽であると仮定します。つまり、a0またはgammaのいずれかまたは両方が0ではありません。このnullを拒否しないと、a0とgamma = 0の両方がつまり、1)ガンマ= 0で単位根が存在し、かつ2)a0 = 0であるため、ドリフト項はありません。こちらがR出力
> summary(ur.df(y=income, type = "drift",lags=1))
>
> ###############################################
> # Augmented Dickey-Fuller Test Unit Root Test #
> ###############################################
>
> Test regression drift
>
>
> Call:
> lm(formula = z.diff ~ z.lag.1 + 1 + z.diff.lag)
>
> Residuals:
> Min 1Q Median 3Q Max
> -0.041910 -0.016484 -0.006994 0.013651 0.074920
>
> Coefficients:
> Estimate Std. Error t value Pr(>|t|)
> (Intercept) 0.43453 0.28995 1.499 0.140
> z.lag.1 -0.07256 0.04873 -1.489 0.143
> z.diff.lag 0.22028 0.13836 1.592 0.118
>
> Residual standard error: 0.0248 on 50 degrees of freedom
> Multiple R-squared: 0.07166, Adjusted R-squared: 0.03452
> F-statistic: 1.93 on 2 and 50 DF, p-value: 0.1559
>
>
> Value of test-statistic is: -1.4891 1.4462
>
> Critical values for test statistics:
> 1pct 5pct 10pct
> tau2 -3.51 -2.89 -2.58
> phi1 6.70 4.71 3.86
Δ Y(t )= a 0 + ga m m a ∗ y(t − 1 )+ a 2 (t )+ e (t )
(ここでa2(t)は時間トレンドの項です)(Enders p。208からの)仮説は次のとおりです:tau:gamma = 0 phi3:gamma = a2 = 0 phi2:a0 = gamma = a2 = 0これはR出力。この場合、テスト統計は-2.4216 2.1927 2.9343です。これらすべてのケースで、これらは「nullを拒否できない」ゾーンに含まれます(以下の重要な値を参照)。上記のように、tau3が意味するのは、ユニットルートのnullを拒否できず、ユニットルートが存在することを意味します。phi3を拒否しないことは、2つのことを意味します。1)ガンマ= 0(単位根)AND 2)時間トレンド項がない、つまり、a2 = 0。このnullを拒否した場合、これらの項の1つまたは両方が0ではなかったことを意味します。phi2を拒否しないことは、3つのことを意味します。1)ガンマ= 0および2)時間トレンド項なしおよび3)ドリフト項なし、つまり= 0、a0 = 0、a2 = 0。
こちらがR出力
> summary(ur.df(y=income, type = "trend",lags=1))
>
> ###############################################
> # Augmented Dickey-Fuller Test Unit Root Test #
> ###############################################
>
> Test regression trend
>
>
> Call:
> lm(formula = z.diff ~ z.lag.1 + 1 + tt + z.diff.lag)
>
> Residuals:
> Min 1Q Median 3Q Max
> -0.036693 -0.016457 -0.000435 0.014344 0.074299
>
> Coefficients:
> Estimate Std. Error t value Pr(>|t|)
> (Intercept) 1.0369478 0.4272693 2.427 0.0190 *
> z.lag.1 -0.1767666 0.0729961 -2.422 0.0192 *
> tt 0.0006299 0.0003348 1.881 0.0659 .
> z.diff.lag 0.2557788 0.1362896 1.877 0.0665 .
> ---
> Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
>
> Residual standard error: 0.02419 on 49 degrees of freedom
> Multiple R-squared: 0.1342, Adjusted R-squared: 0.08117
> F-statistic: 2.531 on 3 and 49 DF, p-value: 0.06785
>
>
> Value of test-statistic is: -2.4216 2.1927 2.9343
>
> Critical values for test statistics:
> 1pct 5pct 10pct
> tau3 -4.04 -3.45 -3.15
> phi2 6.50 4.88 4.16
> phi3 8.73 6.49 5.47
上記の特定の例では、d.Ausデータの場合、両方の検定統計量が「不合格失敗」ゾーン内にあるため、gamma = 0 AND a0 = 0であることを意味し、単位根があることを意味しますが、ドリフト項はありません。