σμ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 スキル›モジュール 3 · 断面と裾›トピック 3.2·~1 時間上級

ファマ・マクベス回帰

スキル上の位置モジュール 3/3 · トピック 2/3
3.2 ファマ・マクベス回帰
—ラボスコア— / 100合計XP—連続日数—

市場が特性に払う対価を推定する

モメンタムが高い、またはボラティリティが低い株は、翌月により多く稼ぐのか、稼ぐならいくらか。全銘柄・全月を1回で回帰すると 18,00018{,}00018,000 件の観測を独立なものとしてまとめますが、同じ月の株は同じ市場ショックを共有するため誤差が相関し、t統計量が大きくなりすぎます。ファマとマクベスは2段階でこれを避けました。まず月ごとに、翌月のリターンを現在の特性に回帰する断面回帰を行い、各月が特性ごとに1つの傾き(その月に市場が払った価格)を与えます。次に、この月次の傾きを時系列として扱います。プレミアムはその平均で、標準誤差はそのばらつきから se=s/T\mathrm{se}=s/\sqrt{T}se=s/T​ と求まり、t統計量は λˉ/(s/T)\bar\lambda/(s/\sqrt T)λˉ/(s/T​) です。

yi,t+1=λ0,t+λ1,txi,t+ui,t+1,λ^=1T∑t=1Tλ^t,t=λ^sλ/Ty_{i,t+1}=\lambda_{0,t}+\lambda_{1,t}x_{i,t}+u_{i,t+1},\qquad \hat\lambda=\frac1T\sum_{t=1}^{T}\hat\lambda_t,\qquad t=\frac{\hat\lambda}{s_\lambda/\sqrt{T}}yi,t+1​=λ0,t​+λ1,t​xi,t​+ui,t+1​,λ^=T1​t=1∑T​λ^t​,t=sλ​/T​λ^​(1)
1か月と4か月ある月に、4銘柄の特性が x=−3,−1,1,3x=-3,-1,1,3x=−3,−1,1,3(すでに平均ゼロ)、翌月リターンが y=−1,0,3,2y=-1,0,3,2y=−1,0,3,2(%)だとします。∑xy=12\sum xy=12∑xy=12、Sxx=20S_{xx}=20Sxx​=20 なので傾きは 0.60.60.6、切片は yˉ=1.0\bar y=1.0yˉ​=1.0 です。4か月分の傾きが 0.6,−0.2,0.5,0.10.6,-0.2,0.5,0.10.6,−0.2,0.5,0.1 だとすると、平均は 0.250.250.25、偏差の2乗和は 0.410.410.41 で s2=0.41/3=0.1367s^2=0.41/3=0.1367s2=0.41/3=0.1367、s=0.370s=0.370s=0.370、se=0.370/4=0.185\mathrm{se}=0.370/\sqrt4=0.185se=0.370/4​=0.185、t=0.25/0.185=1.35t=0.25/0.185=1.35t=0.25/0.185=1.35 です。実データでは、2001年1月から2026年7月までの 307307307 か月、1か月あたり約 565656 銘柄(合計 17,07517{,}07517,075 銘柄月)で、毎月標準化して比べられるようにした3つの特性を使います。直近1か月を除く12か月モメンタム、直近1か月のリターン、12か月ボラティリティです。月次の傾きの平均は、モメンタムが +0.178%+0.178\%+0.178%(t=1.68t=1.68t=1.68)、直近1か月が −0.071%-0.071\%−0.071%(t=−0.83t=-0.83t=−0.83)、ボラティリティが +0.270%+0.270\%+0.270%(t=2.00t=2.00t=2.00、年率 3.2%3.2\%3.2%)です。1回のプール回帰ではボラティリティの t=5.47t=5.47t=5.47 とファマ・マクベスの値の2倍以上になり、月でクラスタリングすると 2.352.352.35 です。ボラティリティのプレミアムは標本の前半・後半ともに正(月 0.27%0.27\%0.27% と 0.27%0.27\%0.27%)ですが、銘柄は現在の大企業なので、生存バイアスにより NVIDIA のようなボラティリティの高い勝者が多く含まれます。

月次の傾き自体が自己相関することもあるので、多くの研究は s/Ts/\sqrt Ts/T​ を数ラグのニューウィー・ウェスト標準誤差に置き換えます。ここでは影響は小さく、ボラティリティの t=2.00t=2.00t=2.00 は4ラグで 1.991.991.99 になります。もう1点、切片 λ^0\hat\lambda_0λ^0​ は平均的な特性を持つ株の平均リターンで、この標本では月 1.00%1.00\%1.00%(t=4.08t=4.08t=4.08)であり、主に市場と生き残った銘柄を映しています。

  • 2段階:月ごとに回帰し、次に月次の傾きの平均とばらつきを見る。
  • 特性を毎月標準化し、傾きを比較可能にして外れ値の影響を減らす。

ラボ:採点される内容

  • 1か月の傾き15点
  • プレミアム10点
  • 傾きのばらつき10点
  • 標準誤差10点
  • t統計量15点
  • プール回帰とファマ・マクベス15点
  • 生存バイアスの影響25点
0/ 7問回答済み
問題
1数値で答える15点未回答

4銘柄の平均ゼロの特性が x=−3,−1,1,3x=-3,-1,1,3x=−3,−1,1,3、翌月リターンが y=−1,0,3,2y=-1,0,3,2y=−1,0,3,2 です。∑xy=12\sum xy=12∑xy=12、Sxx=20S_{xx}=20Sxx​=20 のとき、断面の傾きはいくつですか。小数第1位まで答えてください。

b=b =b=許容:± 0.01
2数値で答える10点未回答

4か月分の傾きは 0.6,−0.2,0.5,0.10.6,-0.2,0.5,0.10.6,−0.2,0.5,0.1 です。平均(推定されたプレミアム)はいくつですか。小数第2位まで答えてください。

m=m =m=許容:± 0.005
3数値で答える10点未回答

平均 0.250.250.25 からの偏差の2乗和は 0.410.410.41 です。T=4T=4T=4 か月のとき、標準偏差 s=0.41/(T−1)s=\sqrt{0.41/(T-1)}s=0.41/(T−1)​ はいくつですか。0.1367=0.370\sqrt{0.1367}=0.3700.1367​=0.370 を使い、小数第3位まで答えてください。

s=s =s=許容:± 0.004
4数値で答える10点未回答

s=0.370s=0.370s=0.370、T=4T=4T=4 か月のとき、se=s/T\mathrm{se}=s/\sqrt Tse=s/T​ はいくつですか。小数第3位まで答えてください。

se=se =se=許容:± 0.002
5数値で答える15点未回答

プレミアムは 0.250.250.25、標準誤差は 0.1850.1850.185 です。t統計量はいくつですか。小数第2位まで答えてください。

t=t =t=許容:± 0.03
6選択15点未回答

プール回帰ではボラティリティ・プレミアムの t = 5.47、ファマ・マクベスでは t = 2.00 です。なぜプールした値は大きすぎるのですか。

7言葉で説明する25点未回答

銘柄リストは現在の大企業で、ボラティリティ・プレミアムは正になりました。生存バイアスがプレミアムを膨らませうる理由を1〜2文で説明してください。

数学モードキーボードまたは LaTeX で入力問題1/7 · Ctrl ⏎ で確認

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

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

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