Diskretna uniformna raspodela

U teoriji verovatnoće i statistici, diskretna uniformna raspodela je simetrična raspodela verovatnoće pri čemu je za konačni broj vrednosti podjednako verovatno da budu uočene; svaka od -{n}- vrednosti ima jednaku verovatnoću 1/-{n}-. Drugim rečima „diskretna uniformna distribucija” je „poznati, konačni broj ishoda koji su podjednako verovatni da se dogode”.

Jednostavan primer diskretne uniformne distribucije je bacanje kocke. Moguće vrednosti su 1, 2, 3, 4, 5, 6, i svaki put kada se baci kocka verovatnoća datih rezultata je 1/6. Ako se bacaju dve kocke i dodaju njihove vrednosti, rezultirajuća raspodjela više nije uniformna, jer svi zbirovi nemaju jednaku verovatnoću.

Sama diskretna uniformna raspodela je inherentno neparametarska. Prikladno je, međutim, prikazati njene vrednosti generalno svim celim brojevima u intervalu [-{a}-,-{b}-], tako da -{a}- i -{b}- postaju glavni parametri distribucije (često se jednostavno razmatra interval [1, -{n}-] sa jednim parametrom -{n}-). Ovim konvencijama može se izraziti funkcija kumulativne distribucije (Шаблон:Jez-eng-lat) diskretne uniformne distribucije za bilo koje -{k ∈ [a,b]}-, kao

$F(k;a,b)=\frac{\lfloor k\rfloor -a+1}{b-a+1}$

Procena maksimuma

Uzorak -{k}- opažanja dobijen je iz uniformne distribucije celih brojeva $1,2,\dots ,N$, s ciljem procene nepoznatog maksimuma -{N}-. Ovaj problem je opšte poznat kao nemački tenkovski problem, po primeni maksimalne procene na procenu nemačke proizvodnje tenkova tokom Drugog svetskog rata.

Nepristrasni procenjivač uniformne minimalne varijanse (Шаблон:Jez-eng-lat) za maksimum je dat sa

$\hat{N}=\frac{k+1}{k}m-1=m+\frac{m}{k}-1$

gde je -{m}- maksimum uzorka i -{k}- je veličina uzorka, uzorkovanog bez zamene. Ovo se može smatrati veoma jednostavnim slučajem procenе maksimalnog razmaka.

To ima varijansu od

$\frac{1}{k}\frac{(N-k)(N+1)}{(k+2)}\approx \frac{N^2}{k^2}\text{ za male uzorke }k\ll N$

tako da je standardna devijacija približno $\frac{N}{k}$, (populaciona) prosečna veličina razmaka između uzoraka; uporedite sa $\frac{m}{k}$ iznad.

Maksimum uzorka je procenjivač maksimalne verovatnoće za populacioni maksimum, mada, kao što je gore diskutovano, on je pristrasan.

Ako uzorci nisu numerisani, već su prepoznatljivi ili markirani, umesto toga se može proceniti veličina populacije metodom označavanja i ponovnog hvatanja.

Derivacija

Za neki ceo broj -{m}- takav da je -{k ≤ m ≤ N}-, verovatnoća da će maksimalni uzorak biti jednak -{m}- može se izračunati na sledeći način. Broj različitih grupa od -{k}- tenkova koji se mogu načiniti od ukupno -{N}- tenkova dat je binomnim koeficijentom $(\genfrac{}{}{0ex}{}{N}{k})$. Pošto se u ovom načinu brojanja permutacija tenkova broji samo jednom, mogu se urediti serijski brojevi i uzeti u obzir maksimum svakog uzorka. Da bi se izračunala verovatnoća, mora se odrediti broj poređanih uzoraka koji se može formirati sa poslednjim elementom jednakim -{m}-, i svi ostali -{k}--1 tenkovi manji ili jednaki -{m}--1. Broj uzoraka sa -{k}--1 tenkova koji se mogu napraviti od ukupno -{m}--1 tenkova dat je binomnim koeficijentom $(\genfrac{}{}{0ex}{}{m-1}{k-1})$, tako da je verovatnoća da je maksimum -{m}- jednaka $P(m)=(\genfrac{}{}{0ex}{}{m-1}{k-1})/(\genfrac{}{}{0ex}{}{N}{k})$.

Ako je dat ukupn broj -{N}- i veličina uzorka -{k}-, očekivana vrednost maksimuma uzorka je

$\begin{array}{cc}\mu =\mathrm{E}[m]& =\sum _{m=k}^{N}m\frac{(\genfrac{}{}{0ex}{}{m-1}{k-1})}{(\genfrac{}{}{0ex}{}{N}{k})}\\ & =\frac{1}{(k-1)!(\genfrac{}{}{0ex}{}{N}{k})}\sum _{m=k}^N\frac{m!}{(m-k)!}\\ & =\frac{k!}{(k-1)!(\genfrac{}{}{0ex}{}{N}{k})}\sum _{m=k}^N(\genfrac{}{}{0ex}{}{m}{k})\\ & =k\frac{(\genfrac{}{}{0ex}{}{N+1}{k+1})}{(\genfrac{}{}{0ex}{}{N}{k})}\\ & =\frac{k(N+1)}{k+1},\end{array}$

gde je korišten identitet hokejaškog štapa $\sum _{m=k}^N(\genfrac{}{}{0ex}{}{m}{k})=(\genfrac{}{}{0ex}{}{N+1}{k+1})$.

Iz ove jednačine, nepoznati kvantitet -{N}- se može izraziti u obliku očekivanja i veličine uzorka kao

$\begin{array}{cc}N& =\mu (1+k^{-1})-1.\end{array}$

Iz linearnosti očekivanja se dobija da je

$\begin{array}{cc}\mu (1+k^{-1})-1& =\mathrm{E}[m(1+k^{-1})-1],\end{array}$

i stoga se nepristrasni procenjivač od -{N}- dobija zamenjivanjem očekivanja opservacijom,

$\begin{array}{cc}\hat{N}& =m(1+k^{-1})-1.\end{array}$

Osim što je nepristrasan ovaj procenjivač takođe dostiže minimalnu varijansu. Da bi se to pokazalo, prvo treba napomenuti da je maksimalan uzorak dovoljna statistička za maksimum populacije, jer je verovatnoća -{P(m;N)}- izražena samo kao funkcija od -{m}-. Zatim se mora pokazati da su statistike za -{m}- takođe kompletne statistike, i posebna vrsta dovoljne statistike. Zatim Lehman-Šefeova teorema implicira da je $\hat{N}$ minimalna procena nepristrasne varijanse od -{N}-.

Varijansa procenjivača izračunava se iz varijance maksimuma uzorka

$\begin{array}{cc}\mathrm{V}\mathrm{a}\mathrm{r}[\hat{N}]& =\frac{(k+1{)}^2}{k^2}\mathrm{V}\mathrm{a}\mathrm{r}[m].\end{array}$

Varijanse maksimuma se zatim izračunava iz očekivanih vrednosti od $m$ i $m^2$. Izračunavanje očekivane vrednosti od $m^2$ je,

$\begin{array}{cc}\mathrm{E}[m^2]& =\sum _{m=k}^N{m}^2\frac{(\genfrac{}{}{0ex}{}{m-1}{k-1})}{(\genfrac{}{}{0ex}{}{N}{k})}\\ & =\frac{1}{(k-1)!(\genfrac{}{}{0ex}{}{N}{k})}\sum _{m=k}^{N}m\frac{m!}{(m-k)!}\\ & =\frac{1}{(k-1)!(\genfrac{}{}{0ex}{}{N}{k})}\sum _{m=k}^N(m+1-1)\frac{m!}{(m-k)!}\\ & =\frac{1}{(k-1)!(\genfrac{}{}{0ex}{}{N}{k})}\sum _{m=k}^N\frac{(m+1)!}{(m-k)!}-\frac{1}{(k-1)!(\genfrac{}{}{0ex}{}{N}{k})}\sum _{m=k}^N\frac{m!}{(m-k)!}\end{array}$

gde je drugi član očekivana vrednost od $m$. Prvi član se može izraziti tako da zavisi od -{k}- i -{N}-,

$\begin{array}{cc}\frac{1}{(k-1)!(\genfrac{}{}{0ex}{}{N}{k})}\sum _{m=k}^N\frac{(m+1)!}{(m-k)!}& =\frac{(k+1)!}{(k-1)!(\genfrac{}{}{0ex}{}{N}{k})}\sum _{m=k}^N(\genfrac{}{}{0ex}{}{m+1}{k+1})\\ & =\frac{k(k+1)}{(\genfrac{}{}{0ex}{}{N}{k})}\sum _{n=k+1}^{N+1}(\genfrac{}{}{0ex}{}{n}{k+1})\\ & =\frac{k(k+1)}{(\genfrac{}{}{0ex}{}{N}{k})}(\genfrac{}{}{0ex}{}{N+2}{k+2})\\ & =\frac{k(N+2)(N+1)}{(k+2)}\end{array}$

pri čemu je napravljena zamena $n=m+1$ i korišten je identitet hokejaškog štapa. Zamenjujući ovaj rezultat i očekivanje od $m$ u jednačini od $E[m^2]$,

$\begin{array}{cc}\mathrm{E}[m^2]& =\frac{k(N+2)(N+1)}{(k+2)}-\frac{k(N+1)}{k+1}\\ & =k(N+1)(\frac{N+2}{k+2}-\frac{1}{k+1})\\ & =\frac{k(N+1)(kN+k+N)}{(k+1)(k+2)}\end{array}$

Variansa od $m$ se zatim dobija,

$\begin{array}{cc}\mathrm{V}\mathrm{a}\mathrm{r}[m]& =\mathrm{E}[m^2]-\mathrm{E}[m{]}^2\\ & =\frac{k(N+1)}{(k+1)}(\frac{kN+k+N}{k+2}-\frac{k(N+1)}{k+1})\\ & =\frac{k(N+1)}{(k+1)}\frac{(N-k)}{(k+2)(k+1)}\\ & =\frac{k(N+1)(N-k)}{(k+1{)}^2(k+2)}\end{array}$

Konačno varijansa procenjivača $\hat{N}$ se može izračunati,

$\begin{array}{cc}\mathrm{V}\mathrm{a}\mathrm{r}[\hat{N}]& =\frac{(k+1{)}^2}{k^2}\mathrm{V}\mathrm{a}\mathrm{r}[m]\\ & =\frac{(k+1{)}^2}{k^2}\frac{k(N+1)(N-k)}{(k+1{)}^2(k+2)}\\ & =\frac{(N+1)(N-k)}{k(k+2)}.\end{array}$

Vidi još

• Dirakova delta funkcija
• Uniformna distribucija (kontinuirana)

Literatura

• Шаблон:Cite journal
• Шаблон:Cite journal
• Шаблон:Cite journal
• Шаблон:Cite journal
• Шаблон:Cite journal
• Шаблон:Cite journal
• Шаблон:Cite conferenceШаблон:Мртва веза
• Шаблон:Cite journal
• Шаблон:Cite journal
• Шаблон:Cite journal
• Шаблон:Cite journal
• Шаблон:Cite book

Spoljašnje veze

• Diskretna slučajna varijabla Шаблон:Wayback