Tutorials › Introduction to Statistics › Choosing Sample Size for a Reliable CLT Approximation

The Central Limit Theorem · Tutorial 545 of 996

Choosing Sample Size for a Reliable CLT Approximation

Use the population and the probability you care about to judge whether a CLT approximation is plausible, and know when to check it more directly.

Intermediate 11 min read

What You'll Learn

  • Explain why no single sample-size cutoff guarantees an accurate CLT approximation.
  • Use population skewness as a rough guide to how quickly a sample mean may become more symmetric.
  • Distinguish central-probability targets from extreme-tail targets.
  • Assess how rare outliers and heavy tails can affect the needed sample size.
  • Recognize why dependent observations cannot be treated as the same amount of information as independent observations.
  • Choose when to use exact calculations or simulation to check a proposed approximation.

Why There Is No Universal Sample-Size Cutoff

The Central Limit Theorem (CLT) explains why the distribution of a sample mean can be approximated by a normal distribution when the sample is sufficiently large and the observations meet appropriate conditions. But “sufficiently large” does not name one sample size that works for every population and every question. A sample mean from a nearly symmetric, well-behaved population may be close to normal at a modest \(n\), while a mean from a strongly skewed population with rare extreme values may need much more data.

Accuracy also depends on what probability is needed. An approximation that is useful for a central probability may be poor for a small upper-tail probability. Choosing \(n\) therefore begins with the population and the specific probability region of interest, not with a rule such as “\(n\) must be at least 30.”

Definition: A CLT approximation for a sample mean is reliable for a stated purpose when its error is acceptably small for the particular probability or range being estimated. To make “acceptably small” precise, specify a target event and, when possible, a tolerable probability error. Reliability for one event does not guarantee reliability for every other event.

For independent, identically distributed observations \(X_1,\ldots,X_n\) with population mean \(\mu\) and finite, positive population variance \(\sigma^2\), the sample mean is \(\bar{x}\). The CLT says that the standardized sample mean approaches a standard normal distribution as \(n\) grows:

$$ \frac{\bar{x}-\mu}{\sigma/\sqrt{n}} \ \text{approaches a standard normal distribution as } n \text{ grows.} $$

This describes a long-run limit, not a guarantee about the accuracy at a particular \(n\). In practical use, judge the approximation using the information available about the population and the event of interest. If the sampling design involves dependence, the independent-observation version of the formula may not apply.

Use Population Shape as a First Guide

A roughly symmetric population without unusually influential observations gives the sample mean a head start: individual observations are not pulling the mean strongly to one side, and the mean is less likely to be dominated by a few extremes. A strongly skewed population typically presents a harder problem. Its sample means can remain skewed at sample sizes that would be adequate for a more symmetric population.

Under the independent, identically distributed model, if the population has finite third absolute moment, its standardized skewness is \(\gamma_1=E[(X-\mu)^3]/\sigma^3\). The standardized skewness of the sample mean is \(\gamma_1/\sqrt{n}\). This relationship is a useful guide to how skewness changes with \(n\), but it is not a formula for the probability error of a normal approximation.

The scaling follows directly from the properties of sums. Write \(Y_i=X_i-\mu\), so that \(E[Y_i]=0\). Independence makes the third central moment of the sum equal to the sum of the individual third moments. Dividing by the standard deviation of that sum, \(\sigma\sqrt{n}\), cubed gives:

$$ \frac{E[(Y_1+\cdots+Y_n)^3]}{(\sigma\sqrt{n})^3} = \frac{nE[Y_1^3]}{n^{3/2}\sigma^3} = \frac{\gamma_1}{\sqrt{n}}. $$

Thus, larger \(n\) reduces the standardized skewness at a rate proportional to \(1/\sqrt{n}\), under the stated assumptions. It does not follow that a particular skewness value guarantees an accurate approximation: skewness alone does not describe every feature of a population’s shape. In particular, a symmetric population can still have rare extreme values or several distinct clusters.

Worked Example: Skewness as a Guide, Not a Guarantee

Suppose observations come independently from an exponential population with mean 10 minutes. For an exponential distribution, the population standard deviation is 10 minutes and the standardized skewness is 2. A researcher wants to use the CLT for the sample mean completion time.

For \(n=25\), the standardized skewness of the sample mean is:

$$ \frac{2}{\sqrt{25}}=0.4. $$

For \(n=100\), it is:

$$ \frac{2}{\sqrt{100}}=0.2. $$

The larger sample has a smaller standardized skewness, so it is a more promising candidate for a normal approximation. The standard error of the sample mean at \(n=100\) is \(10/\sqrt{100}=1\) minute. For example, a cutoff of 12.5 minutes is 2.5 standard errors above the population mean. That is an upper-tail target, so a good central approximation would not by itself establish that the tail probability is accurate.

The model checks matter: the calculation assumes independent observations from the same exponential population, with the stated finite mean, variance, and third moment. If the researcher needs an accurate probability above 12.5 minutes, \(n=100\) should be treated as a candidate, not certified from the skewness calculation. An exact calculation, if available, or a simulation from the specified population can check the target probability directly.

Match the Judgment to the Probability Region

A normal approximation is often used to estimate \(P(\bar{x}>c)\), \(P(\bar{x}<c)\), or a probability that \(\bar{x}\) lies between two values. Standardize the endpoint or endpoints using the model’s mean and standard error. Their positions show whether the requested probability is central or in a tail:

$$ z=\frac{c-\mu}{\sigma/\sqrt{n}}. $$

An endpoint near \(z=0\) is near the center of the sampling distribution. A large positive or negative \(z\) corresponds to a more distant tail. When a target probability is very small, even a modest error in the estimated tail can be important relative to that probability. A normal curve that looks plausible in the center may still give a poor estimate for a far-tail event.

The cutoff also matters in relation to the population’s support. A normal distribution assigns positive probability to every real number, while a population may be bounded or otherwise restricted. If a normal approximation puts substantial probability in a region that the actual sample mean cannot reach, that is a warning sign. Even when the normal probability outside the possible range is small, it can matter if the requested event itself is very rare.

Worked Example: A Bounded Population and an Extreme Tail

Let \(X_1,\ldots,X_4\) be independent observations from a uniform population on the interval from 0 to 1. The population is symmetric about \(\mu=0.5\), and its standard deviation is \(\sigma=1/\sqrt{12}\). Consider the probability that \(\bar{x}>0.9\).

The standard error for \(n=4\) is:

$$ SE=\frac{1/\sqrt{12}}{\sqrt{4}}=\frac{1}{\sqrt{48}}\approx 0.1443. $$

The standardized cutoff is:

$$ z=\frac{0.9-0.5}{1/\sqrt{48}}\approx 2.771. $$

The normal approximation therefore gives \(P(\bar{x}>0.9)\approx P(Z>2.771)\approx 0.00279\), where \(Z\) is standard normal. But the exact probability can be calculated for this small sample. The event \(\bar{x}>0.9\) is equivalent to \(X_1+\cdots+X_4>3.6\). By symmetry of the four uniforms about their sum of 2, this has the same probability as \(X_1+\cdots+X_4<0.4\). For a sum below 0.4, all four observations are below 1. The possible values form a four-dimensional simplex with volume \(0.4^4/4!\), and the joint density is 1. Thus:

$$ P(\bar{x}>0.9)=\frac{0.4^4}{24}\approx 0.001067. $$

The normal estimate is more than twice the exact probability. The observations are independent and come from a symmetric, bounded population, but \(n=4\) is small and the target is an extreme tail. Symmetry by itself did not make this tail approximation accurate. For this known population, the exact probability is preferable; in less tractable settings, simulation can be used to assess the approximation for the event of interest.

Rare Outliers Can Make Large Samples Necessary

A population can have finite variance and still contain rare values far from its center. These observations may occur infrequently but have a large influence on both the shape of the population and the behavior of sample means. The challenge is that a sample can be large yet contain few or none of the rare observations. Increasing \(n\) helps the CLT in the limit, but it does not make every moderate sample automatically representative of such a population.

Worked Example: A Rare Component in a Mixture Population

Imagine an independent population generated as follows: with probability 0.99, an observation is drawn from a normal distribution with mean 0 and standard deviation 1; with probability 0.01, it is drawn from a normal distribution with mean 0 and standard deviation 100. This is a hypothetical mixture population with rare, very large values.

Its population mean is 0. Its variance is the weighted average of the two component variances:

$$ \sigma^2=0.99(1^2)+0.01(100^2)=100.99, \qquad \sigma\approx 10.05. $$

In a sample of 100 observations, the probability that none comes from the rare component is:

$$ P(\text{no rare-component observations})=(0.99)^{100}\approx 0.3660. $$

So, even at \(n=100\), there is about a 36.6% chance that the sample contains no observation from the rare component. At \(n=500\), that chance is \((0.99)^{500}\approx 0.006570\). The larger sample is much more likely to encounter the rare component, but these calculations do not prove that a normal approximation is accurate at \(n=500\). They show why population information about rare values matters to the sample-size decision.

The iid assumption holds by construction, and the mixture has finite variance. If the population model is known, simulate sample means at candidate values of \(n\), then compare simulated probabilities with normal probabilities for the actual target events. If the rare component has been overlooked or is poorly characterized, a large sample from the main component alone may give a misleading picture of the population.

Account for Dependence Before Increasing \(n\)

The basic CLT approximation for \(\bar{x}\) above assumes independent observations. With dependent observations, the variance of a sum includes covariance terms. Consequently, the standard error need not be \(\sigma/\sqrt{n}\), and treating the observations as independent can make the sample mean appear more precise than it is. Increasing the number of measurements does not necessarily add the same amount of information as increasing the number of independent observations.

Worked Example: Repeated Measurements With Positive Correlation

Suppose a hypothetical process produces 300 measurements with common variance \(\sigma^2\). Assume the correlation between observations \(k\) steps apart is \(0.5^k\). For this dependence model, the variance of the mean can be written as:

$$ \operatorname{Var}(\bar{x}) = \frac{\sigma^2}{n} \left[ 1+2\sum_{k=1}^{n-1} \left(1-\frac{k}{n}\right)0.5^k \right]. $$

For large \(n\), the bracketed quantity is approximately \(1+2\sum_{k=1}^{\infty}0.5^k\). Since the geometric series sums to 1, this is approximately 3. Thus, for \(n=300\), the variance of the mean is approximately \(3\sigma^2/300=\sigma^2/100\). Its standard error is approximately \(\sigma/10\), the same as the standard error from 100 independent observations with variance \(\sigma^2\).

This comparison is approximate and uses the stated correlation model; it is not a universal conversion from dependent observations to independent ones. The model’s assumptions must be checked, and a suitable dependent-data analysis may be needed. The practical point is that the nominal sample size of 300 alone does not establish that an independent-observation CLT approximation is appropriate.

A Practical Sample-Size Workflow

Use the following sequence to make a defensible judgment. It does not produce a universal cutoff; it makes clear what evidence supports a proposed \(n\) and where uncertainty remains.

1
State the target probability.
Write down the event to be approximated, such as a central interval or a specified upper-tail probability. A claim that the approximation is “good” is incomplete without a target.
2
Describe the population and sampling design.
Consider symmetry, skewness, boundedness, rare outliers, and whether observations are independent. Use a variance formula appropriate to the design.
3
Assess the sample size as a candidate.
Strong skewness and influential tails suggest that greater \(n\) may be needed. A reduction in standardized skewness as \(1/\sqrt{n}\) is informative, but it cannot certify probability accuracy.
4
Check the relevant region.
Standardize the event boundaries. Treat small tail probabilities and cutoffs near the limits of the population’s support with particular care.
5
Validate when accuracy matters.
Use an exact calculation if available, or simulate from a credible population and sampling model. Compare the probability of the same event under the simulated sample means and the normal approximation.

Common Mistakes and Why This Judgment Matters

The goal is not to find a sample size that guarantees a normal distribution. For a finite sample, the sampling distribution need not be exactly normal. The practical question is whether the approximation is adequate for the intended calculation, given the shape of the population, the dependence structure, and the probability region.

  • Using a universal cutoff. “The sample size is above 30, so the CLT approximation is reliable” ignores population shape, outliers, dependence, and the target event. A cutoff may be a rough classroom convention, but it is not a general guarantee.
  • Treating symmetry as proof of normality. A symmetric population can still be bounded, multimodal, or affected by rare extreme values. Symmetry alone does not settle the accuracy of a finite-sample approximation.
  • Checking only a central probability. A central approximation can look good while a small tail probability is inaccurate. Check the event that the analysis actually needs.
  • Counting dependent observations as independent. The nominal \(n\) does not justify using an independent-sample standard error when observations are correlated.

Sample-size judgment matters because a probability calculation can appear precise while relying on an approximation that is poor for the question being asked. When the population is well understood, its shape and the sampling design guide the choice. When the target is a small tail, the population has rare extremes, or dependence is present, exact methods or simulation can provide a more relevant check. The central lesson is to justify \(n\) in context, rather than treating a single threshold as a rule.

Check Your Understanding

Use the population, sampling design, and probability target to explain each judgment.

  1. Why does the relationship \(\gamma_1/\sqrt{n}\) help assess skewness without guaranteeing an accurate normal approximation?
  2. A researcher has a roughly symmetric population but wants a very small upper-tail probability. What additional concern should shape the sample-size judgment?
  3. In the mixture example, what is the probability that a sample of 100 observations contains no observation from the rare component?
  4. Why can a bounded population make an extreme-tail normal approximation questionable even when the population is symmetric?
  5. What information about dependence is needed before using \(\sigma/\sqrt{n}\) as the standard error of a sample mean?