Tutorials › Statistical Theory I - Probability and Distribution Theory › Centering and Scaling Sums for Normal Limits

Central Limit Theorems · Tutorial 846 of 994

Centering and Scaling Sums for Normal Limits

Build the correctly centered, unit-variance version of an independent sum, and calculate its normalization in iid, weighted, and triangular-array settings.

Advanced 10 min read

What You'll Learn

  • Define the centered, variance-standardized form of an independent sum.
  • Derive its mean and variance from independence and finite second moments.
  • Calculate the normalization for sums of iid variables.
  • Handle negative and unequal weights without losing the correct variance.
  • Normalize a row of a triangular array using row-specific means and variances.
  • Distinguish choosing a normalization from establishing a normal limit.

Why Centering and Scaling Come First

A normal-limit calculation starts by identifying the location and spread of the quantity whose distribution is being studied. For a sum of independent variables, its expectation sets the natural center, while its total standard deviation sets the natural scale. These are determined by the summands’ distributions and the way they enter the sum; they are not chosen simply because a familiar formula uses \(\sqrt{n}\).

This tutorial develops that normalization and calculates it in several settings. The result is a random variable with mean zero and variance one, provided the sum has finite, positive variance. That property is an exact finite-sample calculation. It does not, by itself, say that the resulting distribution is normal or that it converges to a normal distribution. A central limit theorem supplies such a convergence conclusion only under additional conditions.

The Exact Normalization for an Independent Sum

Definition: Let \(X_1,\ldots,X_n\) be independent random variables with finite means \(\mu_i=E[X_i]\) and finite variances \(\sigma_i^2=\operatorname{Var}(X_i)\). Write \(S_n=\sum_{i=1}^n X_i\). Its centered, variance-standardized form is
$$ Z_n=\frac{S_n-\sum_{i=1}^n\mu_i}{\sqrt{\sum_{i=1}^n\sigma_i^2}}, $$ provided \(\sum_{i=1}^n\sigma_i^2>0\).

The numerator subtracts the sum’s expectation. The denominator is the square root of the sum’s variance, which is the sum of the individual variances under independence. A standard deviation is a square root of a variance, so variances—not standard deviations—are added before taking the square root.

Theorem: Under the definition’s assumptions, \(E[Z_n]=0\) and \(\operatorname{Var}(Z_n)=1\).

Proof. Linearity of expectation gives \(E[S_n]=\sum_{i=1}^n\mu_i\), so the numerator has expectation zero. To calculate its variance, center each summand and expand:

$$ \operatorname{Var}(S_n) =E\left[\left(\sum_{i=1}^n(X_i-\mu_i)\right)^2\right] =\sum_{i=1}^n E[(X_i-\mu_i)^2] +2\sum_{1\leq i<j\leq n} E[(X_i-\mu_i)(X_j-\mu_j)]. $$

For \(i\ne j\), independence implies that the centered variables are independent. Their product is integrable, since \(E[|(X_i-\mu_i)(X_j-\mu_j)|]=E[|X_i-\mu_i|]E[|X_j-\mu_j|]<\infty\). Thus the cross-term expectation factors into two zero expectations. It follows that \(\operatorname{Var}(S_n)=\sum_{i=1}^n\sigma_i^2\). Dividing the centered sum by the square root of this positive variance gives \(E[Z_n]=0\) and \(\operatorname{Var}(Z_n)=1\). \(\square\)

Independence is a sufficient condition for the variance sum above; without it, covariances may contribute. In that case, the variance of a finite sum is \(\sum_i\operatorname{Var}(X_i)+2\sum_{i<j}\operatorname{Cov}(X_i,X_j)\), whenever the second moments exist.

Identically Distributed Summands

If \(X_1,\ldots,X_n\) are iid with \(E[X_i]=\mu\) and \(\operatorname{Var}(X_i)=\sigma^2\), where \(0<\sigma^2<\infty\), then the sum has mean \(n\mu\) and variance \(n\sigma^2\). Its standard deviation is \(\sqrt{n}\sigma\), so the normalization becomes

$$ Z_n=\frac{\sum_{i=1}^n(X_i-\mu)}{\sigma\sqrt{n}} =\frac{\sqrt{n}(\bar X_n-\mu)}{\sigma}, \qquad \bar X_n=\frac{1}{n}\sum_{i=1}^n X_i. $$

The two forms are equal because \(\sum_i(X_i-\mu)=n(\bar X_n-\mu)\). This algebra explains the \(\sqrt{n}\) in the normalization of an iid sum or sample mean: the total variance grows as \(n\), and the total standard deviation therefore grows as \(\sqrt{n}\). The formula uses the population standard deviation \(\sigma\), not a sample estimate.

Worked Example: An IID Exponential Sum

Let \(X_1,\ldots,X_{25}\) be iid exponential random variables with rate \(2\), so \(f(x)=2e^{-2x}\) for \(x>0\). Find the centered, variance-standardized sum.

An exponential variable with rate \(2\) has mean \(\mu=1/2\) and variance \(\sigma^2=1/4\). Consequently,

$$ E[S_{25}]=25\left(\frac12\right)=\frac{25}{2}, \qquad \operatorname{Var}(S_{25})=25\left(\frac14\right)=\frac{25}{4}, \qquad \sqrt{\operatorname{Var}(S_{25})}=\frac52. $$

Therefore the required normalization is

$$ Z_{25}=\frac{\sum_{i=1}^{25}X_i-25/2}{5/2}. $$

Equivalently, the numerator is \(\sum_i(X_i-1/2)\), and dividing by \(5/2\) gives \(Z_{25}=\sum_i(X_i-1/2)/(5/2)\). The calculation fixes its mean at zero and its variance at one; it does not assert that this finite-sample variable is exactly normal.

Weighted Sums: Coefficients Affect Both Center and Spread

For a weighted sum \(T=\sum_{i=1}^n a_iX_i\) of independent variables with finite means and variances, its expectation is \(\sum_i a_i\mu_i\), and its variance is \(\sum_i a_i^2\sigma_i^2\). The square on each coefficient matters: a negative weight reverses the contribution to the mean but still contributes nonnegatively to the variance.

Result: If \(\sum_{i=1}^n a_i^2\sigma_i^2>0\), the centered, variance-standardized weighted sum is
$$ \frac{\sum_{i=1}^n a_iX_i-\sum_{i=1}^n a_i\mu_i} {\sqrt{\sum_{i=1}^n a_i^2\sigma_i^2}}. $$

To verify the variance formula, center the weighted summands and expand the square. Each diagonal term has expectation \(a_i^2E[(X_i-\mu_i)^2]=a_i^2\sigma_i^2\). Each cross-term has expectation zero by independence and centering, just as in the theorem above. The weighted sum’s expectation follows directly from linearity. Hence the displayed quotient has mean zero and variance one.

For iid \(X_i\) with mean \(\mu\) and variance \(\sigma^2\), this specializes to

$$ \frac{\sum_{i=1}^n a_iX_i-\mu\sum_{i=1}^n a_i} {\sigma\sqrt{\sum_{i=1}^n a_i^2}}, $$

as long as \(\sigma^2\sum_i a_i^2>0\). In general this is not the same as the equal-weight normalization: the sum of squared weights, not the number of terms alone, determines the variance.

Worked Example: A Weighted IID Bernoulli Sum

Let \(X_1,\ldots,X_4\) be iid Bernoulli variables with \(P(X_i=1)=0.3\) and \(P(X_i=0)=0.7\). Take weights \(a_1=1\), \(a_2=2\), \(a_3=-1\), and \(a_4=3\). Find the normalized weighted sum.

For each summand, \(\mu=0.3\) and \(\sigma^2=0.3(0.7)=0.21\). The sum of weights is \(1+2-1+3=5\), so the weighted sum has expectation \(0.3(5)=1.5\). The sum of squared weights is \(1+4+1+9=15\), giving

$$ \operatorname{Var}(T)=0.21(15)=3.15, \qquad \operatorname{SD}(T)=\sqrt{3.15}. $$

Thus

$$ Z=\frac{X_1+2X_2-X_3+3X_4-1.5}{\sqrt{3.15}}. $$

The negative coefficient on \(X_3\) reduces the weighted mean by \(0.3\), as reflected in the sum of weights. Its variance contribution remains \(1^2(0.21)=0.21\); omitting that contribution or giving it a negative sign would be incorrect.

Triangular Arrays: Normalize Each Row Separately

In a triangular array, the variables in row \(n\) are written \(X_{n,1},\ldots,X_{n,k_n}\). The row length \(k_n\) may change with \(n\), and the distributions may vary across both indices. Suppose the variables within each row are independent, each has a finite mean \(\mu_{n,k}\) and finite variance \(\sigma_{n,k}^2\), and define

$$ S_n=\sum_{k=1}^{k_n}X_{n,k}, \qquad m_n=\sum_{k=1}^{k_n}\mu_{n,k}, \qquad s_n^2=\sum_{k=1}^{k_n}\sigma_{n,k}^2. $$

When \(s_n^2>0\), the row’s normalization is \((S_n-m_n)/s_n\). Both center and scale are indexed by \(n\), because they can change from row to row. Independence is needed only among the variables within each row for this variance calculation; dependence between different rows does not change the normalization of any single row.

Worked Example: Rows of Bernoulli Variables With Changing Probabilities

For row \(n\), let \(X_{n,k}\), \(k=1,\ldots,n\), be independent Bernoulli variables with \(P(X_{n,k}=1)=k/(n+1)\). Find the row normalization, and check it explicitly for \(n=4\).

The mean of the row sum is

$$ m_n=\sum_{k=1}^n\frac{k}{n+1} =\frac{1}{n+1}\frac{n(n+1)}{2} =\frac n2. $$

Each variance is \(p_{n,k}(1-p_{n,k})\), where \(p_{n,k}=k/(n+1)\). Thus

$$ s_n^2 =\sum_{k=1}^n\frac{k}{n+1}\left(1-\frac{k}{n+1}\right) =\frac{(n+1)\sum_{k=1}^n k-\sum_{k=1}^n k^2}{(n+1)^2}. $$

Using \(\sum_{k=1}^n k=n(n+1)/2\) and \(\sum_{k=1}^n k^2=n(n+1)(2n+1)/6\), simplify:

$$ s_n^2 =\frac{n(n+1)\{3(n+1)-(2n+1)\}}{6(n+1)^2} =\frac{n(n+2)}{6(n+1)}. $$

This is positive for every positive integer \(n\). Therefore the normalized row sum is

$$ \frac{\sum_{k=1}^n X_{n,k}-n/2} {\sqrt{n(n+2)/(6(n+1))}}. $$

For \(n=4\), the success probabilities are \(1/5,2/5,3/5,4/5\). Their sum is \(2\), and the variances are \(4/25,6/25,6/25,4/25\), which sum to \(20/25=4/5\). The general formula gives \(4(6)/(6(5))=4/5\) as well. The normalized sum for that row is therefore \((S_4-2)/\sqrt{4/5}\).

What the Normalization Does—and Does Not Do

Subtracting the exact expectation and dividing by the exact standard deviation creates a common scale for sums whose raw magnitudes may differ. It is a bookkeeping step essential to studying normal limits: without it, changing the mean or variance can obscure the behavior of the distribution. But standardizing does not make an arbitrary variable normal, and mean zero with variance one does not uniquely determine a distribution.

For instance, if \(X\) is Bernoulli with success probability \(1/2\), then \(E[X]=1/2\) and \(\operatorname{Var}(X)=1/4\). Its standardized version is \((X-1/2)/(1/2)\), which takes only the values \(-1\) and \(1\). It has mean zero and variance one, but it is not a normal variable. Normal approximation requires a limit result and its hypotheses, not merely this moment calculation.

  • Common mistake: adding standard deviations. Independence gives \(\operatorname{Var}(S_n)=\sum_i\sigma_i^2\), so the denominator is \(\sqrt{\sum_i\sigma_i^2}\), not \(\sum_i\sigma_i\).
  • Common mistake: using \(n\) as the scale in every setting. For iid terms the standard deviation is \(\sigma\sqrt{n}\); for weighted terms it is \(\sigma\sqrt{\sum_i a_i^2}\); for a triangular array it is the square root of that row’s total variance.
  • Common mistake: treating a negative weight as a negative variance contribution. The mean uses \(a_i\), while the variance uses \(a_i^2\). The latter cannot be negative.
  • Common mistake: confusing exact standardization with a normal limit. Unit variance is an identity under the stated assumptions. A normal-limit claim needs additional distributional conditions.
Key takeaway: For an independent sum with finite, positive total variance, subtract the sum of the means and divide by the square root of the sum of the variances. In weighted and triangular-array settings, calculate the coefficients or row-specific moments before choosing the scale.

Check Your Understanding

Use the exact expectations and variances to answer each question.

  1. Independent variables \(X_1,X_2\) have means \(2,-1\) and variances \(4,9\). Write the centered, variance-standardized sum.
  2. For iid variables with variance \(\sigma^2>0\), explain why the standard deviation of their sum is \(\sigma\sqrt{n}\), rather than \(n\sigma\).
  3. If \(T=3X_1-2X_2\), where \(X_1,X_2\) are independent with variances \(1\) and \(5\), calculate \(\operatorname{Var}(T)\).
  4. For a triangular array, why can the center and scale both depend on the row index \(n\)?
  5. Does a mean-zero, variance-one normalization guarantee a normal distribution? Give a brief reason.