Math ClubMath Club
v1 · padrão canônico

Lição 77 — Teorema Central do Limite

A média de n v.a. iid converge à normal independente da distribuição original — a lei mais importante da estatística. Demonstração via função característica, velocidade Berry-Esseen, aplicações de inferência.

Used in: 2.º ano do EM (16-17 anos) · Math B japonês §4.4 · Stochastik LK alemão · H2 Math singapurense cap. 21

XˉndN ⁣(μ,σ2n)(n)\bar X_n \xrightarrow{d} \mathcal{N}\!\left(\mu,\,\frac{\sigma^2}{n}\right) \quad (n \to \infty)

The Central Limit Theorem (CLT) states that the mean of nn independent and identically distributed random variables with mean μ\mu and finite variance σ2\sigma^2 converges in distribution to a normal distribution with parameters μ\mu and σ2/n\sigma^2/n — regardless of the original distribution. This explains why the bell curve appears in natural data, measurement errors, and sample means.

Choose your door

Rigorous notation, full derivation, hypotheses

Formal statement and proof

Lindeberg-Lévy Version

"The central limit theorem is the unofficial sovereign of probability theory." — Grinstead & Snell, Introduction to Probability, §9.1

Version for sums

If Sn=X1++XnS_n = X_1 + \cdots + X_n, then SnN(nμ,nσ2)S_n \approx \mathcal{N}(n\mu,\, n\sigma^2) for large nn.

Zn=SnnμσndN(0,1)Z_n = \frac{S_n - n\mu}{\sigma\sqrt{n}} \xrightarrow{d} \mathcal{N}(0,1)
what this means · Standardization of the sum: same formula, different scaling.

Convergence rate: Berry-Esseen inequality

Proof sketch via characteristic function

Let Yi=(Xiμ)/σY_i = (X_i - \mu)/\sigma (zero mean, unit variance). Taylor expansion of φYi\varphi_{Y_i}:

φYi(t)=1t22+o(t2)(t0).\varphi_{Y_i}(t) = 1 - \frac{t^2}{2} + o(t^2) \quad (t \to 0).

For Zn=(Y1++Yn)/nZ_n = (Y_1 + \cdots + Y_n)/\sqrt{n}:

φZn(t)=[φYi ⁣(tn)]n=[1t22n+o ⁣(1n)]nnet2/2.\varphi_{Z_n}(t) = \left[\varphi_{Y_i}\!\left(\frac{t}{\sqrt{n}}\right)\right]^n = \left[1 - \frac{t^2}{2n} + o\!\left(\frac{1}{n}\right)\right]^n \xrightarrow{n\to\infty} e^{-t^2/2}.

But et2/2e^{-t^2/2} is the characteristic function of N(0,1)\mathcal{N}(0, 1). The Lévy continuity theorem concludes ZndN(0,1)Z_n \xrightarrow{d} \mathcal{N}(0,1). \blacksquare

When the CLT fails

Essential hypotheses

  • Independence (minimally sufficient; relaxable to α\alpha-mixing).
  • Finite variance σ2<\sigma^2 < \infty.
  • nn sufficiently large — practical rule: n30n \geq 30 for non-highly-skewed distributions; n100n \geq 100 for high skewness.

Worked examples

Exercise list

45 exercises · 11 with worked solution (25%)

Application 23Understanding 5Modeling 12Challenge 3Proof 2
  1. Ex. 77.1Application

    XX exponential with μ=1\mu = 1, σ=1\sigma = 1. Write the approximate distribution of Xˉ100\bar X_{100} and calculate σXˉ\sigma_{\bar X}.

    Select the correct option
    Select an option first
    Show solution
    XX exponential with μ=1\mu = 1 and σ=1\sigma = 1. By CLT: Xˉ100N(1,  1/100)=N(1,  0.01)\bar X_{100} \approx \mathcal N(1,\; 1/100) = \mathcal N(1,\; 0.01). Standard deviation of the mean: σXˉ=1/100=0.1\sigma_{\bar X} = 1/\sqrt{100} = 0.1.
  2. Ex. 77.2Application

    XU[0,1]X \sim \mathcal U[0, 1]. Determine μ\mu and σ2\sigma^2 and write the approximate distribution of Xˉ50\bar X_{50} by CLT.

    Select the correct option
    Select an option first
    Show solution
    XU[0,1]X \sim \mathcal U[0, 1]: μ=0.5\mu = 0.5, σ2=1/12\sigma^2 = 1/12. By CLT: Xˉ50N(0.5,  (1/12)/50)=N(0.5,  1/600)\bar X_{50} \approx \mathcal N(0.5,\; (1/12)/50) = \mathcal N(0.5,\; 1/600).
    Show step-by-step (with the why)
    1. Calculate parameters of U[0,1]\mathcal U[0,1]: μ=0.5\mu = 0.5; σ2=(ba)2/12=1/12\sigma^2 = (b-a)^2/12 = 1/12.
    2. Apply CLT with n=50n = 50: Xˉ50N(μ,σ2/n)=N(0.5,  1/600)\bar X_{50} \approx \mathcal N(\mu, \sigma^2/n) = \mathcal N(0.5,\; 1/600).
    3. Standard deviation: σXˉ=1/6000.041\sigma_{\bar X} = \sqrt{1/600} \approx 0.041.
    Trick: For uniform [a,b][a,b]: σ2=(ba)2/12\sigma^2 = (b-a)^2/12.
  3. Ex. 77.3Application

    Roll 100 fair dice. Determine the approximate distribution of the sum S100S_{100}, with E[S]E[S] and Var(S)\text{Var}(S).

    Select the correct option
    Select an option first
    Show solution
    Fair die: μ=3.5\mu = 3.5, σ2=35/12\sigma^2 = 35/12. Sum of 100: S100N(100×3.5,  100×35/12)=N(350,  291.67)S_{100} \approx \mathcal N(100 \times 3.5,\; 100 \times 35/12) = \mathcal N(350,\; 291.67).
  4. Ex. 77.4ApplicationAnswer key

    XBernoulli(0.3)X \sim \text{Bernoulli}(0.3). Write the approximate distribution of Xˉ200\bar X_{200} by CLT and calculate the standard deviation of the sample proportion.

    Select the correct option
    Select an option first
    Show solution
    XBernoulli(0.3)X \sim \text{Bernoulli}(0.3): μ=0.3\mu = 0.3, σ2=0.21\sigma^2 = 0.21. CLT: Xˉ200N(0.3,  0.21/200)=N(0.3,  0.00105)\bar X_{200} \approx \mathcal N(0.3,\; 0.21/200) = \mathcal N(0.3,\; 0.00105). σp^=0.001050.0324\sigma_{\hat p} = \sqrt{0.00105} \approx 0.0324.
  5. Ex. 77.5Application

    μ=50\mu = 50, σ=10\sigma = 10, n=25n = 25. Calculate the standard deviation of Xˉ\bar X.

    Select the correct option
    Select an option first
    Show solution
    σXˉ=σ/n=10/25=2\sigma_{\bar X} = \sigma/\sqrt{n} = 10/\sqrt{25} = 2.
  6. Ex. 77.6Application

    With μ=50\mu = 50, σ=10\sigma = 10, n=25n = 25, calculate P(Xˉ>53)P(\bar X > 53).

    Select the correct option
    Select an option first
    Show solution
    σXˉ=10/25=2\sigma_{\bar X} = 10/\sqrt{25} = 2. z=(5350)/2=1.5z = (53-50)/2 = 1.5. P(Xˉ>53)=P(Z>1.5)0.0668P(\bar X > 53) = P(Z > 1.5) \approx 0.0668.
    Show step-by-step (with the why)
    1. Calculate σXˉ=σ/n=10/5=2\sigma_{\bar X} = \sigma/\sqrt{n} = 10/5 = 2.
    2. Standardize: z=(5350)/2=1.5z = (53-50)/2 = 1.5.
    3. Table: P(Z>1.5)=1Φ(1.5)0.0668P(Z > 1.5) = 1 - \Phi(1.5) \approx 0.0668.
    Note: Always use σ/n\sigma/\sqrt{n} (standard error of the mean), not σ\sigma alone.
  7. Ex. 77.7Application

    With μ=50\mu = 50, σ=10\sigma = 10, n=25n = 25, calculate P(Xˉ<47)P(\bar X < 47).

    Select the correct option
    Select an option first
    Show solution
    z=(4750)/2=1.5z = (47-50)/2 = -1.5. P(Xˉ<47)=Φ(1.5)0.0668P(\bar X < 47) = \Phi(-1.5) \approx 0.0668. Same value as 77.6 by symmetry of the normal.
  8. Ex. 77.8Application

    XX with μ=100\mu = 100, σ=20\sigma = 20, n=100n = 100. Calculate P(98<Xˉ<102)P(98 < \bar X < 102).

    Select the correct option
    Select an option first
    Show solution
    σXˉ=20/100=2\sigma_{\bar X} = 20/\sqrt{100} = 2. z1=(98100)/2=1z_1 = (98-100)/2 = -1, z2=(102100)/2=1z_2 = (102-100)/2 = 1. P(1<Z<1)0.6827P(-1 < Z < 1) \approx 0.6827 — 68% rule.
  9. Ex. 77.9Application

    Sum of 50 iid r.v. with μ=5\mu = 5, σ=2\sigma = 2. Calculate P(S50>270)P(S_{50} > 270).

    Select the correct option
    Select an option first
    Show solution
    S50N(250,  50×4)S_{50} \approx \mathcal N(250,\; 50 \times 4). σS=20014.14\sigma_S = \sqrt{200} \approx 14.14. z=(270250)/14.141.414z = (270-250)/14.14 \approx 1.414. P(S50>270)0.079P(S_{50} > 270) \approx 0.079.
    Show step-by-step (with the why)
    1. Distribution of the sum: S50N(nμ,  nσ2)=N(250,  200)S_{50} \approx \mathcal N(n\mu,\; n\sigma^2) = \mathcal N(250,\; 200).
    2. σS=20014.14\sigma_S = \sqrt{200} \approx 14.14.
    3. Standardize: z=(270250)/14.141.41z = (270-250)/14.14 \approx 1.41.
    4. P(S>270)=1Φ(1.41)0.079P(S > 270) = 1 - \Phi(1.41) \approx 0.079.
    Difference between Xˉ\bar X and SS: E[S]=nμE[S] = n\mu, Var(S)=nσ2\text{Var}(S) = n\sigma^2; for Xˉ\bar X: E[Xˉ]=μE[\bar X] = \mu, Var(Xˉ)=σ2/n\text{Var}(\bar X) = \sigma^2/n.
  10. Ex. 77.10Application

    XX with σ=3\sigma = 3. How many observations for 95% CI with margin of error ±0.5\pm 0.5?

    Select the correct option
    Select an option first
    Show solution
    95% CI half-width = 1.96×σ/n=0.51.96 \times \sigma/\sqrt n = 0.5. n=1.96×3/0.5=11.76\sqrt n = 1.96 \times 3 / 0.5 = 11.76. n=11.762138.3n = 11.76^2 \approx 138.3. Round up: n=139n = 139.
    Show step-by-step (with the why)
    1. Margin of error at 95%: ME=1.96σ/n0.5ME = 1.96 \sigma/\sqrt{n} \leq 0.5.
    2. Isolate n1.96×3/0.5=11.76\sqrt{n} \geq 1.96 \times 3/0.5 = 11.76.
    3. Square: n138.3n \geq 138.3.
    4. Round up: n=139n = 139.
    Trick: To halve the margin of error, quadruple nn.
  11. Ex. 77.11Understanding

    When nn is multiplied by 4, the standard deviation of Xˉ\bar X (=σ/n= \sigma/\sqrt{n}):

    Select the correct option
    Select an option first
    Show solution
    σXˉ=σ/n\sigma_{\bar X} = \sigma/\sqrt{n}. If n4nn \to 4n: σXˉσ/4n=σXˉ/2\sigma_{\bar X} \to \sigma/\sqrt{4n} = \sigma_{\bar X}/2. The standard deviation of the mean is cut in half.
  12. Ex. 77.12UnderstandingAnswer key

    XX has a very skewed distribution (skewness = 3). For what sample size is the CLT reasonable?

    Select the correct option
    Select an option first
    Show solution
    Berry-Esseen: error ρ/(σ3n)\propto \rho/(\sigma^3 \sqrt n), where ρ=E[Xμ3]\rho = E[|X-\mu|^3]. High skewness increases ρ\rho, requiring larger nn for convergence. Rule of thumb: n100n \geq 100 for highly skewed distributions.
  13. Ex. 77.13Application

    Test scores with μ=70\mu = 70, σ=15\sigma = 15. Sample n=36n = 36. Calculate P(Xˉ>75)P(\bar X > 75).

    Select the correct option
    Select an option first
    Show solution
    σXˉ=15/36=2.5\sigma_{\bar X} = 15/\sqrt{36} = 2.5. z=(7570)/2.5=2.0z = (75-70)/2.5 = 2.0. P(Xˉ>75)=1Φ(2)0.0228P(\bar X > 75) = 1 - \Phi(2) \approx 0.0228.
    Show step-by-step (with the why)
    1. Standard error of the mean: σXˉ=15/36=15/6=2.5\sigma_{\bar X} = 15/\sqrt{36} = 15/6 = 2.5.
    2. Standardize the value 75: z=(7570)/2.5=2z = (75 - 70)/2.5 = 2.
    3. Probability: P(Xˉ>75)=P(Z>2)=1Φ(2)0.0228P(\bar X > 75) = P(Z > 2) = 1 - \Phi(2) \approx 0.0228.
    Trick: The 95-5 rule says only ~2.3% lies above μ+2σXˉ\mu + 2\sigma_{\bar X}.
  14. Ex. 77.14Application

    With μ=70\mu = 70, σ=15\sigma = 15, n=36n = 36, calculate P(Xˉ<65)P(\bar X < 65).

    Select the correct option
    Select an option first
    Show solution
    z=(6570)/2.5=2.0z = (65-70)/2.5 = -2.0. P(Xˉ<65)=Φ(2)0.0228P(\bar X < 65) = \Phi(-2) \approx 0.0228. Same value as 77.13 by symmetry of the normal.
  15. Ex. 77.15Application

    With μ=70\mu = 70, σ=15\sigma = 15, n=36n = 36 and Xˉ=72\bar X = 72, construct a 95% CI for μ\mu.

    Select the correct option
    Select an option first
    Show solution
    σXˉ=2.5\sigma_{\bar X} = 2.5. 95% CI: 72±1.96×2.5=72±4.972 \pm 1.96 \times 2.5 = 72 \pm 4.9. Interval: (67.1;  76.9)(67.1;\; 76.9).
  16. Ex. 77.16Application

    Package weight: μ=500\mu = 500 g, σ=50\sigma = 50 g, n=25n = 25. Calculate P(Xˉ>520)P(\bar X > 520).

    Select the correct option
    Select an option first
    Show solution
    σXˉ=50/25=10\sigma_{\bar X} = 50/\sqrt{25} = 10. z=(520500)/10=2z = (520-500)/10 = 2. P(Xˉ>520)=1Φ(2)0.0228P(\bar X > 520) = 1 - \Phi(2) \approx 0.0228.
  17. Ex. 77.17Application

    With μ=500\mu = 500, σ=50\sigma = 50, n=25n = 25, calculate P(485<Xˉ<515)P(485 < \bar X < 515).

    Select the correct option
    Select an option first
    Show solution
    z1=(485500)/10=1.5z_1 = (485-500)/10 = -1.5; z2=(515500)/10=1.5z_2 = (515-500)/10 = 1.5. P(1.5<Z<1.5)0.8664P(-1.5 < Z < 1.5) \approx 0.8664.
  18. Ex. 77.18ApplicationAnswer key

    Response time: μ=50\mu = 50 ms, σ=10\sigma = 10 ms, n=100n = 100. What is the 95% SLA limit?

    Select the correct option
    Select an option first
    Show solution
    σXˉ=10/100=1\sigma_{\bar X} = 10/\sqrt{100} = 1. 95th percentile: x=50+1.645×1=51.645x = 50 + 1.645 \times 1 = 51.645 ms.
  19. Ex. 77.19Application

    Roll a die 1,000 times. Calculate P(Xˉ>3.6)P(\bar X > 3.6).

    Select the correct option
    Select an option first
    Show solution
    Die: μ=3.5\mu = 3.5, σ2=35/12\sigma^2 = 35/12. σXˉ=35/120000.054\sigma_{\bar X} = \sqrt{35/12000} \approx 0.054. z=(3.63.5)/0.0541.85z = (3.6-3.5)/0.054 \approx 1.85. P(Xˉ>3.6)0.032P(\bar X > 3.6) \approx 0.032.
  20. Ex. 77.20Application

    Using the distribution of the sum S1000S_{1000} of 1,000 die rolls, calculate P(S1000>3600)P(S_{1000} > 3600).

    Select the correct option
    Select an option first
    Show solution
    S1000N(3500,  1000×35/12)S_{1000} \approx \mathcal N(3500,\; 1000 \times 35/12). σS=2916.754\sigma_S = \sqrt{2916.7} \approx 54. z=(36003500)/541.85z = (3600-3500)/54 \approx 1.85. P(S>3600)0.032P(S > 3600) \approx 0.032. Same result as 77.19 — confirms equivalence.
  21. Ex. 77.21Application

    XExp(1)X \sim \text{Exp}(1) (μ=1\mu = 1, σ=1\sigma = 1). Calculate P(Xˉ100>1.1)P(\bar X_{100} > 1.1).

    Select the correct option
    Select an option first
    Show solution
    XExp(1)X \sim \text{Exp}(1): μ=1\mu = 1, σ=1\sigma = 1. σXˉ=1/100=0.1\sigma_{\bar X} = 1/\sqrt{100} = 0.1. z=(1.11)/0.1=1z = (1.1-1)/0.1 = 1. P(Xˉ>1.1)=1Φ(1)0.1587P(\bar X > 1.1) = 1 - \Phi(1) \approx 0.1587.
  22. Ex. 77.22ApplicationAnswer key

    Election survey: p=0.40p = 0.40, n=1000n = 1000. Calculate P(p^>0.43)P(\hat p > 0.43).

    Select the correct option
    Select an option first
    Show solution
    σp^=0.40×0.60/10000.01549\sigma_{\hat p} = \sqrt{0.40 \times 0.60/1000} \approx 0.01549. z=(0.430.40)/0.015491.937z = (0.43-0.40)/0.01549 \approx 1.937. P(p^>0.43)0.0263P(\hat p > 0.43) \approx 0.0263.
    Show step-by-step (with the why)
    1. Standard error of the proportion: σp^=p(1p)/n=0.24/10000.01549\sigma_{\hat p} = \sqrt{p(1-p)/n} = \sqrt{0.24/1000} \approx 0.01549.
    2. Standardize: z=(0.430.40)/0.015491.94z = (0.43-0.40)/0.01549 \approx 1.94.
    3. P(p^>0.43)=1Φ(1.94)0.026P(\hat p > 0.43) = 1 - \Phi(1.94) \approx 0.026.
    Trick: For sample proportion, σp^=p(1p)/n\sigma_{\hat p} = \sqrt{p(1-p)/n} — analogous to the standard error of the mean.
  23. Ex. 77.23Modeling

    You hold 50 stocks with daily return μ=0.1%\mu = 0.1\%, σ=2%\sigma = 2\% (independent). What is the distribution of your portfolio's mean daily return?

    Select the correct option
    Select an option first
    Show solution
    50 independent stocks with daily return μ=0.1%\mu = 0.1\%, σ=2%\sigma = 2\%. Rˉ50N(0.001,  (0.02)2/50)=N(0.001,  8×106)\bar R_{50} \approx \mathcal N(0.001,\; (0.02)^2/50) = \mathcal N(0.001,\; 8 \times 10^{-6}). σRˉ=2%/500.28%\sigma_{\bar R} = 2\%/\sqrt{50} \approx 0.28\%. Diversification reduced risk from 2% to 0.28%.
  24. Ex. 77.24Modeling

    ML model: individual error σ=0.5\sigma = 0.5. Calculate the standard deviation of mean error over 1,000 predictions and the 95% CI margin of error.

    Select the correct option
    Select an option first
    Show solution
    σeˉ=0.5/10000.01581\sigma_{\bar e} = 0.5/\sqrt{1000} \approx 0.01581. 95% CI: margin ±1.96×0.01581±0.031\pm 1.96 \times 0.01581 \approx \pm 0.031.
  25. Ex. 77.25Modeling

    Satisfaction survey: margin of error ±3%\pm 3\% at 95%, pp unknown. What is the minimum nn?

    Select the correct option
    Select an option first
    Show solution
    Margin of error 3% at 95%, worst case p=0.5p = 0.5: n=(1.96/0.03)2×0.25=65.332×0.251067n = (1.96/0.03)^2 \times 0.25 = 65.33^2 \times 0.25 \approx 1067. This is why "3% margin of error" Brazilian surveys have ~1000 respondents.
  26. Ex. 77.26Modeling

    Batch of 500 items: μ=100\mu = 100 g, σ=5\sigma = 5 g. Determine the distribution of S500S_{500}.

    Select the correct option
    Select an option first
    Show solution
    S500N(500×100,  500×25)=N(50.000,  12.500)S_{500} \approx \mathcal N(500 \times 100,\; 500 \times 25) = \mathcal N(50.000,\; 12.500). σS=12.500111.8\sigma_{S} = \sqrt{12.500} \approx 111.8 g. Expected total weight: 50 kg with standard deviation ~112 g.
  27. Ex. 77.27Modeling

    Bus wait time: U[0,30]\mathcal U[0, 30] min. Calculate P(Tˉ50>16)P(\bar T_{50} > 16) for 50 passengers.

    Select the correct option
    Select an option first
    Show solution
    U[0,30]\mathcal U[0, 30]: μ=15\mu = 15, σ2=75\sigma^2 = 75. Tˉ50N(15,  75/50=1.5)\bar T_{50} \approx \mathcal N(15,\; 75/50 = 1.5). z=(1615)/1.50.816z = (16-15)/\sqrt{1.5} \approx 0.816. P(Tˉ>16)0.207P(\bar T > 16) \approx 0.207.
  28. Ex. 77.28Modeling

    On an X-bar control chart with n=5n = 5, how are the control limits (UCL and LCL) calculated?

    Select the correct option
    Select an option first
    Show solution
    Shewhart X-bar chart: UCL/LCL = Xˉˉ±3σXˉ=Xˉˉ±3σ/n\bar{\bar X} \pm 3\sigma_{\bar X} = \bar{\bar X} \pm 3\sigma/\sqrt{n}. The 3 comes from P(Z>3)0.27%P(|Z| > 3) \approx 0.27\% — acceptable false-alarm rate.
    Show step-by-step (with the why)
    1. By CLT: Xˉ\bar X of subgroup of size nn is approximately N(μ,σ2/n)\mathcal N(\mu, \sigma^2/n).
    2. 3-sigma limits capture 99.73% of means when the process is in control.
    3. UCL = Xˉˉ+3σ/n\bar{\bar X} + 3\sigma/\sqrt{n}; LCL = Xˉˉ3σ/n\bar{\bar X} - 3\sigma/\sqrt{n}.
    4. For n=5n = 5: total width = 6σ/52.68σ6\sigma/\sqrt{5} \approx 2.68\sigma.
  29. Ex. 77.29ModelingAnswer key

    Survey with margin of error ±5%\pm 5\% at 95% confidence, pp unknown. What is the minimum necessary nn?

    Select the correct option
    Select an option first
    Show solution
    Margin = 1.96×p(1p)/n1.96 \times \sqrt{p(1-p)/n}. For p=0.5p = 0.5 (worst case) and margin = 5%: n=(1.96/0.05)2/4=39.22×0.25384.2n = (1.96/0.05)^2/4 = 39.2^2 \times 0.25 \approx 384.2. Round up: n=385n = 385. Half the margin (3%) requires almost three times more respondents (1067).
  30. Ex. 77.30ModelingAnswer key

    Call time: μ=3\mu = 3 min, σ=1.5\sigma = 1.5 min, 100 calls. Calculate P(total>330 min)P(\text{total} > 330 \text{ min}).

    Select the correct option
    Select an option first
    Show solution
    S100N(300,225)S_{100} \approx \mathcal N(300, 225). σS=15\sigma_S = 15. z=(330300)/15=2z = (330-300)/15 = 2. P(S>330)=1Φ(2)0.023P(S > 330) = 1 - \Phi(2) \approx 0.023.
  31. Ex. 77.31ModelingAnswer key

    A/B test: 10,000 visitors per variant; conversion A = 5%, B = 6%. Is the 1 p.p. lift statistically significant? Calculate zz and the p-value.

    Select the correct option
    Select an option first
    Show solution
    Observed difference: 1%. SE: 2×0.05×0.95/100000.00308\sqrt{2 \times 0.05 \times 0.95/10000} \approx 0.00308. z0.01/0.003083.25z \approx 0.01/0.00308 \approx 3.25. p-value < 0.001. A 1 p.p. difference is statistically significant with 10,000 users per variant.
    Show step-by-step (with the why)
    1. Pooled estimator: p^=(500+600)/20000=0.055\hat p = (500+600)/20000 = 0.055.
    2. SE: 0.055×0.945×(2/10000)0.00322\sqrt{0.055 \times 0.945 \times (2/10000)} \approx 0.00322.
    3. z=0.01/0.003223.11z = 0.01/0.00322 \approx 3.11. p-value 0.002\approx 0.002.
    4. Conclusion: reject H0H_0 — difference is significant.
  32. Ex. 77.32Understanding

    Which option correctly describes the CLT?

    Select the correct option
    Select an option first
    Show solution
    The CLT talks about the distribution of the sample mean (or standardized sum) — not the original distribution. The original distribution stays the same; what converges to normal is the statistic Zn=(Xˉnμ)/(σ/n)Z_n = (\bar X_n - \mu)/(\sigma/\sqrt{n}).
  33. Ex. 77.33Understanding

    Why does the classical CLT not apply to the Cauchy distribution?

    Select the correct option
    Select an option first
    Show solution
    The Cauchy characteristic function is ete^{-|t|} — no second derivative at t=0t = 0 (infinite variance). The Taylor expansion leading to the limit et2/2e^{-t^2/2} doesn't work. The mean of Cauchys is still Cauchy.
  34. Ex. 77.34UnderstandingAnswer key

    When simulating means of die rolls for n{1,5,30}n \in \{1, 5, 30\}, how does the histogram shape evolve?

    Select the correct option
    Select an option first
    Show solution
    Simulating means of fair die rolls: n=1n = 1 = flat uniform; n=5n = 5 = triangular-trapezoidal shape; n=30n = 30 = indistinguishable from normal. Convergence is visible and rapid.
  35. Ex. 77.35Application

    Distance of fly-balls: N(250,502)\mathcal N(250, 50^2) ft. Sample of 49 balls. Calculate P(Xˉ<240)P(\bar X < 240).

    Select the correct option
    Select an option first
    Show solution
    Baseball fly-balls: μ=250\mu = 250 ft, σ=50\sigma = 50. n=49n = 49. σXˉ=50/77.14\sigma_{\bar X} = 50/7 \approx 7.14. z=(240250)/7.141.4z = (240-250)/7.14 \approx -1.4. P(Xˉ<240)=Φ(1.4)0.081P(\bar X < 240) = \Phi(-1.4) \approx 0.081.
  36. Ex. 77.36ProofAnswer key

    Show that the CLT implies the Weak Law of Large Numbers.

    Select the correct option
    Select an option first
    Show solution
    From CLT: XˉnN(μ,σ2/n)\bar X_n \approx \mathcal N(\mu, \sigma^2/n). Chebyshev: P(Xˉnμ>ε)σ2/(nε2)0P(|\bar X_n - \mu| > \varepsilon) \leq \sigma^2/(n\varepsilon^2) \to 0. So XˉnPμ\bar X_n \xrightarrow{P} \mu — Law of Large Numbers. CLT is stronger: it tells how the distribution of the error behaves as it converges.
  37. Ex. 77.37Proof

    Sketch the proof of the CLT via characteristic function, indicating the necessary hypotheses.

    Select the correct option
    Select an option first
    Show solution
    Steps: (1) Define Yi=(Xiμ)/σY_i = (X_i - \mu)/\sigma. (2) φYi(t)=1t2/2+o(t2)\varphi_{Y_i}(t) = 1 - t^2/2 + o(t^2). (3) φZn(t)=[φYi(t/n)]net2/2\varphi_{Z_n}(t) = [\varphi_{Y_i}(t/\sqrt n)]^n \to e^{-t^2/2}. (4) Lévy's theorem concludes convergence in distribution. Requires only finite variance (Lindeberg-Lévy hypothesis).
  38. Ex. 77.38ApplicationAnswer key

    Average time to fill IRS Form 1040: μ=10.53\mu = 10.53 h, σ=2\sigma = 2 h. Sample of 36 taxpayers. Calculate the distribution of Xˉ\bar X and P(Xˉ>12)P(\bar X > 12).

    Select the correct option
    Select an option first
    Show solution
    IRS Form 1040 time: μ=10.53\mu = 10.53 h, σ=2\sigma = 2 h. n=36n = 36. σXˉ=2/6=1/3\sigma_{\bar X} = 2/6 = 1/3. z=(1210.53)/(1/3)=4.41z = (12-10.53)/(1/3) = 4.41. P(Xˉ>12)5×106P(\bar X > 12) \approx 5 \times 10^{-6} — practically impossible.
  39. Ex. 77.39Application

    Marathon times: μ=145\mu = 145 min, σ=14\sigma = 14 min. Sample of 49 runs. Calculate P(142<Xˉ<146)P(142 < \bar X < 146).

    Select the correct option
    Select an option first
    Show solution
    Marathon times: μ=145\mu = 145 min, σ=14\sigma = 14. n=49n = 49. σXˉ=14/7=2\sigma_{\bar X} = 14/7 = 2. z1=(142145)/2=1.5z_1 = (142-145)/2 = -1.5; z2=(146145)/2=0.5z_2 = (146-145)/2 = 0.5. P=Φ(0.5)Φ(1.5)0.69150.0668=0.6247P = \Phi(0.5) - \Phi(-1.5) \approx 0.6915 - 0.0668 = 0.6247.
  40. Ex. 77.40Modeling

    Percentage of fat calories: μ=36%\mu = 36\%, σ=10%\sigma = 10\%. Sample of 16 people. Calculate P(Xˉ>40%)P(\bar X > 40\%).

    Select the correct option
    Select an option first
    Show solution
    % fat calories: μ=36\mu = 36, σ=10\sigma = 10. n=16n = 16. σXˉ=10/4=2.5\sigma_{\bar X} = 10/4 = 2.5. z=(4036)/2.5=1.6z = (40-36)/2.5 = 1.6. P(Xˉ>40%)=1Φ(1.6)0.0548P(\bar X > 40\%) = 1 - \Phi(1.6) \approx 0.0548.
  41. Ex. 77.41Modeling

    Gasoline price (Bay Area): μ=4.59\mu = 4.59 dollars, σ=0.10\sigma = 0.10 dollars. Sample of 16 stations. What is the distribution of Xˉ\bar X and P(Xˉ>4.63)P(\bar X > 4.63)?

    Select the correct option
    Select an option first
    Show solution
    Gasoline Bay Area: μ=4.59\mu = 4.59, σ=0.10\sigma = 0.10. n=16n = 16. σXˉ=0.10/4=0.025\sigma_{\bar X} = 0.10/4 = 0.025. z=(4.634.59)/0.025=1.6z = (4.63-4.59)/0.025 = 1.6. P(Xˉ>4.63)0.05485.5%P(\bar X > 4.63) \approx 0.0548 \approx 5.5\%.
  42. Ex. 77.42ModelingAnswer key

    Air conditioning service: mean time 1 h, σ=1\sigma = 1 h. Budget of 1.1 h per technician, 70 units. Is there sufficient budget with 80% confidence?

    Select the correct option
    Select an option first
    Show solution
    Air conditioning service: μ=1\mu = 1 h, σ=1\sigma = 1 h. n=70n = 70 units. σXˉ=1/700.1195\sigma_{\bar X} = 1/\sqrt{70} \approx 0.1195. z=(1.11)/0.11950.837z = (1.1-1)/0.1195 \approx 0.837. P(Xˉ1.1)Φ(0.84)0.80P(\bar X \leq 1.1) \approx \Phi(0.84) \approx 0.80. 80% chance of not overrunning the budget.
  43. Ex. 77.43Challenge

    Challenge — Lindeberg condition. The Lindeberg-Feller CLT generalizes the CLT to independent but non-i.d. variables. What is the essential condition that replaces the identical-distribution hypothesis, and what does it intuitively require?

    Select the correct option
    Select an option first
    Show solution
    Let X1,X2,X_1, X_2, \ldots be independent (not necessarily i.d.) with E[Xi]=μiE[X_i]=\mu_i, Var(Xi)=σi2\mathrm{Var}(X_i)=\sigma_i^2. Define sn2=i=1nσi2s_n^2 = \sum_{i=1}^n \sigma_i^2. The Lindeberg condition requires that for every ε>0\varepsilon > 0: 1sn2i=1nE[(Xiμi)21Xiμi>εsn]0\frac{1}{s_n^2}\sum_{i=1}^n E\bigl[(X_i-\mu_i)^2 \mathbf{1}_{|X_i-\mu_i|>\varepsilon s_n}\bigr] \to 0. Intuitively: the contribution of any extreme individual value to total variance must be negligible. If satisfied, (Snμi)/sndN(0,1)(S_n - \sum\mu_i)/s_n \xrightarrow{d} N(0,1).
    Show step-by-step (with the why)
    1. Motivation. With i.i.d. variables, classical CLT follows easily from characteristic functions. For non-i.d. variables, we need to prevent a single variable from "dominating" the sum.
    2. Total variance. sn2=σ12+σ22++σn2s_n^2 = \sigma_1^2 + \sigma_2^2 + \cdots + \sigma_n^2. Should grow with nn.
    3. Lindeberg condition. For each ε>0\varepsilon > 0, the fraction of variance explained by large deviations (above εsn\varepsilon s_n) must go to zero.
    4. Practical consequence. If σi2/sn20\sigma_i^2 / s_n^2 \to 0 for each ii (Feller condition), the Lindeberg condition is satisfied for variables with finite fourth moment.
    5. Counterexample. If XnX_n concentrates all variance (e.g. σn2=sn2\sigma_n^2 = s_n^2), Lindeberg fails and CLT does not hold.

    The Lindeberg-Feller CLT is the most general CLT result: it generalizes i.i.d. requiring only that no individual variable "dominates" variance asymptotically.

  44. Ex. 77.44Challenge

    Challenge — Berry-Esseen bound. The Berry-Esseen theorem provides a quantitative bound for the speed of CLT convergence. For i.i.d. variables with third absolute moment ρ=E[Xi3]<\rho = E[|X_i|^3] < \infty, what is the correct form of the bound on the maximum approximation error?

    Select the correct option
    Select an option first
    Show solution
    The Berry-Esseen theorem states: if X1,,XnX_1,\ldots,X_n are i.i.d. with E[Xi]=0E[X_i]=0, σ2=Var(Xi)<\sigma^2 = \mathrm{Var}(X_i) < \infty and ρ=E[Xi3]<\rho = E[|X_i|^3] < \infty, then supxP(Znx)Φ(x)Cρ/(σ3n)\sup_x |P(Z_n \leq x) - \Phi(x)| \leq C\rho/(\sigma^3\sqrt{n}), with universal constant C0.4748C \leq 0.4748. Example: for Bernoulli(1/2)\mathrm{Bernoulli}(1/2): σ2=1/4\sigma^2 = 1/4, ρ=1/8\rho = 1/8, error 0.4748(1/8)/((1/2)3n)=0.4748/n\leq 0.4748 \cdot (1/8)/((1/2)^3 \sqrt{n}) = 0.4748/\sqrt{n}.
    Show step-by-step (with the why)
    1. Statement. Berry-Esseen: supxFn(x)Φ(x)Cρ/(σ3n)\sup_x |F_n(x) - \Phi(x)| \leq C\rho/(\sigma^3\sqrt{n}) with C0.4748C \leq 0.4748 (best constant known in 2024).
    2. Ingredients. ρ=E[Xi3]\rho = E[|X_i|^3] (third absolute moment); σ3\sigma^3 normalizes by scale.
    3. Convergence rate. The error is O(1/n)O(1/\sqrt{n}) — same rate as standard deviation of the sample mean.
    4. Application Bernoulli(1/2). σ=1/2\sigma = 1/2, ρ=E[X3]=1/2\rho = E[|X|^3] = 1/2 (since X3=X|X|^3 = |X|). Error 0.4748/(0.125/0.125n)=0.4748/n\leq 0.4748/(0.125/0.125\cdot\sqrt{n}) = 0.4748/\sqrt{n}. For n=100n = 100: error 0.047\leq 0.047 — reasonable.

    Berry-Esseen quantifies the "n=30 rule": the maximum approximation error decreases as 1/n1/\sqrt{n}, but the constant depends on distribution skewness (via ρ/σ3\rho/\sigma^3). More skewed distributions need larger nn.

  45. Ex. 77.45Challenge

    Challenge — Shewhart control chart. In statistical process control (SPC), limits μ±3σ/n\mu \pm 3\sigma/\sqrt{n} are used to detect out-of-control processes. Based on the CLT, calculate the probability of a false alarm (process in control but point outside limits) and the ARL₀ (average number of subgroups until the first false alarm).

    Select the correct option
    Select an option first
    Show solution
    By Shewhart's 3-sigma rule: control limits are μ±3σ/n\mu \pm 3\sigma/\sqrt{n}. By CLT, XˉN(μ,σ2/n)\bar X \sim N(\mu, \sigma^2/n) when the process is in control. Therefore P(Xˉμ>3σ/n)=P(Z>3)0.0027P(|\bar X - \mu| > 3\sigma/\sqrt{n}) = P(|Z| > 3) \approx 0.0027, or 0.27% — a false alarm occurs about once per 370 subgroups. The average number of subgroups until the first false alarm (ARL₀) is 1/0.00273701/0.0027 \approx 370.
    Show step-by-step (with the why)
    1. Shewhart control chart. Each point is the mean Xˉ\bar X of a subgroup of size nn. Upper control limit (UCL) = μ+3σ/n\mu + 3\sigma/\sqrt{n}; lower control limit (LCL) = μ3σ/n\mu - 3\sigma/\sqrt{n}.
    2. CLT justifies using N(μ,σ2/n)N(\mu, \sigma^2/n). Regardless of individual distribution (when n4n \geq 455 for non-highly-skewed processes).
    3. Probability of false alarm. P(Z>3)=2(1Φ(3))20.00135=0.0027P(|Z| > 3) = 2(1 - \Phi(3)) \approx 2 \cdot 0.00135 = 0.0027.
    4. ARL₀. Average number of subgroups until first false alarm = 1/0.00273701/0.0027 \approx 370. The higher the ARL₀, the less the chart "cries wolf".
    5. Why 3σ? Trade-off: false alarm rare (0.27%) without excessive ARL₀; 2-sigma limits would give 4.5% false alarm, too frequent.

    Shewhart's 3-sigma rule (1924) is the industrial application of CLT. The ARL₀ ≈ 370 is the gold standard in statistical process control (SPC).

Sources

  • OpenIntro Statistics (4th ed) — Diez, Çetinkaya-Rundel, Barr · 2019 · CC-BY-SA. Primary source for exercises 77.2, 77.4, 77.8, 77.11, 77.14–17, 77.22–23, 77.25–26, 77.28, 77.30, 77.33–34.
  • OpenStax Statistics — Illowsky, Dean · 2022 · CC-BY. Source for exercises 77.1, 77.3, 77.5–7, 77.9–10, 77.12–13, 77.18–19, 77.21, 77.24, 77.27, 77.29, 77.31, 77.35 and examples 1–3.
  • Introduction to Probability (Grinstead-Snell) — Grinstead, Snell · 2006 · CC Free · §9 "Central Limit Theorem". Source for challenge exercises 77.43–77.44.
  • Introduction to Probability (Grinstead-Snell) — Grinstead, Snell · Dartmouth · GNU FDL. Source for exercises 77.19–20, 77.26, 77.36–37 and example 5.

Updated on 2026-05-06 · Author(s): Clube da Matemática

Found an error? Open an issue on GitHub or submit a PR — open source forever.