順序統計量(order statistics)とは、確率変数 $X_1,\dots,X_n$ の値を小さい順に並べ直した $X_{(1)}\le\dots\le X_{(n)}$ のことで、最小値・最大値・標本中央値・標本分位点はその例である。$X_i$ が独立で同じ分布関数 $F$ をもつとき、$X_{(k)}\le x$ は「$x$ 以下の値が $k$ 個以上」と同じ事象なので、$X_{(k)}$ の分布関数は二項分布の裾の確率で書ける。この裾はベータ積分に等しく、一様分布の $U_{(k)}$ はベータ分布に従って期待値 $k/(n+1)$ をもち、密度 $f$ をもつ分布では $X_{(k)}$ の密度が $F$ と $f$ で書ける。$\tau$ 分位点がただ 1 つで $k/n\to\tau$ なら、$X_{(k)}$ はその分位点に確率収束する。
$n$ 個の確率変数 $X_1,\dots,X_n$ の値を小さい順に並べ直したものを $X_{(1)}\le X_{(2)}\le\dots\le X_{(n)}$ と書き、$X_{(k)}$ を 第 $k$ 順序統計量 という。最小値 $X_{(1)}$、最大値 $X_{(n)}$、標本中央値、標本分位点(分位点)はどれも順序統計量であり、四分位数(四分位数と箱ひげ図)は順序統計量またはその 2 つの平均である。
順序統計量の分布を調べる鍵は、「$X_{(k)}\le x$」と「$x$ 以下の値が $k$ 個以上ある」が同じ事象だという単純な観察である。$X_i$ が独立で同じ分布に従えば、$x$ 以下の値の個数は二項分布に従うので、$X_{(k)}$ の分布関数は二項分布の裾の確率で書ける(主定理 1)。その裾の確率はベータ関数の不完全積分に等しく、ここから一様分布の順序統計量がベータ分布に従うこと、密度をもつ分布の順序統計量の密度が得られる。さらに、$\tau$ 分位点がただ 1 つのとき、$k/n$ が $\tau$ に近づくなら $X_{(k)}$ がその分位点に確率収束すること(主定理 2)を示す。
この記事を通して、$X_1,\dots,X_n$ は確率空間 $(\Omega,\mathcal F,P)$ の上の実数値の確率変数とし、$x$ 以下の値の個数を
$$
N(x):=\#\{i\mid X_i\le x\}=\sum_{i=1}^n\mathbf 1_{\{X_i\le x\}}
$$
と書く。
各 $\omega\in\Omega$ について、$n$ 個の実数 $X_1(\omega),\dots,X_n(\omega)$ を(同じ値は重複して)小さい順に並べたものを $X_{(1)}(\omega)\le\dots\le X_{(n)}(\omega)$ とする。$X_{(k)}$ を 第 $k$ 順序統計量($k$-th order statistic)、組 $(X_{(1)},\dots,X_{(n)})$ を 順序統計量 という。$X_{(n)}-X_{(1)}$ を 範囲、$n=2m+1$ のときの $X_{(m+1)}$ を 標本中央値 という。
並べ方の定義から、すべての実数 $x$ と $1\le k\le n$ について
$$
\{X_{(k)}\le x\}=\{N(x)\ge k\}
$$
である($k$ 番目に小さい値が $x$ 以下であることと、$x$ 以下の値が $k$ 個以上あることは同じ)。右辺は事象なので、$X_{(k)}$ は確率変数である。また
$$
X_{(k)}=\min_{\lvert S\rvert=k}\ \max_{i\in S}X_i\qquad(S\text{ は }\{1,\dots,n\}\text{ の }k\text{ 元部分集合})
$$
とも書ける。どの $k$ 元集合 $S$ についても $\max_{i\in S}X_i$ は $k$ 番目に小さい値以上であり、小さい方から $k$ 個の添字の集合でちょうどその値になるからである。この表示から、並べ替えの写像 $\mathbb R^n\to\mathbb R^n$ は連続である。
データ $x_1,\dots,x_n$ の経験分布の分位点関数の値は $Q_n(\tau)=x_{(\lceil n\tau\rceil)}$ である(分位点 の例「標本分位点」)。標本分位点は順序統計量そのものである。
$X_1,\dots,X_n$ が互いに独立で、どれも分布関数 $F$ をもつとする。$1\le k\le n$ とすると、すべての実数 $x$ について
$$
P(X_{(k)}\le x)=\sum_{j=k}^n\binom nj F(x)^j\bigl(1-F(x)\bigr)^{n-j}
$$
である。とくに $P(X_{(n)}\le x)=F(x)^n$、$P(X_{(1)}>x)=\bigl(1-F(x)\bigr)^n$ である。
$x$ を固定する。事象 $\{X_1\le x\},\dots,\{X_n\le x\}$ は互いに独立で、どれも確率 $F(x)$ をもつ。よって $N(x)$ は二項分布 $B(n,F(x))$ に従う。実際、$0\le j\le n$ について、事象 $\{N(x)=j\}$ は、$j$ 元集合 $A\subset\{1,\dots,n\}$ ごとの事象 $\bigcap_{i\in A}\{X_i\le x\}\cap\bigcap_{i\notin A}\{X_i>x\}$($\binom nj$ 個あり、互いに交わらない)の和であり、$X_i$ の独立性から各事象の確率は $F(x)^j\bigl(1-F(x)\bigr)^{n-j}$ である(離散確率空間での同じ計算が 二項分布 の命題「成功の回数の分布」である)。定義の後に述べた $\{X_{(k)}\le x\}=\{N(x)\ge k\}$ により
$$
P(X_{(k)}\le x)=P(N(x)\ge k)=\sum_{j=k}^n\binom njF(x)^j\bigl(1-F(x)\bigr)^{n-j}
$$
である。$k=n$ なら和は $j=n$ の項 $F(x)^n$ だけであり、$k=1$ なら $P(X_{(1)}>x)=P(N(x)=0)=(1-F(x))^n$ である。$\square$
$F$ の連続性も密度も使っていないので、この定理は離散分布にもそのまま使える。
$X_1,\dots,X_n$ が互いに独立で、どれも率 $\lambda>0$ の指数分布($F(x)=1-e^{-\lambda x}$、$x\ge0$。一様分布と指数分布)に従うとする。
二項分布の裾の確率は、1 つの積分で書ける。
$1\le k\le n$ とし、
$$
\beta_{k}(t):=\frac{n!}{(k-1)!\,(n-k)!}\,t^{k-1}(1-t)^{n-k}
$$
とおく。$0\le u\le1$ について
$$
\sum_{j=k}^n\binom nju^j(1-u)^{n-j}=\int_0^u\beta_k(t)\,dt
$$
である。
両辺は $u$ の多項式で、$u=0$ ではどちらも $0$ である(左辺は $j\ge k\ge1$ の項だけなので)。導関数が等しいことを示せばよい。$0\le j\le n-1$ について $a_j(u):=n\binom{n-1}ju^j(1-u)^{n-1-j}$、$a_n(u):=0$ とおく。$j\binom nj=n\binom{n-1}{j-1}$、$(n-j)\binom nj=n\binom{n-1}j$ を使うと、$k\le j\le n$ について
$$
\frac d{du}\Bigl(\binom nju^j(1-u)^{n-j}\Bigr)=j\binom nju^{j-1}(1-u)^{n-j}-(n-j)\binom nju^j(1-u)^{n-j-1}=a_{j-1}(u)-a_j(u)
$$
である($j=n$ では第 2 項は $0$ で、$a_n=0$ と合う)。$j=k,\dots,n$ について足すと、隣り合う項が打ち消し合って左辺の導関数は $a_{k-1}(u)-a_n(u)=n\binom{n-1}{k-1}u^{k-1}(1-u)^{n-k}=\beta_k(u)$ になる。これは右辺の導関数である。$\square$
$U_1,\dots,U_n$ が互いに独立で、どれも $(0,1)$ 上の一様分布に従うとする。$1\le k\le n$ について、$U_{(k)}$ は $(0,1)$ の上で $\beta_k$、その外で $0$ となる密度をもつ。すなわち $U_{(k)}$ はパラメータ $(k,\,n-k+1)$ のベータ分布に従う。さらに
$$
E[U_{(k)}]=\frac k{n+1},\qquad V[U_{(k)}]=\frac{k(n-k+1)}{(n+1)^2(n+2)}
$$
である。
$0\le u\le1$ では $F(u)=u$ なので、thm-ord-cdf と lem-ord-beta により $P(U_{(k)}\le u)=\int_0^u\beta_k(t)\,dt$ である。$u< 0$ では両辺 $0$、$u>1$ では $P(U_{(k)}\le u)=1=\int_0^1\beta_k$(lem-ord-beta で $u=1$)である。よって 確率密度関数 の定理「分布関数による特徴づけ」により、上の関数は $U_{(k)}$ の密度である。
期待値は、確率密度関数 の定理「密度による期待値の公式」と Beta関数 の系「正の整数での値」$B(a,b)=\frac{(a-1)!(b-1)!}{(a+b-1)!}$ により
$$
E[U_{(k)}]=\frac{n!}{(k-1)!(n-k)!}B(k+1,n-k+1)=\frac{n!}{(k-1)!(n-k)!}\cdot\frac{k!\,(n-k)!}{(n+1)!}=\frac k{n+1}
$$
であり、同様に $E[U_{(k)}^2]=\frac{n!}{(k-1)!(n-k)!}B(k+2,n-k+1)=\frac{k(k+1)}{(n+1)(n+2)}$ である。よって $V[U_{(k)}]=\frac{k(k+1)}{(n+1)(n+2)}-\frac{k^2}{(n+1)^2}=\frac{k(n-k+1)}{(n+1)^2(n+2)}$ である。$\square$
$E[U_{(k)}]=\frac k{n+1}$ は、$n$ 個の点が区間 $(0,1)$ を $n+1$ 個の部分に分け、どの部分の長さも平均 $\frac1{n+1}$ であることと合っている。$k=n$ なら $E[\max_iU_i]=\frac n{n+1}$ である。
$F$ を分布関数、$Q(\tau)=\inf\{x\mid F(x)\ge\tau\}$ をその分位点関数とし、$U_1,\dots,U_n$ を thm-ord-uniform のとおりとする。$X_1,\dots,X_n$ が互いに独立でどれも分布関数 $F$ をもつなら、組 $(X_{(1)},\dots,X_{(n)})$ は組 $\bigl(Q(U_{(1)}),\dots,Q(U_{(n)})\bigr)$ と同じ同時分布に従う。
$X'_i:=Q(U_i)$ とおく。分位点 の定理「分位点変換と確率積分変換」の 1 により各 $X'_i$ の分布関数は $F$ である。Borel 集合 $B_1,\dots,B_n$ について $\{X'_i\in B_i\}=\{U_i\in Q^{-1}(B_i)\}$ であり($Q$ は単調なので Borel 可測)、$U_i$ の独立性から $P(X'_1\in B_1,\dots,X'_n\in B_n)=\prod_iP(X'_i\in B_i)$ である。よって $(X'_1,\dots,X'_n)$ と $(X_1,\dots,X_n)$ は、直積集合の上で同じ確率を与え、したがって同じ同時分布に従う(直積集合は共通部分について閉じ、$\mathbb R^n$ の Borel 集合族を生成するので、Dynkinのπ-λ定理 の定理「π系の上で一致する測度の一意性」による)。並べ替えの写像は連続なので、それぞれの順序統計量も同じ同時分布に従う。
最後に、$Q$ は単調非減少なので $Q(U_{(1)})\le\dots\le Q(U_{(n)})$ であり、これは値 $Q(U_1),\dots,Q(U_n)$ を並べ替えたものである。よって $X'_{(k)}=Q(U_{(k)})$ である。$\square$
この定理により、たとえば $E[X_{(k)}]=E[Q(U_{(k)})]=\int_0^1Q(u)\beta_k(u)\,du$(右辺が絶対収束するとき)となり、一般の分布の順序統計量の計算が一様分布の場合に帰着する。
$X_1,\dots,X_n$ が互いに独立で、どれも密度 $f$ と分布関数 $F$ をもつとする。$1\le k\le n$ について、$X_{(k)}$ は密度
$$
f_{(k)}(x)=\frac{n!}{(k-1)!\,(n-k)!}\,F(x)^{k-1}\bigl(1-F(x)\bigr)^{n-k}f(x)=\beta_k\bigl(F(x)\bigr)f(x)
$$
をもつ。
要点:thm-ord-cdf と lem-ord-beta から $P(X_{(k)}\le x)=\int_0^{F(x)}\beta_k(u)\,du$ である。これが $\int_{-\infty}^x\beta_k(F(s))f(s)\,ds$ に等しいことを、$F(X)$ が一様分布に従うことを使って示す。
$X:=X_1$ とする。$F$ は密度をもつので連続であり、分位点 の定理「分位点変換と確率積分変換」の 2 により $F(X)$ は $[0,1]$ 上の一様分布に従う。$x$ を固定し、ほとんど確実に $\mathbf 1_{\{X\le x\}}=\mathbf 1_{\{F(X)\le F(x)\}}$ であることを示す。$X\le x$ なら $F(X)\le F(x)$ である。逆向きが崩れるのは事象 $E:=\{X>x,\ F(X)=F(x)\}$ の上だけである。$b:=\sup\{y\ge x\mid F(y)=F(x)\}$ とおく。$b=\infty$ なら $F(x)=1$ で $P(X>x)=0$ である。$b< \infty$ なら $F$ の連続性から $F(b)=F(x)$ で、$E\subset\{x< X\le b\}$ の確率は $F(b)-F(x)=0$ である。
よって、確率密度関数 の定理「密度による期待値の公式」を $X$ と $F(X)$(密度 $\mathbf 1_{[0,1]}$)に使って$$\int_{-\infty}^x\beta_k(F(s))f(s)\,ds=E\bigl[\beta_k(F(X))\mathbf 1_{\{X\le x\}}\bigr]=E\bigl[\beta_k(F(X))\mathbf 1_{\{F(X)\le F(x)\}}\bigr]=\int_0^{F(x)}\beta_k(u)\,du$$である($\beta_k$ は $[0,1]$ 上で有界なので期待値は有限)。これは $P(X_{(k)}\le x)$ に等しいので、確率密度関数 の定理「分布関数による特徴づけ」により $f_{(k)}$ は $X_{(k)}$ の密度である。$\square$
密度の式は次のように読める。$n$ 個のうち 1 個が $x$ のごく近くに落ち($n$ 通り、確率 $\approx f(x)\,dx$)、残りのうち $k-1$ 個が $x$ より左に、$n-k$ 個が右に落ちる($\binom{n-1}{k-1}$ 通り)。同じ考え方で、密度をもつ場合の全体の同時密度は、$x_1< \dots< x_n$ の上で $n!\,f(x_1)\cdots f(x_n)$、その外で $0$ である。
密度をもつので $P(X_i=X_j)=0$($i\ne j$)であり、ほとんど確実に値はすべて異なる。$B\subset\{x_1< \dots< x_n\}$ を Borel 集合とすると、$(X_{(1)},\dots,X_{(n)})\in B$ は、ちょうど 1 つの置換 $\sigma$ について $(X_{\sigma(1)},\dots,X_{\sigma(n)})\in B$ となることと同値である(値がすべて異なるとき)。各 $(X_{\sigma(1)},\dots,X_{\sigma(n)})$ は同時密度 $f(x_1)\cdots f(x_n)$ をもつので、$n!$ 個の置換について足して $P((X_{(1)},\dots,X_{(n)})\in B)=n!\int_Bf(x_1)\cdots f(x_n)\,dx$ を得る。
観測を増やしていくとき、$k/n$ が $\tau$ に近い順序統計量は $\tau$ 分位点に近づく。以下、$X_1,X_2,\dots$ は互いに独立でどれも分布関数 $F$ をもつとし、最初の $n$ 個の第 $k$ 順序統計量を $X_{(k)}^{[n]}$ と書く。
$0< \tau< 1$ とし、実数 $\xi$ がすべての $\varepsilon>0$ について
$$
F(\xi-\varepsilon)< \tau< F(\xi+\varepsilon)
$$
を満たすとする。整数の列 $1\le k_n\le n$ が $k_n/n\to\tau$ を満たすなら、$X_{(k_n)}^{[n]}$ は $\xi$ に確率収束する。すなわち、すべての $\varepsilon>0$ について $P\bigl(\lvert X_{(k_n)}^{[n]}-\xi\rvert>\varepsilon\bigr)\to0$ である。
条件は、$\xi$ が $\tau$ 分位点であり、しかも $\tau$ 分位点がほかにないことを表している(分位点 の定理「分位点の全体」で $Q(\tau)=Q^+(\tau)=\xi$ の場合)。
$\varepsilon>0$ を固定し、$N_n(x):=\#\{i\le n\mid X_i\le x\}$ とおく。$Z_i:=\mathbf 1_{\{X_i\le x\}}-F(x)$ とおくと、$E[Z_i]=0$、$E[Z_i^2]=F(x)\bigl(1-F(x)\bigr)\le\frac14$ であり、$i\ne j$ では独立性から $E[Z_iZ_j]=P(X_i\le x,X_j\le x)-F(x)^2=0$ である。よって $E\bigl[(N_n(x)/n-F(x))^2\bigr]=\frac1{n^2}\sum_{i\le n}E[Z_i^2]\le\frac1{4n}$ であり、$\mathbf 1_{\{\lvert z\rvert\ge\delta\}}\le z^2/\delta^2$ を $z=N_n(x)/n-F(x)$ に当てて期待値をとると(離散確率空間での同じ評価が 大数の法則 の系「相対頻度の法則」である)、すべての $\delta>0$ で
$$
P\Bigl(\Bigl\lvert\frac{N_n(x)}n-F(x)\Bigr\rvert\ge\delta\Bigr)\le\frac1{4n\delta^2}
$$
である。
下側。 $p:=F(\xi-\varepsilon)< \tau$、$\delta:=\frac{\tau-p}2>0$ とおく。$k_n/n\to\tau$ なので、十分大きい $n$ で $k_n/n\ge p+\delta$ である。このとき $\{X_{(k_n)}^{[n]}\le\xi-\varepsilon\}=\{N_n(\xi-\varepsilon)\ge k_n\}\subset\{N_n(\xi-\varepsilon)/n-p\ge\delta\}$ なので、確率は $\frac1{4n\delta^2}$ 以下で、$n\to\infty$ で $0$ に近づく。
上側。 $p':=F(\xi+\varepsilon)>\tau$、$\delta':=\frac{p'-\tau}2>0$ とおく。十分大きい $n$ で $k_n/n\le p'-\delta'$ であり、$\{X_{(k_n)}^{[n]}>\xi+\varepsilon\}=\{N_n(\xi+\varepsilon)< k_n\}\subset\{p'-N_n(\xi+\varepsilon)/n>\delta'\}$ の確率は $\frac1{4n\delta'^2}$ 以下である。
2 つを合わせて $P(\lvert X_{(k_n)}^{[n]}-\xi\rvert>\varepsilon)\to0$ を得る。$\square$
とくに、$\frac12$ 分位点がただ 1 つなら、標本中央値は母集団の中央値に確率収束する。この定理は平均や分散を仮定しないので、期待値をもたない Cauchy分布 でも、標本中央値は位置のパラメータに確率収束する(標本平均がもとの Cauchy 分布に従い続けることは Cauchy分布 の系「1 次結合と標本平均」にある)。
収束の速さについては、$(0,1)$ 上の一様分布から $2m+1$ 個をとったときの標本中央値 $M$ について、$2\sqrt{2m}\,(M-\frac12)$ の分布が $m\to\infty$ で標準正規分布に近づくことが知られている(Dur19 Example 3.2.6(pp. 117–118)。この記事では証明しない)。thm-ord-uniform によれば $V[M]=\frac{(m+1)^2}{(2m+2)^2(2m+3)}=\frac1{4(2m+3)}$ であり、$2\sqrt{2m}$ 倍した分散が $1$ に近づくことと合っている。
$U_1,\dots,U_n$($n\ge2$)が互いに独立で $(0,1)$ 上の一様分布に従うとき、thm-ord-uniform により $E[U_{(n)}-U_{(1)}]=\frac n{n+1}-\frac1{n+1}=\frac{n-1}{n+1}$ である。$n=3$ なら $\frac12$、$n=9$ なら $\frac45$ である。範囲は $n$ とともに $1$ に近づくが、外れ値 1 個で大きく動く量でもある。
$X_1,X_2,X_3$ が互いに独立で、どれも $0,1$ を確率 $\frac12$ ずつとるとする。$F(x)=\frac12$($0\le x< 1$)なので、thm-ord-cdf により標本中央値 $X_{(2)}$ について $P(X_{(2)}\le0)=\binom32\frac18+\binom33\frac18=\frac12$ である。$X_{(2)}$ は $0,1$ を確率 $\frac12$ ずつとる。値が重なる確率は正で($P(X_{(1)}=X_{(2)})\ge P(X_1=X_2=X_3)=\frac14$)、$X_{(k)}$ は密度をもたない。
| 外す条件 | 反例 | 成り立たなくなること |
|---|---|---|
| 互いに独立 | $U$ が $(0,1)$ 上の一様分布で $X_1=X_2=\dots=X_n=U$($n\ge2$) | $P(X_{(n)}\le x)=x\ne x^n$($0< x< 1$)で、thm-ord-cdf の式が成り立たない |
| 密度をもつ | ex-ord-coin | $X_{(k)}$ が正の確率で 1 点に値をとり、cor-ord-density の密度がない |
| $\tau$ 分位点がただ 1 つ | ex-ord-median-oscillate | 標本中央値がどの定数にも確率収束しない(thm-ord-consistency) |
| $k_n/n\to\tau$ | 一様分布で $k_n=1$(最小値) | $X_{(1)}^{[n]}$ は $0$ に確率収束し、$\tau$ 分位点 $\tau$ に近づかない |
$X_1,X_2,\dots$ が互いに独立で、どれも $0,1$ を確率 $\frac12$ ずつとるとし、$\tau=\frac12$ とする。$F$ は $[0,1)$ で $\frac12$ なので、すべての $\varepsilon\in(0,1)$ で $F(\xi-\varepsilon)< \frac12< F(\xi+\varepsilon)$ を満たす $\xi$ はない($\frac12$ 分位点の全体は $[0,1]$)。$n=2m+1$ のとき、標本中央値 $X_{(m+1)}^{[n]}$ が $0$ になるのは $0$ の個数が $m+1$ 以上のときで、$0$ と $1$ の役割を入れ替える対称性から、その確率はちょうど $\frac12$ である。よって標本中央値は、どの奇数 $n$ でも $0,1$ を確率 $\frac12$ ずつとり、どの定数にも確率収束しない。満たさない仮定:$\frac12$ 分位点がただ 1 つであること。
$U$ を $(0,1)$ 上の一様分布に従う確率変数とし、$X_1=U$、$X_2=1-U$ とする。$X_1,X_2$ はどちらも一様分布に従うが独立でない。$X_{(2)}=\max(U,1-U)$ は $[\frac12,1)$ の上の一様分布に従い、$P(X_{(2)}\le x)=2x-1$($\frac12\le x\le1$)である。独立な場合の $x^2$ とは異なり、たとえば $x=\frac34$ では $\frac12$ と $\frac9{16}$ である。満たさない仮定:$X_1,X_2$ が独立であること。
Mathpediaは寄付と、参考文献の書籍リンク(Amazonアソシエイト)の紹介料で運営されています。 支援について / 寄付する