Misplaced Pages

Minimum-variance unbiased estimator

Article snapshot taken from Wikipedia with creative commons attribution-sharealike license. Give it a read and then ask your questions in the chat. We can research this topic together.
(Redirected from Minimum variance unbiased estimator) Unbiased statistical estimator minimizing variance
This article needs additional citations for verification. Please help improve this article by adding citations to reliable sources. Unsourced material may be challenged and removed.
Find sources: "Minimum-variance unbiased estimator" – news · newspapers · books · scholar · JSTOR (November 2009) (Learn how and when to remove this message)

In statistics a minimum-variance unbiased estimator (MVUE) or uniformly minimum-variance unbiased estimator (UMVUE) is an unbiased estimator that has lower variance than any other unbiased estimator for all possible values of the parameter.

For practical statistics problems, it is important to determine the MVUE if one exists, since less-than-optimal procedures would naturally be avoided, other things being equal. This has led to substantial development of statistical theory related to the problem of optimal estimation.

While combining the constraint of unbiasedness with the desirability metric of least variance leads to good results in most practical settings—making MVUE a natural starting point for a broad range of analyses—a targeted specification may perform better for a given problem; thus, MVUE is not always the best stopping point.

Definition

Consider estimation of g ( θ ) {\displaystyle g(\theta )} based on data X 1 , X 2 , , X n {\displaystyle X_{1},X_{2},\ldots ,X_{n}} i.i.d. from some member of a family of densities p θ , θ Ω {\displaystyle p_{\theta },\theta \in \Omega } , where Ω {\displaystyle \Omega } is the parameter space. An unbiased estimator δ ( X 1 , X 2 , , X n ) {\displaystyle \delta (X_{1},X_{2},\ldots ,X_{n})} of g ( θ ) {\displaystyle g(\theta )} is UMVUE if θ Ω {\displaystyle \forall \theta \in \Omega } ,

var ( δ ( X 1 , X 2 , , X n ) ) var ( δ ~ ( X 1 , X 2 , , X n ) ) {\displaystyle \operatorname {var} (\delta (X_{1},X_{2},\ldots ,X_{n}))\leq \operatorname {var} ({\tilde {\delta }}(X_{1},X_{2},\ldots ,X_{n}))}

for any other unbiased estimator δ ~ . {\displaystyle {\tilde {\delta }}.}

If an unbiased estimator of g ( θ ) {\displaystyle g(\theta )} exists, then one can prove there is an essentially unique MVUE. Using the Rao–Blackwell theorem one can also prove that determining the MVUE is simply a matter of finding a complete sufficient statistic for the family p θ , θ Ω {\displaystyle p_{\theta },\theta \in \Omega } and conditioning any unbiased estimator on it.

Further, by the Lehmann–Scheffé theorem, an unbiased estimator that is a function of a complete, sufficient statistic is the UMVUE estimator.

Put formally, suppose δ ( X 1 , X 2 , , X n ) {\displaystyle \delta (X_{1},X_{2},\ldots ,X_{n})} is unbiased for g ( θ ) {\displaystyle g(\theta )} , and that T {\displaystyle T} is a complete sufficient statistic for the family of densities. Then

η ( X 1 , X 2 , , X n ) = E ( δ ( X 1 , X 2 , , X n ) T ) {\displaystyle \eta (X_{1},X_{2},\ldots ,X_{n})=\operatorname {E} (\delta (X_{1},X_{2},\ldots ,X_{n})\mid T)\,}

is the MVUE for g ( θ ) . {\displaystyle g(\theta ).}

A Bayesian analog is a Bayes estimator, particularly with minimum mean square error (MMSE).

Estimator selection

An efficient estimator need not exist, but if it does and if it is unbiased, it is the MVUE. Since the mean squared error (MSE) of an estimator δ is

MSE ( δ ) = var ( δ ) + [ bias ( δ ) ] 2   {\displaystyle \operatorname {MSE} (\delta )=\operatorname {var} (\delta )+^{2}\ }

the MVUE minimizes MSE among unbiased estimators. In some cases biased estimators have lower MSE because they have a smaller variance than does any unbiased estimator; see estimator bias.

Example

Consider the data to be a single observation from an absolutely continuous distribution on R {\displaystyle \mathbb {R} } with density

p θ ( x ) = θ e x ( 1 + e x ) θ + 1 {\displaystyle p_{\theta }(x)={\frac {\theta e^{-x}}{(1+e^{-x})^{\theta +1}}}}

and we wish to find the UMVU estimator of

g ( θ ) = 1 θ 2 {\displaystyle g(\theta )={\frac {1}{\theta ^{2}}}}

First we recognize that the density can be written as

e x 1 + e x exp ( θ log ( 1 + e x ) + log ( θ ) ) {\displaystyle {\frac {e^{-x}}{1+e^{-x}}}\exp(-\theta \log(1+e^{-x})+\log(\theta ))}

Which is an exponential family with sufficient statistic T = log ( 1 + e x ) {\displaystyle T=\log(1+e^{-x})} . In fact this is a full rank exponential family, and therefore T {\displaystyle T} is complete sufficient. See exponential family for a derivation which shows

E ( T ) = 1 θ , var ( T ) = 1 θ 2 {\displaystyle \operatorname {E} (T)={\frac {1}{\theta }},\quad \operatorname {var} (T)={\frac {1}{\theta ^{2}}}}

Therefore,

E ( T 2 ) = 2 θ 2 {\displaystyle \operatorname {E} (T^{2})={\frac {2}{\theta ^{2}}}}

Here we use Lehmann–Scheffé theorem to get the MVUE

Clearly δ ( X ) = T 2 2 {\displaystyle \delta (X)={\frac {T^{2}}{2}}} is unbiased and T = log ( 1 + e x ) {\displaystyle T=\log(1+e^{-x})} is complete sufficient, thus the UMVU estimator is

η ( X ) = E ( δ ( X ) T ) = E ( T 2 2 | T ) = T 2 2 = log ( 1 + e X ) 2 2 {\displaystyle \eta (X)=\operatorname {E} (\delta (X)\mid T)=\operatorname {E} \left(\left.{\frac {T^{2}}{2}}\,\right|\,T\right)={\frac {T^{2}}{2}}={\frac {\log(1+e^{-X})^{2}}{2}}}

This example illustrates that an unbiased function of the complete sufficient statistic will be UMVU, as Lehmann–Scheffé theorem states.

Other examples

k + 1 k m 1 , {\displaystyle {\frac {k+1}{k}}m-1,}
where m is the sample maximum. This is a scaled and shifted (so unbiased) transform of the sample maximum, which is a sufficient and complete statistic. See German tank problem for details.

See also

Bayesian analogs

References

  1. Lee, A. J., 1946- (1990). U-statistics : theory and practice. New York: M. Dekker. ISBN 0824782534. OCLC 21523971.{{cite book}}: CS1 maint: multiple names: authors list (link) CS1 maint: numeric names: authors list (link)
  • Keener, Robert W. (2006). Statistical Theory: Notes for a Course in Theoretical Statistics. Springer. pp. 47–48, 57–58.
  • Keener, Robert W. (2010). Theoretical statistics: Topics for a core course. New York: Springer. DOI 10.1007/978-0-387-93839-4
  • Voinov V. G., Nikulin M.S. (1993). Unbiased estimators and their applications, Vol.1: Univariate case. Kluwer Academic Publishers. pp. 521p.
Statistics
Descriptive statistics
Continuous data
Center
Dispersion
Shape
Count data
Summary tables
Dependence
Graphics
Data collection
Study design
Survey methodology
Controlled experiments
Adaptive designs
Observational studies
Statistical inference
Statistical theory
Frequentist inference
Point estimation
Interval estimation
Testing hypotheses
Parametric tests
Specific tests
Goodness of fit
Rank statistics
Bayesian inference
Correlation
Regression analysis
Linear regression
Non-standard predictors
Generalized linear model
Partition of variance
Categorical / Multivariate / Time-series / Survival analysis
Categorical
Multivariate
Time-series
General
Specific tests
Time domain
Frequency domain
Survival
Survival function
Hazard function
Test
Applications
Biostatistics
Engineering statistics
Social statistics
Spatial statistics
Category: