Standard error



The standard error (SE) of a statistic (usually an estimate of a parameter) is the standard deviation of its sampling distribution or an estimate of that standard deviation. If the statistic is the sample mean, it is called the standard error of the mean (SEM). The standard error is a key ingredient in producing confidence intervals.

The sampling distribution of a mean is generated by repeated sampling from the same population and recording of the sample means obtained. This forms a distribution of different means, and this distribution has its own mean and variance. Mathematically, the variance of the sampling mean distribution obtained is equal to the variance of the population divided by the sample size. This is because as the sample size increases, sample means cluster more closely around the population mean.

Therefore, the relationship between the standard error of the mean and the standard deviation is such that, for a given sample size, the standard error of the mean equals the standard deviation divided by the square root of the sample size. In other words, the standard error of the mean is a measure of the dispersion of sample means around the population mean.

In regression analysis, the term "standard error" refers either to the square root of the reduced chi-squared statistic or the standard error for a particular regression coefficient (as used in, say, confidence intervals).

Exact value
Suppose a statistically independent sample of $$n$$ observations $$ x_1, x_2, \ldots, x_n $$ is taken from a statistical population with a standard deviation of $$\sigma$$. The mean value calculated from the sample, $$\bar{x}$$, will have an associated standard error on the mean, $${\sigma}_\bar{x}$$, given by: $${\sigma}_\bar{x} = \frac{\sigma}{\sqrt{n}}.$$

Practically this tells us that when trying to estimate the value of a population mean, due to the factor $$1/\sqrt{n}$$, reducing the error on the estimate by a factor of two requires acquiring four times as many observations in the sample; reducing it by a factor of ten requires a hundred times as many observations.

Estimate
The standard deviation $$\sigma$$ of the population being sampled is seldom known. Therefore, the standard error of the mean is usually estimated by replacing $$\sigma$$ with the sample standard deviation $$\sigma_{x}$$ instead: $${\sigma}_\bar{x}\ \approx \frac{\sigma_{x}}{\sqrt{n}}.$$

As this is only an estimator for the true "standard error", it is common to see other notations here such as: $$\widehat{\sigma}_{\bar{x}} := \frac{\sigma_{x}}{\sqrt{n}} \qquad \text{ or } \qquad {s}_\bar{x}\ := \frac{s}{\sqrt{n}}.$$

A common source of confusion occurs when failing to distinguish clearly between:


 * the standard deviation of the population ($$\sigma$$),
 * the standard deviation of the sample ($$\sigma_{x}$$),
 * the standard deviation of the mean itself ($$\sigma_{\bar{x}}$$, which is the standard error), and
 * the estimator of the standard deviation of the mean ($$\widehat{\sigma}_{\bar{x}}$$, which is the most often calculated quantity, and is also often colloquially called the standard error).

Accuracy of the estimator
When the sample size is small, using the standard deviation of the sample instead of the true standard deviation of the population will tend to systematically underestimate the population standard deviation, and therefore also the standard error. With n = 2, the underestimate is about 25%, but for n = 6, the underestimate is only 5%. Gurland and Tripathi (1971) provide a correction and equation for this effect. Sokal and Rohlf (1981) give an equation of the correction factor for small samples of n < 20. See unbiased estimation of standard deviation for further discussion.

Derivation
The standard error on the mean may be derived from the variance of a sum of independent random variables, given the definition of variance and some properties thereof. If $$ x_1, x_2, \ldots, x_n $$ is a sample of $$n$$ independent observations from a population with mean $$ \bar{x} $$ and standard deviation $$ \sigma $$, then we can define the total $$ T = (x_1 + x_2 + \cdots + x_n) $$ which due to the Bienaymé formula, will have variance $$ \operatorname{Var}(T) = \big(\operatorname{Var}(x_1) + \operatorname{Var}(x_2) + \cdots + \operatorname{Var}(x_n)\big) = n\sigma^2. $$ where we've approximated the standard deviations, i.e., the uncertainties, of the measurements themselves with the best value for the standard deviation of the population. The mean of these measurements $$\bar{x}$$ is given by $$\bar{x} = T/n .$$

The variance of the mean is then $$\operatorname{Var}(\bar{x}) = \operatorname{Var}\left(\frac{T}{n}\right) = \frac{1}{n^2}\operatorname{Var}(T) = \frac{1}{n^2}n\sigma^2 = \frac{\sigma^2}{n}.$$

The standard error is, by definition, the standard deviation of $$\bar{x}$$ which is the square root of the variance: $$\sigma_{\bar{x}} = \sqrt{\frac{\sigma^2}{n}} = \frac{\sigma}{\sqrt{n}} .$$

For correlated random variables the sample variance needs to be computed according to the Markov chain central limit theorem.

Independent and identically distributed random variables with random sample size
There are cases when a sample is taken without knowing, in advance, how many observations will be acceptable according to some criterion. In such cases, the sample size $$N$$ is a random variable whose variation adds to the variation of $$X$$ such that,$$\operatorname{Var}(T) = \operatorname{E}(N)\operatorname{Var}(X) + \operatorname{Var}(N)\big(\operatorname{E}(X)\big)^2$$ which follows from the law of total variance.

If $$N$$ has a Poisson distribution, then $$\operatorname{E}(N)= \operatorname{Var}(N)$$ with estimator $$n = N$$. Hence the estimator of $$\operatorname{Var}(T)$$ becomes $$nS^2_X + n\bar{X}^2$$, leading the following formula for standard error: $$\operatorname{Standard~Error}(\bar{X})= \sqrt{\frac{S^2_X + \bar{X}^2}{n}}$$ (since the standard deviation is the square root of the variance).

Student approximation when σ value is unknown
In many practical applications, the true value of σ is unknown. As a result, we need to use a distribution that takes into account that spread of possible σ's. When the true underlying distribution is known to be Gaussian, although with unknown σ, then the resulting estimated distribution follows the Student t-distribution. The standard error is the standard deviation of the Student t-distribution. T-distributions are slightly different from Gaussian, and vary depending on the size of the sample. Small samples are somewhat more likely to underestimate the population standard deviation and have a mean that differs from the true population mean, and the Student t-distribution accounts for the probability of these events with somewhat heavier tails compared to a Gaussian. To estimate the standard error of a Student t-distribution it is sufficient to use the sample standard deviation "s" instead of σ, and we could use this value to calculate confidence intervals.

Note: The Student's probability distribution is approximated well by the Gaussian distribution when the sample size is over 100. For such samples one can use the latter distribution, which is much simpler. Also, even though the 'true' distribution of the population is unknown, assuming normality of the sampling distribution makes sense for a reasonable sample size, and under certain sampling conditions, see CLT. If these conditions are not met, then using a Bootstrap distribution to estimate the Standard Error is often a good workaround, but it can be computationally intensive.

Assumptions and usage
An example of how $$\operatorname{SE}$$ is used is to make confidence intervals of the unknown population mean. If the sampling distribution is normally distributed, the sample mean, the standard error, and the quantiles of the normal distribution can be used to calculate confidence intervals for the true population mean. The following expressions can be used to calculate the upper and lower 95% confidence limits, where $$\bar{x}$$ is equal to the sample mean, $$\operatorname{SE}$$ is equal to the standard error for the sample mean, and 1.96 is the approximate value of the 97.5 percentile point of the normal distribution: Upper 95% limit = $\bar{x} + (\operatorname{SE}\times 1.96) $, and

Lower 95% limit = $\bar{x} - (\operatorname{SE}\times 1.96) $. In particular, the standard error of a sample statistic (such as sample mean) is the actual or estimated standard deviation of the sample mean in the process by which it was generated. In other words, it is the actual or estimated standard deviation of the sampling distribution of the sample statistic. The notation for standard error can be any one of SE, SEM (for standard error of measurement or mean), or SE.

Standard errors provide simple measures of uncertainty in a value and are often used because:
 * in many cases, if the standard error of several individual quantities is known then the standard error of some function of the quantities can be easily calculated;
 * when the probability distribution of the value is known, it can be used to calculate an exact confidence interval;
 * when the probability distribution is unknown, Chebyshev's or the Vysochanskiï–Petunin inequalities can be used to calculate a conservative confidence interval; and
 * as the sample size tends to infinity the central limit theorem guarantees that the sampling distribution of the mean is asymptotically normal.

Standard error of mean versus standard deviation
In scientific and technical literature, experimental data are often summarized either using the mean and standard deviation of the sample data or the mean with the standard error. This often leads to confusion about their interchangeability. However, the mean and standard deviation are descriptive statistics, whereas the standard error of the mean is descriptive of the random sampling process. The standard deviation of the sample data is a description of the variation in measurements, while the standard error of the mean is a probabilistic statement about how the sample size will provide a better bound on estimates of the population mean, in light of the central limit theorem.

Put simply, the standard error of the sample mean is an estimate of how far the sample mean is likely to be from the population mean, whereas the standard deviation of the sample is the degree to which individuals within the sample differ from the sample mean. If the population standard deviation is finite, the standard error of the mean of the sample will tend to zero with increasing sample size, because the estimate of the population mean will improve, while the standard deviation of the sample will tend to approximate the population standard deviation as the sample size increases.

Finite population correction (FPC)
The formula given above for the standard error assumes that the population is infinite. Nonetheless, it is often used for finite populations when people are interested in measuring the process that created the existing finite population (this is called an analytic study). Though the above formula is not exactly correct when the population is finite, the difference between the finite- and infinite-population versions will be small when sampling fraction is small (e.g. a small proportion of a finite population is studied). In this case people often do not correct for the finite population, essentially treating it as an "approximately infinite" population.

If one is interested in measuring an existing finite population that will not change over time, then it is necessary to adjust for the population size (called an enumerative study). When the sampling fraction (often termed f) is large (approximately at 5% or more) in an enumerative study, the estimate of the standard error must be corrected by multiplying by a  finite population correction  (a.k.a.: FPC): $$   \operatorname{FPC} = \sqrt{\frac{N-n}{N-1}} $$ which, for large N: $$ \operatorname{FPC} \approx \sqrt{1-\frac{n}{N}} = \sqrt{1-f} $$ to account for the added precision gained by sampling close to a larger percentage of the population. The effect of the FPC is that the error becomes zero when the sample size n is equal to the population size N.

This happens in survey methodology when sampling without replacement. If sampling with replacement, then FPC does not come into play.

Correction for correlation in the sample


If values of the measured quantity A are not statistically independent but have been obtained from known locations in parameter space x, an unbiased estimate of the true standard error of the mean (actually a correction on the standard deviation part) may be obtained by multiplying the calculated standard error of the sample by the factor f: $$f= \sqrt{\frac{1+\rho}{1-\rho}} ,$$ where the sample bias coefficient ρ is the widely used Prais–Winsten estimate of the autocorrelation-coefficient (a quantity between −1 and +1) for all sample point pairs. This approximate formula is for moderate to large sample sizes; the reference gives the exact formulas for any sample size, and can be applied to heavily autocorrelated time series like Wall Street stock quotes. Moreover, this formula works for positive and negative ρ alike. See also unbiased estimation of standard deviation for more discussion.