Rの関数:kpss.test {tseries}

x

検定の対象となる数値ベクトル、または時系列オブジェクトを指定します。

null

帰無仮説の種類を指定する引数です。「データはどのような種類の定常性を持つか」を仮定します。

• "Level"（デフォルト）: Stationarity round a constant level(本ポストでは以降「水準定常性」と呼称します)を帰無仮説とします。データが時間を通じて一定の平均値の周りを変動している状態を「定常」と仮定します。
• "Trend": Stationarity around a deterministic trend）(本ポストでは以降「トレンド定常性」と呼称します)を帰無仮説とします。データが確定的な時間トレンド（直線的な上昇または下降傾向）の周りを変動している状態を「定常」と仮定します。

データの形状に応じて適切な帰無仮説を選択する必要があります。

lshort

検定統計量の計算に使用するラグの切り捨てパラメータの計算方法を制御する論理値です。

• TRUE（デフォルト）: “short”バージョンとして、trunc(4*(length(x)/100)^0.25) という計算式でラグを決定します。
• FALSE: “long”バージョンとして、trunc(12*(length(x)/100)^0.25) という計算式でラグを決定します。

こちらののパラメータは、誤差項の長期的な分散を推定する際に使用されます。

kpss.test検定結果の確認

この出力は、水準定常な過程として意図的に生成された data_level_stationary に対して、null = "Level"（水準定常性）を帰無仮説としてKPSS検定を実行した結果です。

    KPSS Test for Level Stationarity

data:  data_level_stationary
KPSS Level = 0.07606, Truncation lag parameter = 5, p-value = 0.1

結論の導出

KPSS検定は、ADF検定やPP検定とは帰無仮説が逆であるため、結論の解釈には注意が必要です。

• 帰無仮説 (H₀): データは水準定常である。
• 対立仮説 (H₁): データは単位根を持つ（非定常である）。
• 判断ルール: p値 < 有意水準 であれば、帰無仮説を棄却する（非定常であると判断する）。

今回の結果を見てみましょう。

• p値: 0.1
• 設定した有意水準: 0.05

0.1 > 0.05ですので、p値は有意水準よりも大きいです。

ですので、帰無仮説「データは水準定常である」を棄却することはできません。

この検定結果は、データが定常であるという仮説と矛盾しないことを示しています。今回のデータは実際に定常な時系列として生成したものですので、KPSS検定がデータの性質を正しく捉えていることが確認できます。

なお、p値の最大値は0.1として出力される仕様です。

今回の検定過程において、

警告メッセージ:
kpss.test(data_level_stationary, null = "Level") で:
  p-value greater than printed p-value

と出力されていますので、p-value = 0.1は「p値は0.1以上である」と解釈します。

# --- ケース2: 単位根データ ---
cat("--- ケース2：単位根データに null='Level' で検定 ---\n")
# データは非定常ですので、帰無仮説「水準定常である」は棄却される（p値が小さい）と予測されます。
kpss_result2 <- kpss.test(data_unit_root, null = "Level")
print(kpss_result2)

--- ケース2：単位根データに null='Level' で検定 ---

    KPSS Test for Level Stationarity

data:  data_unit_root
KPSS Level = 2.7142, Truncation lag parameter = 5, p-value = 0.01

kpss.test検定結果の確認

この出力は、非定常な単位根過程（ランダムウォーク）として生成された data_unit_root に対して、null = "Level"（水準定常性）を帰無仮説としてKPSS検定を実行した結果です。

    KPSS Test for Level Stationarity

data:  data_unit_root
KPSS Level = 2.7142, Truncation lag parameter = 5, p-value = 0.01

結論の導出

KPSS検定の逆の仮説設定になります。

• 帰無仮説 (H₀): データは水準定常である。
• 対立仮説 (H₁): データは単位根を持つ（非定常である）。
• 判断ルール: p値 < 有意水準 であれば、帰無仮説を棄却する（非定常であると判断する）。

今回の結果を見てみましょう。

• p値: 0.01
• 設定した有意水準: 0.05

0.01 < 0.05ですので、p値は有意水準よりも小さいです。

ですので、帰無仮説「データは水準定常である」を棄却します。

この検定結果は、対立仮説である「データは単位根を持つ（非定常である）」ことを支持します。今回のデータは実際に非定常な単位根過程として生成したものですので、整合が取れていることになります。

なお、p値の最小値は0.01として出力される仕様です。

今回の検定過程において、

警告メッセージ:
kpss.test(data_unit_root, null = "Level") で:
  p-value smaller than printed p-value

と出力されていますので、p-value = 0.01は「p値は0.01以下である」と解釈します。

# --- ケース3: トレンド定常データ ---
cat("--- ケース3：トレンド定常データに null='Trend' で検定 ---\n")
# データはトレンド定常ですので、帰無仮説「トレンド定常である」は棄却されない（p値が大きい）と予測されます。
kpss_result3 <- kpss.test(data_trend_stationary, null = "Trend")
print(kpss_result3)

--- ケース3：トレンド定常データに null='Trend' で検定 ---

    KPSS Test for Trend Stationarity

data:  data_trend_stationary
KPSS Trend = 0.047421, Truncation lag parameter = 5, p-value = 0.1

kpss.test検定結果の確認

この出力は、トレンド定常な過程として意図的に生成された data_trend_stationary に対して、null = "Trend"（トレンド定常性）を帰無仮説としてKPSS検定を実行した結果です。

    KPSS Test for Trend Stationarity

data:  data_trend_stationary
KPSS Trend = 0.047421, Truncation lag parameter = 5, p-value = 0.1

結論の導出

KPSS検定の仮説設定はADF/PP検定とは逆です。

• 帰無仮説 (H₀): データはトレンド定常である。
• 対立仮説 (H₁): データは単位根を持つ（非定常である）。
• 判断ルール: p値 < 有意水準 であれば、帰無仮説を棄却する（非定常であると判断する）。

今回の結果を見てみましょう。

• p値: 0.1
• 設定した有意水準: 0.05

0.1 > 0.05ですので、p値は有意水準よりも大きいです。

ですので、帰無仮説「データはトレンド定常である」を棄却することはできません。

この検定結果は、データがトレンド定常であるという仮説と矛盾しないことを示しています。今回のデータは実際にトレンド定常な時系列として生成したものですので、KPSS検定がデータの性質を正しく捉え、「定常である」という帰無仮説は棄却されないことを確認できます。

# --- ケース4: ドリフト付き単位根データ ---
cat("--- ケース4：ドリフト付き単位根データに null='Trend' で検定 ---\n")
# データは非定常ですので、帰無仮説「トレンド定常である」は棄却される（p値が小さい）と予測されます。
kpss_result4 <- kpss.test(data_unit_root_drift, null = "Trend")
print(kpss_result4)

--- ケース4：ドリフト付き単位根データに null='Trend' で検定 ---

    KPSS Test for Trend Stationarity

data:  data_unit_root_drift
KPSS Trend = 0.12495, Truncation lag parameter = 5, p-value = 0.08899

kpss.test検定結果の確認

この出力は、data_unit_root_drift（ドリフト付き単位根過程）に対して、null = "Trend"（トレンド定常性）を帰無仮説としてKPSS検定を実行した結果です。

    KPSS Test for Trend Stationarity

data:  data_unit_root_drift
KPSS Trend = 0.29387, Truncation lag parameter = 5, p-value = 0.01

結論の導出

KPSS検定の逆の仮説設定を再度確認します。

• 帰無仮説 (H₀): データはトレンド定常である。
• 対立仮説 (H₁): データは単位根を持つ（非定常である）。
• 判断ルール: p値 < 有意水準 であれば、帰無仮説を棄却する（非定常であると判断する）。

今回の結果を見てみましょう。

• p値: 0.01
• 設定した有意水準: 0.05

0.01 < 0.05ですので、p値は有意水準よりも明確に小さいです。

ですので、帰無仮説「データはトレンド定常である」を棄却します。実際に非定常なデータを生成していますので、検定結果は整合が取れていることになります。

以上です。