Stichprobenverteilung

In der Statistik und Wahrscheinlichkeitstheorie wird die Wahrscheinlichkeitsverteilung einer Stichprobenfunktion auch als Stichprobenverteilung der Stichprobenfunktion bezeichnet.

Zur Bezeichnung „Stichprobenverteilung“

Die Bezeichnung „Stichprobenverteilung“ ist zwar gebräuchlich, kann aber missverständlich sein, da damit nicht die Verteilung einer Zufallsstichprobe oder der möglichen Stichprobenwerte bezeichnet wird. Während ‚Verteilung der Stichprobenfunktion‘ unmissverständlich ist, kann die Bezeichnung ‚Stichprobenverteilung der Stichprobenfunktion‘ so missverstanden werden, dass es noch andere Verteilungen einer Stichprobenfunktion gibt. Diese beiden möglichen Missverständnisse können den Grund dafür bilden, dass viele Autoren die Bezeichnung Stichprobenverteilung völlig vermeiden und einfach von der Verteilung der Stichprobenfunktion sprechen.

Der Begriff sampling distribution wurde 1922 von Ronald Aylmer Fisher beiläufig eingeführt und in den Jahren 1928 und 1929 mit der Verwendung im Titel von zwei Aufsätzen etabliert.

Zu unterscheiden ist die Stichprobenverteilungsfunktion, eine seltene Bezeichnung für die empirische Verteilungsfunktion.

Bestimmung der Stichprobenverteilung

Wenn $X_1,\dots ,X_n$ die Stichprobenvariablen einer Zufallsstichprobe vom Umfang $n$ sind und $G=g(X_1,\dots ,X_n)$ eine Stichprobenfunktion ist, dann ist $G$ eine Zufallsvariable, deren Wahrscheinlichkeitsverteilung als Stichprobenverteilung von $G$ bezeichnet wird. Die Stichprobenverteilung der Stichprobenfunktion $G$ hängt über die (messbare) Funktion $g$ von der Wahrscheinlichkeitsverteilung des Stichprobenvektors $𝐗=(X_1,\dots ,X_n)$ ab.

Häufig interessierende Stichprobenfunktionen sind

• die Summenvariable ${\sum }_{i=1}^n{X}_i$, die in diesem Zusammenhang auch Stichprobensumme heißt,
• das arithmetische Mittel ${\overline{X}}_n=\frac{1}{n}\sum _{i=1}^n{X}_i/n$, das in diesem Zusammenhang auch Stichprobenmittel heißt,
• die mittlere quadratische Abweichung vom Stichprobenmittel $S_n^2=\frac{1}{n}\sum _{i=1}^n(X_i-{\overline{X}}_n{)}^2$, die in diesem Zusammenhang auch auch Stichprobenvarianz heißt,
• das Minimum $\min \{X_1,\dots ,X_n\}$ und,
• das Maximum $\max \{X_1,\dots ,X_n\}$.

Allgemeines Vorgehen

Wenn die Wahrscheinlichkeitsverteilung des Stichprobenvektors $𝐗$ und die Funktion $g:{ℝ}^n\to ℝ$ gegeben ist, dann ergibt sich die Verteilungsfunktion – und mit dieser die Stichprobenverteilung – der Stichprobenfunktion $G$ als

$P(G\le x)=P(𝐗\in g^{-1}((-\mathrm{\infty },x])\text{ für alle }x\in ℝ.$

Dabei ist $g^{-1}((-\mathrm{\infty },x])=\{𝐱\in {ℝ}^n\mid g(𝐱)\in (-\mathrm{\infty },x]\}$.

Wenn der häufige Fall vorliegt, dass die Stichprobenvariablen stochastisch unabhängig und identisch verteilt sind und jede Stichprobenvariable die Wahrscheinlichkeitsverteilung $P_X$ hat, dann ist die Wahrscheinlichkeitsverteilung des Zufallsvektor $𝐗$ durch die $n$-fache Produktverteilung $P_X^{\otimes n}$ gegeben.

Bei Anwendungen in der parametrischen Statistik mit stochastisch unabhängigen und identisch verteilten Stichprobenvariablen ist typischerweise die Verteilung einer Stichprobenvariablen nicht exakt bekannt, sondern durch eine parametrische Familie $(P_{X,\theta }{)}_{\theta \in \mathrm{\Theta }}$ eingeschränkt. In diesem Fall ergibt sich auch für $𝐗$ eine parametrische Familie $(P_{X,\theta }^{\otimes n}{)}_{\theta \in \mathrm{\Theta }}$ von Verteilungen und es ergibt sich nicht nur eine Stichprobenverteilung, sondern eine durch den Parameter $\theta$ indizierte Familie von Stichprobenverteilungen der Stichprobenfunktion $G$.

Vorgehen in speziellen Fällen

Für spezielle Verteilungen der Stichprobenvariablen und spezielle Stichprobenfunktionen sind Zusammenhänge bekannt, mit deren Hilfe die Stichprobenverteilung der jeweiligen Stichprobenfunktion angegeben werden kann, ohne den oben angegebenen allgemeinen Weg zu beschreiten.

Unter der Voraussetzung, dass die Stichprobenvariablen $X_1,\dots ,X_n$ stochastisch unabhängig und identisch verteilt sind, gilt z. B.:

• Die Stichprobenverteilung der Summe Bernoulli-verteilter Stichprobenvariablen ist eine Binomialverteilung;

$X_i\sim \mathrm{B}\mathrm{e}\mathrm{r}(p)⟹{\sum }_{i=1}^n{X}_i\sim \mathrm{B}\mathrm{i}\mathrm{n}(n,p)$.

• Die Stichprobenverteilung des Minimums und des Maximums Bernoulli-verteilter Stichprobenvariablen ist eine Bernoulli-Verteilung;

$X_i\sim \mathrm{B}\mathrm{e}\mathrm{r}(p)⟹\min \{X_1,\dots ,X_n\}\sim \mathrm{B}\mathrm{e}\mathrm{r}(p^n),$

$X_i\sim \mathrm{B}\mathrm{e}\mathrm{r}(p)⟹\max \{X_1,\dots ,X_n\}\sim \mathrm{B}\mathrm{e}\mathrm{r}(1-(1-p{)}^n).$

• Die Stichprobenverteilung der Summe normalverteilter Stichprobenvariablen ist eine Normalverteilung;

$X_i\sim \mathrm{N}(\mu ,{\sigma }^2)⟹{\sum }_{i=1}^n{X}_i\sim \mathrm{N}(n\mu ,n{\sigma }^2)$.

• Die Stichprobenverteilung des Stichprobenmittels normalverteilter Stichprobenvariablen ist eine Normalverteilung;

  

  Dichtefunktionen der Stichprobenverteilung des arithmetischen Mittels ${\overline{X}}_n$ stochastisch unabhängiger standardnormalverteilter Zufallsvariablen für unterschiedliche Stichprobenumfänge $n$. Mit größerem Stichprobenumfang $n$ wird die Varianz $\mathrm{V}\mathrm{a}\mathrm{r}[{\overline{X}}_n]=1/n$ kleiner und die Dichtefunktion ist enger um die Stelle $\mu =0$ konzentriert. Die Dichtefunktion der Stichprobenvariable ${\overline{X}}_n$ ist eine Dirac-Folge.

$X_i\sim \mathrm{N}(\mu ,{\sigma }^2)⟹{\overline{X}}_n\sim \mathrm{N}(\mu ,{\sigma }^2/n)$.

• Die Stichprobenverteilung des Produkts lognormalverteilter Stichprobenvariablen ist eine Lognormalverteilung;

$X_i\sim \mathrm{L}\mathrm{N}(\mu ,{\sigma }^2)⟹\prod _{i=1}^n{X}_i\sim \mathrm{L}\mathrm{N}(n\mu ,n{\sigma }^2)$;

• Die Stichprobenverteilung des geometrischen Mittels lognormalverteilter Stichprobenvariablen ist eine Lognormalverteilung;

$X_i\sim \mathrm{L}\mathrm{N}(\mu ,{\sigma }^2)⟹\sqrt[n]{\prod _{i=1}^n{X}_i}\sim \mathrm{L}\mathrm{N}(\mu ,{\sigma }^2/n)$.

• Die Stichprobenverteilung der Stichprobenvarianz $S_n^2$ normalverteilter Stichprobenvariablen ist durch eine Chi-Quadrat-Verteilung mit $n-1$ Freiheitsgraden bestimmt,

$X_i\sim \mathrm{N}(\mu ,{\sigma }^2)⟹\frac{n{S}_n}{{\sigma }^2}\sim {\chi }^2(n-1)$.

Es sind viele weitere ähnliche Zusammenhänge bekannt, die es ermöglichen, die Stichprobenverteilung bestimmter Stichprobenfunktionen unmittelbar anzugeben. Für die Stichprobenverteilung der Summe von Stichprobenvariablen anderer Verteilungen siehe auch Reproduktivitätseigenschaft.

Kennzahlen der Stichprobenverteilung

Wenn die Stichprobenvariablen $X_1,\dots ,X_n$ stochastisch unabhängig und identisch verteilt sind und Kennzahlen der Verteilung der Stichprobenvariablen bekannt sind, können Aussagen über Kennzahlen einer Stichprobenverteilung gemacht werden, z. B. gelten für den Erwartungswert und die Varianz der Summenvariable ${\sum }_{i=1}^n{X}_i$ und des arithmetischen Mittels ${\overline{X}}_n$ die Aussagen:

Aus $\mathrm{E}[X_i]=\mu \in ℝ$ folgt

$\mathrm{E}[{\sum }_{i=1}^n{X}_i]=n$

und

$\mathrm{E}[{\overline{X}}_n]=\mu$.

Aus $\mathrm{V}\mathrm{a}\mathrm{r}[X_i]={\sigma }^2<\mathrm{\infty }$ folgt

$\mathrm{V}\mathrm{a}\mathrm{r}[{\sum }_{i=1}^n{X}_i]=n{\sigma }^2,$

$\mathrm{V}\mathrm{a}\mathrm{r}[{\overline{X}}_n]=\frac{{\sigma }^2}{n}$

und

$\mathrm{E}[S_n^2]=\frac{n-1}{n}{\sigma }^2.$

Approximative Stichprobenverteilung

Wenn die Stichprobenvariablen $X_1,\dots ,X_n$ stochastisch unabhängig und identisch verteilt sind und wenn $\mathrm{E}[X_i]=\mu$ und $0<\mathrm{V}\mathrm{a}\mathrm{r}[X_i]={\sigma }^2<\mathrm{\infty }$ gilt, dann sind für hinreichend großes $n$ die Summe ${\sum }_{i=1}^n{X}_i$ und das arithmetische Mittel ${\overline{X}}_n$ approximativ normalverteilt (siehe Zentraler Grenzwertsatz, Gesetz der großen Zahlen, Gleichung von Bienaymé). Die Stichprobenverteilungen von ${\sum }_{i=1}^n{X}_i$ und ${\overline{X}}_n$ sind also näherungsweise durch Normalverteilungen charakterisierbar;

$\mathrm{V}\mathrm{e}\mathrm{r}\mathrm{t}\mathrm{e}\mathrm{i}\mathrm{l}\mathrm{u}\mathrm{n}\mathrm{g}\left[\sum _{i=1}^n{X}_i\right]\approx \mathrm{N}(n\mu ,n{\sigma }^2),\mathrm{V}\mathrm{e}\mathrm{r}\mathrm{t}\mathrm{e}\mathrm{i}\mathrm{l}\mathrm{u}\mathrm{n}\mathrm{g}[{\overline{X}}_n]\approx \mathrm{N}(\mu ,\frac{{\sigma }^2}{n}).$

Diese Approximationen beruhen auf dem Zentralen Grenzwertsatz der Statistik, der besagt, dass die Folge $(Z_n{)}_{n\in \mathbb{N}}$ der standardisierten Zufallsvariablen

$Z_n:=\frac{\sum _{i=1}^n{X}_i-\mu n}{\sqrt{n{\sigma }^2}}=\sqrt{n}\frac{{\overline{X}}_n-\mu }{\sigma }$

für $n\to \mathrm{\infty }$ in Verteilung gegen eine standardnormalverteilte Zufallsvariable $Z\sim \mathrm{N}(0,1)$ konvergiert. Die Approximation der Verteilung von $Z_n$ durch eine Standardnormalverteilung ist dann äquivalent zu den angegebenen Approximationen für die Stichprobenverteilungen von ${\sum }_{i=1}^n{X}_i$ und ${\overline{X}}_n$.

Statistische Schätzung der Stichprobenverteilung

Wenn die Stichprobenwerte $x_1,\dots ,x_n$ aus einer hinreichend große Zufallsstichprobe vorliegen, kann die empirische Verteilung der Stichprobenwerte als statistische Schätzung der Verteilung der Grundgesamtheit angesehen werden. Die Stichprobenverteilung einer beliebige Stichprobenfunktion kann dann ohne parametrisches Modell mit Hilfe des Bootstrap-Verfahrens geschätzt werden, ohne dass die Verteilung der Stichprobenvariablen bekannt sein muss. Jedoch muss allgemein mathematisch gezeigt werden, dass die Bootstrap-Stichprobenverteilungen mit steigender Zahl der Bootstrap-Stichproben gegen die Stichprobenverteilung konvergieren. Für das Beispiel im Bild ist die Bootstrap-Stichprobenverteilung um $\hat{\mu }=\frac{1}{n}\sum _{i=1}^n{x}_i$ zentriert, und im Allgemeinen nicht um $\mu$, hat jedoch die zu erwartende richtige Streuungsbreite.

Anwendungsbereiche

Die Wahrscheinlichkeitsverteilung einer Stichprobenfunktion dient in der statistischen Schätz- und Testtheorie zur Gewinnung von Aussagen über unbekannte Parameter in der Grundgesamtheit aufgrund einer Stichprobe.

Statistische Schätztheorie

In der statistischen Schätztheorie ist die interessierende Stichprobenverteilung häufig die Wahrscheinlichkeitsverteilung einer Schätzfunktion für einen unbekannten Parameter der Grundgesamtheit. Vorlage:Hauptartikel

Beispiel

Sind die Stichprobenvariablen $X_1,\dots ,X_n$ stochastisch unabhängig und identisch verteilte Bernoulli-Variablen mit dem unbekannten Bernoulli-Parameter $0<p<1$, dann ist die Stichprobenverteilung der Summenvariablen $S_n$ eine Binomialverteilung mit den Parametern $n$ und $p$. Die Stichprobenverteilung der Stichprobenvariablen $S_n$ hängt vom unbekannten Parameter $p$ ab.

Statistische Testtheorie

In der statistischen Testtheorie hängt die Verteilung einer Teststatistik typischerweise von den Stichprobenvariablen $X_1,\dots ,X_n$, dem Stichprobenumfang $n$ und einem spezifizierten Wert eines unbekannten Parameters der Grundgesamtheit ab.

Beispiel

Sind die Stichprobenvariablen stochastisch unabhängig und identisch normalverteilt mit unbekanntem Parameter $\mu$ und bekannter Varianz ${\sigma }^2=1$, dann ist die Zufallsvariable

$T_n=\sqrt{n}({\overline{X}}_n-{\mu }_0)$

die Teststatistik eines Gauß-Tests mit den Hypothesen $H_0:\mu ={\mu }_0$ und $H_1:\mu \ne {\mu }_0$. Die Stichprobenverteilung der Teststatistik $T_n$ ist die Normalverteilung mit dem Erwartungswert $\mu -{\mu }_0$ und der Varianz 1,

$T_n\sim \mathrm{N}(\mu -{\mu }_0,1).$

Die Stichprobenverteilung der Stichprobenfunktion $T_n$ hängt vom unbekannten Parameter $\mu$ ab. Wenn die Nullhypothese richtig ist, also $\mu -{\mu }_0$ gilt, dann ist $T_n$ standardnormalverteilt. Diese spezielle Stichprobenverteilung heißt dann auch die Verteilung der Teststatistik unter $H_0$.

Stichproben aus endlichen Grundgesamtheiten

Im statistischen Methodengebiet der Stichproben aus endlichen Grundgesamtheiten sind die Stichprobenvariablen typischerweise zwar identisch verteilt, aber nicht stochastisch unabhängig, wenn Schemata der Stichprobenziehung berücksichtigt werden, die nicht dem Schema Ziehen mit Zurücklegen entsprechen, das zu stochastisch unabhängigen und identisch verteilten Stichprobenvariablen führt.

Bei vielen Anwendungen erfolgt ein Ziehen ohne Zurücklegen, bei dem die Stichprobenvariablen zwar identisch verteilt, aber nicht stochastisch unabhängig sind. Wenn die Grundgesamtheit aus $N$ unterscheidbaren statistischen Einheiten in der Menge $\{1,2,\dots ,N\}$ besteht, so gibt es beim Ziehen einer Stichprobe vom Umfang $n=2$ ohne Zurücklegen insgesamt $N(N-1)$ mögliche Stichproben, nämlich Paare $\omega =({\omega }_1,{\omega }_2)$ in der Menge

$\mathrm{\Omega }=\{({\omega }_1,{\omega }_2)\mid {\omega }_1,{\omega }_2\in \{1,\dots ,N\},{\omega }_1\ne {\omega }_2\}$

der möglichen Stichproben. Durch dieses Ziehungsschema sind den möglichen Stichproben die Auswahlwahrscheinlichkeiten

$p(\omega )=\frac{1}{N(N-1)}\text{ für alle }\omega \in \mathrm{\Omega }$

zugeordnet. Wenn den Einheiten in $\{1,2,\dots ,N\}$ die Werte $({\xi }_1,\dots ,{\xi }_N)\in {ℝ}^N$ eines statistischen Merkmals zugeordnet sind, so wird in einer Stichprobe $\omega =({\omega }_1,{\omega }_2)\in \mathrm{\Omega }$ das Wertepaar $({\xi }_{{\omega }_1},{\xi }_{{\omega }_2})$ beobachtet. Die Stichprobenvariablen $(X_1,X_2)$, die die zufälligen Werte bei der ersten und zweiten Ziehung beschreiben, haben die gemeinsame zweidimensionale diskrete Wahrscheinlichkeitsverteilung

$P((X_1,X_2)=(x_1,x_2))=\frac{\mathrm{\#}\{\omega \in \mathrm{\Omega }\mid ({\xi }_{{\omega }_1},{\xi }_{{\omega }_2})=(x_1,x_2)\}}{N(N-1)},(x_1,x_2)\in {ℝ}^2.$

Für die Stichprobenfunktion $S=X_1+X_2$ ergibt sich die Stichprobenverteilung als die eindimensionale diskrete Wahrscheinlichkeitsverteilung

$P(S=x)=\frac{\mathrm{\#}\{\omega \in \mathrm{\Omega }\mid {\xi }_{{\omega }_1}+{\xi }_{{\omega }_2}=x\}}{N(N-1)},x\in ℝ.$

Für die Stichprobenfunktion $\overline{X}=(X_1+X_2)/2$ ergibt sich die Stichprobenverteilung als die eindimensionale diskrete Wahrscheinlichkeitsverteilung

$P(\overline{X}=x)=\frac{\mathrm{\#}\{\omega \in \mathrm{\Omega }\mid ({\xi }_{{\omega }_1}+{\xi }_{{\omega }_2})/2=x\}}{N(N-1)},x\in ℝ.$

Beispiel

Für eine Grundgesamt mit $N=3$ Einheiten und ${\xi }_1=10,{\xi }_2=16,{\xi }_3=10$ ergeben sich beim Ziehen ohne Zurücklegen die $N(N-1)=6$ verschiedenen Stichproben $(1,2),(1,3),(2,1),(2,3),(3,1),(3,2)$ vom Umfang $n=2$. Dies haben jeweils die Auswahlwahrscheinlichkeit $1/6$. Die zugehörigen Beobachtungswerte sind $({\xi }_1,{\xi }_2)=({\xi }_3,{\xi }_2)=(10,16)$, $({\xi }_1,{\xi }_3)=({\xi }_3,{\xi }_1)=(10,10)$ und $({\xi }_2,{\xi }_1)=({\xi }_2,{\xi }_3)=(16,10)$. Die gemeinsame diskrete Wahrscheinlichkeitsverteilung von $(X_1,X_2)$ ist

$P((X_1,X_2)=(x_1,x_2))=\{\begin{array}{cc}\frac{2}{6}=\frac{1}{3},& \text{falls }(x_1,x_2)\in \{(10,16),(10,10),(16,10)\}\\ 0& \text{sonst}\end{array}.$

Die Stichprobenverteilung der Stichprobenfunktion $S=X_1+X_2$ ist

$P(S=26)=\frac{4}{6}=\frac{2}{3},P(S=20)=\frac{2}{6}=\frac{1}{3}.$

Die Stichprobenverteilung der Stichprobenfunktion $\overline{X}=(X_1+X_2)/2$ ist

$P(\overline{X}=13)=\frac{2}{3},P(\overline{X}=10)=\frac{1}{3}.$

Bayesianische Inferenzstatistik

Bei der bayesianischen Inferenz wird die A-priori-Verteilung, die auf dem Parameterraum eines zu schätzenden Parameters definiert ist, unter Berücksichtigung der Wahrscheinlichkeitsverteilung des Stichprobenvektors $𝐗=(X_1,\dots ,X_n)$ und unter Verwendung eines realisierten und beobachteten Wertes $𝐱=(x_1,\dots ,x_n)$ des Stichprobenvektors in die A-posteriori-Verteilung transformiert. Dabei ist die A-posteriori-Verteilung proportional zum Produkt aus A-priori-Verteilung und der Likelihoodfunktion. Die Likelihoodfunktion gibt im diskreten Fall die Wahrscheinlichkeit und im stetigen Fall die Wahrscheinlichkeitsdichte des beobachten Wertes $𝐱$ für alternative Parameter an und wird als Funktion auf dem Parameterraum interpretiert.

Aus der A-posteriori-Verteilung können verschiedene Schätzwerte für den zu schätzenden Parameter gewonnen werden, indem z. B. im stetigen Fall der Wert mit maximalen Dichte bzw. im diskreten Fall der Wert mit maximaler Wahrscheinlichkeit, der Median oder der Erwartungswert der A-posteriori-Verteilung als Schätzwert verwendet werden. Ein solcher aus der A-posteriori-Verteilung gewonnener Schätzwert hängt von dem realisierten und beobachteten Wert $𝐱$ des Stichprobenvektors ab. Ersetzt man diesen durch den zufälligen Stichprobenvektor $𝐗$ so ergibt sich die zu dem jeweiligen Schätzwert gehörige Schätzfunktion, deren Wahrscheinlichkeitsverteilung eine Stichprobenverteilung ist.

Falls eine suffiziente Stichprobenfunktion (suffiziente Statistik) existiert, kann die Stichprobenverteilung dieser Stichprobenfunktion an die Stelle der Verteilung des Stichprobenvektors treten, ohne dass sich die resultierende A-posteriori-Verteilung ändert.

Beispiel

Die Stichprobenvariablen $X_1,\dots ,X_n$ seien stochastisch unabhängig und identisch Bernoulli-verteilt mit unbekanntem Bernoulli-Parameter $0<p<1$. Die Wahrscheinlichkeitsverteilung des Stichprobenvektors $𝐗=(X_1,\dots ,X_n)$ ist dann

$P_p(𝐗=𝐱)=\{\begin{array}{cc}{p}^k(1-p{)}^{n-k},& \text{falls }𝐱=(x_1,\dots ,x_n)\in \{0,1{\}}^n\\ 0& \text{sonst}\end{array},$

wobei $k:=\sum _{i=1}^n{x}_i$. Wenn $f_0$ die gegebene Dichtefunktion einer A-Priori-Verteilung auf dem Intervall $(0,1)$ ist und $𝐱$ ein realisierter und beobachteter Wert des Suchprobenvektors $𝐗=(X_1,\dots ,X_n)$ ist, dann ist die Dichtefunktion $f_1$ der A-Posteriori-Verteilung proportional zum Produkt aus $f_0$ und der Likelihoodfunktion

$L_{𝐱}(p):=P_p(𝐗=𝐱),0<p<1.$

Es gilt also

$f_1(p|𝐱)=c\cdot L_{𝐱}(p)\cdot f_0(p),0<p<1$.

Damit ist die A-priori-Verteilung spezifiziert. Dies kann für bestimmte Anwendungen ausreichend sein. Um einen Schätzwert für den Parameter $p$ zu erhalten, muss in einem zweiten Schritt eine Kennzahl der A-Posterior-Verteilung bestimmt werden.

Wenn $f_0$ konstant ist – dies ist die Dichtefunktion einer Beta-Verteilung mit den Parametern $\alpha =\beta =1$ –, dann ist die A-priori-Verteilung proportional zur Likelihoodfunktion und die Stelle mit maximaler A-Priori-Dichte ist der Maximum-Likelihood-Schätzwert $k/n$ für $p$. Die zugehörige Schätzfunktion für $p$ ist $K/n$ mit $K=\sum _{i=1}^n{X}_i$. Die A-Priori-Verteilung ist eine Beta-Verteilung mit den Parametern $\alpha =1+k$ und $\beta =1+n$ und dem Erwartungswert $(k+1)/(n+2)$, so dass $(k+1)/(n+2)$ ein alternativer bayesianischer Schätzwert mit der zugehörigen Schätzfunktion $(K+1)/(n+2)$ ist. Die Stichprobenverteilungen beider Schätzfunktionen sind durch eine Binomialverteilung bestimmt, da $K$ binomialverteilt ist,$K\sim \mathrm{Bin}(n,p)$.

Die Summe $K=\sum _{i=1}^n{X}_i$ der Stichprobenvariablen ist eine suffiziente Stichprobenfunktion für den Parameter $p$ mit der Stichprobenverteilung

$P_p(K=k)=\{\begin{array}{cc}(\genfrac{}{}{0ex}{}{n}{k})p^k(1-p{)}^{n-k},& \text{falls }k\in \{0,1,\dots ,n\}\\ 0& \text{sonst}\end{array}$

und der Likelihoodfunktion

$L_k(p):=P_p(K=k),0<p<1.$

Da sich die Likelihoodfunktionen $L_{𝐱}$ und $L_k$ nur durch einen konstanten Faktor unterscheiden, gilt auch

$f_1(p|𝐱)=c^{\prime }\cdot L_k(p)\cdot f_0(p),0<p<1$.

Die bayessche Inferenz beruhend auf der Verteilung des Stichprobenvektors und auf der Stichprobenverteilung der suffizienten Stichprobenfunktion führt zur selben A-Posteriori-Verteilung.

Zu diesem Beispiel siehe auch Bayessche Statistik#Bayessche Inferenz am Beispiel des Münzwurfes und Suffiziente Statistik#Beispiel: Binomialverteilung.

Anwendung

Anwendung findet die Stichprobenverteilung in der Herleitung von Konfidenzintervallen, siehe dort.

Siehe auch

• Konsistente Schätzfolge

Literatur

• Vorlage:Literatur
• Vorlage:Literatur
• Vorlage:Literatur
• Vorlage:Literatur