σμ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
  • 日本語
  • 中文
Bayesian Inference スキル›モジュール 1 · 信念と証拠›トピック 1.4·~1 時間中級

ベイズ線形回帰

スキル上の位置モジュール 1/3 · トピック 4/4
1.4 ベイズ線形回帰
—ラボスコア— / 100合計XP—連続日数—

正直な不確実性を持つ回帰

通常の回帰は傾きを1つだけ返します。ベイズ回帰は傾きの分布全体を返します。株のリターン yyy と市場のリターン xxx を y=βx+ey=\beta x+ey=βx+e(ノイズ eee の分散は σ2\sigma^{2}σ2)とします。ベータの事前分布として、たとえば m0=1m_0=1m0​=1 を中心とし分散 v0v_0v0​ の正規分布から始めます。データは傾きを最小二乗の値のほうへ押し、事前分布は押し戻します。これは前のトピックと同じ、精度で重み付けした平均です。事前精度は 1/v01/v_01/v0​、データの精度は ∑x2/σ2\sum x^{2}/\sigma^{2}∑x2/σ2 で、後者は市場が大きく動いたときに大きくなります。

Pn=1v0+∑x2σ2,mn=m0/v0+∑xy/σ2Pn,sd(β)=1PnP_n=\frac{1}{v_0}+\frac{\sum x^{2}}{\sigma^{2}},\qquad m_n=\frac{m_0/v_0+\sum xy/\sigma^{2}}{P_n},\qquad \mathrm{sd}(\beta)=\frac{1}{\sqrt{P_n}}Pn​=v0​1​+σ2∑x2​,mn​=Pn​m0​/v0​+∑xy/σ2​,sd(β)=Pn​​1​(1)

注目点が2つあります。第一に、000 を中心とする事前分布はちょうどリッジ回帰になり、傾きは λ=σ2/v0\lambda=\sigma^{2}/v_0λ=σ2/v0​ として ∑xy / (∑x2+λ)\sum xy\,/\,(\sum x^{2}+\lambda)∑xy/(∑x2+λ) です。リッジ回帰は、大きな傾きを疑う事前分布を持つベイズ回帰にほかなりません。第二に、新しい x∗x^{*}x∗ の予測の広がりには2つの源があり、ノイズ σ2\sigma^{2}σ2 と傾きの不確実性(x∗2/Pnx^{*2}/P_nx∗2/Pn​ を加える)です。後者は x∗x^{*}x∗ の二乗で大きくなるので、データから遠い予測ほど不確かです。

数値で見るベータ∑x2=40\sum x^{2}=40∑x2=40、∑xy=48\sum xy=48∑xy=48、σ2=4\sigma^{2}=4σ2=4 とします。最小二乗の傾きは 48/40=1.248/40=1.248/40=1.2 です。事前分布 N(1, 0.52)N(1,\,0.5^{2})N(1,0.52)(つまり v0=0.25v_0=0.25v0​=0.25)では、精度は 4+40/4=144+40/4=144+40/4=14、平均は (4+12)/14=1.143(4+12)/14=1.143(4+12)/14=1.143、標準偏差は 1/14=0.2671/\sqrt{14}=0.2671/14​=0.267 です。000 を中心とするリッジの事前分布なら傾きは 48/(40+16)=0.85748/(40+16)=0.85748/(40+16)=0.857 です。x∗=2x^{*}=2x∗=2 での予測は平均 2.292.292.29、標準偏差 4+4/14=2.07\sqrt{4+4/14}=2.074+4/14​=2.07 です。
  • ベイズの傾き:精度が事前精度とデータ精度の和である分布。
  • リッジ回帰は、000 を中心とする事前分布の場合で λ=σ2/v0\lambda=\sigma^{2}/v_0λ=σ2/v0​。
  • 予測にはノイズと傾きの不確実性が含まれ、後者は x∗2x^{*2}x∗2 で大きくなる。

ラボ:採点される内容

  • 事後精度15点
  • 事後の傾き15点
  • 傾きの不確実性10点
  • リッジの傾き15点
  • 予測の広がり10点
  • 非常に自信のある事前分布10点
  • 用語を結び付ける10点
  • データから遠い所15点
0/ 8問回答済み
問題
1数値で答える15点未回答

∑x2=40\sum x^{2}=40∑x2=40、σ2=4\sigma^{2}=4σ2=4、事前分散 v0=0.25v_0=0.25v0​=0.25 のとき、事後精度 1/v0+∑x2/σ21/v_0+\sum x^{2}/\sigma^{2}1/v0​+∑x2/σ2 はいくつですか。

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

次に事前平均 m0=1m_0=1m0​=1 と ∑xy=48\sum xy=48∑xy=48 を加えます。事後平均 (m0/v0+∑xy/σ2)/14\big(m_0/v_0+\sum xy/\sigma^{2}\big)/14(m0​/v0​+∑xy/σ2)/14 はいくつですか。小数第4位まで答えてください。

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

傾きの事後標準偏差 1/141/\sqrt{14}1/14​ はいくつですか。小数第4位まで答えてください。

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

000 を中心とし v0=0.25v_0=0.25v0​=0.25、σ2=4\sigma^{2}=4σ2=4 の事前分布を使います。つまり λ=σ2/v0=16\lambda=\sigma^{2}/v_0=16λ=σ2/v0​=16 です。∑xy=48\sum xy=48∑xy=48、∑x2=40\sum x^{2}=40∑x2=40 のときのリッジの傾き ∑xy / (∑x2+λ)\sum xy\,/\,(\sum x^{2}+\lambda)∑xy/(∑x2+λ) はいくつですか。小数第4位まで答えてください。

r=r =r=許容:± 0.001
5数値で答える10点未回答

新しい x∗=2x^{*}=2x∗=2 の予測標準偏差は σ2=4\sigma^{2}=4σ2=4 として σ2+x∗2/14\sqrt{\sigma^{2}+x^{*2}/14}σ2+x∗2/14​ です。いくつですか。小数第3位まで答えてください。

p=p =p=許容:± 0.01
6選択10点未回答

事前分散 v0v_0v0​ が非常に小さくなると、事後の傾きはどうなりますか。

7対応づけ10点未回答

各用語をその意味に結び付けてください。

左の各項目に対応するものを選んでください。

最小二乗の傾き
リッジ回帰
事後精度
予測分散
8言葉で説明する15点未回答

1〜2文で、データから遠い xxx の値に対してベイズ回帰の予測範囲が広くなる理由を説明してください。

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

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

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

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