§E11.18指数型確率不等式

最終更新

Markov の不等式を指数関数へ適用すると、確率変数が大きな値を取る確率を積率母関数によって評価することができる。独立な和では積率母関数が積に分解するため、各項の指数積率を個別に評価すればよい。本記事は、この指数 Markov 法から Chernoff の評価を導き、有界確率変数に対する Hoeffding の補題を証明した後、独立な和に対する片側および両側の Hoeffding の不等式を導く。

1 指数 Markov 法と Chernoff の評価

積率母関数§E11.9 定義 4.1はMZ(t)=E[etZ]M_Z(t)=E[e^{tZ}]によって定義される。上側確率を評価するときはt>0t>0を取り、単調増加関数x↦etxx\mapsto e^{tx}に Markov の不等式を適用する。

定理 1.1 (Chernoff の評価).ZZを実確率変数とし、u>0u>0とする。

  1. t>0t>0かつE[etZ]<∞E[e^{tZ}]<\inftyならば、 P(Z≥u)≤e−tuE[etZ]P(Z\geq u)\leq e^{-tu}E[e^{tZ}] が成り立つ。
  2. t>0t>0かつE[e−tZ]<∞E[e^{-tZ}]<\inftyならば、 P(Z≤−u)≤e−tuE[e−tZ]P(Z\leq-u)\leq e^{-tu}E[e^{-tZ}] が成り立つ。

したがって、右辺を有限にするt>0t>0の範囲で下限を取ることができる。

証明.(1)を示す。t>0t>0のときetZe^{tZ}は非負確率変数であり、

{Z≥u}={etZ≥etu}\{Z\geq u\}=\{e^{tZ}\geq e^{tu}\}

である。Markov の不等式§E11.4 定理 3.1をetZe^{tZ}と閾値etue^{tu}へ適用すると、

P(Z≥u)≤E[etZ]etu=e−tuE[etZ]P(Z\geq u) \leq\frac{E[e^{tZ}]}{e^{tu}} =e^{-tu}E[e^{tZ}]

を得る。(2)を示す。下側の評価は、同じ議論を−Z-Zへ適用して得られる。▨

2 区間に入る確率変数の指数積率

Hoeffding の補題の証明では、区間に入る確率変数の分散が区間幅だけで評価されることを用いる。この分散評価を、指数傾斜した確率測度へ適用する。

補題 2.1.WWを実確率変数とし、ある実数a≤ba\leq bについてP(a≤W≤b)=1P(a\leq W\leq b)=1と仮定する。このとき

Var⁡(W)≤(b−a)24\operatorname{Var}(W)\leq\frac{(b-a)^2}{4}

が成り立つ。

証明.m=E[W]m=E[W]と置く。有界性によりmmは有限であり、a≤m≤ba\leq m\leq bである。点ごとに(W−a)(b−W)≥0(W-a)(b-W)\geq0であるから、

W2≤(a+b)W−abW^2\leq(a+b)W-ab

が成り立つ。期待値を取り、m2m^2を引くと、

Var⁡(W)=E[W2]−m2≤(a+b)m−ab−m2=(b−m)(m−a)=(b−a)24−(m−a+b2)2≤(b−a)24\begin{aligned} \operatorname{Var}(W) &=E[W^2]-m^2\\ &\leq(a+b)m-ab-m^2\\ &=(b-m)(m-a)\\ &=\frac{(b-a)^2}{4}-\left(m-\frac{a+b}{2}\right)^2\\ &\leq\frac{(b-a)^2}{4} \end{aligned}

を得る。▨

2.1 証明方針

有界確率変数XXの対数積率母関数をψ(t)=log⁡E[etX]\psi(t)=\log E[e^{tX}]とする。指数傾斜した確率測度の下でψ′′(t)\psi''(t)はXXの分散に等しい。指数傾斜によってXXの値域は変わらないため、補題 2.1がψ′′\psi''を一様に評価する。中心化からψ(0)=ψ′(0)=0\psi(0)=\psi'(0)=0となるので、二階微分の評価を積分して対数積率母関数の二次上界を得る。

補題 2.2 (Hoeffding の補題).XXを実確率変数とし、ある実数a≤ba\leq bについてP(a≤X≤b)=1P(a\leq X\leq b)=1およびE[X]=0E[X]=0を仮定する。このとき、任意のt∈Rt\in\mathbb Rに対して

E[etX]≤exp⁡ ⁣(t2(b−a)28)E[e^{tX}]\leq\exp\!\left(\frac{t^2(b-a)^2}{8}\right)

が成り立つ。

証明.M(t)=E[etX]M(t)=E[e^{tX}]およびψ(t)=log⁡M(t)\psi(t)=\log M(t)と置く。XXは有界であるため、任意の有界なttの区間上で∣X∣jetX|X|^je^{tX}は定数によって抑えられる。優収束定理§E9.7 定理 3.2を差分商へ適用することにより、MMは二回微分可能であり、

M′(t)=E[XetX],M′′(t)=E[X2etX]M'(t)=E[Xe^{tX}], \qquad M''(t)=E[X^2e^{tX}]

となる。

M(t)>0M(t)>0であるから、

ψ′(t)=M′(t)M(t),ψ′′(t)=M′′(t)M(t)−(M′(t)M(t))2.\psi'(t)=\frac{M'(t)}{M(t)}, \qquad \psi''(t)=\frac{M''(t)}{M(t)}-\left(\frac{M'(t)}{M(t)}\right)^2.

確率測度PtP_tを

Pt(A)=E[etX1A]M(t)(A∈F)P_t(A)=\frac{E[e^{tX}\mathbf1_A]}{M(t)} \qquad(A\in\mathcal F)

によって定める。上の二階微分はPtP_tの下でのXXの分散である。Pt(a≤X≤b)=1P_t(a\leq X\leq b)=1であるため、補題 2.1により

ψ′′(t)≤(b−a)24\psi''(t)\leq\frac{(b-a)^2}{4}

がすべてのt∈Rt\in\mathbb Rについて成り立つ。

c=(b−a)2/4c=(b-a)^2/4およびg(t)=ψ(t)−ct2/2g(t)=\psi(t)-ct^2/2と置くと、g′′(t)≤0g''(t)\leq0であるからggは凹関数である。また、M(0)=1M(0)=1とE[X]=0E[X]=0からg(0)=g′(0)=0g(0)=g'(0)=0となる。凹関数は任意の点で接線以下にあるため、

g(t)≤g(0)+g′(0)t=0g(t)\leq g(0)+g'(0)t=0

である。したがって、

ψ(t)≤ct22=t2(b−a)28.\psi(t)\leq\frac{ct^2}{2}=\frac{t^2(b-a)^2}{8}.

両辺の指数を取ると主張を得る。▨

3 独立な有界確率変数の和

各確率変数が入る区間は同一である必要がない。中心化しても区間幅は変わらないため、各幅の二乗和が指数評価を決める。

3.1 証明方針

各XiX_iをYi=Xi−E[Xi]Y_i=X_i-E[X_i]と中心化する。Hoeffding の補題を各YiY_iへ適用し、独立性によって指数積率を積へ分解する。得られた二次指数上界を Chernoff の評価へ代入し、指数の右辺を最小にするttを選ぶ。下側確率には−Yi-Y_iを用い、両側確率には和事象の評価を用いる。

定理 3.1 (Hoeffding の不等式).X1,…,XnX_1,\ldots,X_nを独立な実確率変数とする。各1≤i≤n1\leq i\leq nについて、実数ai≤bia_i\leq b_iが

P(ai≤Xi≤bi)=1P(a_i\leq X_i\leq b_i)=1

を満たすと仮定する。さらに、

V=∑i=1n(bi−ai)2,Tn=∑i=1n(Xi−E[Xi])V=\sum_{i=1}^n(b_i-a_i)^2, \qquad T_n=\sum_{i=1}^n\bigl(X_i-E[X_i]\bigr)

と置く。

V>0V>0ならば、任意のu>0u>0に対して

P(Tn≥u)≤exp⁡ ⁣(−2u2V),P(Tn≤−u)≤exp⁡ ⁣(−2u2V),P(∣Tn∣≥u)≤2exp⁡ ⁣(−2u2V)\begin{aligned} P(T_n\geq u)&\leq\exp\!\left(-\frac{2u^2}{V}\right),\\ P(T_n\leq-u)&\leq\exp\!\left(-\frac{2u^2}{V}\right),\\ P(|T_n|\geq u)&\leq2\exp\!\left(-\frac{2u^2}{V}\right) \end{aligned}

が成り立つ。V=0V=0ならばTn=0T_n=0がほとんど確実に成り立つ。

証明.μi=E[Xi]\mu_i=E[X_i]およびYi=Xi−μiY_i=X_i-\mu_iと置く。有界性により各期待値は有限であり、ai≤μi≤bia_i\leq\mu_i\leq b_iである。したがって、

ai−μi≤Yi≤bi−μia_i-\mu_i\leq Y_i\leq b_i-\mu_i

がほとんど確実に成り立ち、この区間の幅はbi−aib_i-a_iである。さらにE[Yi]=0E[Y_i]=0であるから、補題 2.2により、任意のt∈Rt\in\mathbb Rに対して

E[etYi]≤exp⁡ ⁣(t2(bi−ai)28)E[e^{tY_i}]\leq\exp\!\left(\frac{t^2(b_i-a_i)^2}{8}\right)

となる。

Y1,…,YnY_1,\ldots,Y_nは独立であり、各etYie^{tY_i}は有限区間上の指数関数として有界である。したがって各因子は可積分であり、有限積∏i=1netYi=etTn\prod_{i=1}^ne^{tY_i}=e^{tT_n}も有界であって可積分である。2≤k≤n2\leq k\leq nに対して、∏i=1k−1etYi\prod_{i=1}^{k-1}e^{tY_i}はY1,…,Yk−1Y_1,\ldots,Y_{k-1}が生成するシグマ加法族に関して可測であり、etYke^{tY_k}と独立である。積の期待値の分解§E11.7 定理 3.1をk=2k=2から順に適用すると、

E[etTn]=∏i=1nE[etYi]≤exp⁡ ⁣(t28∑i=1n(bi−ai)2)=exp⁡ ⁣(t2V8)\begin{aligned} E[e^{tT_n}] &=\prod_{i=1}^nE[e^{tY_i}]\\ &\leq\exp\!\left(\frac{t^2}{8}\sum_{i=1}^n(b_i-a_i)^2\right)\\ &=\exp\!\left(\frac{t^2V}{8}\right) \end{aligned}

を得る。

V>0V>0とu>0u>0を仮定する。t>0t>0に対して定理 1.1を用いると、

P(Tn≥u)≤exp⁡ ⁣(−tu+t2V8).P(T_n\geq u) \leq\exp\!\left(-tu+\frac{t^2V}{8}\right).

右辺の指数はt=4u/Vt=4u/Vで最小となり、最小値は−2u2/V-2u^2/Vである。この指数最小化から上側の評価を得る。−Y1,…,−Yn-Y_1,\ldots,-Y_nも独立であり、それぞれ幅bi−aib_i-a_iの区間に入るため、同じ議論を−Tn-T_nへ適用すると下側の評価を得る。最後に、

{∣Tn∣≥u}={Tn≥u}∪{Tn≤−u}\{|T_n|\geq u\} =\{T_n\geq u\}\cup\{T_n\leq-u\}

へ和事象の評価を適用すると両側の評価を得る。

V=0V=0ならば、すべてのiiについてai=bia_i=b_iである。各XiX_iは定数aia_iにほとんど確実に等しいため、Yi=0Y_i=0、したがってTn=0T_n=0がほとんど確実に成り立つ。▨

4 具体例

例 4.1 (独立な Bernoulli 確率変数の平均).X1,…,XnX_1,\ldots,X_nを独立とし、XiX_iが母数pi∈[0,1]p_i\in[0,1]の Bernoulli 分布に従うとする。同一分布は仮定しない。各XiX_iは区間[0,1][0,1]に入るためV=nV=nであり、pˉ=n−1∑i=1npi\bar p=n^{-1}\sum_{i=1}^np_iと置くと、任意のε>0\varepsilon>0に対して

P ⁣(∣1n∑i=1nXi−pˉ∣≥ε)≤2e−2nε2P\!\left(\left|\frac1n\sum_{i=1}^nX_i-\bar p\right|\geq\varepsilon\right) \leq2e^{-2n\varepsilon^2}

が成り立つ。期待値が同一でなくても、独立性と共通の有界区間だけで指数的な評価を得ることができる。

5 演習

問題 5.1.

  1. 補題 2.1で等号が成り立つための分布を一つ与え、証明中の二つの不等式がともに等号になることを確認せよ。
  2. 補題 2.2の証明で、指数傾斜後もXXが区間[a,b][a,b]にほとんど確実に入ることを、PtP_tの定義から示せ。
  3. XiX_iが区間[−ci,ci][-c_i,c_i]に入り、E[Xi]=0E[X_i]=0を満たす場合について、P(∣∑iXi∣≥u)P(|\sum_iX_i|\geq u)の上界をc1,…,cnc_1,\ldots,c_nによって書き下せ。
  4. Hoeffding の不等式の証明で、各etYie^{tY_i}の可積分性だけでなく積etTne^{tT_n}の可積分性も確認する必要がある理由を、積率母関数の積公式の仮定と対応させて説明せよ。

6 扱った範囲と後続単元との境界

本記事は、指数 Markov 法、Chernoff の評価、Hoeffding の補題、および独立な有界確率変数の有限和に対する片側・両側の Hoeffding の不等式を扱った。各確率変数の区間は異なってよい。Bernstein・Bennett 型不等式、一般の劣 Gauss 確率変数は扱っていない。マルチンゲール差分列に対する Azuma–Hoeffding の不等式は、フィルトレーションとマルチンゲールを扱う「確率過程」へ委ねる。

参考文献

  1. Wassily Hoeffding, Probability Inequalities for Sums of Bounded Random Variables, Journal of the American Statistical Association 58 (1963), no. 301, 13–30.Hoeffding の補題と独立な有界確率変数の和に対する評価を参考にした。
  2. Stéphane Boucheron, Gábor Lugosi, and Pascal Massart, Concentration Inequalities: A Nonasymptotic Theory of Independence, Oxford University Press, Oxford, 2013.指数 Markov 法と集中不等式の定式化を参考にした。

前提記事