σμ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
  • 日本語
  • 中文
Optimization スキル›モジュール 4 · 制約と正則化›トピック 4.3·~50 分中級

最適化としてのL1・L2正則化

スキル上の位置モジュール 4/5 · トピック 3/4
4.3 最適化としてのL1・L2正則化
—ラボスコア— / 100合計XP—連続日数—

縮める罰則

正則化は、大きな係数にコストがかかるように、当てはめの損失にペナルティを加えます:min⁡w ∥y−Xw∥2+λ Ω(w)\min_{\mathbf w}\ \lVert\mathbf y-X\mathbf w\rVert^2+\lambda\,\Omega(\mathbf w)minw​ ∥y−Xw∥2+λΩ(w)。リッジは Ω=∥w∥22\Omega=\lVert\mathbf w\rVert_2^2Ω=∥w∥22​、ラッソは Ω=∥w∥1\Omega=\lVert\mathbf w\rVert_1Ω=∥w∥1​ を使います。これも普通の最適化問題で、Ω(w)≤t\Omega(\mathbf w)\le tΩ(w)≤t の下で損失を最小化する等価な制約付きの形があり、そこでは λ\lambdaλ がまさに予算 ttt のKKT乗数です。λ\lambdaλ が大きいほど予算は小さくなります。

min⁡w 12(w−z)2+λ∣w∣ ⇒ w=sign(z)max⁡(∣z∣−λ, 0),min⁡w 12(w−z)2+λ2w2 ⇒ w=z1+λ\min_{w}\ \tfrac12(w-z)^{2}+\lambda\lvert w\rvert\ \Rightarrow\ w=\mathrm{sign}(z)\max(\lvert z\rvert-\lambda,\,0),\qquad \min_{w}\ \tfrac12(w-z)^{2}+\tfrac{\lambda}{2}w^{2}\ \Rightarrow\ w=\frac{z}{1+\lambda}wmin​ 21​(w−z)2+λ∣w∣ ⇒ w=sign(z)max(∣z∣−λ,0),wmin​ 21​(w−z)2+2λ​w2 ⇒ w=1+λz​(1)

1変数の問題で違いがわかります。リッジは罰則なしの値 zzz を 1+λ1+\lambda1+λ で割るので、すべての係数が同じ割合で縮みますがゼロにはなりません。ラッソは ∣z∣|z|∣z∣ から λ\lambdaλ を引いてゼロで切り捨てるので、∣z∣≤λ|z|\le\lambda∣z∣≤λ の係数はちょうどゼロになります。これがラッソが変数を選択する理由です。幾何的には、ℓ1\ell_1ℓ1​ の球は軸上に角を持ち、損失の等高線はたいていその角で接しますが、丸い ℓ2\ell_2ℓ2​ の球には角がありません。解析的には、∣w∣|w|∣w∣ はゼロで傾きが −λ-\lambda−λ から +λ+\lambda+λ へ跳ぶ折れ点を持ち、最小点がその折れ点に位置することがあります。

2つの罰則の数値例λ=1\lambda=1λ=1 とします。z=3z=3z=3 ではラッソは max⁡(3−1,0)=2\max(3-1,0)=2max(3−1,0)=2、リッジは 3/2=1.53/2=1.53/2=1.5 です。z=0.8z=0.8z=0.8 ではラッソは 0.8≤10.8\le10.8≤1 なので 000、リッジはまだ 0.40.40.4 です。回帰の例:x=(1,2)x=(1,2)x=(1,2)、y=(1,3)y=(1,3)y=(1,3) では xTx=5x^{T}x=5xTx=5、xTy=7x^{T}y=7xTy=7 で、リッジは w=7/(5+λ)w=7/(5+\lambda)w=7/(5+λ) となり、λ=0\lambda=0λ=0 で 1.41.41.4、λ=2\lambda=2λ=2 でちょうど 111 です。

金融では、リッジのペナルティは相関する説明変数への回帰を安定させ(共分散の縮小推定の親戚です)、ラッソは疎なインデックス・トラッキング・ポートフォリオを作ります。落とし穴が1つあります。ロングオンリーで全額投資のポートフォリオでは ∑i∣wi∣=∑iwi=1\sum_i|w_i|=\sum_i w_i=1∑i​∣wi​∣=∑i​wi​=1 なので、ペナルティ λ∥w∥1\lambda\lVert\mathbf w\rVert_1λ∥w∥1​ は定数 λ\lambdaλ にすぎず、何も変えません。非ゼロの持ち高や取引を少なくしたいなら、取引量 ∥w−w′∥1\lVert\mathbf w-\mathbf w'\rVert_1∥w−w′∥1​ に罰則を課します。これが総仕上げの回転率とコストの上限がまさにしていることです。ベイズ的には、リッジは係数のガウス事前分布に、ラッソはラプラス事前分布に対応します。

  • リッジ:1変数では w=z/(1+λ)w=z/(1+\lambda)w=z/(1+λ)。滑らかに縮み、ちょうどゼロにはなりません。
  • ラッソ:ソフト閾値 sign(z)max⁡(∣z∣−λ,0)\mathrm{sign}(z)\max(|z|-\lambda,0)sign(z)max(∣z∣−λ,0)。小さな係数をちょうどゼロにします。
  • 制約付きの形:重み λ\lambdaλ は予算 Ω(w)≤t\Omega(\mathbf w)\le tΩ(w)≤t のKKT乗数です。
λは標本外で選ぶλ\lambdaλ を大きくすると分散とバイアスのトレードオフになり、係数は小さく安定しますが、系統的にゼロへ引き寄せられます。λ\lambdaλ は標本内の当てはまり(常に λ=0\lambda=0λ=0 を好みます)ではなく、交差検証かウォークフォワード検証で選びます。時系列では時間順を守ったブロックを使います。

ラボ:採点される内容

  • ラッソの係数10点
  • リッジの係数10点
  • ゼロになる係数15点
  • 1変数のリッジ回帰15点
  • ラッソが正確なゼロを作る理由15点
  • 意味のない罰則15点
  • 事前分布と効果10点
  • λを大きくすると何が起きるか10点
0/ 8問回答済み
問題
1数値で答える10点未回答

罰則なしの値は z=3z=3z=3、λ=1\lambda=1λ=1 です。ラッソ(ソフト閾値)の値はいくつですか。

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

同じ z=3z=3z=3、λ=1\lambda=1λ=1 で、リッジ z/(1+λ)z/(1+\lambda)z/(1+λ) の値はいくつですか。

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

今度は罰則なしの値が z=0.8z=0.8z=0.8、λ=1\lambda=1λ=1 です。ラッソの値はいくつですか。

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

x=(1,2)x=(1,2)x=(1,2)、y=(1,3)y=(1,3)y=(1,3) のとき、リッジは w=xTy/(xTx+λ)w=x^{T}y/(x^{T}x+\lambda)w=xTy/(xTx+λ) です。λ=2\lambda=2λ=2 での www はいくつですか。

w=w =w=許容:± 0.01
5選択15点未回答

ラッソが一部の係数をちょうどゼロにし、リッジはそうしない理由を説明しているのはどれですか。

6選択15点未回答

ポートフォリオはロングオンリー(wi≥0w_i\ge0wi​≥0)で全額投資(∑iwi=1\sum_i w_i=1∑i​wi​=1)です。λ∥w∥1\lambda\lVert\mathbf w\rVert_1λ∥w∥1​ を加えても何も変わらないのはなぜですか。

7対応づけ10点未回答

各項目をその説明に結び付けてください。

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

リッジ
ラッソ
大きい λ\lambdaλ
小さい λ\lambdaλ
8言葉で説明する10点未回答

1〜2文で、ペナルティの重み λ が大きくなると、係数と当てはめのバイアス・分散に何が起きるかを述べてください。

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

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

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

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