Kurs:Statistik für Anwender/Hypergeometrischverteilte Zufallsvariable

Hypergeometrisch verteilte ZV

Situation

In einer Menge von ${\textstyle N\in \mathbb {N} }$ Objekten sind ${\textstyle K\in \{1,...,N\}}$ Objekte mit einer bestimmten Eigenschaft ausgezeichnet. Nun werden daraus ${\textstyle n\in \{1,...,n\}}$ Objekte zufällig ausgewählt (gezogen). Wichtig ist dabei, dass die Ziehung zufällig und unabhängig von der Eigenschaft ist, d.h. die ausgezeichneten Objekte haben dieselbe Chance gezogen zu werden, wie die anderen Objekte.

Hypergeometrisch verteilte ZV

Die ZV ${\textstyle A}$ beschreibt die Zahl ${\textstyle k\in \{1,...,n\}}$ der ausgezeichneten Objekte unter den Gezogenen.

Man sagt: ${\textstyle A}$ ist hypergeometrisch verteilt mit ${\textstyle K}$ Ausgezeichneten bei ${\textstyle N}$ Objekten insgesamt (bzw. mit ${\textstyle N-K}$ Nicht-Ausgezeichneten) und ${\textstyle n}$ Gezogenen.
Die möglichen Werte von ${\textstyle A}$ sind dann ${\textstyle 0,\ldots ,n}$ und es gilt: $$P(A=k)={\frac {{K \choose k}{N-K \choose n-k}}{N \choose n}}\quad {\text{für }}k\in \{0,\ldots ,n\}$$

Beispiel 1

Für ${\textstyle N=70,\ K=12}$ und ${\textstyle n=20}$ ist beispielsweise:
$${\begin{array}{rcccl}P(A=4)&=&0.244497\\P(A=10)&=&0.000021\\P(A=15)&=&0\end{array}}$$

Beispiel 2.1

Für ${\textstyle N=12,\ K=7,\ n=8}$ ist:
$${\begin{array}{c||c|c|c|c|c|c|c|c|c}k&0&1&2&3&4&5&6&7&8\\\hline &&&&&&&&&\\P(A=k)&&&&&&&&&\\\parallel &\parallel &\parallel &\parallel &\parallel &\parallel &\parallel &\parallel &\parallel \\&0&0&0&0.071&0.354&0.424&0.141&0.010&0\end{array}}$$

Beispiel 2.2

Beispiel 3

Hier einige weitere Beispiele:

Interkatives Beispiel

Interaktive Shiny-App zur Hypergeometrischen Verteilung:
Download und Link

Wahrscheinlichkeiten

Es folgt:
$$P(A\leq k)=\sum \limits _{j=0}^{k}{\frac {{K \choose j}{N-K \choose n-j}}{N \choose n}}$$
$$P(A\geq k)=\sum \limits _{j=k}^{n}{\frac {{K \choose j}{N-K \choose n-j}}{N \choose n}}$$
$$P(k\leq A\leq \ell )=\sum \limits _{j=k}^{\ell }{\frac {{K \choose j}{N-K \choose n-j}}{N \choose n}}$$

Beispiel 4

Für ${\textstyle N=70,\ K=12}$ und ${\textstyle n=20}$ ist beispielsweise:
$${\begin{array}{rclcl}P(3\leq A\leq 5)&=&0.658525\\&&&&\\P(A\leq 4)&=&0.778049\\&&&&\\P(A\geq 4)&=&0.466448\end{array}}$$

Berechnung in R

In R:

$${\begin{array}{|c|c|c|c|}\hline \quad \color {blue}{{\text{dhyper(}}k,K,N-K,n)}&{\text{ergibt:}}&P(A=k)&=&{\frac {{K \choose k}{N-K \choose n-k}}{N \choose n}}\\\hline \quad \color {blue}{{\text{phyper(}}k,K,N-K,n)}&{\text{ergibt:}}&P(A\leq k)&=&\sum \limits _{j=0}^{k}{\frac {{K \choose j}{N-K \choose n-j}}{N \choose n}}\\\hline \quad \color {blue}{1-{\text{phyper(}}k-1,K,N-K,n)}&{\text{ergibt:}}&P(A\geq k))&=&\sum \limits _{j=k}^{n}{\frac {{K \choose j}{N-K \choose n-j}}{N \choose n}}\\\hline \quad \color {blue}{{\text{phyper(}}\ell ,K,N-K,n)-{\text{phyper(}}k-1,K,N-K,n)}&{\text{ergibt:}}&P(k\leq A\leq \ell )&=&\sum \limits _{j=k}^{\ell }{\frac {{K \choose j}{N-K \choose n-j}}{N \choose n}}\\\hline \end{array}}$$

Aufgabe 1

Berechnen Sie für eine hypergeometrisch verteilte ZV ${\textstyle A}$ mit den jeweils angegebenen Werten für ${\textstyle N,K}$ und ${\textstyle n}$ die angegebenen Wahrscheinlichkeiten:

• Für ${\textstyle N=14,\ K=4}$ und ${\textstyle n=7}$: ${\textstyle \quad P(A=k)}$ für alle ${\textstyle k=0,\ldots ,7}$
• Für ${\textstyle N=25,\ K=8}$ und ${\textstyle n=8}$: ${\textstyle \quad P(A\leq 2),\ P(A\geq 4),\ P(1\leq A\leq 3)}$
• Für ${\textstyle N=140,\ K=30}$ und ${\textstyle n=20}$: ${\textstyle \quad P(A\leq 7),\ P(A\geq 5),\ P(3\leq A\leq 6)}$

Beispiel 5

• (Ziehen ohne Zurücklegen) Aus einer Lostrommel, die ${\textstyle N}$ Kugeln enthält, von denen ${\textstyle K}$ rot sind, werden ohne Zurücklegen ${\textstyle n}$ Kugeln gezogen. Die ZV für die Anzahl der roten Kugeln unter den Gezogenen ist hypergeometrisch verteilt.

Beispiel 6

• In einem Teich befinden sich ${\textstyle 350}$ Fische einer Art, von denen ${\textstyle 80}$ markiert sind. Nun werden ${\textstyle 45}$ Fische gefangen. Die ZV für die Zahl der markierten Fische unter den Gefangenen ist hypergeometrisch verteilt mit ${\textstyle N=350,\ K=80}$ und ${\textstyle n=45}$. (Voraussetzung: Die markierten Fische sind über den See gleichmäßig verteilt und lassen sich genauso leicht fangen, wie die Übrigen.)

Beispiel 7

• In einer Klasse befinden sich ${\textstyle 12}$ Jungen und ${\textstyle 15}$ Mädchen. Es werden ${\textstyle 10}$ Schüler/innen für ein Projekt ausgelost. Die ZV, die die Zahl der Jungen unter den Ausgelosten angibt, ist hypergeometrisch verteilt mit ${\textstyle N=27,\ K=12}$ und ${\textstyle n=10}$.

Aufgabe 2

Wie groß ist beim Lotto (6 aus 49) die Wahrscheinlichkeit, genau ${\textstyle k}$ Richtige zu haben (${\textstyle k=0,\ldots ,6}$).

Aufgabe 3

Bei einem Multiple-Choice Test gibt es 20 Aussagen, von denen genau 10 richtig sind. Ein unvorbereiteter Teilnehmer kreuzt willkürlich genau 10 Aussagen an. Wie groß ist die Wahrscheinlichkeit, dass dabei $${\text{(i)}}\ {\text{mindestens 6}}\quad \quad {\text{(ii)}}\ {\text{weniger als 4}}\quad \quad {\text{(iii)}}\ {\text{zwischen 1 und 5}}$$ der angekreuzten Aussagen richtig sind?

Aufgabe 4

Unter 500 Glühbirnen in einem Karton befinden sich 35 defekte. Bei einer Qualitätskontrolle werden 50 Birnen getestet. Wie groß${\textstyle \;}$ist die Wahrscheinlichkeit, dass $${\text{(i)}}\ {\text{mindestens 4}}\quad \quad {\text{(ii)}}\ {\text{zwischen 1 und 3}}\quad \quad {\text{(iii)}}\ {\text{keine}}$$ der Birnen defekt ist?

Erwartungswert und Varianz

Für eine hypergeometrisch verteilte ZV ${\textstyle A}$ mit ${\textstyle N,K,n}$ wie bisher gilt:
$$E(A)=n\cdot {\frac {K}{N}}\quad {\text{und}}\quad V(A)=n\cdot {\frac {K}{N}}\cdot \left(1-{\frac {K}{N}}\right)\cdot {\frac {N-n}{N-1}}$$

Beispiel 8

Für ${\textstyle N=9,\ K=4}$ und ${\textstyle n=6}$ haben wir oben bereits die Wahrscheinlichkeitsverteilung bestimmt. Daraus ergibt sich: $${\begin{array}{rcccl}E(A)&=&2.667\\V(A)&=&0.5556\end{array}}$$

Beispiel 9

Für ${\textstyle N=30,\ K=18}$ und ${\textstyle n=12}$ berechnen wir zunächst
${\textstyle P(A=k)={\frac {{18 \choose k}{12 \choose 12-k}}{30 \choose 12}}}$
für alle möglichen Werte ${\textstyle k=0,\ldots ,12}$:
$${\begin{array}{|c||c|c|c|c|c|c|c|}\hline \hline k&0&1&2&3&4&5&6\\P(T=k)&<0.001&<0.001&0.0001&0.0021&0.0175&0.0785&0.1983\\\hline \hline k&7&8&9&10&11&12&\\P(A=k)&0.2914&0.2504&0.1237&0.0334&0.0044&0.0002&\\\hline \hline \end{array}}$$
Daraus ergibt sich:
$${\begin{array}{rclclccclcl}E(A)&=&7.2\\V(A)&=&1.7876\end{array}}$$

Aufgabe 5

Bestimmen Sie für eine hypergeometrsich verteilte ZV ${\textstyle Z}$ mit ${\textstyle N=100}$, ${\textstyle K=30}$ und ${\textstyle n=20}$ die nachfolgenden Werte:

1. ${\textstyle P(Z=10)}$
2. ${\textstyle P(Z\leq 17)}$
3. ${\textstyle P(Z\geq 11)}$
4. ${\textstyle P(Z=20)}$

Bestimmen Sie Erwartungswert und Varianz.

Aufgabe 6

Wenn für die Vorlesung 79 Menschen angemeldet sind, von denen ca. 30 Personen auch regelmäßig in die Vorlesung kommen. Die Klausur wird von 40 Personen geschrieben. Wie hoch ist die Wahrscheinlichkeit, dass

1. alle 30 Personen aus der Vorlesung die Klausur mitschreiben?
2. mindestens 20 Personen aus der Vorlesung die Klausur mitschreiben?
3. maximal 20 Personen aus der Vorlesung die Klausur mitschreiben?

Bestimmen Sie auch Erwartungswert und Varianz.

Schätzung der Zahl der ausgezeichneten Objekte K

Beispiel 1

Von ${\textstyle N=1000}$ Glübirnen einer Lieferung sind eine unbekannte Anzahl ${\textstyle K}$ defekt. Man testet ${\textstyle n=15}$ zufällig ausgewählte Birnen und stellt fest, dass ${\textstyle k=2}$ davon defekt sind. Wie kann man daraus auf die Zahl ${\textstyle K}$ schließen?

Situation

Situation:

$${\text{Es sind }}N,n{\text{ und }}k{\text{ bekannt, aber nicht }}K.{\text{Wie kann man }}K{\text{ sinnvoll schätzen?}}$$

Genauer:

• ${\textstyle N\in \mathbb {N} }$ und ${\textstyle n\in \{1,\ldots ,N\}}$ sind feststehend und bekannt. Oft kann man ${\textstyle n}$ selbst festlegen.
• ${\textstyle k\in \{0,\ldots ,n\}}$ entsteht zufällig, ist dann aber bekannt.
• ${\textstyle K\in \{0,\ldots ,N\}}$ steht fest, ist aber nicht bekannt.

$${\text{System mit Parameter }}K{\stackrel {\text{zufällig}}{\longrightarrow }}{\text{Daten }}k{\stackrel {\text{methodisch}}{\longrightarrow }}{\text{Schätzung }}{\hat {K}}{\text{ für }}K$$ Wiederum ist dabei folglich die Schätzung zufällig.

Erwartungstreue Punktschätzungen für K

Durch ${\textstyle {\hat {K}}={\frac {k\cdot N}{n}}}$ erhält man eine erwartungstreue Schätzung für ${\textstyle K}$.

Genauer: Die Zahl ${\textstyle k}$ hängt vom Zufall ab und wird (vor der Datenerhebung) durch die ZV ${\textstyle A}$ beschrieben. Da die Schätzung für (die feste aber unbekannte Zahl) ${\textstyle K}$ von ${\textstyle k}$ abhängt, ist sie ebenfalls vom Zufall abhängig. Die Schätzung ${\textstyle {\hat {K}}={\frac {A\cdot N}{n}}}$ kann somit als ZV beschrieben werden. Dabei gilt dann (unabhängig vom unbekannten Wert
${\textstyle K\in \{0,\ldots ,N\}}$) stets ${\textstyle E({\hat {K}})=E\left({\frac {A\cdot N}{n}}\right)=K}$.

Beispiel 2

$${\begin{array}{ccccrcccl}N=15&n=8&k=3&\Rightarrow &{\hat {K}}&=&5.625\\N=40&n=12&k=3&\Rightarrow &{\hat {K}}&=&33.33\\N=1000&n=15&k=2&\Rightarrow &{\hat {K}}&=&133.33\end{array}}$$

Das Maximum-Likelihood-Prinzip 1

Mit der Maximum-Likelihood-Methode wird ${\textstyle K\in \{0,\ldots ,N\}}$ (basierend auf der zufälligen Zahl ${\textstyle k}$) so geschätzt, dass die Wahrscheinlichkeit ${\textstyle P(A=k)}$ maximal wird. Wir suchen also die Maximumstelle der Likelihood-Funktion
$$L:\{0,\ldots ,N\}\to [0,1],\ L(K)={\frac {{K \choose k}\cdot {N-K \choose n-k}}{N \choose n}}$$

Das Maximum-Likelihood-Prinzip 2

Man stellt fest:
Die Maximumstelle(n) von ${\textstyle L}$ ist/sind:
$$\left\{{\begin{array}{lcc}{\hat {K}}_{1}={\frac {k\cdot (N+1)}{n}}-1\quad {\text{und}}\quad {\hat {K}}_{2}={\frac {k\cdot (N+1)}{n}}&,&{\text{falls}}\ {\frac {k\cdot (N+1)}{n}}\in \mathbb {N} \ {\text{ist.}}\\{\hat {K}}=\left\lfloor {\frac {k\cdot (N+1)}{n}}\right\rfloor &,&{\text{falls}}\ {\frac {k\cdot (N+1)}{n}}\notin \mathbb {N} \ {\text{ist.}}\end{array}}\right\}$$
(dabei bezeichnet ${\textstyle \lfloor x\rfloor }$ die größte ganze Zahl, die kleiner oder gleich ${\textstyle x}$ ist)

Beispiel 3.1

Beispiel 3.2

$N=15,n=8,k=3:{\hat {K_{1}}}={\frac {3\cdot 16}{8}}-1=5$ und $\;{\hat {K_{2}}}={\frac {3\cdot 16}{8}}=6$
$N=40,n=12,k=10:{\hat {K}}=\left\lfloor {\frac {10\cdot 41}{12}}\right\rfloor =34$
$N=1000,n=15,k=2:{\hat {K}}=\left\lfloor {\frac {2\cdot 1001}{15}}\right\rfloor =133$

Intervallschätzungen für K

Ziel

Gib eine Methode an, mit der man aus ${\textstyle k}$ ein Intervall ${\textstyle [K_{U},K_{O}]}$ bestimmen kann, so dass die Wahrscheinlichkeit, dass sich ein Intervall ergibt, das ${\textstyle K}$ enthält, garantiert (also für jeden denkbaren Wert von ${\textstyle K}$) größer oder gleich einem vorgegebenen Konfidenzniveau ${\textstyle \delta }$ ist.
$${\text{System mit Parameter }}K{\stackrel {\text{zufällig}}{\longrightarrow }}{\text{Daten }}k$$

$${\stackrel {\text{methodisch}}{\longrightarrow }}{\text{(Intervall-)Schätzung }}[K_{U},K_{O}]{\text{ für }}K$$ $${\textbf {Ziel:}}{\text{ Für jeden möglichen Wert von }}K:\quad P\left([K_{U},K_{O}]\ni K\right)\geq \delta$$ Eine sinnvolle Möglichkeit wird im Folgenden beschrieben:

IVS für K, die ein gegebenes Konfidenzniveau einhält

Gegeben seien ${\textstyle N\in \mathbb {N} }$ und ${\textstyle n\in \{1,\ldots ,n\}}$. Unbekannt sei ${\textstyle K\in \{0,\ldots ,N\}}$. Weiter sei ein Konfidenzniveau ${\textstyle \delta \in ]0,1[}$ vorgegeben.
Basierend auf der zufälligen Zahl ${\textstyle k}$ geht man nun wie folgt vor:

Bestimmung von Ku

• Man bestimmt ${\textstyle K_{U}}$ als die kleinstmögliche Zahl mit $$\color {blue}{{\text{phyper(}}k-1,K_{U},N-K_{U},n)}$$${\textstyle =\sum \limits _{j=0}^{k-1}{\frac {{K_{U} \choose j}\cdot {N-K_{U} \choose n-j}}{N \choose n}}<{\frac {1+\delta }{2}}}$

Bestimmung von Ko

• Man bestimmt ${\textstyle K_{O}}$ als die größtmögliche Zahl mit $$\color {blue}{{\text{phyper(}}k-1,K_{O},N-K_{O},n)}$$ ${\textstyle =\sum \limits _{j=0}^{k}{\frac {{K_{O} \choose j}\cdot {N-K_{O} \choose n-j}}{N \choose n}}>{\frac {1-\delta }{2}}}$

Konfidenzniveau

Ohne weiter in die mathematischen Hintergründe einzusteigen, halten wir fest, dass die folgende (bei Intervallschätzungen immer zu erreichende) Bedingung bei diesem Verfahren garantiert erfüllt ist:
$$P\left([K_{U},K_{O}]\ni K\right)\geq \delta \quad {\text{(D.h. das Konfidenzniveau }}\delta {\text{ wird eingehalten.)}}$$
Man beachte, dass der Aussage "${\textstyle K\in [K_{U},K_{O}]}$${\textstyle \;}$" eine Wahrscheinlichkeit zugeschrieben werden kann, weil die Intervallgrenzen ${\textstyle K_{U}}$ und ${\textstyle K_{O}}$ zufällig sind (und nicht etwa der unbekannte Wert ${\textstyle K}$).

Beispiel 1.1

Wir betrachten erneut den Fall ${\textstyle n=1000,n=15,k=2}$ und führen eine Intervallschätzung zum Niveau ${\textstyle \delta =0.95}$ durch.

• Wir suchen also zunächst die kleinstmögliche Zahl ${\textstyle K_{U}}$ mit $$\quad \color {blue}{{\text{phyper(}}1,K_{U},1000-K_{U},15)}$$${\textstyle <0.975}$
  Durch Ausprobieren findet man:
  ${\textstyle \quad K_{U}=17}$

Beispiel 1.2

• Analog suchen wir die größtmögliche Zahl mit $$\quad \color {blue}{{\text{phyper(}}2,K_{O},1000-K_{O},15)}$$${\textstyle >0.025}$
  Durch Ausprobieren findet man: ${\textstyle \quad K_{0}=402}$
  

Damit ist ${\textstyle [K_{U},K_{O}]=[17,402]}$ das gesuchte Konfidenzintervall zu ${\textstyle \delta =0.95}$.

Beispiel 2.1

Für ${\textstyle N=200}$ und ${\textstyle n=18}$ berechnet man abhängig von ${\textstyle k}$ die folgenden ML-Schätzungen und Intervallschätzungen zum Vertrauensniveau ${\textstyle \delta =0.7}$:
$${\begin{array}{|c||c|c|c|c|c|c|c|c|c|c|}\hline k&0&1&2&3&4&5&6&7&8&9\\\hline {\hat {K}}&0&11&22&33&44&55&67&78&89&100\\\hline K\in &[0,19]&[2,34]&[8,47]&[16,60]&[24,72]&[33,84]&[43,95]&[52,106]&[62,117]&[73,127]\\\hline \end{array}}$$
$${\begin{array}{|c||c|c|c|c|c|c|c|c|c|}\hline k&10&11&12&13&14&15&16&17&18\\\hline {\hat {K}}&111&122&134&145&156&167&178&189&201\\\hline K\in &[83,138]&[94,148]&[105,157]&[116,167]&[128,176]&[140,184]&[153,192]&[166,198]&[181,200]\\\hline \end{array}}$$

Beispiel 2.2

• Angenommen, es ist ${\textstyle K=120}$. Dann ist die Intervallschätzung für ${\textstyle k\in \{9,10,12,12,13\}}$ korrekt. Die Wahrscheinlichkeit dafür ist:
  $$P(9\leq A\leq 13)=0.792$$
• Angenommen, es ist ${\textstyle K=48}$. Dann ist die Intervallschätzung für ${\textstyle k\in \{3,4,5,6\}}$ korrekt. Die Wahrscheinlichkeit dafür ist:
  $$P(3\leq A\leq 6)=0.749$$

Beispiel 2.3

• Angenommen, es ist ${\textstyle K=199}$. Dann ist die Intervallschätzung nur für ${\textstyle k=18}$ korrekt. Die Wahrscheinlichkeit dafür ist:
  

$$P(A=18)=0.91$$

Korrektheit der IVS

Das mathematische Modell garantiert, dass die Intervallschätzung bei beliebigem ${\textstyle K}$ immer mindestens mit der Wahrscheinlichkeit ${\textstyle \delta }$ korrekt ist.

Aufgabe 1

In Ihrem Wohnort stehen ${\textstyle N=153}$ Wohngebaude. Sie wissen, dass nur Eines von Vieren einen Keller hat.

1. Führen Sie für ${\textstyle K}$ eine Punktschätzung durch, stellen Sie die Maximum-Likelihood-Funktion auf und plotten Sie diese in R.
2. Geben Sie die Formeln für die Intervallschätzung mit ${\textstyle \delta =0.9}$ für ${\textstyle K}$ an.

Schätzung der Gesamtzahl der Objekte N

Beispiel (Capture-Recapture)

In einem See befindet sich eine unbekannte Anzahl ${\textstyle N}$ von Fischen einer Art. Man möchte wissen, wie groß ${\textstyle N}$ in etwa ist. Dazu fängt man eine (kleinere) Anzahl ${\textstyle K}$ von Fischen und markiert sie. Dann setzt man sie wieder aus und wartet einen angemessenen Zeitraum. Dann fängt man in einem zweiten Fischzug ${\textstyle n}$ Fische und bestimmt die Anzahl ${\textstyle k}$ der markierten Fische unter ihnen.
Beispielsweise hat man ${\textstyle K=100}$ Fische markiert und unter ${\textstyle n=50}$ gefangenen Fischen ${\textstyle k=14}$ markierte Fische wiedergefunden.

Wie kann man daraus eine sinnvolle Schätzung für ${\textstyle N}$ abgeben?

Situation

Situation:

$${\text{ Es sind }}K,n{\text{ und }}k{\text{bekannt, aber nicht }}N.{\text{Wie kann man }}N{\text{ sinnvoll schätzen?}}$$ Genauer:

• ${\textstyle K\in \mathbb {N} }$ und ${\textstyle n\in \mathbb {N} }$ sind fest und bekannt. Manchmal kann man ${\textstyle K}$ und ${\textstyle n}$ selbst festlegen.
• ${\textstyle k\in \{0,\ldots ,n\}}$ entsteht zufällig, ist dann aber bekannt.
• ${\textstyle N\in \mathbb {N} }$ mit ${\textstyle N\geq \max(n,K)}$ steht fest, ist aber nicht bekannt.

$${\text{System mit Parameter }}N{\stackrel {\text{zufällig}}{\longrightarrow }}{\text{Daten }}k{\stackrel {\text{methodisch}}{\longrightarrow }}{\text{Schätzung }}{\hat {N}}{\text{ für }}N$$ Wiederum ist damit die Schätzung zufällig.

Erwartungstreue Punktschätzungen für N

Durch ${\textstyle {\hat {N}}={\frac {K\cdot n}{k}}}$ erhält man eine Schätzung für ${\textstyle N}$.

Dabei gilt: ${\textstyle {\frac {1}{\hat {N}}}={\frac {k}{K\cdot n}}}$ ist erwartungstreu für ${\textstyle {\frac {1}{N}}}$.

Genauer: Die Zahl ${\textstyle k}$ hängt vom Zufall ab und wird (vor der Datenerhebung) durch die ZV ${\textstyle A}$ beschrieben. Da die Schätzung für (die feste aber unbekannte Zahl) ${\textstyle N}$ von ${\textstyle k}$ abhängt, ist sie ebenfalls vom Zufall abhängig. Die Schätzung ${\textstyle {\hat {N}}={\frac {K\cdot n}{A}}}$ kann somit als ZV beschrieben werden. Dabei gilt dann (unabhängig vom unbekannten Wert ${\textstyle N\in \mathbb {N} }$) stets ${\textstyle E\left({\frac {1}{\hat {N}}}\right)=E\left({\frac {A}{K\cdot n}}\right)={\frac {1}{N}}}$.

Beispiel 1

$${\begin{array}{ccccrcccl}K=25&n=16&k=10&\Rightarrow &{\hat {N}}&=&40\\K=25&n=16&k=10&\Rightarrow &{\hat {N}}&=&133.33\\K=100&n=50&k=14&\Rightarrow &{\hat {N}}&=&357.14\end{array}}$$

Schätzung für N mit der Maximum-Likelihood-Methode

Mit der Maximum-Likelihood-Methode wird ${\textstyle N\in \{\max(K,n),\ldots \}}$ (basierend auf der zufälligen Zahl ${\textstyle k}$) so geschätzt, dass die Wahrscheinlichkeit ${\textstyle P(A=k)}$ maximal wird. Wir suchen also die Maximumstelle der Likelihood-Funktion
$$L:\{\max(K,n),\ldots \}\to [0,1],\ L(N)={\frac {{K \choose k}\cdot {N-K \choose n-k}}{N \choose n}}$$

Maximumstellen

Man stellt fest:
Die Maximumstelle(n) ist/sind von ${\textstyle L}$: $$\left\{{\begin{array}{lcc}{\hat {N}}_{1}={\frac {K\cdot n}{k}}-1\quad {\text{und}}\quad {\hat {N}}_{2}={\frac {K\cdot n}{k}}&,&{\text{falls}}\ {\frac {K\cdot n}{k}}\in \mathbb {N} \ {\text{ist.}}\\{\hat {N}}=\left\lfloor {\frac {K\cdot n}{k}}\right\rfloor &,&{\text{falls}}\ {\frac {K\cdot n}{k}}\notin \mathbb {N} \ {\text{ist.}}\end{array}}\right\}$$ (dabei bezeichnet ${\textstyle \lfloor x\rfloor }$ die größte ganze Zahl, die kleiner oder gleich ${\textstyle x}$ ist)

Beispiel 2.1

Beispiel 2.2

$K=25,n=16,k=10:{\hat {N}}=\left\lfloor {\frac {25\cdot 16}{10}}\right\rfloor =40$
$K=25,n=16,k=3:{\hat {N}}=\left\lfloor {\frac {25\cdot 16}{3}}\right\rfloor =133$
$K=100,n=50,k=14:{\hat {N}}=\left\lfloor {\frac {100\cdot 50}{14}}\right\rfloor =357$

Intervallschätzungen für N

Ziel

Gib eine Methode an, mit der man aus ${\textstyle k}$ ein Intervall ${\textstyle [N_{U},N_{O}]}$ bestimmen kann, so dass die Wahrscheinlichkeit, dass sich ein Intervall ergibt, das ${\textstyle N}$ enthält, auf jeden Fall (also für jeden denkbaren Wert von ${\textstyle N}$) mindestens ein vorgegebenes Konfidenzniveau ${\textstyle \delta }$ ist.
$${\text{System mit Parameter }}N{\stackrel {\text{zufällig}}{\longrightarrow }}{\text{Daten }}k$$$${\stackrel {\text{methodisch}}{\longrightarrow }}{\text{(Intervall-)Schätzung }}[N_{U},N_{O}]{\text{ für }}N$$
$${\textbf {Ziel:}}{\text{ Für jeden möglichen Wert von }}N:\quad P\left([N_{U},N_{O}]\ni N\right)\geq \delta$$
Eine sinnvolle Möglichkeit wird im Folgenden beschrieben.

IVS für N, die ein gegebenes Konfidenzniveau einhält

Gegeben seien ${\textstyle K,n\in \{1,\ldots ,n\}}$. Unbekannt sei ${\textstyle N\in \{\max(K,n),\ldots ,N\}}$. Weiter sei ein Konfidenzniveau ${\textstyle \delta \in ]0,1[}$ vorgegeben.
Basierend auf der zufälligen Zahl ${\textstyle k}$ geht man nun wie folgt vor:

Bestimmung von Nu

• Man bestimmt ${\textstyle N_{U}}$ als die kleinstmögliche Zahl mit
  $$\quad \color {blue}{{\text{phyper(}}k,K,N_{U}-K,n)}$$${\textstyle =\sum \limits _{j=0}^{k}{\frac {{K \choose j}\cdot {N_{U}-K \choose n-j}}{N_{U} \choose n}}>{\frac {1-\delta }{2}}}$

Bestimmung von No

• Man bestimmt ${\textstyle N_{O}}$ als die größtmögliche Zahl mit
  $$\quad \color {blue}{{\text{phyper(}}k-1,K,N_{O}-K,n)}$$${\textstyle =\sum \limits _{j=0}^{k-1}{\frac {{K \choose j}\cdot {N_{O}-K \choose n-j}}{N_{O} \choose n}}<{\frac {1+\delta }{2}}}$

Konfidenzniveau

Wir halten fest, dass die folgende (bei Intervallschätzungen immer zu erreichende) Bedingung bei diesem Verfahren garantiert erfüllt ist: $$P\left([N_{U},N_{O}]\ni N\right)\geq \delta \quad {\text{(D.h. das Konfidenzniveau }}\delta {\text{ wird eingehalten.)}}$$ Man beachte, dass der Aussage "${\textstyle N\in [N_{U},N_{O}]}$"${\textstyle \;}$ eine Wahrscheinlichkeit zugeschrieben werden kann, weil die Intervallgrenzen ${\textstyle N_{U}}$ und ${\textstyle N_{O}}$ zufällig sind (und nicht etwa der unbekannte Wert ${\textstyle N}$).

Beispiel 1.1

Wir betrachten erneut den Fall ${\textstyle K=100,n=50,k=14}$ und führen eine Intervallschätzung zum Niveau ${\textstyle \delta =0.8}$ durch.

• Wir suchen also zunächst die kleinstmögliche Zahl ${\textstyle N_{U}}$ mit $$\quad \color {blue}{{\text{phyper(}}14,100,N_{U}-100,50)}$$${\textstyle >0.1}$ Durch Ausprobieren findet man: ${\textstyle \quad N_{U}=272}$

Beispiel 1.2

• Analog suchen wir die größtmögliche Zahl mit$$\quad \color {blue}{{\text{phyper(}}13,100,N_{O}-100,50)}$$ ${\textstyle <0.9}$ Durch Ausprobieren findet man: ${\textstyle \quad N_{0}=499}$
  

Damit ist ${\textstyle [N_{U},N_{O}]=[272,499]}$ das gesuchte Konfidenzintervall zu ${\textstyle \delta =0.8}$.

Aufgabe 1

Einige Zeit nach einem Wiederansiedlungsversuch einer Spezies wollen Sie wissen, ob dieser geglückt ist und die Spezies sich vermehrt hat.

1. Geben Sie an, wie Sie zur Überprüfung vorgehen würden (Sie können nicht alle Exemplare zählen).
2. Seien nun ${\textstyle K=18,n=17\;{\text{und}}\;k=7}$. Bestimmen Sie ${\textstyle N}$ mittels der einfachen Punktschätzung, stellen Sie die Maximum-Likelihood-Funktion auf, plotten diese in R und führen eine Intervallschätzung mit ${\textstyle \delta =0.9}$ durch.
3. Was fällt Ihnen hinsichtlich der verschiedenen Schätzungen auf? Vergleichen Sie Ihre Ergebnisse.

Seiteninformation

Diese Lernresource können Sie als Wiki2Reveal-Foliensatz darstellen.

Wiki2Reveal

Dieser Wiki2Reveal Foliensatz wurde für den Lerneinheit Kurs:Statistik für Anwender' erstellt der Link für die Wiki2Reveal-Folien wurde mit dem Wiki2Reveal-Linkgenerator erstellt.

• Die Seite wurde als Dokumententyp PanDocElectron-SLIDE erstellt.
• Link zur Quelle in Wikiversity: https://de.wikiversity.org/wiki/Kurs:Statistik%20f%C3%BCr%20Anwender/Hypergeometrischverteilte%20Zufallsvariable
• siehe auch weitere Informationen zu Wiki2Reveal und unter Wiki2Reveal-Linkgenerator.