Weighted arithmetic mean

Last updated

The weighted arithmetic mean is similar to an ordinary arithmetic mean (the most common type of average), except that instead of each of the data points contributing equally to the final average, some data points contribute more than others. The notion of weighted mean plays a role in descriptive statistics and also occurs in a more general form in several other areas of mathematics.

Contents

If all the weights are equal, then the weighted mean is the same as the arithmetic mean. While weighted means generally behave in a similar fashion to arithmetic means, they do have a few counterintuitive properties, as captured for instance in Simpson's paradox.

Examples

Basic example

Given two school classes, one with 20 students, and one with 30 students, the grades in each class on a test were:

Morning class = 62, 67, 71, 74, 76, 77, 78, 79, 79, 80, 80, 81, 81, 82, 83, 84, 86, 89, 93, 98
Afternoon class = 81, 82, 83, 84, 85, 86, 87, 87, 88, 88, 89, 89, 89, 90, 90, 90, 90, 91, 91, 91, 92, 92, 93, 93, 94, 95, 96, 97, 98, 99

The mean for the morning class is 80 and the mean of the afternoon class is 90. The unweighted mean of the two means is 85. However, this does not account for the difference in number of students in each class (20 versus 30); hence the value of 85 does not reflect the average student grade (independent of class). The average student grade can be obtained by averaging all the grades, without regard to classes (add all the grades up and divide by the total number of students):

Or, this can be accomplished by weighting the class means by the number of students in each class. The larger class is given more "weight":

Thus, the weighted mean makes it possible to find the mean average student grade without knowing each student's score. Only the class means and the number of students in each class are needed.

Convex combination example

Since only the relative weights are relevant, any weighted mean can be expressed using coefficients that sum to one. Such a linear combination is called a convex combination.

Using the previous example, we would get the following weights:

Then, apply the weights like this:

Mathematical definition

Formally, the weighted mean of a non-empty finite multiset of data with corresponding non-negative weights is

which expands to:

Therefore, data elements with a high weight contribute more to the weighted mean than do elements with a low weight. The weights cannot be negative. Some may be zero, but not all of them (since division by zero is not allowed).

The formulas are simplified when the weights are normalized such that they sum up to , i.e.:

.

For such normalized weights the weighted mean is then:

.

Note that one can always normalize the weights by making the following transformation on the original weights:

.

Using the normalized weight yields the same results as when using the original weights:

The ordinary mean is a special case of the weighted mean where all data have equal weights.

The standard error of the weighted mean (unit input variances), can be shown via uncertainty propagation to be:

Statistical properties

The weighted sample mean, , is itself a random variable. Its expected value and standard deviation are related to the expected values and standard deviations of the observations, as follows. For simplicity, we assume normalized weights (weights summing to one).

If the observations have expected values

then the weighted sample mean has expectation

In particular, if the means are equal, , then the expectation of the weighted sample mean will be that value,

For uncorrelated observations with variances , the variance of the weighted sample mean is[ citation needed ]

whose square root can be called the standard error of the weighted mean (general case).[ citation needed ]

Consequently, if all the observations have equal variance, , the weighted sample mean will have variance

where . The variance attains its maximum value, , when all weights except one are zero. Its minimum value is found when all weights are equal (i.e., unweighted mean), in which case we have , i.e., it degenerates into the standard error of the mean, squared.

Note that because one can always transform non-normalized weights to normalized weights all formula in this section can be adapted to non-normalized weights by replacing all .

Variance weights

For the weighted mean of a list of data for which each element potentially comes from a different probability distribution with known variance , one possible choice for the weights is given by the reciprocal of variance:

The weighted mean in this case is:

and the standard error of the weighted mean (with variance weights) is:

Note this reduces to when all . It is a special case of the general formula in previous section,

The equations above can be combined to obtain:

The significance of this choice is that this weighted mean is the maximum likelihood estimator of the mean of the probability distributions under the assumption that they are independent and normally distributed with the same mean.

Correcting for over- or under-dispersion

Weighted means are typically used to find the weighted mean of historical data, rather than theoretically generated data. In this case, there will be some error in the variance of each data point. Typically experimental errors may be underestimated due to the experimenter not taking into account all sources of error in calculating the variance of each data point. In this event, the variance in the weighted mean must be corrected to account for the fact that is too large. The correction that must be made is

where is the reduced chi-squared:

The square root can be called the standard error of the weighted mean (variance weights, scale corrected).

When all data variances are equal, , they cancel out in the weighted mean variance, , which again reduces to the standard error of the mean (squared), , formulated in terms of the sample standard deviation (squared),

Bootstrapping validation

It has been shown by bootstrapping methods that the following is an accurate estimation for the square of the standard error of the mean (general case): [1]

where . Further simplification leads to

Weighted sample variance

Typically when a mean is calculated it is important to know the variance and standard deviation about that mean. When a weighted mean is used, the variance of the weighted sample is different from the variance of the unweighted sample.

The biased weighted sample variance is defined similarly to the normal biased sample variance :

where for normalized weights. If the weights are frequency weights (and thus are random variables), it can be shown that is the maximum likelihood estimator of for iid Gaussian observations.

For small samples, it is customary to use an unbiased estimator for the population variance. In normal unweighted samples, the N in the denominator (corresponding to the sample size) is changed to N  1 (see Bessel's correction). In the weighted setting, there are actually two different unbiased estimators, one for the case of frequency weights and another for the case of reliability weights.

Frequency weights

If the weights are frequency weights (where a weight equals the number of occurrences), then the unbiased estimator is:

This effectively applies Bessel's correction for frequency weights.

For example, if values are drawn from the same distribution, then we can treat this set as an unweighted sample, or we can treat it as the weighted sample with corresponding weights , and we get the same result either way.

If the frequency weights are normalized to 1, then the correct expression after Bessel's correction becomes

where the total number of samples is (not ). In any case, the information on total number of samples is necessary in order to obtain an unbiased correction, even if has a different meaning other than frequency weight.

Note that the estimator can be unbiased only if the weights are not standardized nor normalized, these processes changing the data's mean and variance and thus leading to a loss of the base rate (the population count, which is a requirement for Bessel's correction).

Reliability weights

If the weights are instead non-random (reliability weights[ definition needed ]), we can determine a correction factor to yield an unbiased estimator. Assuming each random variable is sampled from the same distribution with mean and actual variance , taking expectations we have,

where and . Therefore, the bias in our estimator is , analogous to the bias in the unweighted estimator (also notice that is the effective sample size). This means that to unbias our estimator we need to pre-divide by , ensuring that the expected value of the estimated variance equals the actual variance of the sampling distribution.

The final unbiased estimate of sample variance is:

, [2]

where .

The degrees of freedom of the weighted, unbiased sample variance vary accordingly from N  1 down to 0.

The standard deviation is simply the square root of the variance above.

As a side note, other approaches have been described to compute the weighted sample variance. [3]

Weighted sample covariance

In a weighted sample, each row vector (each set of single observations on each of the K random variables) is assigned a weight .

Then the weighted mean vector is given by

And the weighted covariance matrix is given by: [4]

Similarly to weighted sample variance, there are two different unbiased estimators depending on the type of the weights.

Frequency weights

If the weights are frequency weights, the unbiased weighted estimate of the covariance matrix , with Bessel's correction, is given by: [4]

Note that this estimator can be unbiased only if the weights are not standardized nor normalized, these processes changing the data's mean and variance and thus leading to a loss of the base rate (the population count, which is a requirement for Bessel's correction).

Reliability weights

In the case of reliability weights, the weights are normalized:

(If they are not, divide the weights by their sum to normalize prior to calculating :

Then the weighted mean vector can be simplified to

and the unbiased weighted estimate of the covariance matrix is: [5]

The reasoning here is the same as in the previous section.

Since we are assuming the weights are normalized, then and this reduces to:

If all weights are the same, i.e. , then the weighted mean and covariance reduce to the unweighted sample mean and covariance above.

Vector-valued estimates

The above generalizes easily to the case of taking the mean of vector-valued estimates. For example, estimates of position on a plane may have less certainty in one direction than another. As in the scalar case, the weighted mean of multiple estimates can provide a maximum likelihood estimate. We simply replace the variance by the covariance matrix and the arithmetic inverse by the matrix inverse (both denoted in the same way, via superscripts); the weight matrix then reads: [6]

The weighted mean in this case is:

(where the order of the matrix-vector product is not commutative), in terms of the covariance of the weighted mean:

For example, consider the weighted mean of the point [1 0] with high variance in the second component and [0 1] with high variance in the first component. Then

then the weighted mean is:

which makes sense: the [1 0] estimate is "compliant" in the second component and the [0 1] estimate is compliant in the first component, so the weighted mean is nearly [1 1].

Accounting for correlations

In the general case, suppose that , is the covariance matrix relating the quantities , is the common mean to be estimated, and is a design matrix equal to a vector of ones (of length ). The Gauss–Markov theorem states that the estimate of the mean having minimum variance is given by:

and

where:

Decreasing strength of interactions

Consider the time series of an independent variable and a dependent variable , with observations sampled at discrete times . In many common situations, the value of at time depends not only on but also on its past values. Commonly, the strength of this dependence decreases as the separation of observations in time increases. To model this situation, one may replace the independent variable by its sliding mean for a window size .

Exponentially decreasing weights

In the scenario described in the previous section, most frequently the decrease in interaction strength obeys a negative exponential law. If the observations are sampled at equidistant times, then exponential decrease is equivalent to decrease by a constant fraction at each time step. Setting we can define normalized weights by

where is the sum of the unnormalized weights. In this case is simply

approaching for large values of .

The damping constant must correspond to the actual decrease of interaction strength. If this cannot be determined from theoretical considerations, then the following properties of exponentially decreasing weights are useful in making a suitable choice: at step , the weight approximately equals , the tail area the value , the head area . The tail area at step is . Where primarily the closest observations matter and the effect of the remaining observations can be ignored safely, then choose such that the tail area is sufficiently small.

Weighted averages of functions

The concept of weighted average can be extended to functions. [7] Weighted averages of functions play an important role in the systems of weighted differential and integral calculus. [8]

See also

Related Research Articles

Normal distribution Probability distribution

In probability theory, a normaldistribution is a type of continuous probability distribution for a real-valued random variable. The general form of its probability density function is

Standard deviation Measure of the amount of variation or dispersion of a set of values

In statistics, the standard deviation is a measure of the amount of variation or dispersion of a set of values. A low standard deviation indicates that the values tend to be close to the mean of the set, while a high standard deviation indicates that the values are spread out over a wider range.

Variance Statistical measure of how far values spread from their average

In probability theory and statistics, variance is the expectation of the squared deviation of a random variable from its mean. In other words, it measures how far a set of numbers is spread out from their average value. Variance has a central role in statistics, where some ideas that use it include descriptive statistics, statistical inference, hypothesis testing, goodness of fit, and Monte Carlo sampling. Variance is an important tool in the sciences, where statistical analysis of data is common. The variance is the square of the standard deviation, the second central moment of a distribution, and the covariance of the random variable with itself, and it is often represented by , , or .

In probability theory, the central limit theorem (CLT) establishes that, in many situations, when independent random variables are added, their properly normalized sum tends toward a normal distribution even if the original variables themselves are not normally distributed. The theorem is a key concept in probability theory because it implies that probabilistic and statistical methods that work for normal distributions can be applicable to many problems involving other types of distributions.

Multivariate normal distribution Generalization of the one-dimensional normal distribution to higher dimensions

In probability theory and statistics, the multivariate normal distribution, multivariate Gaussian distribution, or joint normal distribution is a generalization of the one-dimensional (univariate) normal distribution to higher dimensions. One definition is that a random vector is said to be k-variate normally distributed if every linear combination of its k components has a univariate normal distribution. Its importance derives mainly from the multivariate central limit theorem. The multivariate normal distribution is often used to describe, at least approximately, any set of (possibly) correlated real-valued random variables each of which clusters around a mean value.

Students <i>t</i>-distribution Probability distribution

In probability and statistics, Student's t-distribution is any member of a family of continuous probability distributions that arise when estimating the mean of a normally-distributed population in situations where the sample size is small and the population's standard deviation is unknown. It was developed by English statistician William Sealy Gosset under the pseudonym "Student".

In probability theory and statistics, covariance is a measure of the joint variability of two random variables. If the greater values of one variable mainly correspond with the greater values of the other variable, and the same holds for the lesser values, the covariance is positive. In the opposite case, when the greater values of one variable mainly correspond to the lesser values of the other,, the covariance is negative. The sign of the covariance therefore shows the tendency in the linear relationship between the variables. The magnitude of the covariance is not easy to interpret because it is not normalized and hence depends on the magnitudes of the variables. The normalized version of the covariance, the correlation coefficient, however, shows by its magnitude the strength of the linear relation.

Covariance matrix measure of covariance of components of a random vector

In probability theory and statistics, a covariance matrix is a square matrix giving the covariance between each pair of elements of a given random vector. Any covariance matrix is symmetric and positive semi-definite and its main diagonal contains variances.

Cross-correlation

In signal processing, cross-correlation is a measure of similarity of two series as a function of the displacement of one relative to the other. This is also known as a sliding dot product or sliding inner-product. It is commonly used for searching a long signal for a shorter, known feature. It has applications in pattern recognition, single particle analysis, electron tomography, averaging, cryptanalysis, and neurophysiology. The cross-correlation is similar in nature to the convolution of two functions. In an autocorrelation, which is the cross-correlation of a signal with itself, there will always be a peak at a lag of zero, and its size will be the signal energy.

In statistics, sometimes the covariance matrix of a multivariate random variable is not known but has to be estimated. Estimation of covariance matrices then deals with the question of how to approximate the actual covariance matrix on the basis of a sample from the multivariate distribution. Simple cases, where observations are complete, can be dealt with by using the sample covariance matrix. The sample covariance matrix (SCM) is an unbiased and efficient estimator of the covariance matrix if the space of covariance matrices is viewed as an extrinsic convex cone in Rp×p; however, measured using the intrinsic geometry of positive-definite matrices, the SCM is a biased and inefficient estimator. In addition, if the random variable has normal distribution, the sample covariance matrix has Wishart distribution and a slightly differently scaled version of it is the maximum likelihood estimate. Cases involving missing data require deeper considerations. Another issue is the robustness to outliers, to which sample covariance matrices are highly sensitive.

In statistics, the bias of an estimator is the difference between this estimator's expected value and the true value of the parameter being estimated. An estimator or decision rule with zero bias is called unbiased. In statistics, "bias" is an objective property of an estimator. Bias can also be measured with respect to the median, rather than the mean, in which case one distinguishes median-unbiased from the usual mean-unbiasedness property. Bias is a distinct concept from consistency. Consistent estimators converge in probability to the true value of the parameter, but may be biased or unbiased; see bias versus consistency for more.

A ratio distribution is a probability distribution constructed as the distribution of the ratio of random variables having two other known distributions. Given two random variables X and Y, the distribution of the random variable Z that is formed as the ratio Z = X/Y is a ratio distribution.

The sample mean and the sample covariance are statistics computed from a sample of data on one or more random variables.

In probability theory and statistics, the normal-gamma distribution is a bivariate four-parameter family of continuous probability distributions. It is the conjugate prior of a normal distribution with unknown mean and precision.

In statistics, pooled variance is a method for estimating variance of several different populations when the mean of each population may be different, but one may assume that the variance of each population is the same. The numerical estimate resulting from the use of this method is also called the pooled variance.

In statistics, Bessel's correction is the use of n − 1 instead of n in the formula for the sample variance and sample standard deviation, where n is the number of observations in a sample. This method corrects the bias in the estimation of the population variance. It also partially corrects the bias in the estimation of the population standard deviation. However, the correction often increases the mean squared error in these estimations. This technique is named after Friedrich Bessel.

In probability and statistics, the class of exponential dispersion models (EDM) is a set of probability distributions that represents a generalisation of the natural exponential family. Exponential dispersion models play an important role in statistical theory, in particular in generalized linear models because they have a special structure which enables deductions to be made about appropriate statistical inference.

In probability theory and statistics, the negative multinomial distribution is a generalization of the negative binomial distribution to more than two outcomes.

In statistics, inverse-variance weighting is a method of aggregating two or more random variables to minimize the variance of the weighted average. Each random variable is weighted in inverse proportion to its variance, i.e. proportional to its precision.

In statistics, effective sample size is a notion defined for a sample from a distribution when the observations in the sample are correlated or weighted.

References

  1. Gatz, Donald F.; Smith, Luther (June 1995). "The standard error of a weighted mean concentration—I. Bootstrapping vs other methods". Atmospheric Environment. 29 (11): 1185–1193. doi:10.1016/1352-2310(94)00210-C.
  2. "GNU Scientific Library – Reference Manual: Weighted Samples". Gnu.org. Retrieved 22 December 2017.
  3. "Weighted Standard Error and its Impact on Significance Testing (WinCross vs. Quantum & SPSS), Dr. Albert Madansky" (PDF). Analyticalgroup.com. Retrieved 22 December 2017.
  4. 1 2 Price, George R. (April 1972). "Extension of covariance selection mathematics" (PDF). Annals of Human Genetics. 35 (4): 485–490. doi:10.1111/j.1469-1809.1957.tb01874.x.
  5. Mark Galassi, Jim Davies, James Theiler, Brian Gough, Gerard Jungman, Michael Booth, and Fabrice Rossi. GNU Scientific Library - Reference manual, Version 1.15, 2011. Sec. 21.7 Weighted Samples
  6. James, Frederick (2006). Statistical Methods in Experimental Physics (2nd ed.). Singapore: World Scientific. p. 324. ISBN   981-270-527-9.
  7. G. H. Hardy, J. E. Littlewood, and G. Pólya. Inequalities (2nd ed.), Cambridge University Press, ISBN   978-0-521-35880-4, 1988.
  8. Jane Grossman, Michael Grossman, Robert Katz. The First Systems of Weighted Differential and Integral Calculus, ISBN   0-9771170-1-4, 1980.

Further reading