σμimplied vol · KKmax(S−K,0)udβρᵢⱼfat tailsPnL · drawdownyield · TdSₜ = μSₜ dt + σSₜ dWₜC = S·N(d₁) − K·e⁻ʳᵀ·N(d₂)σₚ² = wᵀΣwd₁ = [ln(S/K) + (r + σ²⁄2)T] ⁄ σ√T∂V/∂t + ½σ²S²·∂²V/∂S² + rS·∂V/∂S − rV = 0Sharpe = (Rₚ − R_f) ⁄ σₚVaR₉₉ = μ − 2.33σρ = Cov(X,Y) ⁄ σₓσᵧE[Rₚ] = Σ wᵢ·E[Rᵢ]φ(x) = e^(−x²/2) ⁄ √2πE[Rᵢ] = R_f + βᵢ(E[Rₘ] − R_f)dXₜ = θ(μ − Xₜ)dt + σ dWₜσ²ₜ = ω + α·ε²ₜ₋₁ + β·σ²ₜ₋₁C − P = S − K·e⁻ʳᵀΔ = ∂V/∂S Γ = ∂²V/∂S²p = (e^(rΔt) − d) ⁄ (u − d)rₜ = ln(Pₜ ⁄ Pₜ₋₁)E[Xₜ₊₁ | ℱₜ] = Xₜf* = (bp − q) ⁄ bz = (x − μ) ⁄ σP(N=k) = λᵏe^(−λ) ⁄ k!Θ = −∂V/∂Tν = S·φ(d₁)·√Tmin wᵀΣw s.t. wᵀμ = μ*∫Σ∂∞λ
uantCraft Academyロードマップスキル
Ctrl K
AO
+1.2k 人の学習者
ログイン参加する ↗
ロードマップスキル
AO
+1.2k 人の学習者
ログイン
  • English
  • Français
  • Español
  • العربية
  • Português
  • Deutsch
  • 日本語
  • 中文
σμimplied vol · KKmax(S−K,0)udβρᵢⱼfat tailsPnL · drawdownyield · TdSₜ = μSₜ dt + σSₜ dWₜC = S·N(d₁) − K·e⁻ʳᵀ·N(d₂)σₚ² = wᵀΣwd₁ = [ln(S/K) + (r + σ²⁄2)T] ⁄ σ√T∂V/∂t + ½σ²S²·∂²V/∂S² + rS·∂V/∂S − rV = 0Sharpe = (Rₚ − R_f) ⁄ σₚVaR₉₉ = μ − 2.33σρ = Cov(X,Y) ⁄ σₓσᵧE[Rₚ] = Σ wᵢ·E[Rᵢ]φ(x) = e^(−x²/2) ⁄ √2πE[Rᵢ] = R_f + βᵢ(E[Rₘ] − R_f)dXₜ = θ(μ − Xₜ)dt + σ dWₜσ²ₜ = ω + α·ε²ₜ₋₁ + β·σ²ₜ₋₁C − P = S − K·e⁻ʳᵀΔ = ∂V/∂S Γ = ∂²V/∂S²p = (e^(rΔt) − d) ⁄ (u − d)rₜ = ln(Pₜ ⁄ Pₜ₋₁)E[Xₜ₊₁ | ℱₜ] = Xₜf* = (bp − q) ⁄ bz = (x − μ) ⁄ σP(N=k) = λᵏe^(−λ) ⁄ k!Θ = −∂V/∂Tν = S·φ(d₁)·√Tmin wᵀΣw s.t. wᵀμ = μ*∫Σ∂∞λ
uantCraft Academyロードマップスキル
Ctrl K
AO
+1.2k 人の学習者
ログイン参加する ↗
ロードマップスキル
AO
+1.2k 人の学習者
ログイン
  • English
  • Français
  • Español
  • العربية
  • Português
  • Deutsch
  • 日本語
  • 中文
statsmodels スキル›モジュール 1 · 回帰とその検証›トピック 1.7·~50 分中級

ニューウェイ–ウェスト標準誤差を手計算

スキル上の位置モジュール 1/3 · トピック 7/10
1.7 ニューウェイ–ウェスト標準誤差を手計算
—ラボスコア— / 100合計XP—連続日数—

誤差が過去を覚えているとき

通常の標準誤差は、各誤差が他と独立だと仮定します。重なり合うリターンはこれを壊します。1〜3月のリターンと2〜4月のリターンは2か月を共有するので、誤差は一緒に動きます。係数は正しいままですが、古典的な標準誤差は小さすぎ、t統計量は大きすぎます。ニューウェイ–ウェスト(HAC)共分散がこれを直します。通常の「サンドイッチ」(XᵀX)⁻¹ S (XᵀX)⁻¹ から出発し、中身を各行と残差の積 xₜeₜ で作ります:S=Σₜeₜ²xₜxₜᵀ+ Σₖ₌₁^LwₖΣₜeₜeₜ₋ₖ (xₜxₜ₋ₖᵀ+ xₜ₋ₖxₜᵀ)。バートレット重み wₖ=1-k/(L+1) は古い自己共分散をゼロに近づけ、行列を有効に保ちます。

statsmodels では sm.OLS(y, X).fit(cov_type='HAC', cov_kwds={'maxlags': L}) で同じ行列が得られ、res.tvalues はそのときロバスト標準誤差を使います。一度は手で書く価値があります。ラグ0ではホワイトの不均一分散に頑健な共分散に一致すること、正の自己相関があるとロバスト分散が古典的な分散より大きくなること、ラグ数の選択は判断であること(よくある目安は重なりの長さ程度)が見えてきます。SPYの3か月リターンでは、2005年5月から2026年8月までの重なり合う256個の観測の平均が 2.78%、古典的な t は 6.10 です。ニューウェイ–ウェストでは2ラグで 4.18、6ラグで 3.65 に下がります。なお有意ですが、最初の数字が示すより証拠は弱いのです。

  • 重なりは記憶を意味する:重なるリターンは誤差が自己相関をもつので、古典的なt統計量は証拠を過大に見せる。
  • サンドイッチ:ニューウェイ–ウェストは (X'X)^-1 S (X'X)^-1 で、S は x_t e_t とバートレット重み付き自己共分散から作る。
  • ライブラリを自分の式で確かめる:手計算の行列がライブラリのものと一致すれば、そのオプションの意味が分かったということ。

ラボ:採点される内容

  • 重なるリターン10点
  • 古典的なt統計量10点
  • バートレット重み10点
  • 手計算のニューウェイ–ウェスト25点
  • statsmodelsでのHAC t統計量15点
  • SPYの研究30点
エクスプローラー · lab-sx-13
▾ lab-sx-13
▾ data
…
bash · ブラウザ内
learner@quantcraft:~/lab-sx-13$
Pythonは未読み込みPythonLn 1, Col 1Spaces: 4UTF-8

このトピックをリセットしますか?

lab-sx-13 の内容はすべてスターターファイルに戻り、ターミナルも消去されます。現在の解答は失われ、元に戻せません。

リセットされるもの
  • main.py → スターターコード
  • ターミナル履歴とテスト出力
  • 未提出の変更
保持されるもの
  • ラボの最高スコア
  • XPと連続日数
  • 講義の進捗とメモ

コードを実行するにはログインしてください。

レクチャーはすべて無料で読めます。コードの実行、テスト、ラボの提出にはアカウントの作成かログインが必要です。

  • ブラウザ内IDEとbashターミナル。インストール不要
  • 即時採点、XP、スキルでの現在地
アカウントを作成ログイン
または
GitHubで続ける