σμ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
  • 日本語
  • 中文
Statistics スキル›モジュール 3 · モデルと回帰›トピック 3.1·~55 分中級

最尤法

スキル上の位置モジュール 3/4 · トピック 1/4
3.1 最尤法
—ラボスコア— / 100合計XP—連続日数—

データを最も起こりやすくするパラメータ

データが未知のパラメータ θ\thetaθ を持つ分布族から来ていると仮定します(正規分布の平均と分散、ポアソン分布の率など)。尤度 L(θ)=∏if(xi;θ)L(\theta)=\prod_i f(x_i;\theta)L(θ)=∏i​f(xi​;θ) は、実際に観測したデータの確率を θ\thetaθ の関数として見たものです。最尤推定は、観測したデータを最も起こりやすくする θ^\hat\thetaθ^ を選びます。回帰、ボラティリティモデル、機械学習の大半を支える主力の手法です。

ℓ(θ)=∑i=1nln⁡f(xi;θ),θ^=arg⁡max⁡θ ℓ(θ),∂ℓ∂θ∣θ^=0\ell(\theta)=\sum_{i=1}^{n}\ln f(x_i;\theta),\qquad \hat\theta=\arg\max_\theta\,\ell(\theta),\qquad \frac{\partial\ell}{\partial\theta}\Big|_{\hat\theta}=0ℓ(θ)=i=1∑n​lnf(xi​;θ),θ^=argθmax​ℓ(θ),∂θ∂ℓ​​θ^​=0(1)

実際には対数尤度 ℓ\ellℓ を最大化します。小さな数の多数の積はアンダーフローしやすく微分もしにくい一方、対数の和は安定で、対数が増加関数なので最大点は同じです。微分を0とおくと、標準的な場合に閉じた式が得られます。ポアソンでは λ^=xˉ\hat\lambda=\bar xλ^=xˉ、ベルヌーイでは p^\hat pp^​ が成功の割合、正規分布では μ^=xˉ\hat\mu=\bar xμ^​=xˉ、σ^2=1n∑(xi−xˉ)2\hat\sigma^2=\frac1n\sum(x_i-\bar x)^2σ^2=n1​∑(xi​−xˉ)2 です。n−1n-1n−1 でなく nnn で割るので、小標本では分散がやや小さめに推定されます。

最尤推定の3つの簡単な例1分間の約定数 2,4,3,32,4,3,32,4,3,3 からポアソン率は λ^=12/4=3\hat\lambda=12/4=3λ^=12/4=3 です。40回のうち12回勝てば p^=12/40=0.3\hat p=12/40=0.3p^​=12/40=0.3 で、その対数尤度は 12ln⁡0.3+28ln⁡0.7≈−24.43512\ln0.3+28\ln0.7\approx-24.43512ln0.3+28ln0.7≈−24.435 で、他のどの ppp でもこれより小さくなります。1,2,3,4,51,2,3,4,51,2,3,4,5 では正規分布の最尤推定は μ^=3\hat\mu=3μ^​=3、σ^2=10/5=2\hat\sigma^2=10/5=2σ^2=10/5=2 で、不偏標本分散は 10/4=2.510/4=2.510/4=2.5 です。最尤推定量のほうがやや小さく、その差は nnn が大きくなると消えます。
  • 尤度:観測データの確率をパラメータの関数として見たもの。
  • 最尤推定量:尤度を最大にするパラメータ。安定のため対数尤度を最大化する。
  • 性質:緩い条件のもとで一致性と漸近正規性を持つが、小標本で常に不偏とは限らない。
誤ったモデルへのよい当てはめは誤りのまま最尤法が見つけるのは、選んだ分布族の中での最良のパラメータです。日次リターンの裾が重いのに正規分布を当てはめると、最良の正規分布が得られますが、その分散は暴落リスクを依然として過小評価します。信頼する前に、当てはめたモデルをデータ(テールの件数やヒストグラム)と比べ、スチューデントtのような裾の重い分布族も検討しましょう。

ラボ:採点される内容

  • ポアソン率10点
  • 正規分布の分散15点
  • 勝率10点
  • 対数尤度15点
  • 最尤推定が行うこと15点
  • 最尤推定量の性質15点
  • なぜ対数尤度を最大化するか20点
0/ 7問回答済み
問題
1数値で答える10点未回答

4分間の約定数は 2, 4, 3, 3 です。ポアソン率の最尤推定値はいくつですか。

λ^=\hat\lambda =λ^=許容:± 0.001
2数値で答える15点未回答

標本 1, 2, 3, 4, 5 について、正規分布の分散の最尤推定値(n で割る)はいくつですか。

σ^2=\hat\sigma^{2} =σ^2=許容:± 0.001
3数値で答える10点未回答

40回の取引のうち12回が勝ちでした。勝つ確率の最尤推定値はいくつですか。

p^=\hat p =p^​=許容:± 0.001
4数値で答える15点未回答

対数尤度 12 ln(p) + 28 ln(1 - p) を p = 0.3 で評価してください(自然対数)。小数第3位まで答えてください。

ℓ=\ell =ℓ=許容:± 0.01
5選択15点未回答

最尤推定は何を選びますか。

6選択15点未回答

正しい記述をすべて選んでください。

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

尤度ではなく対数尤度を最大化する理由を2文で述べてください。

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

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

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

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