実数値関数の複雑さの尺度
計算学習理論 ( 機械学習 と 計算理論 ) において、 ハンス・ラーデマッハー にちなんで名付けられた ラーデマッハー複雑度は、 確率分布 に関する集合のクラスの豊かさを測定します 。この概念は、実数値関数にも拡張できます。
定義
集合のラデマッハ複雑度
集合Aが与えられたとき 、 A のラデマッハ複雑性は 次のように定義される: [1] [2] : 326
あ
⊆
R
メートル
{\displaystyle A\subseteq \mathbb {R} ^{m}}
ラッド
(
あ
)
:=
1
メートル
え
σ
[
すする
1つの
∈
あ
∑
私
=
1
メートル
σ
私
1つの
私
]
{\displaystyle \operatorname {Rad} (A):={\frac {1}{m}}\mathbb {E} _{\sigma }\left[\sup _{a\in A}\sum _{i=1}^{m}\sigma _{i}a_{i}\right]}
ここで、 は 、 ラデマッハ分布 から抽出された独立したランダム変数です 。 つまり、 、 の場合です。 著者によっては、最大値を取る前に合計の絶対値を取る人もいますが、が 対称で ある場合は 違いはありません。
σ
1
、
σ
2
、
…
、
σ
メートル
{\displaystyle \sigma _{1},\sigma _{2},\dots ,\sigma _{m}}
広報
(
σ
私
=
+
1
)
=
広報
(
σ
私
=
−
1
)
=
1
/
2
{\displaystyle \Pr(\sigma _{i}=+1)=\Pr(\sigma _{i}=-1)=1/2}
私
=
1
、
2
、
…
、
メートル
{\displaystyle i=1,2,\dots ,m}
1つの
=
(
1つの
1
、
…
、
1つの
メートル
)
{\displaystyle a=(a_{1},\ldots ,a_{m})}
あ
{\displaystyle A}
関数クラスのラデマッハ複雑度
を点のサンプルとし、 上の実数値関数の関数クラスを考えます 。 すると 、 の 経験的ラデマッハ複雑度は 次 のように定義されます。
S
=
{
ず
1
、
ず
2
、
…
、
ず
メートル
}
⊂
ず
{\displaystyle S=\{z_{1},z_{2},\dots ,z_{m}\}\subset Z}
ふ
{\displaystyle {\mathcal {F}}}
ず
{\displaystyle Z}
ふ
{\displaystyle {\mathcal {F}}}
S
{\displaystyle S}
ラッド
S
(
ふ
)
=
1
メートル
え
σ
[
すする
ふ
∈
ふ
∑
私
=
1
メートル
σ
私
ふ
(
ず
私
)
]
{\displaystyle \operatorname {Rad} _{S}({\mathcal {F}})={\frac {1}{m}}\mathbb {E} _{\sigma }\left[\sup _{f\in {\mathcal {F}}}\sum _{i=1}^{m}\sigma _{i}f(z_{i})\right]}
これは、前の定義を使って次のように書くこともできる: [2] : 326
ラッド
S
(
ふ
)
=
ラッド
(
ふ
∘
S
)
{\displaystyle \operatorname {Rad} _{S}({\mathcal {F}})=\operatorname {Rad} ({\mathcal {F}}\circ S)}
ここで は 関数の合成 を表します 。つまり、
ふ
∘
S
{\displaystyle {\mathcal {F}}\circ S}
ふ
∘
S
:=
{
(
ふ
(
ず
1
)
、
…
、
ふ
(
ず
メートル
)
)
∣
ふ
∈
ふ
}
{\displaystyle {\mathcal {F}}\circ S:=\{(f(z_{1}),\ldots ,f(z_{m}))\mid f\in {\mathcal {F}}\}}
を 上の確率分布と します 。 サンプルサイズ に対する 関数クラスの ラーデマッハ複雑度 は次のようになります。
ポ
{\displaystyle P}
ず
{\displaystyle Z}
ふ
{\displaystyle {\mathcal {F}}}
ポ
{\displaystyle P}
メートル
{\displaystyle m}
ラッド
ポ
、
メートル
(
ふ
)
:=
え
S
〜
ポ
メートル
[
ラッド
S
(
ふ
)
]
{\displaystyle \operatorname {Rad} _{P,m}({\mathcal {F}}):=\mathbb {E} _{S\sim P^{m}}\left[\operatorname {Rad} _{S}({\mathcal {F}})\right]}
ここで、上記の期待値は、 に従って生成された 同一独立分布 (iid)サンプルに対して取られます 。
S
=
(
ず
1
、
ず
2
、
…
、
ず
メートル
)
{\displaystyle S=(z_{1},z_{2},\dots ,z_{m})}
ポ
{\displaystyle P}
直感
ラデマッハ複雑度は、通常、分類に使用されるモデルの関数クラスに適用され、任意のラベル付けの下で確率空間から抽出された点を分類する能力を測定することを目的としています。関数クラスが十分に豊富な場合、 期待値の下でのランダムな抽出によってシミュレートされたラベルの各配置に適切に適応できる関数が含まれ、合計のこの量が最大化されます。
σ
私
{\displaystyle \sigma_{i}}
例
1.には 単一のベクトル(例 : )が含まれます。
あ
{\displaystyle A}
あ
=
{
(
1つの
、
b
)
}
⊂
R
2
{\displaystyle A=\{(a,b)\}\subset \mathbb {R} ^{2}}
ラッド
(
あ
)
=
1
2
⋅
(
1
4
⋅
(
1つの
+
b
)
+
1
4
⋅
(
1つの
−
b
)
+
1
4
⋅
(
−
1つの
+
b
)
+
1
4
⋅
(
−
1つの
−
b
)
)
=
0
{\displaystyle \operatorname {Rad} (A)={1 \over 2}\cdot \left({1 \over 4}\cdot (a+b)+{1 \over 4}\cdot (ab)+{1 \over 4}\cdot (-a+b)+{1 \over 4}\cdot (-ab)\right)=0}
同じことがすべてのシングルトン仮説クラスにも当てはまります。 [3] : 56
2.には 2 つのベクトル (例 : ) が含まれます。
あ
{\displaystyle A}
あ
=
{
(
1
、
1
)
、
(
1
、
2
)
}
⊂
R
2
{\displaystyle A=\{(1,1),(1,2)\}\subset \mathbb {R} ^{2}}
ラッド
(
あ
)
=
1
2
⋅
(
1
4
⋅
最大
(
1
+
1
、
1
+
2
)
+
1
4
⋅
最大
(
1
−
1
、
1
−
2
)
+
1
4
⋅
最大
(
−
1
+
1
、
−
1
+
2
)
+
1
4
⋅
最大
(
−
1
−
1
、
−
1
−
2
)
)
=
1
8
(
3
+
0
+
1
−
2
)
=
1
4
{\displaystyle {\begin{aligned}\operatorname {Rad} (A)&={1 \over 2}\cdot \left({1 \over 4}\cdot \max(1+1,1+2)+{1 \over 4}\cdot \max(1-1,1-2)+{1 \over 4}\cdot \max(-1+1,-1+2)+{1 \over 4}\cdot \max(-1-1,-1-2)\right)\\[5pt]&={1 \over 8}(3+0+1-2)={1 \over 4}\end{aligned}}}
ラデマッハ複雑性の使用
Rademacher 複雑度は、関数クラスの学習可能性 に関するデータ依存の上限を導出するために使用できます 。直感的には、Rademacher 複雑度が小さい関数クラスの方が学習が容易です。
代表性の制限
機械学習 では、 サンプルデータの真の分布を表す トレーニングセット が必要です。これは、 代表性 の概念を使用して定量化できます 。 サンプルが抽出される確率分布をで表します。仮説(潜在的な分類器)の集合をで表し 、 対応 する 誤差関数の集合をで表します。つまり、すべての仮説に対して、 各トレーニングサンプル(特徴、ラベル)を分類器の誤差にマッピングする 関数が存在します (この場合、仮説と分類器は同じ意味で使用されます)。たとえば、が バイナリ分類器を表す場合、誤差関数は0–1損失関数です。つまり、誤差関数は、サンプルを正しく分類する 場合は0を返し 、そうでない場合は1を返します。基礎となる仮説が無関係な場合は、インデックスを省略して、 代わりにを記述します 。定義:
S
{\displaystyle S}
ポ
{\displaystyle P}
H
{\displaystyle H}
ふ
{\displaystyle F}
h
∈
H
{\displaystyle h\in H}
ふ
h
∈
ふ
{\displaystyle f_{h}\in F}
h
{\displaystyle h}
h
{\displaystyle h}
ふ
h
{\displaystyle f_{h}}
h
{\displaystyle h}
ふ
{\displaystyle f}
ふ
h
{\displaystyle f_{h}}
ら
ポ
(
ふ
)
:=
え
ず
〜
ポ
[
ふ
(
ず
)
]
{\displaystyle L_{P}(f):=\mathbb {E} _{z\sim P}[f(z)]}
–実分布上の 何らかの誤差関数の期待誤差 。
ふ
∈
ふ
{\displaystyle f\in F}
P
{\displaystyle P}
L
S
(
f
)
:=
1
m
∑
i
=
1
m
f
(
z
i
)
{\displaystyle L_{S}(f):={1 \over m}\sum _{i=1}^{m}f(z_{i})}
–サンプル上の 何らかの誤差関数の推定誤差 。
f
∈
F
{\displaystyle f\in F}
S
{\displaystyle S}
および に関する サンプル の代表性は 次のように定義されます。
S
{\displaystyle S}
P
{\displaystyle P}
F
{\displaystyle F}
Rep
P
(
F
,
S
)
:=
sup
f
∈
F
(
L
P
(
f
)
−
L
S
(
f
)
)
{\displaystyle \operatorname {Rep} _{P}(F,S):=\sup _{f\in F}(L_{P}(f)-L_{S}(f))}
代表性は小さいほど良いです。これは、過剰適合を 回避する方法を提供するためです 。つまり、分類器の実際の誤差が推定誤差よりそれほど高くないことを意味し、推定誤差が低い分類器を選択すると、実際の誤差も低くなります。ただし、代表性の概念は相対的であるため、異なるサンプル間で比較することはできません。
サンプルの代表性の期待値は、関数クラスのラデマッハ複雑度によって制限される: [2] :326
E
S
∼
P
m
[
Rep
P
(
F
,
S
)
]
≤
2
⋅
E
S
∼
P
m
[
Rad
(
F
∘
S
)
]
{\displaystyle \mathbb {E} _{S\sim P^{m}}[\operatorname {Rep} _{P}(F,S)]\leq 2\cdot \mathbb {E} _{S\sim P^{m}}[\operatorname {Rad} (F\circ S)]}
一般化誤差の制限
ラデマッハ複雑度が小さい場合、経験的リスク最小化を 使用して仮説クラス H を学習することが可能です 。
例えば、(バイナリ誤差関数の場合) [2] :328 あらゆる に対して 、少なくとも の確率で 、あらゆる仮説に対して :
δ
>
0
{\displaystyle \delta >0}
1
−
δ
{\displaystyle 1-\delta }
h
∈
H
{\displaystyle h\in H}
L
P
(
h
)
−
L
S
(
h
)
≤
2
Rad
(
F
∘
S
)
+
4
2
ln
(
4
/
δ
)
m
{\displaystyle L_{P}(h)-L_{S}(h)\leq 2\operatorname {Rad} (F\circ S)+4{\sqrt {2\ln(4/\delta ) \over m}}}
ラデマッハ複雑性の限界
ラデマッハ複雑度は小さいほど良いので、様々な関数集合のラデマッハ複雑度の上限を持つことは有用である。集合のラデマッハ複雑度の上限を決めるには、以下の規則を使用することができる 。 [2] : 329–330
A
⊂
R
m
{\displaystyle A\subset \mathbb {R} ^{m}}
1. 内のすべてのベクトル が定数ベクトルによって変換される場合 、Rad( A )は変化しません。
A
{\displaystyle A}
a
0
∈
R
m
{\displaystyle a_{0}\in \mathbb {R} ^{m}}
2. 内のすべてのベクトル にスカラーを掛けると 、Rad( A )に が掛けられます 。
A
{\displaystyle A}
c
∈
R
{\displaystyle c\in \mathbb {R} }
|
c
|
{\displaystyle |c|}
3. . [3] : 56
Rad
(
A
+
B
)
=
Rad
(
A
)
+
Rad
(
B
)
{\displaystyle \operatorname {Rad} (A+B)=\operatorname {Rad} (A)+\operatorname {Rad} (B)}
4. (Kakade & Tewari Lemma) 内のすべてのベクトルが Lipschitz 関数 によって操作される場合 、 Rad( A ) は (最大で)関数の Lipschitz 定数 で乗算されます。特に、内のすべてのベクトルが 収縮マッピング によって操作される場合 、 Rad( A ) は厳密に減少します。
A
{\displaystyle A}
A
{\displaystyle A}
5.の 凸包 のラデマッハ複雑度は Rad( A )に等しい。
A
{\displaystyle A}
6. (マサートの補題) 有限集合のラーデマッハ複雑度は集合の大きさとともに対数的に増大する。正式には、 を のベクトル の集合とし 、 を のベクトルの平均とすると 、次のようになる。
A
{\displaystyle A}
N
{\displaystyle N}
R
m
{\displaystyle \mathbb {R} ^{m}}
a
¯
{\displaystyle {\bar {a}}}
A
{\displaystyle A}
Rad
(
A
)
≤
max
a
∈
A
‖
a
−
a
¯
‖
⋅
2
log
N
m
{\displaystyle \operatorname {Rad} (A)\leq \max _{a\in A}\|a-{\bar {a}}\|\cdot {{\sqrt {2\log N}} \over m}}
特に、が バイナリベクトルの集合である場合、ノルムは最大で となるため 、次のようになります。
A
{\displaystyle A}
m
{\displaystyle {\sqrt {m}}}
Rad
(
A
)
≤
2
log
N
m
{\displaystyle \operatorname {Rad} (A)\leq {\sqrt {2\log N \over m}}}
VC 次元 が で ある 集合族 を 仮定します。 の 成長関数は 次のように制限されることが知られています 。
H
{\displaystyle H}
d
{\displaystyle d}
H
{\displaystyle H}
全員 :
m
>
d
+
1
{\displaystyle m>d+1}
Growth
(
H
,
m
)
≤
(
e
m
/
d
)
d
{\displaystyle \operatorname {Growth} (H,m)\leq (em/d)^{d}}
これは、 最大で個の 要素を持つすべての集合に対して、が成り立つことを意味します 。集合族は 上の 2 値ベクトルの集合として考えることができます 。これを Massart の補題に代入すると、次のようになります。
h
{\displaystyle h}
m
{\displaystyle m}
|
H
∩
h
|
≤
(
e
m
/
d
)
d
{\displaystyle |H\cap h|\leq (em/d)^{d}}
H
∩
h
{\displaystyle H\cap h}
R
m
{\displaystyle \mathbb {R} ^{m}}
Rad
(
H
∩
h
)
≤
2
d
log
(
e
m
/
d
)
m
{\displaystyle \operatorname {Rad} (H\cap h)\leq {\sqrt {2d\log(em/d) \over m}}}
より高度な技術( ダドリーのエントロピー境界 とハウスラーの上限 [4] )を用いると、例えば、定数が存在し、 ヴァプニク・チェルヴォネンキス次元 を持つ任意の クラスの ラデマッハ複雑度の上限が であることを示すことができます 。
C
{\displaystyle C}
{
0
,
1
}
{\displaystyle \{0,1\}}
d
{\displaystyle d}
C
d
m
{\displaystyle C{\sqrt {\frac {d}{m}}}}
以下の境界は、 内のベクトル の定数集合に対する線形演算に関連している 。 [2] : 332–333
S
{\displaystyle S}
m
{\displaystyle m}
R
n
{\displaystyle \mathbb {R} ^{n}}
1. 単位球 内のベクトルと 内のベクトルのドット積の集合を定義します 。次に、
A
2
=
{
(
w
⋅
x
1
,
…
,
w
⋅
x
m
)
∣
‖
w
‖
2
≤
1
}
=
{\displaystyle A_{2}=\{(w\cdot x_{1},\ldots ,w\cdot x_{m})\mid \|w\|_{2}\leq 1\}=}
S
{\displaystyle S}
Rad
(
A
2
)
≤
max
i
‖
x
i
‖
2
m
{\displaystyle \operatorname {Rad} (A_{2})\leq {\max _{i}\|x_{i}\|_{2} \over {\sqrt {m}}}}
2. 1 ノルムの単位球内のベクトルと
内のベクトルのドット積の集合を定義します。次に、次の操作を行います。
A
1
=
{
(
w
⋅
x
1
,
…
,
w
⋅
x
m
)
∣
‖
w
‖
1
≤
1
}
=
{\displaystyle A_{1}=\{(w\cdot x_{1},\ldots ,w\cdot x_{m})\mid \|w\|_{1}\leq 1\}=}
S
{\displaystyle S}
Rad
(
A
1
)
≤
max
i
‖
x
i
‖
∞
⋅
2
log
(
2
n
)
m
{\displaystyle \operatorname {Rad} (A_{1})\leq \max _{i}\|x_{i}\|_{\infty }\cdot {\sqrt {2\log(2n) \over m}}}
次の境界は、集合のラデマッハ複雑性 とその外部 被覆数 (和集合に が含まれるよう な半径のボールの数)を関連付けるものである 。この境界はダドリーによるものである。 [2] : 338
A
{\displaystyle A}
r
{\displaystyle r}
A
{\displaystyle A}
が長さ(ノルム)が最大 のベクトルの集合であると します 。このとき、すべての整数 に対して次の式が成り立ちます 。
A
⊂
R
m
{\displaystyle A\subset \mathbb {R} ^{m}}
c
{\displaystyle c}
M
>
0
{\displaystyle M>0}
Rad
(
A
)
≤
c
⋅
2
−
M
m
+
6
c
m
⋅
∑
i
=
1
M
2
−
i
log
(
N
c
⋅
2
−
i
ext
(
A
)
)
{\displaystyle \operatorname {Rad} (A)\leq {c\cdot 2^{-M} \over {\sqrt {m}}}+{6c \over m}\cdot \sum _{i=1}^{M}2^{-i}{\sqrt {\log \left(N_{c\cdot 2^{-i}}^{\text{ext}}(A)\right)}}}
特に、が の d 次元部分空間 内にある場合 、次のようになります。
A
{\displaystyle A}
R
m
{\displaystyle \mathbb {R} ^{m}}
∀
r
>
0
:
N
r
ext
(
A
)
≤
(
2
c
d
/
r
)
d
{\displaystyle \forall r>0:N_{r}^{\text{ext}}(A)\leq (2c{\sqrt {d}}/r)^{d}}
これを前の境界に代入すると、Rademacher 複雑度の境界は次のようになります。
Rad
(
A
)
≤
6
c
m
⋅
(
d
log
(
2
d
)
+
2
d
)
=
O
(
c
d
log
(
d
)
m
)
{\displaystyle \operatorname {Rad} (A)\leq {6c \over m}\cdot {\bigg (}{\sqrt {d\log(2{\sqrt {d}})}}+2{\sqrt {d}}{\bigg )}=O{\bigg (}{c{\sqrt {d\log(d)}} \over m}{\bigg )}}
ガウス複雑度
ガウス複雑度 は、同様の物理的意味を持つ同様の複雑度であり、 の代わりにランダム変数を使用してラーデマッハ複雑度から得ることができます。 ここで、 は平均がゼロで分散が 1 で ある ガウス i.id ランダム変数、つまり です 。ガウス複雑度とラーデマッハ複雑度は、対数因子まで同等であることが知られています。
g
i
{\displaystyle g_{i}}
σ
i
{\displaystyle \sigma _{i}}
g
i
{\displaystyle g_{i}}
g
i
∼
N
(
0
,
1
)
{\displaystyle g_{i}\sim {\mathcal {N}}(0,1)}
ラデマッハ複雑性とガウス複雑性の同値性
集合が与えられた場合、 [5] が成り立ちます 。
ここで は A のガウス複雑度です。例として、L1 球のラデマッハ複雑度とガウス複雑度を考えてみましょう。ラデマッハ複雑度はちょうど 1 で与えられますが、ガウス複雑度は のオーダーです(これは、 サブガウス ランダム変数の集合の上限の既知の特性を適用することで示せます )。 [5]
A
⊆
R
n
{\displaystyle A\subseteq \mathbb {R} ^{n}}
G
(
A
)
2
log
n
≤
Rad
(
A
)
≤
π
2
G
(
A
)
{\displaystyle {\frac {G(A)}{2{\sqrt {\log {n}}}}}\leq {\text{Rad}}(A)\leq {\sqrt {\frac {\pi }{2}}}G(A)}
G
(
A
)
{\displaystyle G(A)}
log
d
{\displaystyle {\sqrt {\log d}}}
参考文献
^ Balcan, Maria-Florina (2011年11月15日〜17日). 「機械学習理論 - Rademacher Complexity」 (PDF) . 2016年 12月10日 閲覧 。
^ abcdefg 第 26 章、 Shalev-Shwartz, Shai、Ben-David, Shai (2014)。 機械学習を理解する - 理論からアルゴリズムまで 。ケンブリッジ大学出版局 。ISBN 9781107057135 。
^ ab 毛利、メリヤル ;ロスタミザデ、アフシン。タルウォーカー、アメート (2012)。 機械学習の基礎 。米国、マサチューセッツ州:MIT Press。 ISBN 9780262018258 。
^ Bousquet, O. (2004). 統計学習理論入門. 生物サイバネティクス , 3176 (1), 169–207. doi :10.1007/978-3-540-28650-9_8
^ ab Wainwright, Martin (2019). 高次元統計:非漸近的観点。ケンブリッジ、イギリス。pp. 演習 5.5。ISBN 978-1-108-62777-1 . OCLC 1089254580. {{cite book}}: CS1 maint: location missing publisher (link)
Peter L. Bartlett、Shahar Mendelson (2002) Rademacher と Gaussian Complexities: Risk Bounds and Structural Results . Journal of Machine Learning Research 3 463–482 ページ
Giorgio Gnecco、Marcello Sanguineti (2008) Rademacher の計算量による近似誤差の境界 。応用数学科学、第 2 巻、2008 年、第 4 号、153–176 ページ