Lecture 14

Concentration inequalities

Course edition
Fall 2026
Last updated
October 19, 2026

Part V · High dimensions: tools and phenomena · Lecture 1 of 4

Student PDFReport an issue

A guided route

Lecture map

High and mid topics are both examined. Follow the topic links, and open optional branches when they help.

These notes support STAT 709 at the University of Wisconsin–Madison. Use them alongside the cited textbooks; corrections and questions are welcome.

High priorityOverview

  • Limits. A limit theorem describes a distribution as the sample size grows.

  • Finite samples. A concentration inequality instead controls a probability at the sample size in hand.

  • Tools. The central tools are the exponential Markov bound, the subgaussian moment-generating-function bound, and Bernstein’s two deviation scales.

  • Use. We track which assumptions justify each bound and how to turn a tail probability into a confidence radius.

This begins Part V. The part overview connects these lectures and includes a final self-check.

Recap from earlier lecturesFrom Lectures 1, 4, 5 and 10: moments and tails · 1 min

Reminder only: this block introduces no new exam requirements.

Learning goals

High priorityIntroduction

In this lecture, we will introduce some basics for understanding high-dimensional phenomena.

What are concentration inequalities?

From the LLN, we know that an empirical mean \(\hat \mu\)μ̂\hat \mu “concentrates” around the population mean \(\mu\)μ\mu. This is characterized by the asymptotic notion of convergence in probability. When we zoom in on the difference \(\hat \mu - \mu\)μ̂−μ\hat \mu - \mu by multiplying it with a factor of \(\sqrt{n}\)n\sqrt{n}, then, under the i.i.d. finite-variance CLT assumptions, \(\sqrt{n}(\hat\mu-\mu)\)n(μ̂−μ)\sqrt{n}(\hat\mu-\mu) has a normal distribution asymptotically.

Concentration inequalities provide non-asymptotic bounds on the difference as an alternative way of characterizing the concentration. Usually we will see statements like “With very high probability, a certain event (e.g. inequality) happens.” Often, the exception probability decays exponentially in the number of variables. Compared with classical asymptotic theory, concentration inequalities often give us refined quantitative control of extreme/exception probabilities.

In this lecture, we use \(C\)CC and \(c\)cc for positive absolute constants, which may change from line to line. An unsubscripted vector norm \(\|a\|\)∥a∥\|a\| means \(\|a\|_2\)∥a∥2\|a\|_2.

1. Motivation for concentration inequalities

Mid priorityMotivation for concentration inequalities

Recall the central limit theorem: if \(X_1, X_2, \ldots, X_n \stackrel{\iid}{\sim} X\)X1,X2,…,Xn∼i.i.d.XX_1, X_2, \ldots, X_n \stackrel{\mathrm{i.i.d.}}{\sim} X satisfy \(\mathbb{E} X = 0\)𝔼X=0\mathbb{E} X = 0 and \(\mathbb{E} X^2 = 1\)𝔼X2=1\mathbb{E} X^2 = 1, then \[\frac{1}{\sqrt{n}}\sum_{i=1}^n X_i \stackrel{d}{\longrightarrow} \mathcal{N}(0, 1).\]1n∑i=1nXi→d𝒩(0,1).\frac{1}{\sqrt{n}}\sum_{i=1}^n X_i \stackrel{d}{\longrightarrow} \mathcal{N}(0, 1). Recall the Gaussian tail probability bounds from Lecture 1.

Proposition 1Gaussian tail probability

Let \(G \sim \mathcal{N}(0, 1)\)G∼𝒩(0,1)G \sim \mathcal{N}(0, 1), then for any \(t > 0\)t>0t > 0, \[\left(\frac{1}{t} - \frac{1}{t^3}\right) \frac{1}{\sqrt{2\pi}} e^{-t^2/2} \le \mathbb{P}(G \ge t) \le \frac{1}{t} \frac{1}{\sqrt{2\pi}} e^{-t^2/2}.\](1t−1t3)12πe−t2/2≤ℙ(G≥t)≤1t12πe−t2/2.\left(\frac{1}{t} - \frac{1}{t^3}\right) \frac{1}{\sqrt{2\pi}} e^{-t^2/2} \le \mathbb{P}(G \ge t) \le \frac{1}{t} \frac{1}{\sqrt{2\pi}} e^{-t^2/2}.

Let \(\Phi\)Φ\Phi and \(\phi\)ϕ\phi be the standard Gaussian CDF and density, respectively. From the above proposition, if \(t\to\infty\)t→∞t\to\infty, then \[1-\Phi(t)=(1+o(1))\frac{\phi(t)}{t}.\]1−Φ(t)=(1+o(1))ϕ(t)t.1-\Phi(t)=(1+o(1))\frac{\phi(t)}{t}. If we pretend that the following approximation is valid, \[\mathbb{P}\big(\frac{1}{n}\sum_{i=1}^n X_i > t \big) = \mathbb{P}\big(\frac{1}{\sqrt{n}}\sum_{i=1}^n X_i > \sqrt{n}\,t \big) \approx \mathbb{P}(G > \sqrt{n}t),\]ℙ(1n∑i=1nXi>t)=ℙ(1n∑i=1nXi>nt)≈ℙ(G>nt),\mathbb{P}\big(\frac{1}{n}\sum_{i=1}^n X_i > t \big) = \mathbb{P}\big(\frac{1}{\sqrt{n}}\sum_{i=1}^n X_i > \sqrt{n}\,t \big) \approx \mathbb{P}(G > \sqrt{n}t), then the tail probability of empirical mean of \(n\)nn i.i.d. random variables is \(O\left(e^{-{nt^2/2}}\right)\)O(e−nt2/2)O\left(e^{-{nt^2/2}}\right). The CLT alone does not justify this approximation: its convergence at fixed thresholds gives no relative-error control at the moving threshold \(\sqrt n\,t\)nt\sqrt n\,t. Exponentially decaying bounds do hold under stronger tail assumptions, which we develop below.

Proposition 2Berry–Esseen inequality

Suppose that \(X_1,X_2,\ldots, X_n\)X1,X2,…,XnX_1,X_2,\ldots, X_n are i.i.d. random variables with zero mean, variance one and finite absolute third moment. Then, for all \(t \in \R\)t∈ℝt \in \mathbb{R}, \[\left| \mathbb{P}\left(\frac{1}{\sqrt{n}}\sum_{i=1}^n X_i \ge t\right) - \mathbb{P}(G \ge t)\right| \le \frac{C\mathbb{E}|X_1|^3}{\sqrt{n}},\]|ℙ(1n∑i=1nXi≥t)−ℙ(G≥t)|≤C𝔼|X1|3n,\left| \mathbb{P}\left(\frac{1}{\sqrt{n}}\sum_{i=1}^n X_i \ge t\right) - \mathbb{P}(G \ge t)\right| \le \frac{C\mathbb{E}|X_1|^3}{\sqrt{n}}, where \(C\)CC is a positive constant.

The error term is \(O\left(\frac{1}{\sqrt{n}}\right)\)O(1n)O\left(\frac{1}{\sqrt{n}}\right), much larger than \(O\left(e^{-nt^2/2}\right)\)O(e−nt2/2)O\left(e^{-nt^2/2}\right) for a fixed \(t>0\)t>0t>0 as \(n\to\infty\)n→∞n\to\infty. Moreover, the error term cannot be made smaller in general due to the next simple example.

Thus Berry–Esseen gives a valid additive bound for \(\PR(n^{-1}\sum_iX_i\ge t)\)ℙ(n−1∑iXi≥t)\mathbb{P}(n^{-1}\sum_iX_i\ge t), but its \(n^{-1/2}\)n−1/2n^{-1/2} error can be much larger than the exponentially small probability we want to bound.

Example 1

Let \(X_1, X_2, \ldots, X_n \stackrel{\iid}{\sim} \mathrm{Bernoulli}(\frac{1}{2})\)X1,X2,…,Xn∼i.i.d.Bernoulli(12)X_1, X_2, \ldots, X_n \stackrel{\mathrm{i.i.d.}}{\sim} \mathrm{Bernoulli}(\frac{1}{2}), where \(n\)nn is an even integer, then \[\mathbb{P}\left(\sum_{i=1}^{n}X_i = \frac{n}{2}\right) = 2^{-n} \binom{n}{n/2} = 2^{-n} \frac{n!}{[(n/2)!]^2}.\]ℙ(∑i=1nXi=n2)=2−n(nn/2)=2−nn![(n/2)!]2.\mathbb{P}\left(\sum_{i=1}^{n}X_i = \frac{n}{2}\right) = 2^{-n} \binom{n}{n/2} = 2^{-n} \frac{n!}{[(n/2)!]^2}. Applying Stirling’s approximation: \[n! = (1+ o(1)) \sqrt{2\pi n} \left(\frac{n}{e}\right)^n,\]n!=(1+o(1))2πn(ne)n,n! = (1+ o(1)) \sqrt{2\pi n} \left(\frac{n}{e}\right)^n, we have \[\mathbb{P}\left(\sum_{i=1}^{n}X_i = \frac{n}{2}\right) = (1 + o(1))\sqrt{\frac{2}{\pi}}\frac{1}{\sqrt{n}}\]ℙ(∑i=1nXi=n2)=(1+o(1))2π1n\mathbb{P}\left(\sum_{i=1}^{n}X_i = \frac{n}{2}\right) = (1 + o(1))\sqrt{\frac{2}{\pi}}\frac{1}{\sqrt{n}} such that \[\mathbb{P}\left(\frac{1}{\sqrt{n}}\sum_{i=1}^{n}(X_i-\mathbb E X_i) = 0\right) = (1+o(1))\sqrt{\frac{2}{\pi}}\frac{1}{\sqrt{n}}.\]ℙ(1n∑i=1n(Xi−𝔼Xi)=0)=(1+o(1))2π1n.\mathbb{P}\left(\frac{1}{\sqrt{n}}\sum_{i=1}^{n}(X_i-\mathbb E X_i) = 0\right) = (1+o(1))\sqrt{\frac{2}{\pi}}\frac{1}{\sqrt{n}}. Note that \(\mathbb{P}(G = 0) = 0\)ℙ(G=0)=0\mathbb{P}(G = 0) = 0. In this example, the jump has order \(n^{-1/2}\)n−1/2n^{-1/2}, so any continuous CDF differs from this discrete CDF on at least one side of the jump by half that size. After standardizing the Bernoulli variance, Berry–Esseen gives a matching upper order.

Mid priorityPractice: A continuous CDF cannot match a jump

See Exercise 1 in Section 7 for the full statement and complete solution.

2. Subgaussian variables and Hoeffding’s inequality

High prioritySubgaussian variables and Hoeffding’s inequality

To achieve our goal of bounding tail probabilities, we will need to build the machinery from scratch. First, we introduce an important definition.

Theorem 1Subgaussian random variables

Let \(X\)XX be a random variable with \(\mathbb{E}X = 0\)𝔼X=0\mathbb{E}X = 0. Consider three conditions.

  • (i) \(\mathbb{P}(|X| \ge t) \le 2e^{-t^2/K_1^2},~\forall\;t \ge 0\)ℙ(|X|≥t)≤2e−t2/K12,∀t≥0\mathbb{P}(|X| \ge t) \le 2e^{-t^2/K_1^2},~\forall\;t \ge 0.

  • (ii) \(\sup_{p\ge 1}\frac{1}{\sqrt{p}}(\mathbb{E}|X|^p)^{1/p} \le K_2\)supp≥11p(𝔼|X|p)1/p≤K2\sup_{p\ge 1}\frac{1}{\sqrt{p}}(\mathbb{E}|X|^p)^{1/p} \le K_2.

  • (iii) \(\mathbb{E}[e^{\lambda X}] \le e^{K_3^2 \lambda^2},~\forall\;\lambda \in \mathbb{R}\)𝔼[eλX]≤eK32λ2,∀λ∈ℝ\mathbb{E}[e^{\lambda X}] \le e^{K_3^2 \lambda^2},~\forall\;\lambda \in \mathbb{R}.1

The three conditions are equivalent up to absolute constants: if one holds with constant \(K_j\)KjK_j, each other condition holds with a suitable \(K_i\le CK_j\)Ki≤CKjK_i\le CK_j. We say \(X\)XX is a subgaussian random variable if any of the three conditions holds.

Proof. (i)\(\Rightarrow\)⇒\Rightarrow(ii): Without loss of generality, assume \(K_1 = 1\)K1=1K_1 = 1. Otherwise set \(X := X / K_1\)X:=X/K1X := X / K_1. Observe \[\begin{aligned} \mathbb{E}|X|^p &= \int_0^{\infty}\mathbb{P}(|X| \ge t) pt^{p-1}dt \qquad \text{(see Lecture 4, or use Fubini's theorem)}\\ &\le \int_0^\infty 2p t^{p-1}e^{-t^2}dt\\ & = p \int_0^\infty u^{(p-2)/2} e^{-u}\; du \qquad \text{(change of variables $u = t^2$)} \\ &= p\, \Gamma(p/2) \qquad \text{(Gamma function)} \\ &\le 3p(\frac{p}{2})^{p/2}\qquad \text{(by property of Gamma functions).} \end{aligned}\]𝔼|X|p=∫0∞ℙ(|X|≥t)ptp−1dt(see Lecture 4, or use Fubini’s theorem)≤∫0∞2ptp−1e−t2dt=p∫0∞u(p−2)/2e−udu(change of variables u=t2)=pΓ(p/2)(Gamma function)≤3p(p2)p/2(by property of Gamma functions).\begin{aligned} \mathbb{E}|X|^p &= \int_0^{\infty}\mathbb{P}(|X| \ge t) pt^{p-1}dt \qquad \text{(see Lecture 4, or use Fubini's theorem)}\\ &\le \int_0^\infty 2p t^{p-1}e^{-t^2}dt\\ & = p \int_0^\infty u^{(p-2)/2} e^{-u}\; du \qquad \text{(change of variables $u = t^2$)} \\ &= p\, \Gamma(p/2) \qquad \text{(Gamma function)} \\ &\le 3p(\frac{p}{2})^{p/2}\qquad \text{(by property of Gamma functions).} \end{aligned} (ii)\(\Rightarrow\)⇒\Rightarrow(iii): Let \(\widetilde X\)X̃\widetilde X be an independent copy, and \(Z=X-\widetilde X\)Z=X−X̃Z=X-\widetilde X. Exchanging the copies shows that \(Z\)ZZ is symmetric. Jensen and independence give \(\E e^{\lambda X}\le\E e^{\lambda Z}\)𝔼eλX≤𝔼eλZ\mathbb{E}e^{\lambda X}\le\mathbb{E}e^{\lambda Z}. From (ii) and Minkowski, \((\E|Z|^{2m})^{1/(2m)}\le2K_2\sqrt{2m}\)(𝔼|Z|2m)1/(2m)≤2K22m(\mathbb{E}|Z|^{2m})^{1/(2m)}\le2K_2\sqrt{2m}. Symmetry, the nonnegative series for \(\cosh\)cosh\cosh, and Tonelli give \[\E e^{\lambda Z}=\E\cosh(\lambda Z)=1+\sum_{m\ge1}\frac{\lambda^{2m}\E Z^{2m}}{(2m)!}.\]𝔼eλZ=𝔼cosh(λZ)=1+∑m≥1λ2m𝔼Z2m(2m)!.\mathbb{E}e^{\lambda Z}=\mathbb{E}\cosh(\lambda Z)=1+\sum_{m\ge1}\frac{\lambda^{2m}\mathbb{E}Z^{2m}}{(2m)!}. Using \((2m)!\ge(2m/e)^{2m}\)(2m)!≥(2m/e)2m(2m)!\ge(2m/e)^{2m} and \(m!\le m^m\)m!≤mmm!\le m^m, the \(m\)mm-th coefficient is at most \((2e^2K_2^2\lambda^2)^m/m!\)(2e2K22λ2)m/m!(2e^2K_2^2\lambda^2)^m/m!. Therefore \(\E e^{\lambda X}\le\exp(2e^2K_2^2\lambda^2)\)𝔼eλX≤exp(2e2K22λ2)\mathbb{E}e^{\lambda X}\le\exp(2e^2K_2^2\lambda^2). The finite bound also justifies exponential integrability; Tonelli applies to nonnegative terms.

(iii)\(\Rightarrow\)⇒\Rightarrow(i): Without loss of generality, assume \(K_3 = 1\)K3=1K_3 = 1. For any \(\lambda > 0\)λ>0\lambda > 0, \(t > 0\)t>0t > 0, \[\begin{aligned} \mathbb{P}(X \ge t) &\le \mathbb{P}(e^{\lambda X} \ge e^{\lambda t}) \\ &\le e^{-\lambda t} \mathbb{E} e^{\lambda X} \qquad \text{(by Markov's inequality)}\\ &\le e^{-\lambda t + \lambda^2}\\ &\le e^{-t^2/4} \qquad \text{(choose $\lambda = t/2$ to minimize $-\lambda t + \lambda^2$).} \end{aligned}\]ℙ(X≥t)≤ℙ(eλX≥eλt)≤e−λt𝔼eλX(by Markov’s inequality)≤e−λt+λ2≤e−t2/4(choose λ=t/2 to minimize −λt+λ2).\begin{aligned} \mathbb{P}(X \ge t) &\le \mathbb{P}(e^{\lambda X} \ge e^{\lambda t}) \\ &\le e^{-\lambda t} \mathbb{E} e^{\lambda X} \qquad \text{(by Markov's inequality)}\\ &\le e^{-\lambda t + \lambda^2}\\ &\le e^{-t^2/4} \qquad \text{(choose $\lambda = t/2$ to minimize $-\lambda t + \lambda^2$).} \end{aligned} Repeating the same for \(-X\)−X-X, we have \[\mathbb{P}(-X \ge t) \le e^{-t^2/4}.\]ℙ(−X≥t)≤e−t2/4.\mathbb{P}(-X \ge t) \le e^{-t^2/4}. Therefore, \[\mathbb{P}(|X| \ge t) \le 2e^{-t^2/4},\qquad \forall\;t > 0.\]ℙ(|X|≥t)≤2e−t2/4,∀t>0.\mathbb{P}(|X| \ge t) \le 2e^{-t^2/4},\qquad \forall\;t > 0. ◻

  • Moment growth. Taking \(p\)ppth roots in (i)\(\Rightarrow\)⇒\Rightarrow(ii) gives \(\|X\|_{L^p}\le C\sqrt p\)∥X∥Lp≤Cp\|X\|_{L^p}\le C\sqrt p, since \((3p)^{1/p}\)(3p)1/p(3p)^{1/p} is uniformly bounded for \(p\ge1\)p≥1p\ge1.

  • Definition. The smallest \(K_2\)K2K_2 in condition (ii) is called the subgaussian norm of \(X\)XX, denoted by \(\|X\|_{\psi_2}\)∥X∥ψ2\|X\|_{\psi_2}. By Minkowski inequality, it is easy to check that it satisfies the requirements of a norm (see Exercise 5).

  • Interpretation. Like standard deviation, it measures the spread of a centered variable, but finite subgaussian norm additionally controls the tails.

  • Convention. We use the moment definition for variables with any mean; only the uncentered MGF bound (iii) requires mean zero. Some textbooks use an exponential-moment definition of \(\|X\|_{\psi_2}\)∥X∥ψ2\|X\|_{\psi_2}; these conventions agree up to absolute constants.

Example 2

If \(X \sim \mathcal{N}(0, \sigma^2)\)X∼𝒩(0,σ2)X \sim \mathcal{N}(0, \sigma^2), then \(\|X\|_{\psi_2} \le C \sigma\)∥X∥ψ2≤Cσ\|X\|_{\psi_2} \le C \sigma.

Example 3

If \(|X| \le M\)|X|≤M|X| \le M a.s., then \(\|X\|_{\psi_2} \le CM\)∥X∥ψ2≤CM\|X\|_{\psi_2} \le CM.

For a sum of independent subgaussian random variables, we have the desired Gaussian tail bound.

Theorem 2Hoeffding’s inequality

Let \(X_1, X_2, \ldots, X_n\)X1,X2,…,XnX_1, X_2, \ldots, X_n be independent random variables, such that \(\mathbb{E} X_i = 0\)𝔼Xi=0\mathbb{E} X_i = 0 and each \(X_i\)XiX_i is subgaussian. Let \(a = (a_1,\ldots, a_n)\)a=(a1,…,an)a = (a_1,\ldots, a_n) be a nonrandom vector. Denote \(K := \max_{1\le i \le n} \|X_i\|_{\psi_2}\)K:=max1≤i≤n∥Xi∥ψ2K := \max_{1\le i \le n} \|X_i\|_{\psi_2}. Then

\[\mathbb{P}\left(\left|\sum_{i=1}^{n}a_i X_i\right|\ge t\right) \le 2\exp\left(-\frac{ct^2}{K^2\|a\|^2}\right)\]ℙ(|∑i=1naiXi|≥t)≤2exp(−ct2K2∥a∥2)\mathbb{P}\left(\left|\sum_{i=1}^{n}a_i X_i\right|\ge t\right) \le 2\exp\left(-\frac{ct^2}{K^2\|a\|^2}\right)

(1)

holds for \(t>0\)t>0t>0 and a positive absolute constant \(c\)cc. If \(K\|a\|_2=0\)K∥a∥2=0K\|a\|_2=0, the sum is zero almost surely, so this degenerate case is handled directly.

In particular, if we choose \(a_i = 1/n\)ai=1/na_i = 1/n for all \(i\)ii then Hoeffding’s inequality says that we have the Gaussian tail bound.

Apart from the constant, this bound cannot be improved. In the case \(X_i \stackrel{\iid}{\sim} N(0,1)\)Xi∼i.i.d.N(0,1)X_i \stackrel{\mathrm{i.i.d.}}{\sim} N(0,1), we have \(\sum_{i=1}^{n}a_i X_i \sim N(0, \| a \|^2)\)∑i=1naiXi∼N(0,∥a∥2)\sum_{i=1}^{n}a_i X_i \sim N(0, \| a \|^2), so by Proposition 1 its tail has exponent \(-t^2/(2\|a\|^2)\)−t2/(2∥a∥2)-t^2/(2\|a\|^2) up to the stated polynomial factor.

Proof. Of Theorem 2.

\[\begin{aligned} \mathbb{E} \exp\left(\lambda \sum_{i=1}^{n}a_i X_i \right) &= \prod_{i=1}^{n} \mathbb{E} \exp(\lambda a_i X_i) \qquad \text{(by independence)}\\ &\le \prod_{i=1}^{n} \exp(C\lambda^2a_i^2\|X_i\|_{\psi_2}^2)\qquad\text{(by Theorem~1~(iii))}\\ &\le \exp(C\lambda^2 \sum_{i=1}^{n} a_i^2 K^2)\\ &= \exp(C\lambda^2 \|a\|^2 K^2). \end{aligned}\]𝔼exp(λ∑i=1naiXi)=∏i=1n𝔼exp(λaiXi)(by independence)≤∏i=1nexp(Cλ2ai2∥Xi∥ψ22)(by Theorem 1 (iii))≤exp(Cλ2∑i=1nai2K2)=exp(Cλ2∥a∥2K2).\begin{aligned} \mathbb{E} \exp\left(\lambda \sum_{i=1}^{n}a_i X_i \right) &= \prod_{i=1}^{n} \mathbb{E} \exp(\lambda a_i X_i) \qquad \text{(by independence)}\\ &\le \prod_{i=1}^{n} \exp(C\lambda^2a_i^2\|X_i\|_{\psi_2}^2)\qquad\text{(by Theorem~1~(iii))}\\ &\le \exp(C\lambda^2 \sum_{i=1}^{n} a_i^2 K^2)\\ &= \exp(C\lambda^2 \|a\|^2 K^2). \end{aligned} Then use Theorem 1 (i). ◻

Subgaussian scaling.

A centered subgaussian variable has a Gaussian-type MGF bound. For independent weighted sums, the squared coefficient scales add; a Chernoff argument converts that sum into a tail bound.

High priorityPractice: Unequal weights and concentration

3. Subexponential variables and Bernstein’s inequality

High prioritySubexponential variables and Bernstein’s inequality
Quick Explanation
Up to a constant

We introduce another notion that characterizes the fast decay of the tail probability, though not as fast as subgaussian tails.

Theorem 3Subexponential random variables

Let \(X\)XX be a random variable with \(\mathbb{E}X = 0\)𝔼X=0\mathbb{E}X = 0.

  • (i) \(\mathbb{P}(|X| \ge t) \le 2e^{-t/K_1},~\forall\;t \ge 0\)ℙ(|X|≥t)≤2e−t/K1,∀t≥0\mathbb{P}(|X| \ge t) \le 2e^{-t/K_1},~\forall\;t \ge 0.

  • (ii) \(\sup_{p\ge 1}\frac{1}{p}(\mathbb{E}|X|^p)^{1/p} \le K_2\)supp≥11p(𝔼|X|p)1/p≤K2\sup_{p\ge 1}\frac{1}{p}(\mathbb{E}|X|^p)^{1/p} \le K_2.

  • (iii) \(\mathbb{E}[e^{\lambda X}] \le e^{K_3^2\lambda^2},~\forall\;|\lambda| \le 1 / K_3\)𝔼[eλX]≤eK32λ2,∀|λ|≤1/K3\mathbb{E}[e^{\lambda X}] \le e^{K_3^2\lambda^2},~\forall\;|\lambda| \le 1 / K_3.

The three conditions are equivalent up to absolute constants: each condition with \(K_j\)KjK_j implies the others with suitable \(K_i\le CK_j\)Ki≤CKjK_i\le CK_j. See [1] (Proposition 2.7.1) for the parallel proof.

The smallest \(K_2\)K2K_2 in condition (ii) is called the subexponential norm of \(X\)XX, denoted by \(\|X\|_{\psi_1}\)∥X∥ψ1\|X\|_{\psi_1}. By definition, subgaussian random variables are necessarily subexponential random variables with \(\|X\|_{\psi_1} \le \|X\|_{\psi_2}\)∥X∥ψ1≤∥X∥ψ2\|X\|_{\psi_1} \le \|X\|_{\psi_2}; this is because \[\frac{1}{p}(\mathbb{E}|X|^p)^{1/p} \le \frac{1}{\sqrt{p}}(\mathbb{E}|X|^p)^{1/p}, \qquad \text{for any}~p \ge 1\]1p(𝔼|X|p)1/p≤1p(𝔼|X|p)1/p,for anyp≥1\frac{1}{p}(\mathbb{E}|X|^p)^{1/p} \le \frac{1}{\sqrt{p}}(\mathbb{E}|X|^p)^{1/p}, \qquad \text{for any}~p \ge 1 so taking \(\sup_{p\ge 1}\)supp≥1\sup_{p\ge 1} yields \(\|X\|_{\psi_1} \le \|X\|_{\psi_2}\)∥X∥ψ1≤∥X∥ψ2\|X\|_{\psi_1} \le \|X\|_{\psi_2}. In fact, intuitively, a subexponential random variable can be viewed as the square of a subgaussian random variable. This intuition is captured by the following lemma.

Centering.

For random variables with nonzero mean, conditions (i) and (ii) in Theorem 3 (and similarly in Theorem 1) are still equivalent. Moreover, we have the bound \[\| X - \E X \|_{\psi_2} \le C\|X \|_{\psi_2}.\]∥X−𝔼X∥ψ2≤C∥X∥ψ2.\| X - \mathbb{E}X \|_{\psi_2} \le C\|X \|_{\psi_2}. See Lemma 2.6.8 of the book [1].

Lemma 1

Suppose that \(X\)XX, \(Y\)YY are subgaussian (not necessarily independent). Then \(XY\)XYXY is subexponential with \[\|XY\|_{\psi_1}\le2\|X\|_{\psi_2}\|Y\|_{\psi_2}.\]∥XY∥ψ1≤2∥X∥ψ2∥Y∥ψ2.\|XY\|_{\psi_1}\le2\|X\|_{\psi_2}\|Y\|_{\psi_2}.

Proof. We use the second condition of subgaussian and subexponential random variables. \[\begin{aligned} \left[\mathbb{E}(|X|^p|Y|^p)\right]^{1/p} &\le \left[\mathbb{E}|X|^{2p} \mathbb{E}|Y|^{2p}\right]^{1/(2p)} \qquad \text{ (Cauchy--Schwarz Inequality)} \\ &\le \|X\|_{\psi_2}\sqrt{2p}\,\|Y\|_{\psi_2}\sqrt{2p} =2p\|X\|_{\psi_2}\|Y\|_{\psi_2}. \end{aligned}\][𝔼(|X|p|Y|p)]1/p≤[𝔼|X|2p𝔼|Y|2p]1/(2p) (Cauchy–Schwarz Inequality)≤∥X∥ψ22p∥Y∥ψ22p=2p∥X∥ψ2∥Y∥ψ2.\begin{aligned} \left[\mathbb{E}(|X|^p|Y|^p)\right]^{1/p} &\le \left[\mathbb{E}|X|^{2p} \mathbb{E}|Y|^{2p}\right]^{1/(2p)} \qquad \text{ (Cauchy--Schwarz Inequality)} \\ &\le \|X\|_{\psi_2}\sqrt{2p}\,\|Y\|_{\psi_2}\sqrt{2p} =2p\|X\|_{\psi_2}\|Y\|_{\psi_2}. \end{aligned} ◻

Example 4

Let \(G \sim \mathcal{N}(0, 1)\)G∼𝒩(0,1)G \sim \mathcal{N}(0, 1), then \(G^2\)G2G^2 (namely chi-square \(\chi_1^2\)χ12\chi_1^2) is subexponential, \(\|G^2\|_{\psi_1} \le C\)∥G2∥ψ1≤C\|G^2\|_{\psi_1} \le C.

Example 5

Let \(X \sim \text{Exp}(1)\)X∼Exp(1)X \sim \text{Exp}(1) be an exponential random variable, then \(\|X\|_{\psi_1} \le C\)∥X∥ψ1≤C\|X\|_{\psi_1} \le C.

Example 6

Let \(X\)XX be a random variable such that \(|X| \le M\)|X|≤M|X| \le M is bounded for \(M \ge 0\)M≥0M \ge 0. Then, \(\|X\|_{\psi_1} \le CM\)∥X∥ψ1≤CM\|X\|_{\psi_1} \le CM.

Now we are in a position to state the concentration inequality for a sum of subexponential random variables.

Theorem 4Bernstein’s inequality

Let \(X_1, \ldots, X_n\)X1,…,XnX_1, \ldots, X_n be independent random variables, such that \(\mathbb{E} X_i = 0\)𝔼Xi=0\mathbb{E} X_i = 0 and each \(X_i\)XiX_i is subexponential. Let \(a = (a_1,\ldots,a_n)\)a=(a1,…,an)a = (a_1,\ldots,a_n) be a nonrandom vector. Denote \(K := \max_{1\le i \le n} \|X_i\|_{\psi_1}\)K:=max1≤i≤n∥Xi∥ψ1K := \max_{1\le i \le n} \|X_i\|_{\psi_1}. Then

\[\mathbb{P}\left(\big|\sum_{i=1}^n a_iX_i\big| \ge t\right) \le 2\exp\left(-c\min\left(\frac{t^2}{K^2\|a\|^2}, \frac{t}{K\|a\|_{\infty}}\right)\right),\]ℙ(|∑i=1naiXi|≥t)≤2exp(−cmin(t2K2∥a∥2,tK∥a∥∞)),\mathbb{P}\left(\big|\sum_{i=1}^n a_iX_i\big| \ge t\right) \le 2\exp\left(-c\min\left(\frac{t^2}{K^2\|a\|^2}, \frac{t}{K\|a\|_{\infty}}\right)\right),

(2)

for \(t>0\)t>0t>0, with a positive absolute constant \(c\)cc; a zero scale again means the sum vanishes almost surely.

Subexponential random variables have “heavier” tails compared with subgaussian random variables. As expected, the upper bound in Equation 2 is worse than the bound in Hoeffding’s theorem, since it involves an additional term \(t / (K \| a \|_{\infty})\)t/(K∥a∥∞)t / (K \| a \|_{\infty}) in the exponent.

Proof. We apply a similar strategy as in the proof for Hoeffding’s inequality. \[\mathbb{P}\left(\sum_{i=1}^n a_iX_i \ge t\right) \le e^{-\lambda t} \prod_{i=1}^{n} \mathbb{E}e^{a_i \lambda X_i}.\]ℙ(∑i=1naiXi≥t)≤e−λt∏i=1n𝔼eaiλXi.\mathbb{P}\left(\sum_{i=1}^n a_iX_i \ge t\right) \le e^{-\lambda t} \prod_{i=1}^{n} \mathbb{E}e^{a_i \lambda X_i}. If \(|\lambda| \le \frac{C}{\|a\|_{\infty}K}\)|λ|≤C∥a∥∞K|\lambda| \le \frac{C}{\|a\|_{\infty}K}, then by (iii) and the same argument for the negative sum, \[\mathbb{P}\left(\big|\sum_{i=1}^n a_iX_i\big| \ge t\right) \le 2\exp\left(-\lambda t + C\lambda^2K^2\|a\|^2\right).\]ℙ(|∑i=1naiXi|≥t)≤2exp(−λt+Cλ2K2∥a∥2).\mathbb{P}\left(\big|\sum_{i=1}^n a_iX_i\big| \ge t\right) \le 2\exp\left(-\lambda t + C\lambda^2K^2\|a\|^2\right).

  • Choose a legal parameter. Write the exponent as \(-\lambda t+A\lambda^2\)−λt+Aλ2-\lambda t+A\lambda^2, where \(A=CK^2\|a\|_2^2\)A=CK2∥a∥22A=CK^2\|a\|_2^2, and let \(L=c_0/(K\|a\|_\infty)\)L=c0/(K∥a∥∞)L=c_0/(K\|a\|_\infty) be the permitted MGF range. For a nonzero scale, choose \[\lambda=\min\{t/(2A),L/2\}.\]λ=min{t/(2A),L/2}.\lambda=\min\{t/(2A),L/2\}.

  • Quadratic regime. If \(t\le AL\)t≤ALt\le AL, the exponent is \(-t^2/(4A)\)−t2/(4A)-t^2/(4A).

  • Linear regime. If \(t>AL\)t>ALt>AL, it is at most \(-Lt/2+AL^2/4\le-Lt/4\)−Lt/2+AL2/4≤−Lt/4-Lt/2+AL^2/4\le-Lt/4. These are the quadratic and linear terms in the claimed bound. If the sum is identically zero, handle it directly.

 ◻

Theorem 5Bernstein’s inequality for bounded random variables

Let \(X_1, \ldots, X_n\)X1,…,XnX_1, \ldots, X_n be independent random variables, such that \(\mathbb{E} X_i = 0\)𝔼Xi=0\mathbb{E} X_i = 0 and \(\max_{1\le i \le n}|X_i| \le M\)max1≤i≤n|Xi|≤M\max_{1\le i \le n}|X_i| \le M a.s. Then \[\mathbb{P}\left(\big|\sum_{i=1}^n X_i\big| \ge t\right) \le 2\exp\left(-\frac{t^2/2}{\sigma^2 + Mt/3}\right),\]ℙ(|∑i=1nXi|≥t)≤2exp(−t2/2σ2+Mt/3),\mathbb{P}\left(\big|\sum_{i=1}^n X_i\big| \ge t\right) \le 2\exp\left(-\frac{t^2/2}{\sigma^2 + Mt/3}\right), where \(\sigma^2 = \sum_{i=1}^n \mathbb{E}X_i^2\)σ2=∑i=1n𝔼Xi2\sigma^2 = \sum_{i=1}^n \mathbb{E}X_i^2 is total variance.

Note that \[\max\{\sigma^2, Mt/3\} \le \sigma^2 + Mt/3 \le 2 \max\{\sigma^2, Mt/3\},\]max{σ2,Mt/3}≤σ2+Mt/3≤2max{σ2,Mt/3},\max\{\sigma^2, Mt/3\} \le \sigma^2 + Mt/3 \le 2 \max\{\sigma^2, Mt/3\}, so the exponent in the upper bound of Theorem 5 is roughly (up to a constant2) \[-\frac{t^2}{\max\{\sigma^2, Mt\}} = - \min\left\{ \frac{t^2}{\sigma^2}, \frac{t}{M} \right\}\]−t2max{σ2,Mt}=−min{t2σ2,tM}-\frac{t^2}{\max\{\sigma^2, Mt\}} = - \min\left\{ \frac{t^2}{\sigma^2}, \frac{t}{M} \right\} This is similar to the exponent in the bound in (2)

Remark 1

There is a popular variant of Bernstein’s inequality.

\[\mathbb{P}\left(\big|\sum_{i=1}^{n} X_i\big| \ge \frac{C(\sigma + M)}{\delta}\right) \le 2\exp\left(-\min\left(\frac{1}{\delta}, \frac{1}{\delta^2}\right)\right), \qquad \text{for any}~\delta > 0,\]ℙ(|∑i=1nXi|≥C(σ+M)δ)≤2exp(−min(1δ,1δ2)),for anyδ>0,\mathbb{P}\left(\big|\sum_{i=1}^{n} X_i\big| \ge \frac{C(\sigma + M)}{\delta}\right) \le 2\exp\left(-\min\left(\frac{1}{\delta}, \frac{1}{\delta^2}\right)\right), \qquad \text{for any}~\delta > 0,

(3)

where \(C\)CC is a certain constant. This inequality is a consequence of Theorem 5. In fact, note that \[\frac{-t^2/2}{\sigma^2+Mt/3}\le\frac{-t^2}{4\max\{\sigma^2,Mt/3\}} = -\min\left\{\frac{t^2}{4\sigma^2}, \frac{3t}{4M}\right\}.\]−t2/2σ2+Mt/3≤−t24max{σ2,Mt/3}=−min{t24σ2,3t4M}.\frac{-t^2/2}{\sigma^2+Mt/3}\le\frac{-t^2}{4\max\{\sigma^2,Mt/3\}} = -\min\left\{\frac{t^2}{4\sigma^2}, \frac{3t}{4M}\right\}. So choosing \(t = \frac{C(\sigma + M)}{\delta}\)t=C(σ+M)δt = \frac{C(\sigma + M)}{\delta} with a sufficiently large \(C\)CC leads to \[\frac{t^2}{4\sigma^2} \ge \frac{1}{\delta^2}, \qquad \frac{3t}{4M} \ge \frac{1}{\delta}\]t24σ2≥1δ2,3t4M≥1δ\frac{t^2}{4\sigma^2} \ge \frac{1}{\delta^2}, \qquad \frac{3t}{4M} \ge \frac{1}{\delta} and therefore leads to the upper bound in Equation 3.

Instructor’s NoteTwo costs of a large deviation

A moderate fluctuation can be built from many small contributions, so a variance scale is natural. A very large fluctuation can be driven by a large individual contribution, so the largest weight enters too. Bernstein’s minimum in the exponent keeps both possibilities in view.

Two Bernstein regimes.

Subexponential sums have a quadratic exponent for moderate deviations and a linear exponent for large deviations. Products of subgaussians naturally have this weaker tail behavior.

High priorityPractice: Use a Bernstein radius

See Exercise 3 in Section 7 for the full statement and complete solution.

4. Understanding the two tail probabilities

Mid priorityUnderstanding the two tail probabilities

For the link between tails and expectation, revisit Lecture 4: area under the survival curve.

Vershynin’s Bernstein illustration: Gaussian-type small deviations and exponential-type large deviations.

Figure excerpt: Vershynin, High-Dimensional Probability, Figure 2.3.

Let \(X_1, X_2,\ldots,X_n\)X1,X2,…,XnX_1, X_2,\ldots,X_n be independent random variables, and \(a = (a_1,\ldots,a_n)\)a=(a1,…,an)a = (a_1,\ldots,a_n) a nonrandom vector. We are interested in the tail probability of the sum \(S_n=\sum_{i=1}^n a_iX_i\)Sn=∑i=1naiXiS_n=\sum_{i=1}^n a_iX_i. When each \(X_i\)XiX_i has mean zero and variance one, \(\Var(S_n)=\|a\|^2\)Var(Sn)=∥a∥2\mathrm{Var}(S_n)=\|a\|^2; use this normalization in the following comparison, with uniformly bounded tail norms.

  1. Gaussian distribution. In the canonical case \(X_i \stackrel{d}{=} N(0,1)\)Xi=dN(0,1)X_i \stackrel{d}{=} N(0,1), we have \(S_n \stackrel{d}{=} N(0, \|a\|^2)\)Sn=dN(0,∥a∥2)S_n \stackrel{d}{=} N(0, \|a\|^2), so the distribution of \(S_n\)SnS_n is completely determined by its variance. According to Proposition 1, the tail probability is mainly determined by \(e^{-t^2 / (2\|a\|^2)}\)e−t2/(2∥a∥2)e^{-t^2 / (2\|a\|^2)} (ignoring poly-factor).

  2. Asymptotic normal distribution. Under certain conditions we expect the asymptotic normality \(S_n / \| a\| \stackrel{d}{\to} N(0, 1)\)Sn/∥a∥→dN(0,1)S_n / \| a\| \stackrel{d}{\to} N(0, 1), so approximately \(S_n \stackrel{d}{\approx} N(0, \|a\|^2)\)Sn≈dN(0,∥a∥2)S_n \stackrel{d}{\approx} N(0, \|a\|^2). However, as we explained in Section 1, the CLT controls fixed standardized thresholds; it does not by itself give a relative approximation for increasingly rare tails.

  3. Subgaussian distribution. If the tail probability of each individual \(X_i\)XiX_i is “lighter” (in the sense of Theorem 1(i)) than the standard Gaussian tail, then we still get a “light” tail controlled by \(e^{-ct^2 / (2\|a\|^2)}\)e−ct2/(2∥a∥2)e^{-ct^2 / (2\|a\|^2)} (up to constant \(c\)cc in the exponent).

  4. Subexponential distribution. If the tail of \(X_i\)XiX_i is slightly heavier than Gaussian, then \(S_n\)SnS_n has an upper bound with two deviation regimes in Equation 2—one is Gaussian \(e^{-ct^2 / (2\|a\|^2)}\)e−ct2/(2∥a∥2)e^{-ct^2 / (2\|a\|^2)} and another is exponential \(e^{-ct / \|a\|_{\infty}}\)e−ct/∥a∥∞e^{-ct / \|a\|_{\infty}}. The tail probability is controlled by the Gaussian tail when \(t\)tt is moderately large and by the exponential tail when \(t\)tt is very large. A Gaussian bound is unavailable in general: consider \(a = (1,0,\ldots,0)\)a=(1,0,…,0)a = (1,0,\ldots,0) and \(X_1\stackrel d=\mathrm{Exp}(1)-1\)X1=dExp(1)−1X_1\stackrel d=\mathrm{Exp}(1)-1, where the tail probability is purely exponential.

  • Two tail scales. Bernstein’s bound has two scales: the accumulated tail scale \(K\|a\|_2\)K∥a∥2K\|a\|_2 and the largest individual scale \(K\|a\|_\infty\)K∥a∥∞K\|a\|_\infty.

  • Two regimes. The former gives a quadratic exponent for moderate deviations; the latter gives a linear exponent for large deviations.

  • Variance versus a proxy. For independent unit-variance summands, \(\|a\|_2\)∥a∥2\|a\|_2 is the actual standard deviation, but in general \(K\|a\|_2\)K∥a∥2K\|a\|_2 is only a tail proxy. The accompanying figure illustrates the two regimes.

Theorem 5 controls the tail probability in a slightly different way: first using a variance quantity \(\sigma^2\)σ2\sigma^2 similar to \(\|a\|^2\)∥a∥2\|a\|^2 above, second using a crude bound \(M\)MM on \(X_i\)XiX_i (implies subexponential norm at most \(M\)MM) to control the potentially “heavier” tail of \(X_i\)XiX_i.

High-probability bound.

Finally, to get a high-probability bound on \(S_n\)SnS_n, we conclude that \(|S_n|\)|Sn||S_n| is bounded by the following quantity with probability at least \(1-\delta\)1−δ1-\delta.

  1. Gaussian case: \(C\sigma\sqrt{\log(2/\delta)}\)Cσlog(2/δ)C\sigma\sqrt{\log(2/\delta)}, where \(\sigma^2=\Var(S_n)\)σ2=Var(Sn)\sigma^2=\mathrm{Var}(S_n). For subgaussian summands, replace \(\sigma\)σ\sigma by the tail scale \(K\|a\|_2\)K∥a∥2K\|a\|_2; this need not be comparable to the actual standard deviation.

  2. Subexponential case: \(CK[\|a\|_2\sqrt{\log(2/\delta)}+\|a\|_\infty\log(2/\delta)]\)CK[∥a∥2log(2/δ)+∥a∥∞log(2/δ)]CK[\|a\|_2\sqrt{\log(2/\delta)}+\|a\|_\infty\log(2/\delta)]. For bounded summands, the variance-sensitive radius below is sharper than the coarser \(C(\sigma+M)\log(2/\delta)\)C(σ+M)log(2/δ)C(\sigma+M)\log(2/\delta) form when \(0<\delta\le1/2\)0<δ≤1/20<\delta\le1/2.

A sharper explicit radius.

  • Failure probability. For bounded centered summands, let \(u=\log(2/\alpha)\)u=log(2/α)u=\log(2/\alpha) with \(0<\alpha<1\)0<α<10<\alpha<1. A useful Bernstein radius is \[r_B(\alpha)=\sqrt{2\sigma^2u}+\frac{2M}{3}u,\qquad \PR(|S_n|\ge r_B(\alpha))\le\alpha.\]rB(α)=2σ2u+2M3u,ℙ(|Sn|≥rB(α))≤α.r_B(\alpha)=\sqrt{2\sigma^2u}+\frac{2M}{3}u,\qquad \mathbb{P}(|S_n|\ge r_B(\alpha))\le\alpha.

  • Interpretation. The square-root term represents accumulated variance; the linear term records the size of one bounded summand.

  • Example. For \(\sigma^2=1\)σ2=1\sigma^2=1, \(M=1\)M=1M=1, \(u=1\)u=1u=1 and \(u=4\)u=4u=4, the radii are \(\sqrt2+2/3\)2+2/3\sqrt2+2/3 and \(\sqrt8+8/3\)8+8/3\sqrt8+8/3.

Turn a tail bound into an error bar.

For centered independent bounded summands, a sufficient Bernstein radius at failure probability \(\alpha\)α\alpha is \(\sqrt{2\sigma^2\log(2/\alpha)}+(2M/3)\log(2/\alpha)\)2σ2log(2/α)+(2M/3)log(2/α)\sqrt{2\sigma^2\log(2/\alpha)}+(2M/3)\log(2/\alpha), with total variance \(\sigma^2\)σ2\sigma^2 and bound \(M\)MM.

Mid priorityPractice: Invert the bounded Bernstein bound

See Exercise 4 in Section 7 for the full statement and complete solution.

5. Further reading

Further reading

The recent book [1] on high-dimensional probability provides an excellent introduction to this topic. The concentration material and the deviation sketch are adapted from Chapter 2 of this book. The course uses the moment-norm convention stated above, so numerical constants may differ from the book’s exponential-moment norm. The author’s first-edition text contains the cited numbering.

Key ideas

  • Subgaussian and subexponential scales describe different tail behaviors.

  • Chernoff bounds combine Markov with control of exponential moments.

  • Hoeffding and Bernstein bound independent sums at finite sample sizes.

  • Bernstein has quadratic and linear deviation regimes.

Detailed reference

Detailed reference

High prioritySubgaussian scaling

A centered subgaussian variable has a Gaussian-type MGF bound. For independent weighted sums, the squared coefficient scales add; a Chernoff argument converts that sum into a tail bound. Return to the earlier explanation.

High priorityTwo Bernstein regimes

Subexponential sums have a quadratic exponent for moderate deviations and a linear exponent for large deviations. Products of subgaussians naturally have this weaker tail behavior. Return to the earlier explanation.

Mid priorityTurn a tail bound into an error bar

For centered independent bounded summands, a sufficient Bernstein radius at failure probability \(\alpha\)α\alpha is \(\sqrt{2\sigma^2\log(2/\alpha)}+(2M/3)\log(2/\alpha)\)2σ2log(2/α)+(2M/3)log(2/α)\sqrt{2\sigma^2\log(2/\alpha)}+(2M/3)\log(2/\alpha), with total variance \(\sigma^2\)σ2\sigma^2 and bound \(M\)MM. Return to the earlier explanation.

Looking ahead

Lecture 15 applies these tail bounds to random vectors and their norms. Independence and union bounds turn one-dimensional control into geometric statements in many dimensions. Continue to Lecture 15.

7. Exercises and complete solutions8 exercises · complete solutions

Mid priorityA continuous CDF cannot match a jump
Exercise 1A continuous CDF cannot match a jump#

If a CDF \(F\)FF has a jump of size \(q\)qq at \(x_0\)x0x_0 and \(G\)GG is continuous there, prove \(\sup_x|F(x)-G(x)|\ge q/2\)supx|F(x)−G(x)|≥q/2\sup_x|F(x)-G(x)|\ge q/2.

Reveal solution

Let \(a=F(x_0-)\)a=F(x0−)a=F(x_0-) and \(b=F(x_0)\)b=F(x0)b=F(x_0), so \(b-a=q\)b−a=qb-a=q. Continuity gives \(G(x_0-)=G(x_0)=g\)G(x0−)=G(x0)=gG(x_0-)=G(x_0)=g. Since \(q\le|b-g|+|g-a|\)q≤|b−g|+|g−a|q\le|b-g|+|g-a|, one of the two errors is at least \(q/2\)q/2q/2. The left error is attained as a one-sided limit and hence is bounded by the supremum.

High priorityUnequal weights and concentration
Exercise 2Unequal weights and concentration#

For independent mean-zero variables with subgaussian norm at most \(K\)KK, compare equal weights \(a_i=1/n\)ai=1/na_i=1/n with weights \(a_1=1\)a1=1a_1=1, \(a_i=0\)ai=0a_i=0 for \(i>1\)i>1i>1 in Hoeffding’s inequality.

Reveal solution

Equal weights give \(\|a\|^2=1/n\)∥a∥2=1/n\|a\|^2=1/n and the bound \(2e^{-cnt^2/K^2}\)2e−cnt2/K22e^{-cnt^2/K^2}. A single nonzero weight gives \(\|a\|^2=1\)∥a∥2=1\|a\|^2=1 and \(2e^{-ct^2/K^2}\)2e−ct2/K22e^{-ct^2/K^2}. Averaging independent terms produces the sample-size improvement; merely having \(n\)nn observed variables does not.

High priorityUse a Bernstein radius
Exercise 3Use a Bernstein radius#

Let independent centered \(Y_i\)YiY_i satisfy \(|Y_i|\le2\)|Yi|≤2|Y_i|\le2 and \(\sum_i\Var(Y_i)=9\)∑iVar(Yi)=9\sum_i\mathrm{Var}(Y_i)=9. Use Bernstein to give a radius with two-sided failure probability at most \(2e^{-4}\)2e−42e^{-4}.

Reveal solution

Take \(u=4\)u=4u=4, \(\sigma^2=9\)σ2=9\sigma^2=9 and \(M=2\)M=2M=2. The sufficient radius is \(\sqrt{2(9)(4)}+2(2)(4)/3=6\sqrt2+16/3\)2(9)(4)+2(2)(4)/3=62+16/3\sqrt{2(9)(4)}+2(2)(4)/3=6\sqrt2+16/3. Indeed this radius satisfies \(t^2\ge2u(\sigma^2+Mt/3)\)t2≥2u(σ2+Mt/3)t^2\ge2u(\sigma^2+Mt/3), so Bernstein gives the required bound. The exact positive quadratic root would be slightly smaller.

Mid priorityInvert the bounded Bernstein bound
Exercise 4Invert the bounded Bernstein bound#

Prove the explicit radius \(r_B(\alpha)\)rB(α)r_B(\alpha) displayed above, assuming \(\sigma^2,M>0\)σ2,M>0\sigma^2,M>0.

Reveal solution

Let \(a=\sqrt{2\sigma^2u}\)a=2σ2ua=\sqrt{2\sigma^2u} and \(b=2Mu/3\)b=2Mu/3b=2Mu/3. For \(t=a+b\)t=a+bt=a+b, \(t^2-(2uM/3)t-2u\sigma^2=(a+b)^2-b(a+b)-a^2=ab\ge0\)t2−(2uM/3)t−2uσ2=(a+b)2−b(a+b)−a2=ab≥0t^2-(2uM/3)t-2u\sigma^2=(a+b)^2-b(a+b)-a^2=ab\ge0. Therefore \(t^2/[2(\sigma^2+Mt/3)]\ge u\)t2/[2(σ2+Mt/3)]≥ut^2/[2(\sigma^2+Mt/3)]\ge u and the tail bound is at most \(2e^{-u}=\alpha\)2e−u=α2e^{-u}=\alpha. If all summands vanish the assertion is immediate.

High priorityThe moment subgaussian norm
Exercise 5The moment subgaussian norm#

Verify homogeneity, definiteness and the triangle inequality for \(\|X\|_{\psi_2}=\sup_{p\ge1}p^{-1/2}\|X\|_{L^p}\)∥X∥ψ2=supp≥1p−1/2∥X∥Lp\|X\|_{\psi_2}=\sup_{p\ge1}p^{-1/2}\|X\|_{L^p}.

Reveal solution

Homogeneity follows from that of each \(L^p\)LpL^p norm. If the supremum is zero, its \(p=1\)p=1p=1 term gives \(\E|X|=0\)𝔼|X|=0\mathbb{E}|X|=0, so \(X=0\)X=0X=0 a.s.; the converse is immediate. Minkowski gives \(p^{-1/2}\|X+Y\|_{L^p}\le p^{-1/2}\|X\|_{L^p}+p^{-1/2}\|Y\|_{L^p}\)p−1/2∥X+Y∥Lp≤p−1/2∥X∥Lp+p−1/2∥Y∥Lpp^{-1/2}\|X+Y\|_{L^p}\le p^{-1/2}\|X\|_{L^p}+p^{-1/2}\|Y\|_{L^p}, and taking suprema proves the triangle inequality on the space where the norm is finite.

High priorityOptimize a quadratic MGF bound
Exercise 6Optimize a quadratic MGF bound#

If \(\E e^{\lambda Z}\le e^{v\lambda^2/2}\)𝔼eλZ≤evλ2/2\mathbb{E}e^{\lambda Z}\le e^{v\lambda^2/2} for every real \(\lambda\)λ\lambda and \(v>0\)v>0v>0, prove \(\PR(|Z|\ge t)\le2e^{-t^2/(2v)}\)ℙ(|Z|≥t)≤2e−t2/(2v)\mathbb{P}(|Z|\ge t)\le2e^{-t^2/(2v)}.

Reveal solution

For \(\lambda>0\)λ>0\lambda>0, Markov gives \(\PR(Z\ge t)\le e^{-\lambda t+v\lambda^2/2}\)ℙ(Z≥t)≤e−λt+vλ2/2\mathbb{P}(Z\ge t)\le e^{-\lambda t+v\lambda^2/2}. The minimum occurs at \(\lambda=t/v\)λ=t/v\lambda=t/v, giving \(e^{-t^2/(2v)}\)e−t2/(2v)e^{-t^2/(2v)}. Repeat for \(-Z\)−Z-Z and use the union bound. At \(t=0\)t=0t=0, the bound is trivially valid.

High priorityCentering a tail norm
Exercise 7Centering a tail norm#

For either moment norm \(\psi_1\)ψ1\psi_1 or \(\psi_2\)ψ2\psi_2, show \(\|X-\E X\|_{\psi_j}\le2\|X\|_{\psi_j}\)∥X−𝔼X∥ψj≤2∥X∥ψj\|X-\mathbb{E}X\|_{\psi_j}\le2\|X\|_{\psi_j}.

Reveal solution

The norm of the constant \(\E X\)𝔼X\mathbb{E}X is \(|\E X|\)|𝔼X||\mathbb{E}X|, because the supremum of \(p^{-1/j}\)p−1/jp^{-1/j} is 1. Also \(|\E X|\le\E|X|\le\|X\|_{\psi_j}\)|𝔼X|≤𝔼|X|≤∥X∥ψj|\mathbb{E}X|\le\mathbb{E}|X|\le\|X\|_{\psi_j}. The triangle inequality gives the claim.

Mid priorityLocate the two-scale crossover
Exercise 8Locate the two-scale crossover#

For a nonzero weight vector, find where \(t^2/(K^2\|a\|_2^2)=t/(K\|a\|_\infty)\)t2/(K2∥a∥22)=t/(K∥a∥∞)t^2/(K^2\|a\|_2^2)=t/(K\|a\|_\infty) for \(t>0\)t>0t>0. Which exponent is smaller on either side?

Reveal solution

Cancel \(t\)tt to obtain \(t_*=K\|a\|_2^2/\|a\|_\infty\)t*=K∥a∥22/∥a∥∞t_*=K\|a\|_2^2/\|a\|_\infty. For \(0<t<t_*\)0<t<t*0<t<t_* the quadratic term is smaller; for \(t>t_*\)t>t*t>t_* the linear term is smaller. Equal weights \(1/n\)1/n1/n give \(t_*=K\)t*=Kt_*=K, while the exponent in the quadratic range grows proportionally to \(n\)nn.

References

[1]
R. Vershynin. High-dimensional probability: An introduction with applications in data science. Cambridge university press. 2018.

Terminology notes


  1. Note: \(\mathbb{E}[e^{\lambda X}]\)𝔼[eλX]\mathbb{E}[e^{\lambda X}] is the moment-generating function (MGF) of \(X\)XX.↩︎

  2. When we say two quantities \(s\)ss and \(t\)tt are up to a constant, we mean \(ct \le s \le Ct\)ct≤s≤Ctct \le s \le Ct for absolute constants \(c,C>0\)c,C>0c,C>0.↩︎

What changed

October 19, 2026: Corrected notation, wording and cross-references after the course-wide errata review; refreshed the linked HTML and PDF editions.

October 14, 2026. Added issue-reporting links to the webpage and PDF. The mathematics and priority assignments are unchanged in this update.

October 14, 2026 (figure follow-up). Restored the original two-regime tail illustration. The previously removed numbered Bernstein figure and explorer remain removed.

Comments

Discuss an idea or ask a public question. Include a section or exercise number when relevant. For errors in the notes or website, Report an issue.

Public comments use Disqus. Enable JavaScript to load them.

Email the instructor for a private question or if comments are unavailable.