縮める罰則
正則化は、大きな係数にコストがかかるように、当てはめの損失にペナルティを加えます:minw ∥y−Xw∥2+λΩ(w)。リッジは Ω=∥w∥22、ラッソは Ω=∥w∥1 を使います。これも普通の最適化問題で、Ω(w)≤t の下で損失を最小化する等価な制約付きの形があり、そこでは λ がまさに予算 t のKKT乗数です。λ が大きいほど予算は小さくなります。
wmin 21(w−z)2+λ∣w∣ ⇒ w=sign(z)max(∣z∣−λ,0),wmin 21(w−z)2+2λw2 ⇒ w=1+λz(1) 1変数の問題で違いがわかります。リッジは罰則なしの値 z を 1+λ で割るので、すべての係数が同じ割合で縮みますがゼロにはなりません。ラッソは ∣z∣ から λ を引いてゼロで切り捨てるので、∣z∣≤λ の係数はちょうどゼロになります。これがラッソが変数を選択する理由です。幾何的には、ℓ1 の球は軸上に角を持ち、損失の等高線はたいていその角で接しますが、丸い ℓ2 の球には角がありません。解析的には、∣w∣ はゼロで傾きが −λ から +λ へ跳ぶ折れ点を持ち、最小点がその折れ点に位置することがあります。
2つの罰則の数値例λ=1 とします。
z=3 ではラッソは
max(3−1,0)=2、リッジは
3/2=1.5 です。
z=0.8 ではラッソは
0.8≤1 なので
0、リッジはまだ
0.4 です。回帰の例:
x=(1,2)、
y=(1,3) では
xTx=5、
xTy=7 で、リッジは
w=7/(5+λ) となり、
λ=0 で
1.4、
λ=2 でちょうど
1 です。
金融では、リッジのペナルティは相関する説明変数への回帰を安定させ(共分散の縮小推定の親戚です)、ラッソは疎なインデックス・トラッキング・ポートフォリオを作ります。落とし穴が1つあります。ロングオンリーで全額投資のポートフォリオでは ∑i∣wi∣=∑iwi=1 なので、ペナルティ λ∥w∥1 は定数 λ にすぎず、何も変えません。非ゼロの持ち高や取引を少なくしたいなら、取引量 ∥w−w′∥1 に罰則を課します。これが総仕上げの回転率とコストの上限がまさにしていることです。ベイズ的には、リッジは係数のガウス事前分布に、ラッソはラプラス事前分布に対応します。
- リッジ:1変数では w=z/(1+λ)。滑らかに縮み、ちょうどゼロにはなりません。
- ラッソ:ソフト閾値 sign(z)max(∣z∣−λ,0)。小さな係数をちょうどゼロにします。
- 制約付きの形:重み λ は予算 Ω(w)≤t のKKT乗数です。
λは標本外で選ぶλ を大きくすると分散とバイアスのトレードオフになり、係数は小さく安定しますが、系統的にゼロへ引き寄せられます。λ は標本内の当てはまり(常に λ=0 を好みます)ではなく、交差検証かウォークフォワード検証で選びます。時系列では時間順を守ったブロックを使います。 ラボ:採点される内容
- ラッソの係数10点
- リッジの係数10点
- ゼロになる係数15点
- 1変数のリッジ回帰15点
- ラッソが正確なゼロを作る理由15点
- 意味のない罰則15点
- 事前分布と効果10点
- λを大きくすると何が起きるか10点