Theorem (Sub-exponential properties). Let XX be a random variable. Then the following properties are equivalent; the parameters Ki>0K_i > 0 appearing in these properties differ from each other by at most an absolute constant factor.

  1. The tails of XX satisfyPr⁡{∣X∣≥t}≤2exp⁡(−t/K1),∀t≥0. \Pr \{ |X| \ge t \} \le 2 \exp( − t / K_1 ), \quad \forall t \ge 0.
  2. The moments of XX satisfy∥X∥Lp=(E∣X∣p)1/p≤K1⋅2p⋅Γ(p)p=K2p,∀p≥1. \Vert X \Vert_{L_p} = ( \mathbb{E} |X|^p )^{1/p} \le K_1 \cdot \sqrt[p]{2 p \cdot \Gamma (p)} = K_2 p, \quad \forall p \ge 1.
  3. The MGF of ∣X∣|X| satisfiesEexp⁡(λ∣X∣)≤exp⁡(K3λ),∀λ,  s.t.,  0≤λ≤1/K3. \mathbb{E} \exp( \lambda |X| ) \le \exp( K_3 \lambda ), \quad \forall \lambda, \, \text{ s.t., } \, 0 \le \lambda \le 1 / K_3.
  4. The MGF of ∣X∣|X| is bounded at some point, namelyEexp⁡(∣X∣/K4)≤2. \mathbb{E} \exp( |X| / K_4 ) \le 2.
  5. Moreover, if EX=0\mathbb{E} X = 0 then properties 1 - 4 are also equivalent to the following property: The MGF of XX satisfiesEexp⁡(λX)≤exp⁡(K52λ2),∀λ,  s.t.,  ∣λ∣≤1/K5. \mathbb{E} \exp( \lambda X ) \le \exp( K_5^2 \lambda^2 ), \quad \forall \lambda, \, \text{ s.t., } \, | \lambda | \le 1/K_5.

Proof.

1→2:1 \rightarrow 2:

∥X∥Lpp=E∣X∣p=∫0∞Pr⁡[∣X∣p≥t] dt=∫0∞Pr⁡[∣X∣≥t1/p] dt≤∫0∞2exp⁡(−t1/p/K1) dtu≐t1/p/K1≤∫0∞2K1p⋅p⋅up−1exp⁡(−u) du=2K1p⋅p⋅Γ(p). \begin{aligned} \Vert X \Vert_{L_p}^p &= \mathbb{E}|X|^p \\ &= \int_{0}^\infty \Pr[ |X|^p \ge t ] \, dt \\ &= \int_{0}^\infty \Pr[ |X| \ge t^{1 / p} ] \, dt \\ &\le \int_{0}^\infty 2 \exp(− t^{1 / p} / K_1) \, dt & u \doteq t^{1 / p} / K_1 \\ &\le \int_{0}^\infty 2 K_1^p \cdot p \cdot u^{p - 1} \exp(− u ) \, du \\ &= 2 K_1^p \cdot p \cdot \Gamma (p). \end{aligned}

Noting that Γ(p)≤pp\Gamma(p) \le p^p proves property 2.

2→3:2 \rightarrow 3: By monotone convergence theorem, we have

Eexp⁡(λ∣X∣)=E(∑n∈Nλn∣X∣nn!)=1+∑n∈N+λnE∣X∣nn!≤1+∑n∈N+λn⋅2K1n⋅n⋅Γ(n)n! since  E∣X∣n≤2K1n⋅n⋅Γ(n)≤1+2⋅∑n∈N+λnK1n=1+2λK11−λK1 provided that  λK1<1=1+λK11−λK1≤e3K1λ provided that  λK1∈[0,4/5]. \begin{aligned} \mathbb{E} \exp( \lambda |X|) &= \mathbb{E} \left( \sum_{n \in \mathbb{N}} \frac{\lambda^n |X|^n}{n!}\right) \\ &= 1 + \sum_{n \in \mathbb{N}^+} \frac{\lambda^n \mathbb{E} |X|^n}{n!} \\ &\le 1 + \sum_{n \in \mathbb{N}^+} \frac{\lambda^n \cdot 2 K_1^n \cdot n \cdot \Gamma(n)}{n!} & \text{ since }\, \mathbb{E} |X|^n \le 2 K_1^n \cdot n \cdot \Gamma(n) \\ &\le 1 + 2 \cdot \sum_{n \in \mathbb{N}^+} \lambda^n K_1^n \\ &= 1 + \frac{2 \lambda K_1}{1 - \lambda K_1} & \text{ provided that }\, \lambda K_1 < 1 \\ &= \frac{1 + \lambda K_1}{1 - \lambda K_1} \\ &\le e^{ 3 K_1 \lambda } & \text{ provided that }\, \lambda K_1 \in [0, 4 / 5]. \end{aligned}

3→4:3 \rightarrow 4:

It is trivial to prove property 4 by setting λ=1/K3\lambda = 1 / K_3 in property 3.

4→1:4 \rightarrow 1:

Pr⁡{∣X∣≥t}=Pr⁡{e∣X∣/K4≥et/K4}≤Ee∣X∣/K4et/K4≤2exp⁡(−t/K4),∀t≥0. \begin{aligned} \Pr \{ |X| \ge t \} &= \Pr \{ e^{|X| / K_4} \ge e^{t / K_4} \} \\ &\le \frac{\mathbb{E} e^{|X| / K_4}}{e^{t / K_4}} \\ &\le 2 \exp(−t / K_4), &\forall t \ge 0. \end{aligned}

2→5:2 \rightarrow 5: By monotone convergence theorem, we have

Eexp⁡(λX)=E(∑n∈NλnXnn!)=1+E(λX)+E(∑n∈NλnXnn!)=1+E(∑n≥2λnXnn!)≤1+E(∑n≥2λn∣X∣nn!)=1+∑n≥2λnE∣X∣nn!≤1+∑n≥2λn⋅2K1n⋅n⋅Γ(n)n! since  E∣X∣n≤2K1n⋅n⋅Γ(n)≤1+2⋅∑n≥2λnK1n=1+2λ2K121−λK1 provided that  λK1<1=1+4λ2K12 provided that  λK1≤1/2≤e4K12λ2. \begin{aligned} \mathbb{E} \exp( \lambda X) &= \mathbb{E} \left( \sum_{n \in \mathbb{N}} \frac{\lambda^n X^n}{n!}\right) \\ &= 1 + \mathbb{E} (\lambda X) + \mathbb{E} \left( \sum_{n \in \mathbb{N}} \frac{\lambda^n X^n}{n!}\right) \\ &= 1 + \mathbb{E} \left( \sum_{n \ge 2} \frac{\lambda^n X^n}{n!}\right) \\ &\le 1 + \mathbb{E} \left( \sum_{n \ge 2} \frac{\lambda^n |X|^n}{n!}\right) \\ &= 1 + \sum_{n \ge 2} \frac{\lambda^n \mathbb{E} |X|^n}{n!} \\ &\le 1 + \sum_{n \ge 2} \frac{\lambda^n \cdot 2 K_1^n \cdot n \cdot \Gamma(n)}{n!} & \text{ since }\, \mathbb{E} |X|^n \le 2 K_1^n \cdot n \cdot \Gamma(n) \\ &\le 1 + 2 \cdot \sum_{n \ge 2} \lambda^n K_1^n \\ &= 1 + \frac{2 \lambda^2 K_1^2}{1 - \lambda K_1} & \text{ provided that }\, \lambda K_1 < 1 \\ &= 1 + 4 \lambda^2 K_1^2 & \text{ provided that }\, \lambda K_1 \le 1 / 2 \\ &\le e^{ 4 K_1^2 \lambda^2 }. \end{aligned}

5→2:5 \rightarrow 2: Taking λ=1/K5\lambda = 1 / K_5 and λ=−1/K5\lambda = - 1 / K_5 gives

Eexp⁡(X/K5)≤exp⁡(1),Eexp⁡(−X/K5)≤exp⁡(1), \begin{aligned} \mathbb{E}\exp( X / K_5 ) \le \exp( 1 ), \qquad \mathbb{E}\exp( - X / K_5 ) \le \exp( 1 ), \end{aligned}

Further, for each x∈Rx \in \mathbb{R}, and p>0p > 0, since ∣x∣/p≤1+∣x/p∣≤e∣x∣/p|x| / p \le 1 + |x / p| \le e^{|x| / p}, it holds that

∣x∣p/pp≤e∣x∣≤ex+e−x,⟺∣x/K5∣p/pp≤e∣x/K5∣≤ex/K5+e−x/K5. |x|^p / p^p \le e^{|x|} \le e^x + e^{-x}, \quad \Longleftrightarrow \quad |x / K_5|^p / p^p \le e^{|x / K_5|} \le e^{x / K_5} + e^{- x / K_5}.

It follows that

E[∣X/K5∣p/pp]≤E[eX/K5]+E[e−X/K5]≤2, \mathbb{E} \left[ |X / K_5|^p / p^p \right] \le \mathbb{E} \left[ e^{X / K_5} \right] + \mathbb{E} \left[ e^{- X / K_5} \right] \le 2,

i.e.,

E[∣X∣p]≤2⋅pp⋅K5p. \mathbb{E} \left[ |X|^p \right] \le 2 \cdot p^p \cdot K_5^p.

□\square

A random variable XX that satisfies one of the equivalent properties 1 - 5 in previous theorem is called a sub-exponential random variable.

Definition (Sub-Exponential Norm). The sub-exponential norm of XX, denoted ∥X∥ψ1\Vert X \Vert_{\psi_1}, is defined to be the smallest K3K_3 in property 3. In other words,

∥X∥ψ1=inf⁡{t>0:Eexp⁡(∣X∣/t)≤2}. \Vert X \Vert_{\psi_1} = \inf \{t > 0: \mathbb{E}\exp(|X|/t) \le 2 \}.

Lemma (Sub-gaussian squared is sub-exponential). A random variable XX is sub-gaussian if and only if X2X^2 is sub-exponential. Moreover,

∥X2∥ψ1=∥X∥ψ22. \Vert X^2 \Vert_{\psi_1} = \Vert X \Vert_{\psi_2}^2.

Proof. By definition we have

∥X2∥ψ1≐inf⁡{t>0:Eexp⁡(∣X∣2/t)≤2}.∥X∥ψ2≐inf⁡{t>0:Eexp⁡(∣X∣2/t2)≤2}. \begin{align*} \Vert X^2 \Vert_{\psi_1} &\doteq \inf \{t > 0: \mathbb{E}\exp(|X|^2/t) \le 2 \}. \\ \Vert X \Vert_{\psi_2} &\doteq \inf \{t > 0: \mathbb{E}\exp(|X|^2/t^2) \le 2 \}. \end{align*}

Hence the claim follows.

□\square

Lemma (The product of sub-gaussians is sub-exponential). Let XX and YY be sub-gaussian random variables. Then XYXY is sub-exponential. Moreover,

∥XY∥ψ1≤∥X∥ψ2∥Y∥ψ2. \Vert XY \Vert_{\psi_1} \le \Vert X \Vert_{\psi_2} \Vert Y \Vert_{\psi_2}.

Proof. Assume that ∥X∥ψ2=a\Vert X \Vert_{\psi_2} = a, and ∥Y∥ψ2=b.\Vert Y \Vert_{\psi_2} = b. Then by definition, we have

Eexp⁡(∣X∣2/a2)≤2,Eexp⁡(∣Y∣2/b2)≤2. \begin{aligned} \mathbb{E}\exp( |X|^2 / a^2 ) \le 2, \\ \mathbb{E}\exp( |Y|^2 / b^2 ) \le 2. \end{aligned}

Further

E[exp⁡(∣XY∣/(a⋅b))]≤E[exp⁡(12⋅(X2/a2+Y2/b2))]≤E[exp⁡(12⋅(X2/a2))⋅exp⁡(12⋅(Y2/b2))]≤E[12⋅(exp⁡(12⋅(X2/a2)))2+12⋅(exp⁡(12⋅(Y2/b2)))2]≤12⋅(2+2)=2. \begin{aligned} \mathbb{E} \left[ \exp ( | X Y | / (a \cdot b) ) \right] &\le \mathbb{E} \left[ \exp \bigg( \frac{1}{2} \cdot (X^2 / a^2 + Y^2 / b^2) \bigg) \right] \\ &\le \mathbb{E} \left[ \exp \bigg( \frac{1}{2} \cdot (X^2 / a^2 ) \bigg) \cdot \exp \bigg( \frac{1}{2} \cdot (Y^2 / b^2) \bigg) \right] \\ &\le \mathbb{E} \left[ \frac{1}{2} \cdot \left( \exp \bigg( \frac{1}{2} \cdot (X^2 / a^2 ) \bigg) \right)^2 + \frac{1}{2} \cdot \left( \exp \bigg( \frac{1}{2} \cdot (Y^2 / b^2) \bigg) \right)^2 \right] \\ &\le \frac{1}{2} \cdot \left( 2 + 2 \right) = 2. \end{aligned}

□\square

Theorem (Centering). Given sub-exponential random variables XX, it holds that

∥X−EX∥ψ1≤C∥X∥ψ1. \Vert X − \mathbb{E} X \Vert_{\psi_1} \le C \Vert X \Vert_{\psi_1}.

Proof. Since ∥⋅∥ψ1\Vert \cdot \Vert_{\psi_1} is a norm, it holds that

∥X−EX∥ψ1≤∥X∥ψ1+∥EX∥ψ1, \begin{aligned} \Vert X − \mathbb{E} X \Vert_{\psi_1} &\le \Vert X \Vert_{\psi_1} + \Vert \mathbb{E} X \Vert_{\psi_1}, \end{aligned}

it suffices to bound ∥EX∥ψ1\Vert \mathbb{E} X \Vert_{\psi_1}. Define t=∣EX∣/ln⁡2t = |\mathbb{E} X| / \ln 2. Then

E[exp⁡(∣EX∣/t)]=E[exp⁡(ln⁡2)]=2. \mathbb{E} \left[ \exp \left( |\mathbb{E} X| / t \right) \right] = \mathbb{E} \left[ \exp \left( \ln 2 \right) \right] = 2.

Therefore,

∥EX∥ψ1=1ln⁡2⋅∣EX∣≤1ln⁡2⋅E∣X∣=∥X∥ψ1ln⁡2⋅E[∣X∣/∥X∥ψ1]≤∥X∥ψ1ln⁡2⋅(1+E[∣X∣∥X∥ψ1]+12!⋅E[(∣X∣∥X∥ψ1)2]+13!⋅E[(∣X∣∥X∥ψ1)3]+⋯ )=(a)1ln⁡2⋅E[exp⁡(∣X∣∥X∥ψ1)]=2⋅∥X∥ψ1ln⁡2, \begin{aligned} \Vert \mathbb{E} X \Vert_{\psi_1} &= \frac{1}{\ln 2} \cdot |\mathbb{E} X | \\ &\le \frac{1}{\ln 2} \cdot \mathbb{E} |X| \\ &= \frac{\Vert X \Vert_{\psi_1}}{\ln 2} \cdot \mathbb{E} \left[ |X| / \Vert X \Vert_{\psi_1} \right] \\ &\le \frac{\Vert X \Vert_{\psi_1}}{\ln 2} \cdot \left( 1 + \mathbb{E} \left[ \frac{|X|}{\Vert X \Vert_{\psi_1}} \right] + \frac{1}{2!} \cdot \mathbb{E} \left[ \left( \frac{|X|}{\Vert X \Vert_{\psi_1}} \right)^2 \right] + \frac{1}{3!} \cdot \mathbb{E} \left[ \left( \frac{|X|}{\Vert X \Vert_{\psi_1}} \right)^3 \right] + \cdots \right) \\ &\stackrel{(a)}{=} \frac{1}{\ln 2} \cdot \mathbb{E} \left[ \exp \left( \frac{|X|}{\Vert X \Vert_{\psi_1}} \right) \right] \\ &= 2 \cdot \frac{ \Vert X \Vert_{\psi_1} }{\ln 2}, \end{aligned}

where (a)(a) follows from monotone convergence theorem.

□\square

Theorem (Bernstein’s Inequality). Let X1,…,XnX_1, \ldots, X_n be independent mean-zero sub-exponential random variables. Then, for every t≥0t \ge 0, we have

Pr⁡[∑i∈[n]Xi≥t]≤2exp⁡(−cmin⁡(t2∑i∈[n]∥Xi∥ψ12,tmax⁡i∈[n]∥Xi∥ψ1)), \Pr \left[ \sum_{i \in [n]} X_i \ge t \right] \le 2 \exp \left( − c \min \left( \frac{t^2}{\sum_{i \in [n]} \Vert X_i \Vert_{\psi_1}^2}, \frac{t}{\max_{i \in [n]} \Vert X_i \Vert_{\psi_1}} \right) \right),

where c>0c > 0 is an absolute constant.

Proof. Let λ>0\lambda > 0. Since the KiK_i in each property in Theorem 1 differs by at most an absolute constant, and by property (5) , when $ \lambda \le \frac{ c' }{ \max_{i \in [n]} \Vert X_i \Vert_{\psi_1} } $ for some absolute constant cc, it holds that

Pr⁡[∑i∈[n]Xi≥t]≤∏i∈[n]Eexp⁡(λ⋅Xi)exp⁡(λ⋅t)≤exp⁡(∑i∈[n]λ2⋅c∗⋅∥Xi∥ψ12−λ⋅t), \begin{aligned} \Pr \left[ \sum_{i \in [n]} X_i \ge t \right] \le \frac{ \prod_{i \in [n]} \mathbb{E} \exp \left( \lambda \cdot X_i \right) }{ \exp \left( \lambda \cdot t \right) } \le \exp \left( \sum_{i \in [n]} \lambda^2 \cdot c^* \cdot \Vert X_i \Vert_{\psi_1}^2 - \lambda \cdot t \right), \end{aligned}

where c∗=1/c′c^* = 1 / c'. The RHS is minimized for

λ=min⁡{t2⋅c∗⋅∑i∈[n]∥Xi∥ψ12,c′max⁡i∈[n]∥Xi∥ψ1}, \lambda = \min \left\lbrace \frac{t}{ 2 \cdot c^* \cdot \sum_{i \in [n]} \Vert X_i \Vert_{\psi_1}^2 } ,\quad \frac{ c' }{ \max_{i \in [n]} \Vert X_i \Vert_{\psi_1} } \right\rbrace,

which gives

Pr⁡[∑i∈[n]Xi≥t]≤exp⁡(∑i∈[n]λ2⋅c∗⋅∥Xi∥ψ12−λ⋅t). \begin{aligned} \Pr \left[ \sum_{i \in [n]} X_i \ge t \right] \le \exp \left( \sum_{i \in [n]} \lambda^2 \cdot c^* \cdot \Vert X_i \Vert_{\psi_1}^2 - \lambda \cdot t \right). \end{aligned}

When λ=t2⋅c∗⋅∑i∈[n]∥Xi∥ψ12≤c′max⁡i∈[n]∥Xi∥ψ1,\lambda = \frac{t}{ 2 \cdot c^* \cdot \sum_{i \in [n]} \Vert X_i \Vert_{\psi_1}^2 } \le \frac{ c' }{ \max_{i \in [n]} \Vert X_i \Vert_{\psi_1} }, we have

∑i∈[n]λ2⋅c∗⋅∥Xi∥ψ12−λ⋅t=t24⋅c∗⋅∑i∈[n]∥Xi∥ψ12−t22⋅c∗⋅∑i∈[n]∥Xi∥ψ12=−t24⋅c∗⋅∑i∈[n]∥Xi∥ψ12. \sum_{i \in [n]} \lambda^2 \cdot c^* \cdot \Vert X_i \Vert_{\psi_1}^2 - \lambda \cdot t = \frac{t^2}{ 4 \cdot c^* \cdot \sum_{i \in [n]} \Vert X_i \Vert_{\psi_1}^2 } - \frac{t^2}{ 2 \cdot c^* \cdot \sum_{i \in [n]} \Vert X_i \Vert_{\psi_1}^2 } = - \frac{t^2}{ 4 \cdot c^* \cdot \sum_{i \in [n]} \Vert X_i \Vert_{\psi_1}^2 }.

When λ=c′max⁡i∈[n]∥Xi∥ψ1≤t2⋅c∗⋅∑i∈[n]∥Xi∥ψ12,\lambda = \frac{ c' }{ \max_{i \in [n]} \Vert X_i \Vert_{\psi_1} } \le \frac{t}{ 2 \cdot c^* \cdot \sum_{i \in [n]} \Vert X_i \Vert_{\psi_1}^2 }, we have

∑i∈[n]λ2⋅c∗⋅∥Xi∥ψ12−λ⋅t=(c′max⁡i∈[n]∥Xi∥ψ1)2⋅∑i∈[n]c∗⋅∥Xi∥ψ12−(c′max⁡i∈[n]∥Xi∥ψ1)t≤(c′max⁡i∈[n]∥Xi∥ψ1)⋅t2⋅c∗⋅∑i∈[n]∥Xi∥ψ12⋅∑i∈[n]c∗⋅∥Xi∥ψ12−(c′max⁡i∈[n]∥Xi∥ψ1)t=−c′t2⋅max⁡i∈[n]∥Xi∥ψ1. \begin{aligned} \sum_{i \in [n]} \lambda^2 \cdot c^* \cdot \Vert X_i \Vert_{\psi_1}^2 - \lambda \cdot t &= \left( \frac{ c' }{ \max_{i \in [n]} \Vert X_i \Vert_{\psi_1} } \right)^2 \cdot \sum_{i \in [n]} c^* \cdot \Vert X_i \Vert_{\psi_1}^2 - \left( \frac{ c' }{ \max_{i \in [n]} \Vert X_i \Vert_{\psi_1} } \right) t \\ &\le \left( \frac{ c' }{ \max_{i \in [n]} \Vert X_i \Vert_{\psi_1} } \right) \cdot \frac{t}{ 2 \cdot c^* \cdot \sum_{i \in [n]} \Vert X_i \Vert_{\psi_1}^2 } \cdot \sum_{i \in [n]} c^* \cdot \Vert X_i \Vert_{\psi_1}^2 - \left( \frac{ c' }{ \max_{i \in [n]} \Vert X_i \Vert_{\psi_1} } \right) t \\ &= - \frac{ c' t }{ 2 \cdot \max_{i \in [n]} \Vert X_i \Vert_{\psi_1} }. \end{aligned}

It concludes that

Pr⁡[∑i∈[n]Xi≥t]≤exp⁡(max⁡{−t24⋅c∗⋅∑i∈[n]∥Xi∥ψ12, −c′t2⋅max⁡i∈[n]∥Xi∥ψ1})=exp⁡(−min⁡{t24⋅c∗⋅∑i∈[n]∥Xi∥ψ12, c′t2⋅max⁡i∈[n]∥Xi∥ψ1}). \begin{aligned} \Pr \left[ \sum_{i \in [n]} X_i \ge t \right] &\le \exp \left( \max \left\lbrace - \frac{t^2}{ 4 \cdot c^* \cdot \sum_{i \in [n]} \Vert X_i \Vert_{\psi_1}^2 }, \, - \frac{ c' t }{ 2 \cdot \max_{i \in [n]} \Vert X_i \Vert_{\psi_1} } \right\rbrace \right) \\ &= \exp \left( - \min \left\lbrace \frac{t^2}{ 4 \cdot c^* \cdot \sum_{i \in [n]} \Vert X_i \Vert_{\psi_1}^2 }, \, \frac{ c' t }{ 2 \cdot \max_{i \in [n]} \Vert X_i \Vert_{\psi_1} } \right\rbrace \right). \end{aligned}

□\square

Reference

[1] R. Vershynin, "High-Dimensional Probability: An Introduction with Applications in Data Science." in Cambridge Series in Statistical and Probabilistic Mathematics. Cambridge: Cambridge University Press, 2018. doi: 10.1017/9781108231596.