next up previous contents
Nächste Seite: Suffizienz Aufwärts: Güteeigenschaften von Punktschätzern Vorherige Seite: Erwartungstreue; mittlerer quadratischer Fehler   Inhalt


Ungleichung von Cramér-Rao

Definition
 


Beispiel
$ \;$ Poisson-verteilte Stichprobenvariablen


Zunächst leiten wir eine allgemeine untere Schranke für die Varianz von Schätzern mit gewissen Regularitätseigenschaften her, die in der Literatur Ungleichung von Cramér-Rao genannt wird.


Das parametrische Modell % latex2html id marker 27401
$ \{P_\theta,\theta\in\Theta\}$ genüge den folgenden Regularitätsbedingungen:

  1. Die Familie % latex2html id marker 27403
$ \{P_\theta,\theta\in\Theta\}$ bestehe entweder nur aus diskreten Verteilungen oder nur aus absolutstetigen Verteilungen, wobei % latex2html id marker 27405
$ \Theta\subset\mathbb{R}$ ein offenes Intervall sei.
  2. Die Menge $ B=\{x\in\mathbb{R}:L(x;\theta)>0\}$ hänge nicht von % latex2html id marker 27409
$ \theta\in\Theta$ ab, wobei die Likelihood-Funktion $ L(x;\theta)$ gegeben ist durch

    $\displaystyle L(x;\theta)=\left\{\begin{array}{ll} p(x;\theta) & \mbox{im diskr...
... Fall,}\\
f(x;\theta) & \mbox{im absolutstetigen Fall}
\end{array}\right.
$

    und $ p(x;\theta)$ bzw. $ f(x;\theta)$ die Wahrscheinlichkeitsfunktion bzw. Dichte von $ P_\theta$ ist.
  3. Die Ableitung $ \frac{\partial}{\partial\theta}L(x;\theta)$ existiere für beliebige % latex2html id marker 27423
$ \theta\in\Theta$ und $ x\in B$.
  4. Vertauschbarkeit von Ableitung und Summe/Integral: Für jedes % latex2html id marker 27427
$ \theta\in\Theta$ gelte

    $\displaystyle \sum\limits_{x\in B}
 \frac{\partial}{\partial\theta}L(x;\theta)=
 \frac{d}{d\theta}\sum\limits_{x\in B} L(x;\theta)=0$   im diskreten Fall (28)

    bzw.

    $\displaystyle \int\limits_B \frac{\partial}{\partial\theta}L(x;\theta)\,dx
 =\frac{d}{d\theta}\int\limits_B L(x;\theta)\,dx =0$   im absolutstetigen Fall. (29)


Theorem 2.2   Sei % latex2html id marker 27434
$ \{P_\theta,\,\theta\in\Theta\}$ eine Familie von Verteilungen, die den Regularitätsbedingungen $ 1$-$ 4$ genügt, und sei % latex2html id marker 27440
$ \widehat\theta:\mathbb{R}^n\to\Theta$ eine Stichprobenfunktion, so dass für jedes % latex2html id marker 27442
$ \theta\in\Theta$ Dann gilt für jedes % latex2html id marker 27454
$ \theta\in\Theta$

$\displaystyle {\rm Var\,}_\theta\,\widehat\theta(X_1,\ldots,X_n)\ge
 \frac{\dis...
...(\Bigl(
 \frac{\partial}{\partial\theta}\,\log
 L(X_1;\theta)\Bigr)^2\Bigr)}\;.$ (31)

Beweis
 

Korollar 2.1   Sei % latex2html id marker 27530
$ \{P_\theta,\,\theta\in\Theta\}$ eine Familie von Verteilungen und % latex2html id marker 27532
$ \widehat\theta:\mathbb{R}^n\to\Theta$ eine Stichprobenfunktion, die den Bedingungen von Theorem  % latex2html id marker 27534
$ \ref{the.ung.rao}$ genügen. Falls $ \widehat\theta(X_1,\ldots,X_n)$ ein erwartungstreuer Schätzer für $ \theta$ ist, dann gilt für jedes % latex2html id marker 27540
$ \theta\in\Theta$

$\displaystyle {\rm Var\,}_\theta\,\widehat\theta(X_1,\ldots,X_n)\ge
 \frac{1}{\...
...(\Bigl(
 \frac{\partial}{\partial\theta}\,\log
 L(X_1;\theta)\Bigr)^2\Bigr)}\;.$ (32)

Beweis
 
Wir diskutieren nun zwei Beispiele von parametrischen Verteilungsfamilien, für die die Bedingungen von Theorem 2.2 bzw. Korollar 2.1 erfüllt sind. Außerdem zeigen wir, dass für diese beiden Verteilungsfamilien das Stichprobenmittel $ \overline X_n$ ein Schätzer ist, der ,,optimal'' im Sinne der Ungleichungen (31) bzw. (32) ist.
Beispiele
 
  1. Normalverteilte Stichprobenvariablen
    • Falls $ X_i\sim$ N $ (\mu,\sigma^2)$, wobei die Varianz $ \sigma^2$ bekannt sei, dann gilt % latex2html id marker 27564
$ \Theta=B=\mathbb{R}$ und für jedes $ x\in\mathbb{R}$

      $\displaystyle L(x;\mu)=\frac{1}{\sqrt{2\pi}\sigma}\;\exp\Bigl(-\frac{(x-\mu)^2}{2\sigma^2}\Bigr)$ (33)

      bzw.

      $\displaystyle \frac{d }{d \mu}\,L(x;\mu)=
 \frac{1}{\sqrt{2\pi}\sigma}\;\frac{x-\mu}{\sigma^2}\;
 \exp\Bigl(-\frac{(x-\mu)^2}{2\sigma^2}\Bigr)\;.$ (34)

    • Somit gilt

      $\displaystyle \int\limits_{\mathbb{R}} \frac{d }{d \mu}L(x;\mu)\,dx
 =\frac{d}{d\mu}\int\limits_{\mathbb{R}} L(x;\mu)\,dx =0 \,.$ (35)

    • Die Regularitätsbedingungen 1-4, die unmittelbar vor Theorem 2.2 formuliert wurden, sind also erfüllt. Wir zeigen nun, dass der Schätzer $ \overline X_n$ für $ \mu$ die Bedingungen von Theorem 2.2 erfüllt.
    • Weil $ \overline X_n$ normalverteilt und ein erwartungstreuer Schätzer für $ \mu$ ist, gilt offenbar $ {\mathbb{E}\,}_\mu(\overline
X_n^2)<\infty$, und die Ableitung $ (d/d\mu)\, {\mathbb{E}\,}_\mu\overline
X_n=1$ ist wohldefiniert.
    • Die Gültigkeit von (30) mittels vollständiger Induktion,
      • denn aus (33) und (34) folgt, dass für $ n=1$

        $\displaystyle \int\limits_{\mathbb{R}}x\;\frac{d }{d \mu}\,L(x;\mu)\,dx =
 \int...
...\mu)\,dx=\frac{1}{\sigma^2}\;
 \bigl({\mathbb{E}\,}_\mu(X_i^2)-\mu^2\bigr)=1\,.$ (36)

      • Für ein $ n\in\mathbb{N}$ gelte nun

        $\displaystyle \int\limits_{\mathbb{R}^n}(x_1+\ldots+x_n)\;\frac{d }{d \mu}\,
 L(x_1,\ldots,x_n;\mu)\,d(x_1,\ldots,x_n) =n\,.$ (37)

      • Weil
        $\displaystyle { \frac{d }{d \mu}\,
L(x_1,\ldots,x_{n+1};\mu)\,d(x_1,\ldots,x_{n+1})}$
          $\displaystyle =$ $\displaystyle \frac{d }{d \mu}\,
L(x_1,\ldots,x_n;\mu)\,L(x_{n+1};\mu)+L(x_1,\ldots,x_n;\mu)\;\frac{d
}{d \mu} \,L(x_{n+1};\mu)$  

        und weil

        $\displaystyle \int\limits_{\mathbb{R}}\;\frac{d }{d \mu}\,L(x_{i};\mu)\,dx_{i}
=0\,,\qquad\forall\,i=1,\ldots,n+1\,,
$

        ergibt sich nun aus (36) und (37), dass

        $\displaystyle \int\limits_{\mathbb{R}^{n+1}}(x_1+\ldots+x_{n+1})\;\frac{d }{d \mu}\,
L(x_1,\ldots,x_{n+1};\mu)\,d(x_1,\ldots,x_{n+1}) =n+1\,.
$

    • Außerdem gilt

      $\displaystyle \frac{d }{d \mu}\,\log L(x,\mu)=\;\frac{d }{d
\mu}\;\Bigl(-\;\fr...
...(2\pi\sigma^2)-\;
\frac{(x-\mu)^2}{2\sigma^2}\Bigr)=\;\frac{x-\mu}{\sigma^2}
$

      und somit

      $\displaystyle {\rm Var\,}_\mu\,\overline X_n =\;\frac{\sigma^2}{n}\;=\;
\frac{...
...bb{E}\,}_\mu\Bigl(\Bigl( \frac{d }{d \mu}\,\log
L(X_1;\mu)\Bigr)^2\Bigr)}\;.
$

  2. Poisson-verteilte Stichprobenvariablen
    • Falls $ X_i\sim$ Poi$ (\lambda)$, dann gilt $ B=\mathbb{N}$ und $ L(x;\lambda)=(\lambda^x/x!)\,e^{-\lambda}$ für jedes $ x\in\mathbb{N}$ bzw.

      $\displaystyle \frac{d }{d \lambda}\,L(x;\lambda)=\left\{\begin{array}{ll}
 -e^{...
...ambda^{x-1}}{(x-1)!}\;e^{-\lambda}\,,&
 \mbox{falls $x>0$.}
 \end{array}\right.$ (38)

    • Somit gilt

      $\displaystyle \sum\limits_{x\in\mathbb{N}}\frac{d }{d \lambda}\,L(x;\lambda)=
 ...
...gr)\frac{\lambda^{k-1}}{(k-1)!}\;e^{-\lambda}
 =-e^{-\lambda}+e^{-\lambda}=0\,.$ (39)

    • Die Regularitätsbedingungen 1-4 von Theorem 2.2 sind also erfüllt.
    • Außerdem genügt der Schätzer $ \overline X_n$ für $ \lambda$ den Bedingungen von Theorem 2.2, denn für jedes $ \lambda>0$ gilt

      $\displaystyle {\mathbb{E}\,}_\lambda\vert\overline X_n\vert={\mathbb{E}\,}_\lam...
...line X_n
=\lambda\,,\qquad{\rm Var\,}_\lambda\overline X_n=\frac{\lambda}{n}
$

      und somit $ (d/d\lambda)\,{\mathbb{E}\,}_\lambda\overline X_n=1$.
    • Hieraus und aus (38) ergibt sich die Gültigkeit der Bedingung (30) erneut mittels vollständiger Induktion:
      • Für $ n=1$ gilt

        $\displaystyle \sum\limits_{x\in\mathbb{N}} x\,\frac{d }{d \lambda}\,L(x;\lambda...
...da}{k}\Bigr)\frac{\lambda^{k-1}}{(k-1)!}\;e^{-\lambda}
 =\lambda+1-\lambda=1\,.$ (40)

      • Wir nehmen nun an, dass

        $\displaystyle \sum\limits_{x_1,\ldots,x_{n-1}\in\mathbb{N}}
 (x_1+\ldots+x_{n-1})\,\frac{d }{d
 \lambda}\,L(x_1,\ldots,x_{n-1};\lambda)=n-1\,.$ (41)

      • Dann gilt
        $\displaystyle { \sum\limits_{x_1,\ldots,x_n\in\mathbb{N}}
(x_1+\ldots+x_n)\,\frac{d }{d \lambda}\,L(x_1,\ldots,x_n;\lambda)}$
          $\displaystyle =$ $\displaystyle \sum\limits_{x_1,\ldots,x_{n-1}\in\mathbb{N}}\sum\limits_{x_n\in\...
...(\Bigl(\frac{d }{d
\lambda}\,L(x_1,\ldots,x_{n-1};\lambda)\Bigr)
L(x_n;\lambda)$  
            $\displaystyle +L(x_1,\ldots,x_{n-1};\lambda)
\frac{d }{d \lambda}\,L(x_n;\lambda)\Bigr)$  
          % latex2html id marker 27646
$\displaystyle \stackrel{(\ref{poi.bed.rao}),(\ref{ind.poi.rao})}{=}$ $\displaystyle \sum\limits_{x_1,\ldots,x_{n-1}\in\mathbb{N}}\Bigl((x_1+\ldots+x_{n-1})
\frac{d }{d \lambda}\,L(x_1,\ldots,x_{n-1};\lambda)$  
            $\displaystyle + \lambda\frac{d }{d \lambda}\,L(x_1,\ldots,x_{n-1};\lambda)+
L(x_1,\ldots,x_{n-1};\lambda)\Bigr)$  
          % latex2html id marker 27652
$\displaystyle \stackrel{(\ref{ind.ann.rao}),(\ref{poi.bed.rao})}{=}$ $\displaystyle (n-1)+0+1=n\,.$  

    • Außerdem gilt für jedes $ x\in\mathbb{N}$

      $\displaystyle \frac{d }{d \lambda}\,\log
L(x;\lambda)=\frac{x}{\lambda}-1=\frac{1}{\lambda}(x-\lambda)\,.
$

    • Hieraus folgt, dass

      $\displaystyle \frac{1}{\displaystyle n\,{\mathbb{E}\,}_\lambda\Bigl(\Bigl( \fra...
...
n\,{\mathbb{E}\,}_\lambda\bigl((X_1-\lambda)^2\bigr)} =\frac{\lambda}{n}\,.
$

    • Weil $ {\rm Var\,}_\lambda\overline X_n=\lambda/n$ gilt, ist damit gezeigt, dass in der Klasse derjenigen Schätzer, die die Bedingungen von Theorem 2.2 erfüllen, das Stichprobenmittel $ \overline X_n$ bester erwartungstreuer Schätzer für $ \lambda$ ist.


Beachte
$ \;$ Es gibt jedoch Familien % latex2html id marker 27670
$ \{P_\theta,\theta\in\Theta\}$ von Verteilungen,
Beispiel
$ \;$ Gleichverteilte Stichprobenvariablen

next up previous contents
Nächste Seite: Suffizienz Aufwärts: Güteeigenschaften von Punktschätzern Vorherige Seite: Erwartungstreue; mittlerer quadratischer Fehler   Inhalt
Ursa Pantle 2004-07-14