σμ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
  • 日本語
  • 中文
Matrix Decompositions スキル›モジュール 4 · 共分散の推定›トピック 4.3·~55 分上級

レドイト・ウルフの縮小推定

スキル上の位置モジュール 4/5 · トピック 3/3
4.3 レドイト・ウルフの縮小推定
—ラボスコア— / 100合計XP—連続日数—

推定を単純なものに近づける

クリッピングは固有値ごとにシグナルかどうかを決めます。縮小推定はもっと滑らかな道を取り、ノイズの多い標本行列 SSS を、単純で条件の良い目標 FFF と混ぜます。Σ∗=(1−δ)S+δF\Sigma^{*}=(1-\delta)S+\delta FΣ∗=(1−δ)S+δF で、強度 δ∈[0,1]\delta\in[0,1]δ∈[0,1] が目標をどれだけ信頼するかを表します。古典的な目標は平均分散 μ=tr⁡(S)/N\mu=\operatorname{tr}(S)/Nμ=tr(S)/N を使ったスケール付き単位行列 F=μIF=\mu IF=μI で、全資産の分散が同じで無相関だと主張します。これは確実に間違いですが安定していて、少し混ぜるだけで多くのノイズを打ち消せます。

Σ∗=(1−δ)S+δμI,μ=tr⁡SN,δ=min⁡ ⁣(1,b2d2),d2=∥S−μI∥F2N\Sigma^{*}=(1-\delta)S+\delta\mu I,\qquad \mu=\frac{\operatorname{tr}S}{N},\qquad \delta=\min\!\Big(1,\frac{b^{2}}{d^{2}}\Big),\qquad d^{2}=\frac{\lVert S-\mu I\rVert_F^{2}}{N}Σ∗=(1−δ)S+δμI,μ=NtrS​,δ=min(1,d2b2​),d2=N∥S−μI∥F2​​(1)

強度は自由なパラメータではありません。レドイトとウルフは、この目標に対する最適な δ\deltaδ に閉じた式があることを示しました。分子 b2b^2b2 は SSS のノイズの大きさ(1日分の行列 xtxtT\mathbf x_t\mathbf x_t^{T}xt​xtT​ の SSS からの平均2乗距離を TTT で割り、d2d^2d2 で頭打ちにしたもの)を表し、分母 d2d^2d2 は SSS が目標からどれだけ離れているかを表します。ノイズが大きく目標に近い SSS には大きな δ\deltaδ が与えられます。縮小には3つの整った性質があります。トレースを保つこと、各固有値を平均に近づける(λ↦(1−δ)λ+δμ\lambda\mapsto(1-\delta)\lambda+\delta\muλ↦(1−δ)λ+δμ)ので大きい方は縮み小さい方は増えること、そして SSS が特異でも δ>0\delta>0δ>0 なら Σ∗\Sigma^{*}Σ∗ が正定値になることです。TTT が増えるとノイズ項が減り δ→0\delta\to0δ→0 になります。

2×2の例S=(4112)S=\begin{pmatrix}4&1\\1&2\end{pmatrix}S=(41​12​) とすると、μ=(4+2)/2=3\mu=(4+2)/2=3μ=(4+2)/2=3、S−μI=(111−1)S-\mu I=\begin{pmatrix}1&1\\1&-1\end{pmatrix}S−μI=(11​1−1​) のフロベニウスノルムの2乗は4なので d2=4/2=2d^2=4/2=2d2=4/2=2 です。ノイズ項が b2=0.5b^2=0.5b2=0.5 なら強度は δ=0.5/2=0.25\delta=0.5/2=0.25δ=0.5/2=0.25。代わりに δ=0.2\delta=0.2δ=0.2 とすると Σ∗=0.8S+0.2⋅3I=(3.80.80.82.2)\Sigma^{*}=0.8S+0.2\cdot3I=\begin{pmatrix}3.8&0.8\\0.8&2.2\end{pmatrix}Σ∗=0.8S+0.2⋅3I=(3.80.8​0.82.2​) です。SSS の固有値は 3±2=4.4143\pm\sqrt2=4.4143±2​=4.414 と 1.5861.5861.586、Σ∗\Sigma^{*}Σ∗ の固有値は 0.8(4.414)+0.6=4.1310.8(4.414)+0.6=4.1310.8(4.414)+0.6=4.131 と 0.8(1.586)+0.6=1.8690.8(1.586)+0.6=1.8690.8(1.586)+0.6=1.869。条件数は 2.782.782.78 から 2.212.212.21 に改善します。
  • 混合:目標 F=μIF=\mu IF=μI(μ=tr⁡S/N\mu=\operatorname{tr}S/Nμ=trS/N)で Σ∗=(1−δ)S+δF\Sigma^{*}=(1-\delta)S+\delta FΣ∗=(1−δ)S+δF。
  • 強度:δ=min⁡(1,b2/d2)\delta=\min(1,b^2/d^2)δ=min(1,b2/d2)。ノイズが大きく目標に近い行列ほど強く縮小される。
  • 効果:トレース保存、固有値が平均へ近づく、δ>0\delta>0δ>0 なら正定値、TTT が増えると δ→0\delta\to0δ→0。
意図的なバイアス縮小推定は目標が間違っているために意図的にバイアスを加え、その代わりに分散を大きく減らします。推定の期待2乗誤差は下がり、多数の資産の共分散推定ではたいてい非常に有利な取引になります。リッジ回帰と同じバイアスと分散のトレードオフで、これがレドイト・ウルフが多くのライブラリで共分散推定器の既定になっている理由です。

ラボ:採点される内容

  • 目標の水準10点
  • 縮小した行列20点
  • 強度15点
  • 縮小後の固有値15点
  • 新しい条件数15点
  • データが増えれば縮小は弱まる15点
  • バイアスと分散のトレードオフ10点
0/ 7問回答済み
問題
1数値で答える10点未回答

S=(4112)S=\begin{pmatrix}4&1\\1&2\end{pmatrix}S=(41​12​) の目標は μI\mu IμI(μ=tr⁡S/N\mu=\operatorname{tr}S/Nμ=trS/N)です。μ\muμ はいくつですか。

mu=mu =mu=許容:± 0.01
2行列・ベクトル20点未回答

δ=0.2\delta=0.2δ=0.2、目標 3I3I3I のとき、S=(4112)S=\begin{pmatrix}4&1\\1&2\end{pmatrix}S=(41​12​) について Σ∗=0.8S+0.2⋅3I\Sigma^{*}=0.8S+0.2\cdot3IΣ∗=0.8S+0.2⋅3I を計算してください。

Σ∗=\Sigma^{*} =Σ∗=
3数値で答える15点未回答

ノイズ項 b2=0.5b^2=0.5b2=0.5、距離 d2=2d^2=2d2=2 のとき、レドイト・ウルフの強度 δ=min⁡(1,b2/d2)\delta=\min(1,b^2/d^2)δ=min(1,b2/d2) はいくつですか。

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

SSS の最大固有値は 3+2≈4.4143+\sqrt2\approx4.4143+2​≈4.414 です。δ=0.2\delta=0.2δ=0.2、μ=3\mu=3μ=3 のとき、(1−δ)λ+δμ(1-\delta)\lambda+\delta\mu(1−δ)λ+δμ を使って Σ∗\Sigma^{*}Σ∗ の最大固有値を小数第3位まで求めてください。

lambdamax=lambda_max =lambdam​ax=許容:± 0.005
5数値で答える15点未回答

縮小後の行列の固有値は4.131と1.869です。条件数(最大÷最小)を小数第2位まで求めてください。

kappa=kappa =kappa=許容:± 0.01
6選択15点未回答

N を固定したとき、観測数 T が増えるとレドイト・ウルフの最適強度はどうなりますか。

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

誤った目標に向けて縮小しても共分散の推定が改善されうる理由を1〜2文で述べてください。

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

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

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

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