Definition (Sub-Gaussian Random Variable). A random variable XX is called a sub-gaussian random variable, if it exhibits one of the following equivalent properties: let Ki>0,i∈[5]K_i > 0, i \in [5] be absolute constants.

  1. The tails of XX satisfyPr⁡[∣X∣≥t]≤2exp⁡(−t2K12),∀t≥0. \Pr[|X| \ge t] \le 2 \exp \left(− \frac{t^2}{K_1^2} \right), \quad \forall t \ge 0.
  2. The moments of XX satisfy∥X∥Lp=(E∣X∣p)1/p≤K1p⋅Γ(p/2)p≤K2p,∀p≥1. \Vert X \Vert_{L^p} = ( \mathbb{E}|X|^p )^{1/p} \le K_1 \sqrt[p]{p \cdot \Gamma(p / 2)} \le K_2 \sqrt{p}, \quad \forall p \ge 1.
  3. The MGF of X2X^2 satisfiesEexp⁡(λ2X2)≤exp⁡(K32λ2)∀λ  s.t. , ∣λ∣≤1/K3. \mathbb{E} \exp(\lambda^2 X^2) \le \exp(K_3^2 \lambda^2) \quad \forall \lambda \, \text{ s.t. }, \, |\lambda| \le 1 / K_3 .
  4. The MGF of X2X^2 is bounded at some point, namelyEexp⁡(X2/K42)≤2. \mathbb{E} \exp(X^2 / K_4^2) \le 2.
  5. Moreover, if E[X]=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)∀λ∈R. \mathbb{E} \exp(\lambda X) \le \exp(K_5^2 \lambda^2) \quad \forall \lambda \in \mathbb{R}.

Proof.

1→2:1 \rightarrow 2:

E[∣X∣p]=∫0∞Pr⁡[∣X∣p≥t] dt=∫0∞Pr⁡[∣X∣≥t1/p] dt≤∫0∞2⋅exp⁡(−t2/pK12) dtu=t2/p/K12→t=K1pup/2=∫0∞2⋅K1p⋅(p/2)⋅up/2−1e−u du=K1p⋅p⋅Γ(p/2) \begin{aligned} \mathbb{E} \big[ |X|^p \big] &= \int_0^\infty \Pr\big[ |X|^p \ge t \big] \, d t \\ &= \int_0^\infty \Pr\big[ |X| \ge t^{1 / p} \big] \, d t \\ &\le \int_0^\infty 2 \cdot \exp \left(− \frac{t^{2/ p}}{K_1^2} \right) \, d t & u = t^{2/p} / K_1^2 \rightarrow t = K_1^p u^{p / 2} \\ &= \int_0^\infty 2 \cdot K_1^p \cdot (p / 2) \cdot u^{p / 2 - 1} e^{-u} \, d u \\ &= K_1^p \cdot p \cdot \Gamma(p / 2) \end{aligned}

Applying the Stirling's approximation of Γ(x)≤xx\Gamma(x) \le x^x, and taking the pp root of both sides given

E[∣X∣p]≤K1⋅pp⋅p/2. \mathbb{E} \big[ |X|^p \big] \le K_1 \cdot \sqrt[p]{p} \cdot \sqrt{p / 2}.

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

Eexp⁡(λ2X2)=E(1+∑n∈N+λ2nX2nn!)=1+∑n∈N+E(λ2nX2nn!)≤1+∑n∈N+(λ2K12)n⋅2n⋅Γ(n)n!since E∣X∣2n≤K12n⋅2n⋅Γ(2n/2)≤1+2⋅∑n∈N+(λ2K12)n=1+2K12λ21−K12λ2provided that K12λ2<1=1+K12λ21−K12λ2=exp⁡(3⋅K12λ2),1+x1−x≤e3x, ∀x∈[0,4/5]. \begin{aligned} \mathbb{E} \exp(\lambda^2 X^2) &= \mathbb{E} \left( 1 + \sum_{n \in \mathbb{N}^+} \frac{\lambda^{2n} X^{2n}}{ n! } \right) \\ &= 1 + \sum_{n \in \mathbb{N}^+} \mathbb{E} \left( \frac{\lambda^{2n} X^{2n}}{ n! } \right) \\ &\le 1 + \sum_{n \in \mathbb{N}^+} \frac{ ( \lambda^2 K_1^2 )^n \cdot 2n \cdot \Gamma(n) }{ n! } & \text{since } \mathbb{E} |X|^{2n} \le K_1^{2n} \cdot 2 n \cdot \Gamma(2n / 2) \\ &\le 1 + 2 \cdot \sum_{n \in \mathbb{N}^+} (\lambda^2 K_1^2 )^n \\ &= 1 + \frac{2 K_1^2 \lambda^2}{1 - K_1^2 \lambda^2} & \text{provided that } K_1^2 \lambda^2 < 1 \\ &= \frac{1 + K_1^2 \lambda^2}{1 - K_1^2 \lambda^2} \\ &= \exp(3 \cdot K_1^2 \lambda^2), & \frac{1 + x}{1 - x} \le e^{3x}, \, \forall x \in [0, 4 / 5]. \end{aligned}

3→4:3 \rightarrow 4: Taking λ=1/K3\lambda = 1 / K_3 in property 3 proves property 4.

4→1:4 \rightarrow 1: By Markov's inequality, we have

Pr⁡[∣X∣≥t]=Pr⁡[X2/K42≥t2/K42]=Pr⁡[exp⁡(X2/K42)≥exp⁡(t2/K42)]≤Eexp⁡(X2/K42)exp⁡(t2/K42)≤2exp⁡(−t2/K42). \begin{aligned} \Pr[|X| \ge t] &= \Pr[X^2 / K_4^2 \ge t^2 / K_4^2] \\ &= \Pr[\exp (X^2 / K_4^2 ) \ge \exp( t^2 / K_4^2 )] \\ &\le \frac{ \mathbb{E} \exp(X^2 / K_4^2) }{ \exp( t^2 / K_4^2 ) } \\ &\le 2 \exp( - t^2 / K_4^2 ). \end{aligned}

3→53 \rightarrow 5 (when ∣λ∣⋅K3≤1|\lambda| \cdot K_3 \le 1): Since ex≤x+ex2,∀x∈Re^x \le x + e^{x^2}, \forall x \in R, we have

Eexp⁡(λX)≤E(λX+exp⁡(λ2X2))=E(exp⁡(λ2X2))≤exp⁡(K32λ2). \begin{aligned} \mathbb{E} \exp(\lambda X) &\le \mathbb{E} \left( \lambda X + \exp(\lambda^2 X^2) \right) \\ &= \mathbb{E} \left( \exp(\lambda^2 X^2) \right) \\ &\le \exp(K_3^2 \lambda^2). \end{aligned}

4→54 \rightarrow 5 (when ∣λ∣⋅K3≥1|\lambda| \cdot K_3 \ge 1): Since λX≤(X2/K32+λ2⋅K32)/2\lambda X \le (X^2 / K_3^2 + \lambda^2 \cdot K_3^2) / 2, we have

Eexp⁡(λX)≤Eexp⁡(X2/K32+λ2⋅K322)=exp⁡(λ2⋅K322)⋅Eexp⁡(X2/K322)≤exp⁡(λ2⋅K322)exp⁡(12)exp⁡(λ2⋅K32), \begin{aligned} \mathbb{E} \exp(\lambda X) &\le \mathbb{E} \exp \left( \frac{X^2 / K_3^2 + \lambda^2 \cdot K_3^2}{2} \right) \\ &= \exp \left( \frac{\lambda^2 \cdot K_3^2}{2} \right) \cdot \mathbb{E} \exp \left( \frac{X^2 / K_3^2 }{2} \right) \\ &\le \exp \left( \frac{\lambda^2 \cdot K_3^2}{2} \right) \exp \left( \frac{1}{2} \right) \\ & \exp \left( \lambda^2 \cdot K_3^2 \right), \end{aligned}

where the final inequality follows as λ2K32≥1\lambda^2 K_3^2 \ge 1.

5→1:5 \rightarrow 1: When λ>0\lambda > 0,

Pr⁡[X≥t]=Pr⁡[exp⁡(λX)≥exp⁡(λt)]≤exp⁡(−λt)⋅Eexp⁡(λX)≤exp⁡(−λt+K52λ2)=exp⁡(−t24K52) choosing λ=1/(2K52). \begin{aligned} \Pr[X \ge t] &= \Pr[ \exp(\lambda X) \ge \exp(\lambda t) ] \\ &\le \exp(- \lambda t) \cdot \mathbb{E} \exp(\lambda X) \\ &\le \exp( - \lambda t + K_5^2 \lambda^2 ) \\ &= \exp \left( - \frac{t^2}{4 K_5^2} \right) & \text{ choosing } \lambda = 1 / (2 K_5^2). \end{aligned}

When λ<0\lambda < 0,

Pr⁡[X≤−t]=Pr⁡[exp⁡(λX)≥exp⁡(λt)]≤exp⁡(−λt)⋅Eexp⁡(λX)≤exp⁡(−λt+K52λ2)=exp⁡(−t24K52) choosing λ=1/(2K52). \begin{aligned} \Pr[X \le -t] &= \Pr[ \exp(\lambda X) \ge \exp(\lambda t) ] \\ &\le \exp(- \lambda t) \cdot \mathbb{E} \exp(\lambda X) \\ &\le \exp( - \lambda t + K_5^2 \lambda^2 ) \\ &= \exp \left( - \frac{t^2}{4 K_5^2} \right) & \text{ choosing } \lambda = 1 / (2 K_5^2). \end{aligned}

□\square

Remark. If

Eexp⁡(λX)≤exp⁡(K52λ2)∀λ∈R, \mathbb{E} \exp(\lambda X) \le \exp(K_5^2 \lambda^2) \quad \forall \lambda \in \mathbb{R},

then it holds that

E[X]=0, and E[∣X∣k]∈O(k!), ∀k∈N+. \mathbb{E} [ X ] = 0, \quad \text{ and } \quad \mathbb{E} [ |X|^k ] \in O(k!), \, \forall k \in \mathbb{N}^+.

Proof. Since exp⁡(∣λX∣)≤exp⁡(−λX)+exp⁡(λX)\exp(|\lambda X|) \le \exp(-\lambda X) + \exp(\lambda X), it holds that exp⁡(∣λX∣)\exp(|\lambda X|) is integrable and

Eexp⁡(∣λX∣)≤2⋅exp⁡(K52λ2). \mathbb{E} \exp(|\lambda X|) \le 2 \cdot \exp(K_5^2 \lambda^2).

It follows that Eexp⁡(∣X∣)∈O(1)\mathbb{E} \exp(|X|) \in O(1). Since ∣X∣k≤k!⋅exp⁡(∣X∣)|X|^k \le k! \cdot \exp(|X|), it is integrable and

E∣X∣k∈O(k!). \mathbb{E} |X|^k \in O(k!).

Further, as for each n∈Nn \in \mathbb{N}, it holds that

∑k=0nλkXkk!≤∑k=0n∣λ∣k∣X∣kk!≤exp⁡(∣λX∣), \sum_{k = 0}^n \frac{\lambda^k X^k}{k!} \le \sum_{k = 0}^n \frac{|\lambda|^k |X|^k}{k!} \le \exp(|\lambda X|),

by dominated convergence theorem, we have

1+λE[X]+λ22E[X2]+∑k=3∞λkE[Xk]k!=Eexp⁡(λX)≤exp⁡(K52λ2)=1+K52λ2+∑k=2∞(K5λ)2kk!. 1 + \lambda \mathbb{E}[X] + \frac{\lambda^2}{2} \mathbb{E}[X^2] + \sum_{k = 3}^\infty \frac{\lambda^k \mathbb{E}[X^k] }{k!} = \mathbb{E} \exp(\lambda X) \le \exp(K_5^2 \lambda^2) = 1 + K_5^2 \lambda^2 + \sum_{k = 2}^\infty \frac{ (K_5\lambda)^{2k} }{k!}.

When ∣λ∣<1|\lambda| < 1 and λ→0\lambda \rightarrow 0, clearly

∑k=2∞(K5λ)2kk!∈o(λ2). \sum_{k = 2}^\infty \frac{ (K_5\lambda)^{2k} }{k!} \in o(\lambda^2).

Observe that

λkE[Xk]k!≤λkE[∣X∣k]k!∈O(λk). \frac{\lambda^k \mathbb{E}[X^k] }{k!} \le \frac{\lambda^k \mathbb{E}[|X|^k] }{k!} \in O(\lambda^k).

Therefore,

∣∑k=3∞λkE[Xk]k!∣∈O(λ3)∈o(λ2). \left| \sum_{k = 3}^\infty \frac{\lambda^k \mathbb{E}[X^k] }{k!} \right| \in O(\lambda^3) \in o(\lambda^2).

Hence,

λE[X]+λ22E[X2]+o(λ2)≤K52λ2+o(λ2). \lambda \mathbb{E}[X] + \frac{\lambda^2}{2} \mathbb{E}[X^2] + o(\lambda^2) \le K_5^2 \lambda^2 + o(\lambda^2).

Taking λ>0\lambda > 0 and letting λ→0\lambda \rightarrow 0 gives E[X]≤0\mathbb{E}[X] \le 0.

Taking λ<0\lambda < 0 and letting λ→0\lambda \rightarrow 0 gives E[X]≥0\mathbb{E}[X] \ge 0.

It concludes that E[X]=0\mathbb{E} [ X ] = 0, and therefore E[X2]≤2K52\mathbb{E}[X^2] \le 2 K_5^2.

□\square

Definition. The sub-gaussian norm of XX, denoted ∥X∥ψ2\Vert X \Vert_{\psi_2}, is defined to be the smallest K4K_4 in property 4, i.e.,

∥X∥ψ2=inf⁡{t>0:Eexp⁡(X2/t2)≤2}. \Vert X \Vert_{\psi_2} = \inf \{ t > 0: \mathbb{E} \exp( X^2 / t^2) \le 2 \}.

We first show that Eexp⁡(X2/∥X∥ψ22)≤2\mathbb{E} \exp( X^2 / \Vert X \Vert_{\psi_2}^2) \le 2. Let t1,t2,t3,…t_1, t_2, t_3, \ldots be a decreasing sequence such that ti−∥X∥ψ2≤1/kt_i - \Vert X \Vert_{\psi_2} \le 1 / k. Then by monotone convergence theorem, we have

Eexp⁡(X2/∥X∥ψ22)=Elim⁡i→∞exp⁡(X2/ti2)=lim⁡i→∞Eexp⁡(X2/ti2)≤2. \mathbb{E} \exp( X^2 / \Vert X \Vert_{\psi_2}^2) = \mathbb{E} \lim_{i \rightarrow \infty} \exp( X^2 / t_i^2 ) = \lim_{i \rightarrow \infty} \mathbb{E} \exp( X^2 / t_i^2 ) \le 2.

We restate the previous theorem in terms of the sub-gaussian norm.

  1. Pr⁡(∣X∣≥t)≤2exp⁡(−ct2/∥X∥ψ2),∀t≥0\Pr (|X| \ge t) \le 2 \exp( −c t^2 / \Vert X \Vert_{\psi_2} ), \quad \forall t \ge 0;
  2. ∥X∥Lp≤Cp∥X∥ψ2,∀p≥1\Vert X \Vert_{L_p} \le C \sqrt{p} \Vert X \Vert_{\psi_2}, \quad \forall p \ge 1
  3. Eexp⁡(X2/∥X∥ψ2)≤2\mathbb{E} \exp(X^2 / \Vert X \Vert_{\psi_2} ) \le 2;
  4. if E[X]=0\mathbb{E} [ X ] = 0 then Eexp⁡(λX)≤exp⁡(Cλ2∥X∥ψ2),∀λ∈R\mathbb{E} \exp( \lambda X) \le \exp( C \lambda^2 \Vert X \Vert_{\psi_2} ) , \quad \forall \lambda \in \mathbb{R}.

Here C,c>0C, c > 0 are absolute constants.

Theorem. ∥⋅∥ψ2\Vert \cdot \Vert_{\psi_2} is a norm on the space of sub-gaussian random variables.

Proof. Let XX and YY be sub-gaussian random variables, and a=∥X∥ψ2,b=∥Y∥ψ2a = \Vert X \Vert_{\psi_2}, b = \Vert Y \Vert_{\psi_2}. Then

Eexp⁡(X+Ya+b)2=Eexp⁡(Xa⋅aa+b+Yb⋅ba+b)2≤E[aa+bexp⁡(Xa)2+ba+b⋅exp⁡(Yb)2]≤2. \begin{aligned} \mathbb{E} \exp \left( \frac{X + Y}{a + b} \right)^2 &= \mathbb{E} \exp \left( \frac{X}{a} \cdot \frac{a}{a + b} + \frac{Y}{b} \cdot \frac{b}{a + b} \right)^2 \\ &\le \mathbb{E} \left[ \frac{a}{a + b} \exp \left( \frac{X}{a} \right)^2 + \frac{b}{a + b} \cdot \exp \left( \frac{Y}{b} \right)^2 \right] \le 2. \end{aligned}

Hence,

∥X+Y∥ψ2≤a+b=∥X∥ψ2+∥Y∥ψ2. \Vert X + Y \Vert_{\psi_2} \le a + b = \Vert X \Vert_{\psi_2} + \Vert Y \Vert_{\psi_2}.

□\square

Theorem (Sums of independent sub-gaussians) Let X1,…,XNX_1, \ldots, X_N be independent mean-zero sub-gaussian random variables. Then ∑i∈[N]Xi\sum_{i \in [N]} X_i is also a sub-gaussian random variable, and

∥∑i∈[N]Xi∥ψ22≤C⋅∑i∈[N]∥X∥ψ22 \left\Vert \sum_{i \in [N]} X_i \right\Vert_{\psi_2}^2 \le C \cdot \sum_{i \in [N]} \Vert X \Vert_{\psi_2}^2

where CC is an absolute constant.

Proof.

Eexp⁡(λ∑i∈[N]Xi)=∏i∈[N]Eexp⁡(λXi)≤∏i∈[N]exp⁡(Cλ2∥Xi∥ψ22)=exp⁡(∑i∈[N]Cλ2∥Xi∥ψ22). \begin{aligned} \mathbb{E} \exp \left( \lambda \sum_{i \in [N]} X_i \right) &= \prod_{i \in [N]} \mathbb{E} \exp \left( \lambda X_i \right) \\ &\le \prod_{i \in [N]} \exp \left( C \lambda^2 \Vert X_i \Vert_{\psi_2}^2 \right) = \exp \left( \sum_{i \in [N]} C \lambda^2 \Vert X_i \Vert_{\psi_2}^2 \right). \end{aligned}

Hence, ∑i∈[N]Xi\sum_{i \in [N]} X_i is sug-gaussian, and ∥Xi∥ψ22∈O(∑i∈[N]∥Xi∥ψ22)\Vert X_i \Vert_{\psi_2}^2 \in O( \sum_{i \in [N]} \Vert X_i \Vert_{\psi_2}^2 ).

□\square

Theorem (General Hoeffding Inequality). Let X1,…,XNX_1, \ldots, X_N be independent mean-zero sub-gaussian random variables. Then,

Pr⁡(∣∑i∈[N]Xi∣≥t)≤2exp⁡(−ct2∑i∈[N]∥X∥ψ22),t≥0. \Pr \left( \left| \sum_{i \in [N]} X_i \right| \ge t \right) \le 2 \exp \left( − \frac{c t^2}{\sum_{i \in [N]} \Vert X \Vert_{\psi_2}^2} \right), \quad t \ge 0.

Proof. This can be viewed as a corollary as the previous theorem.

□\square

Theorem (Khintchine’s inequality) Let X1,…,XNX_1, \ldots, X_N be independent sub-gaussian random variables with zero means and unit variances, and let a=(a1,…,aN)∈RNa = (a_1,\ldots,a_N) \in \mathbb{R}^N. Then for every p∈[2,∞)p \in [2, \infty) we have

∥a∥2≤∥∑i∈[N]aiXi∥Lp≤CKp∥a∥2 \Vert a \Vert_2 \le \left\Vert \sum_{i \in [N]} a_i X_i \right\Vert_{L_p} \le C K \sqrt{p} \Vert a \Vert_2

where K=max⁡i∥Xi∥ψ2K = \max_i \Vert X_i \Vert_{\psi_2} and CC is an absolute constant.

Proof. First, since p≥2p \ge 2, we have

∥∑i∈[N]aiXi∥Lp≥∥∑i∈[N]aiXi∥L2=(E[(∑i∈[N]aiXi)2])1/2=(∑i,j∈[N]E[aiajXiXj])1/2=∥a∥2. \begin{aligned} \left\Vert \sum_{i \in [N]} a_i X_i \right\Vert_{L_p} &\ge \left\Vert \sum_{i \in [N]} a_i X_i \right\Vert_{L_2} \\ &= \left( \mathbb{E} \left[ \left( \sum_{i \in [N]} a_i X_i \right)^2 \right] \right)^{1 / 2} \\ &= \left( \sum_{i, j \in [N]} \mathbb{E} \left[ a_i a_j X_i X_j \right] \right)^{1 / 2} = \Vert a \Vert_2. \end{aligned}

Let X≐∑i∈[N]aiXiX \doteq \sum_{i \in [N]} a_i X_i. Then

∥X∥ψ22≤∑i∈[n]∥aiXi∥ψ22=∑i∈[n]ai2⋅∥Xi∥ψ22≤K2⋅∥a∥22. \Vert X \Vert_{\psi_2}^2 \le \sum_{i \in [n]} \Vert a_i X_i \Vert_{\psi_2}^2 = \sum_{i \in [n]} a_i^2 \cdot \Vert X_i \Vert_{\psi_2}^2 \le K^2 \cdot \Vert a \Vert_2^2.

It follows that

E[∣X∣p]∈O(p⋅∥X∥ψ2)⊆O(p⋅K⋅∥a∥2). \mathbb{E} [ |X|^p ] \in O( \sqrt{p} \cdot \Vert X \Vert_{\psi_2} ) \subseteq O( \sqrt{p} \cdot K \cdot \Vert a \Vert_2 ).

□\square

Theorem (Centering). If XX is a sub-gaussian random variable then X−E[X]X − \mathbb{E} [ X ] is sub-gaussian too and

∥X−E[X]∥ψ2≤C∥X∥ψ2, \Vert X − \mathbb{E} [ X ] \Vert_{\psi_2} \le C \Vert X \Vert_{\psi_2},

where CC is an absolute constant.

Proof. Since ∥⋅∥ψ2\Vert \cdot \Vert_{\psi_2} is a norm,

∥X−E[X]∥ψ2≤∥X∥ψ2+∥E[X]∥ψ2. \Vert X − \mathbb{E} [ X ] \Vert_{\psi_2} \le \Vert X \Vert_{\psi_2} + \Vert \mathbb{E} [ X ] \Vert_{\psi_2}.

It remains to bound the second term. By definition, we have

E[exp⁡((E[X])2((1/ln⁡2)1/2⋅E[X])2)]=2, \mathbb{E} \left[ \exp \left( \frac{(\mathbb{E} [ X ])^2}{ \big( (1 / \ln 2)^{1 / 2} \cdot\mathbb{E} [ X ] \big)^2} \right) \right] = 2,

Therefore,

∥E[X]∥ψ2=1ln⁡2⋅E[X]≤1ln⁡2⋅E∣X∣∈O(∥X∥ψ2). \begin{aligned} \Vert \mathbb{E} [ X ] \Vert_{\psi_2} &= \frac{1}{\sqrt{\ln 2}} \cdot \mathbb{E} [ X ] \le \frac{1}{\sqrt{\ln 2}} \cdot \mathbb{E} | X | \in O( \Vert X \Vert_{\psi_2} ). \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.