0% found this document useful (0 votes)
26 views14 pages

Minimum Variance Unbiased Estimation

This document discusses minimum variance unbiased estimation and the Cramer-Rao lower bound (CRLB). It begins with an example of estimating the DC level of a noisy signal. It then defines minimum variance unbiased estimation and provides the intuition behind the CRLB. The CRLB gives a lower bound on the variance of any unbiased estimator. Examples are provided to illustrate the CRLB for different cases including Gaussian signals. The document concludes with an example of using the CRLB for range estimation from noisy measurements.

Uploaded by

Sushruth
Copyright
© All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
0% found this document useful (0 votes)
26 views14 pages

Minimum Variance Unbiased Estimation

This document discusses minimum variance unbiased estimation and the Cramer-Rao lower bound (CRLB). It begins with an example of estimating the DC level of a noisy signal. It then defines minimum variance unbiased estimation and provides the intuition behind the CRLB. The CRLB gives a lower bound on the variance of any unbiased estimator. Examples are provided to illustrate the CRLB for different cases including Gaussian signals. The document concludes with an example of using the CRLB for range estimation from noisy measurements.

Uploaded by

Sushruth
Copyright
© All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd

Estimation: chapter 2+3

Minimum variance unbiased estimation + the CRLB


Natasha Devroye
devroye@[Link]
[Link]
Spring 2011

Estimation: a first example


Estimate the DC level, A, of a signal given noisy measurements x[0], x[1], ... x
[N-1] where

x[n] are samples of this!

Find a few estimators


Compare their performance

mean?
variance?
pdf?

Estimation: a first example


Estimators of the DC level, A

Estimation: definitions

Estimation: definitions

How would you pick a ``good estimator?


Vector versions....

Minimum variance unbiased estimation

Why?

So?

Minimum variance unbiased estimation

Give a counter-example! (b1pg.20)

The Cramer-Rao Lower Bound


the CRLB give a lower bound on the variance of ANY UNBIASED estimator
does NOT guarantee bound can be obtained
IF find an estimator whose variance = CRLB then its MVUE
otherwise can use Ch.5 tools (Rao-Blackwell-Lehmann-Scheffe Theorem and
Neyman-Fisher Factorization Theorem) to construct a better estimator from
any unbiased one - possibly the MVUE if conditions are met

The Cramer-Rao Lower Bound (CRLB)


Use?
Intuition?

The Cramer-Rao Lower Bound (CRLB)

CRLB examples

CRLB proof

CRLB T or F

What is I()?

Why information?
non-negative
additive for independent observations

Vector form of the CRLB

Vector form of the CRLB

Vector form of the CRLB examples

What can we conclude?

CRLB for transformations

Vector CRLB for transformations

Example of vector CRLB with transformation

CRLB for General Gaussian Case


When observations are Gaussian and one knows the dependence of the
mean and covariance matrix on the unknown parameters, we know the
closed form of the CRLB (or Fisher information matrix):

CRLB for Gaussians examples

Example: range estimation

678$9$:)('.$6;<&=)<&,($/>,?*.=
678$9$:)('.$6;<&=)<&,($/>,?*.=

-.&/0%12345607833
"+#)333333/9/:7.939;7.3#&<,($"=
-.&/0%12345607833
"+#)333333/9/:7.939;7.3#&<,($"=
>7?71;73>7@67?219/8333 "+#'( !!)
>7?71;73>7@67?219/8333 "+#'( !!)
A7&05.73-1%73B76&C8 !!
A7&05.73-1%73B76&C8 !!
!"#$%&'()*$+,-.*
Continuous time signal model!"#$%&'()*$+,-.*
& + # ) # " + # % ! ) " %+ # )
, $ # $ $ # $ "!
!
$
!#!
"
& +# ) #" +"#+*!# %) ! ! ) " %+# )
$
!!#!
"

"

" + # *! ! )

"+#)
"+#)

D&/E61%127E
FG1273H&5001&/
D&/E61%127E
FG1273H&5001&/
/%0$,1$!2"3

D4I &+#)
J3K%L
D4I &+#)
J3K%L

$"

#$45
/%0$,1$!2"3
#$45
(*

(*

"+#''( $! " )

"+#''( !!)

)'

!(%&'

, $ # $ $ # $" " ! !(%&'

)'

!"#$

!"#$

Taken from [Link]

Example: range estimation


/0123&45678067'1&-9:&,7210;&<'=3;
Discrete time
signal model
/0123&45678067'1&-9:&,7210;&<'=3;
#$%&'(&$)%*
#$%&'(&$)%*
"! " !"#
"! " !"#

+,-&'(&$)%*
+,-&'(&$)%*
"#.!
"#.!

*!
*!

%&'()*+,-*./+! 0+#"1!"2*3
%&'()*+,-*./+!
0+#"1!"2*3
#4$5+0+#6$
!7
#4$5+0+#6$!7

&

)"

&
#"1! #"! !"1!
#"1! #"! !"1!

)"

'4$ 5 $ & 6 $! # & ( 7 % #4$ 5


'4$ 5 $ & 6 $! # & ( 7 % #4$ 5

89+:;<=*+
89+:;<=*+
>&?22<&@+AB<2*
>&?22<&@+AB<2*
C&.+"! " !"#
C&.+"! " !"#

$ $ ED#D ! D % # #
$ $ ED#D ! D % # #

&4$F&(5G ;&2++ @B@HI*.B+2&'()*2+2=&.=<@J+&=+$(


&4$F&(5G ;&2++ @B@HI*.B+2&'()*2+2=&.=<@J+&=+$(

E + $ + $( # #
E + $ + $( # #
$( + $ + $ ( % + # #
$( + $ + $ ( % + # #
$( % + + $ + % # #
$( % + + $ + % # #

* #4$ 5
*'#4$ 5
'
'4$ 5 $'') & 6 $! # & ( 7 % #4$ 5
'4$ 5 $ )'& 6 $! # & ( 7 % #4$ 5
'(' #4$ 5
'( #4$ 5

!"#$

!"#$

Taken from [Link]

Example: range estimation


!"#$%&'()*+")*,#&-!./
%&'()**+,(-.)/0)10(2345(16-7+.(8&1(-9:/)+(;(<=%>
H+7:(9/I()/0(J66*(
/&/KL61&(.61M-

G)1B/F$ A -

.@
% "#

+ ,&E"D/ C (
! )) ,/ $ $ &&
'
" # ?*

.@

"$ % ! "#

" # "$

.@
"$ % ! "# +

" # "$

) ,& B# A
) ,#
*

NO*+&9.(2)+P7+7-QQQ

(
&
&
# # "$ "/ $ '

+ ,& B "$ " / $ A (


))
&&
,/ $
*
'

.@
! "# +

,& B# A
! )) ,#
" #? *

(
&
&
# # "$ '

R-6()**1&O9M).9&/>(/$'('$ "$
ST6/(0&(PT)/:6(&8(G)19)U+6-QQ
!"#$
Taken from [Link]

Example: range estimation

/'%0"&1+2$*'2$3%&4/5(&673%289

-../012.+03412.3+567826.2.0+449 +33,:;<2./02=>267?18,+49
*+,(-)# ' +

,&
&

# "! ' (! (& ' $


%
" %&
* !% & (& #

'# " &


"! ' (! ( & ' $

!%

&

%
" %&
& (& #

"! ' (! ( & ' $

$!
'# " &

!%

&

%
" %&
& (& #
$!
"!

*+,(-)# ' +
$!
'# " &

"
&
&
!% .&2( / ) ( ( ' %(
!

$!
0

$!
'# " &

$ ! ) ! ! & (& '%&


%

!"#$%"&'&()&*"'+,-".
0

@A2AB1:,10
C2D+,.1*+4

!10 .&2( / ) ( ( ' %(


0
&
!10 ) ( ( ' %&
&

&

D+,.1*+4

*+,! )

!10 .&2( / ) ( ( ' %(


0
&
!10 ) ( ( ' %&
&

&

*+,! 6.2EFGH2IJK22(LM'
!"#$

-2?>312:N2EHOFK

Taken from [Link]

Example: range estimation


&"#$%'()*+,"*+-#'.&/0'12-#*34
&"#$%'()*+,"*+-#'.&/0'12-#*34
%&'()*+,-&-*'.-/&*0-*/11'2-*/+*+,-*3456*7(*+,-*.-8/9:

%&'()*+,-&-*'.-/&*0-*/11'2-*/+*+,-*3456*7(*+,-*.-8/9:
#
2/1<#>$ ; "
<&-= ? ;
#?
&'(
!
%
!"#
2/1<#>$ ; "
<&-= ? ;
?
&'( ! %!"#

@7*)-+*+,-*3456*7(*+,-*!"#$%A*B&-*C+1/(&DE*7D*F/1G&H*1-&B8+:*
@7*)-+*+,-*3456*7(*+,-*!"#$%A*B&-*C+1/(&DE*7D*F/1G&H*1-&B8+:*
?

) *( &
$ *(+%
*(+%(> + ''
0'+,* (
?
#>$
*#) $*$( &
*(+%(> (+ '' % $$ *(+%#>$ 0'+,*
( *# $ %

J*)#$,"*?
( J*)#$,"*?

3456*'&*'(2-1&-89*F17F71+'7(/8*+7:
3456*'&*'(2-1&-89*F17F71+'7(/8*+7:
3456*'&*'(2-1&-89*F17F71+'7(/8*+7:
KK LM4*N-/&B13456*'&*'(2-1&-89*F17F71+'7(/8*+7:
LM4*N-/&B1LM4*N-/&B1KK KK4NL*6O*N-/&B1LM4*N-/&B14NL*6O*N-/&B1KK 4NL*6O*N-/&B14NL*6O*N-/&B1-

)? " I
2/1< (> ; "
<G ? ;
) ?%"?I
&'(
!
>
!"#
2/1< ( ; "
<G ? ;
?
&'( ! %!"#

L7*+,-*3456*+-88&*B&A
K 3,77&-*&')(/8*0'+,*8/1)-*%
L7*+,-*3456*+-88&*B&A !"#
K P(&B1-*+,/+*LM4*'&*8/1)K 3,77&-*&')(/8*0'+,*8/1)-*%!"#
K 6-++-1*7(*M-/1Q9"8/1)-*+/1)-+&
K P(&B1-*+,/+*LM4*'&*8/1)K O,'=,*'&*Q-++-1R**
K 6-++-1*7(*M-/1Q9"8/1)-*+/1)-+&
S T7BQ8-*+1/(&G'++-.*-(-1)9R
K O,'=,*'&*Q-++-1R**
S T7BQ8-*4NL*Q/(.0'.+,R
S T7BQ8-*+1/(&G'++-.*-(-1)9R
S T7BQ8-*4NL*Q/(.0'.+,R

!"#$

!"#$
Taken from [Link]

Common questions

Powered by AI

The Fisher Information Matrix is critical because it quantifies the amount of information that an observable random variable carries about an unknown parameter upon which the likelihood depends. In the context of the Cramer-Rao Lower Bound, the Fisher Information is directly used to formulate the bound itself, as the inverse of the Fisher Information Matrix provides the lower bound on the variance of unbiased estimators. It is non-negative and additive for independent observations, making it a foundational tool for deriving and analyzing the CRLB .

The attainability of the CRLB in practical scenarios is determined by several factors including the choice of estimator, the nature and quality of data, and the validity of underlying assumptions of the model. For instance, estimators need to be unbiased and sample data must exhibit properties like independence and identical distribution. Additionally, models should comply with regularity conditions, like the differentiability of probability density functions, which are necessary for the CRLB to hold and be achievable. Violations of these conditions or assumptions can make attaining the CRLB impractical .

Unbiasedness is a central assumption in the context of the Cramer-Rao Lower Bound because the bound applies only to unbiased estimators. An estimator is unbiased if its expected value equals the true value of the parameter across all possible samples. This ensures that the estimator does not systematically over or underestimate the parameter. The CRLB essentially sets the benchmark of variance for these unbiased estimators, serving as a measure against which the precision of any unbiased estimator is judged .

The primary purpose of the Cramer-Rao Lower Bound (CRLB) is to provide a lower bound on the variance of any unbiased estimator. It essentially defines the theoretical limit of precision an estimator can achieve. However, while CRLB specifies the minimum variance, it does not guarantee that this bound can be attained by any particular estimator. The CRLB is significant because if an estimator can achieve variance equal to the CRLB, it is considered the minimum variance unbiased estimator (MVUE).

When applying CRLB to estimators of transformation parameters, challenges include properly accounting for the influence of the transformation on both the estimator and the information contained in the data. The variance bound needs adjustment for transformations, which can be complex depending on the relationship between the transformed and original parameters. Moreover, ensuring the conditions of independence and differentiability required for CRLB are satisfied increases complexity, especially for non-linear transformations .

Vector forms of the CRLB offer advantages when modeling with transformations by providing a comprehensive framework that can capture multi-dimensional parameter relationships and their mutual influences on variance bounds. They facilitate the consideration of dependencies and covariances among parameters, essential for accurately understanding the transformed parameter space. Such a framework aids in deriving correct variance bounds and assists in enhancing the estimation procedures in multidimensional and transformation-rich contexts .

In vector estimation scenarios, the CRLB is derived using the Fisher Information Matrix instead of a single Fisher Information value used in scalar cases. The CRLB in these cases is specified as a matrix that provides variance bounds on each component of the vector parameter being estimated. Unlike scalar cases where variance is a single value, applying the CRLB involves analyzing the entire matrix to understand the precision and limits of multiple correlated estimates simultaneously. This complexity requires additional mathematical tools and understanding of matrix operations .

Achieving equality in the Cramer-Rao Lower Bound for a particular estimator signifies that the estimator has the minimum possible variance among all unbiased estimators, thus it is the Minimum Variance Unbiased Estimator (MVUE). This demonstrates optimal estimation precision given the available data, and affirms that further unbiased improvements in variance are not possible with respect to the unknown parameter being estimated. It validates the efficiency and effectiveness of the estimator within the given statistical model .

The Rao-Blackwell-Lehmann-Scheffe theorem helps in constructing improved estimators by converting any unbiased estimator into an unbiased estimator with a variance that is less than or equal to the original estimator. This is achieved by conditioning the original estimator on a sufficient statistic of the data. The theorem is used to derive the minimum variance unbiased estimator (MVUE) when conditions are met, thus systematically enhancing estimation accuracy and efficiency .

The Gaussian assumption impacts the application of the CRLB by simplifying the derivation of the bound. When observations are Gaussian and the dependence of mean and covariance matrix on unknown parameters is known, the CRLB or the Fisher Information Matrix can be provided in a closed form. This assumption allows for straightforward application and understanding of the CRLB since the statistical properties of Gaussian distributions are well-characterized, making it easier to compute and apply in practice .

You might also like