Geordende steekproef

In de statistiek vormen de naar grootte gerangschikte elementen van een steekproef $X_1,\dots ,X_n$ van continue stochastische variabelen, die onderling onafhankelijk zijn, maar niet noodzakelijk gelijkverdeeld, de geordende steekproef, meestal genoteerd als

$X_{(1)}\le \dots \le X_{(n)}$.

Met $X_{(k)}$ wordt het steekproefelement aangeduid met het rangnummer $k$. De notatie $X_{k:n}$ wordt ook gebruikt, waaraan tevens de steekproefomvang is te zien. Als er geen knopen zijn, geldt dus:

$X_{(k)}=X_i⟺r_i=k$

Als $x_1,\dots ,x_n$ de uitkomst van de steekproef is, worden de geordende resultaten genoteerd als:

$x_{(1)}\le \dots \le x_{(n)}$

De elementen in de geordende steekproef zijn stochastisch afhankelijk en elk van de elementen is een steekproeffunctie van de oorspronkelijke steekproef. In het bijzonder is

$X_{(1)}=\min (X_1,\dots ,X_n)$

en

$X_{(n)}=\max (X_1,\dots ,X_n)$

Verdeling

In de meeste gevallen worden gelijkverdeelde variabelen beschouwd, die dus een aselecte steekproef vormen.

In het algemene geval is de verdeling gecompliceerder en wordt deze gegeven door de stelling van Bapat–Beg, die in 1989 gepubliceerd werd door Bapat en Beg. De auteurs gaven geen bewijs, maar in 1994 gaf Hande een eenvoudig bewijs van de stelling.

Aselecte steekproef

Voor een aselecte steekproef van $X,$ dus voor onderling onafhankelijke en gelijkverdeelde $X$-en, is de simultane verdeling voor $y_1\le y_2\le \dots \le y_n$ gegeven door de kansdichtheid:

$f_{X_{(1)},\dots ,X_{(n)}}(y_1,\dots ,y_n)=n!\prod _k{f}_X(y_k)$

De verdelingsfunctie van $X_{(k)}$ wordt gegeven door:

$F_{X_{(k)}}(x)=P(X_{(k)}\le x)=$

$=\sum _{m=k}^{n}P(X_{(1)}\le x,\dots ,X_{(m)}\le x,X_{(m+1)}>x,\dots ,X_{(n)}>x)$

$=\sum _{m=k}^n(\genfrac{}{}{0ex}{}{n}{m})P(X_1\le x,\dots ,X_m\le x,X_{m+1}>x,\dots ,X_n>x)$

$=\sum _{m=k}^n(\genfrac{}{}{0ex}{}{n}{m}){(F_X(x))}^m{(1-F_X(x))}^{n-m}$,

want elk van de $(\genfrac{}{}{0ex}{}{n}{m})$ gebeurtenissen

$\{X_{i_1}\le x,\dots ,X_{i_m}\le x,X_{i_{m+1}}>x,\dots ,X_{i_n}>x\}$

heeft dezelfde kans als

$\{X_1\le x,\dots ,X_m\le x,X_{m+1}>x,\dots ,X_n>x\}$

De dichtheid van $X_{(k)}$ is:

$f_{X_{(k)}}(x)=n(\genfrac{}{}{0ex}{}{n-1}{k-1}){(F_X(x))}^{k-1}{f}_X(x)(1-(F_X(x)){)}^{n-k}$

Immers:

$f_{X_{(k)}}(x)\mathrm{d}x\approx P(X_{(k)}\in [x,x+\mathrm{d}x))=$

$=n(\genfrac{}{}{0ex}{}{n-1}{k-1})P(X_1\le x,\dots ,X_{k-1},X_k\in [x,x+\mathrm{d}x)),X_{k+1}>x,\dots ,X_n>x)$

$=n(\genfrac{}{}{0ex}{}{n-1}{k-1})F_X(x)\cdot \dots \cdot F_X(x)f_X(x)\mathrm{d}x(1-F_X(x))\cdot \dots \cdot (1-F_X(x))$

Minimum en maximum

Voor het minimum geldt dus:

$F_{X_{(1)}}(x)=1-{(1-F_X(x))}^n$ en $f_{X_{(1)}}(x)=n{(1-F_X(x))}^{n-1}{f}_X(x)$,

en voor het maximum:

$F_{X_{(n)}}(x)={(F_X(x))}^n$ en $f_{X_{(n)}}(x)=n{(F_X(x))}^{n-1}{f}_X(x)$

Uniforme verdeling op (0,1)

Voor een aselecte steekproef $U_1,\dots ,U_n$ uit de uniforme verdeling op het interval (0,1) is:

$f_{U_{(k)}}(u)=n(\genfrac{}{}{0ex}{}{n-1}{m-1})u^{k-1}(1-u{)}^{n-k},(u\in (0,1))$

Dit betekent dat $U_{(k)}$ een bètaverdeling heeft met parameters $k$ en $n+1-k$:

$U_{(k)}\sim B(k,n+1-k)$

Stelling van Bapat-Beg

De stochastische variabelen $X_1,\dots ,X_n$ zijn onderling onafhankelijk en hebben verdelingsfuncties $F_i=F_{X_i},i=1,\dots ,n$. De simultane verdelingsfunctie van de elementen $X_{(r_1)}\le X_{(r_2)}\le \dots \le X_{(n_k)}$ van de geordende steekproef wordt voor $x_1<x_2<\dots <x_k$ gegeven door:

$F_{X_{(r_1)},X_{(r_2)},\dots ,X_{(r_k)}}(x_1,x_2,\dots ,x_k)=P(X_{(r_1)}\le x_1,X_{(r_2)}\le x_2,\dots ,X_{(r_k)}\le x_k)=$

$=\sum _{n_k=r_k}^n\dots \sum _{n_2=r_2}^{n_3}\sum _{n_1=r_1}^{n_2}\frac{P_{n_1,\dots ,n_k}(x_1,\dots ,x_k)}{n_1!(n_2-n_1)!\dots (n-n_k)!}$,

waarin

$P_{n_1,\dots ,n_k}(x_1,\dots ,x_k)=$

$=\mathrm{per}\left(\left[\begin{array}{ccccc}{F}_{11}\dots F_{11}& F_{12}-F_{11}\dots F_{12}-F_{11}& \dots & F_{1k}-F_{1,k-1}\dots F_{1k}-F_{1,k-1}& 1-F_{1k}\dots 1-F_{1k}\\ F_{21}\dots F_{21}& F_{22}-F_{21}\dots F_{22}-F_{21}& \dots & F_{2k}-F_{2,k-1}\dots F_{2k}-F_{2,k-1}& 1-F_{2k}\dots 1-F_{2k}\\ ⋮& ⋮& & ⋮\\ \underset{n_1}{\underbrace{F_{n1}\dots F_{n1}}}& \underset{n_2-n_1}{\underbrace{F_{n2}-F_{n1}\dots F_{n2}-F_{n1}}}& \dots & \underset{n_k-n_{k-1}}{\underbrace{F_{nk}-F_{n,k-1}\dots F_{nk}-F_{n,k-1}}}& \underset{n-n_k}{\underbrace{1-F_{nk}\dots 1-F_{nk}}}\end{array}\right]\right)$

de permanent is van de genoemde matrix met $F_{rm}=F_r(x_m)$ en onder de accolades de getallen staan, die het aantal kolommen aangeven.

Toepassing

De geordende steekproef en de rangnummers spelen een belangrijke rol in de verdelingsvrije statistiek.

Als de verdelingsfunctie van de verdeling waaruit de steekproef getrokken is, bekend is, kan de geordende steekproef herleid worden tot de geordende steekproef uit de uniforme verdeling, en de eigenschappen aan de hand hiervan bestudeerd worden.

Literatuur

• Bapat, R. B.; Beg, M. I. (1989). "Order Statistics for Nonidentically Distributed Variables and Permanents". Sankhyā: The Indian Journal of Statistics, Series A (1961-2002) 51 (1): 79–93. JSTOR 25050725. MR 1065561.
• David, H. A. Order Statistics, 2nd ed. New York: Wiley, 1981.
• Gibbons, J. D. and Chakraborti, S. (Eds.). Nonparametric Statistic Inference, 3rd ed. exp. rev. New York: Dekker, 1992.
• Hande, Sayaji (1994). "A Note on Order Statistics for Nondentically Distributed Variables". Sankhyā: The Indian Journal of Statistics, Series A (1961-2002) 56 (2): 365–368. JSTOR 25050995. MR 1664921.
• Hogg, R. V. and Craig, A. T. Introduction to Mathematical Statistics, 3rd ed. New York: Macmillan, 1970.
• Rose, C. and Smith, M. D. "Order Statistics." §9.4 in Mathematical Statistics with Mathematica. New York: Springer-Verlag, pp. 311-322, 2002.

Websites

• MathWorld. Order Statistic.
• M. Güngör: "On Joint Distributions of Order Statistics from innid Variables" Sjabloon:Pdf