順序統計量とは|定義・公式とアクチュアリー試験の関連問題 | acpass
確率・用語・変数変換・順序統計量
ひとことで言うと 「n個の標本を小さい順に並べ替えたときの、k番目の値」を表します。最小値も最大値も、この特別な場合(k=1またはk=n)です。
こんなデータが従う n個の部品のうち最初に壊れる部品の寿命(最小値) n人の試験の最高点(最大値) ある期間の最高気温・最低気温 保険金支払いの中で最も大きい1件の金額 n個のサンプルの中央値
標本そのものではなく、「並べ替えた後の順番」に注目した統計量です。最大値・最小値・中央値はいずれも順序統計量の特別な場合です。
n個の標本を並べ替えると順序統計量X(1)以下…以下X(n)になる(最小値と最大値を強調) X3 X1 X5 X2 X4 X(1) X(2) X(3) X(4) X(5) 観測順 大きさの順
拡大 5個の標本 X 1 X_1 X 1
数式で表すと F X ( n ) ( t ) = F ( t ) n , F X ( 1 ) ( t ) = 1 − ( 1 − F ( t ) ) n F_{X_{(n)}}(t)=F(t)^n,\ F_{X_{(1)}}(t)=1-(1-F(t))^n F X ( n ) ( t ) = F ( t ) n , F X ( 1 ) ( t ) = 1 − ( 1 − F ( t ) ) n 記法は X ( 1 ) ≤ X ( 2 ) ≤ ⋯ ≤ X ( n ) X_{(1)}\le X_{(2)}\le\dots\le X_{(n)} X ( 1 ) ≤ X ( 2 ) ≤ ⋯ ≤ X ( n )
試験に出る性質 最大値の分布
F X ( n ) ( x ) = [ F ( x ) ] n F_{X_{(n)}}(x)=[F(x)]^n F X ( n ) ( x ) = [ F ( x ) ] n
例で見る 独立に5個の U ( 0 , 1 ) U(0,1) U ( 0 , 1 ) をとったとき、最大値が0.9以下である確率は
P ( X ( 5 ) ≤ 0.9 ) = [ F ( 0.9 ) ] 5 = 0.9 5 ≈ 0.590 P(X_{(5)}\le 0.9)=[F(0.9)]^5=0.9^5\approx 0.590 P ( X ( 5 ) ≤ 0.9 ) = [ F ( 0.9 ) ] 5 = 0. 9 5 ≈ 0.590
つまずきポイント 最小値の式に最大値の式 [ F ( x ) ] n [F(x)]^n [ F ( x ) ] n をそのまま使ってしまう(最小値は 1 − [ 1 − F ( x ) ] n 1-[1-F(x)]^n 1 − [ 1 − F ( x ) ] n )定着クイズ Q1 X ( n ) X_{(n)} X ( n ) (最大値)の分布関数は?
A. 1 − [ 1 − F ( x ) ] n 1-[1-F(x)]^n 1 − [ 1 − F ( x ) ] n B. [ F ( x ) ] n [F(x)]^n [ F ( x ) ] n C. n F ( x ) nF(x) n F ( x ) Q2 X ( 1 ) X_{(1)} X ( 1 ) (最小値)の分布関数は?
A. [ F ( x ) ] n [F(x)]^n [ F ( x ) ] n B. F ( x ) / n F(x)/n F ( x ) / n C. 1 − [ 1 − F ( x ) ] n 1-[1-F(x)]^n 1 − [ 1 − F ( x ) ] n Q3 独立な 指数分布 E x p ( λ ) \mathrm{Exp}(\lambda) Exp ( λ ) をn個とったときの最小値が従う分布は?
A. 指数分布 E x p ( n λ ) \mathrm{Exp}(n\lambda) Exp ( nλ ) B. 指数分布 E x p ( λ / n ) \mathrm{Exp}(\lambda/n) Exp ( λ / n ) C. ガンマ分布 G a m m a ( n , λ ) \mathrm{Gamma}(n,\lambda) Gamma ( n , λ ) A. F ( t ) n F(t)^n F ( t ) n B. n F ( t ) nF(t) n F ( t ) C. 1 − ( 1 − F ( t ) ) n 1-(1-F(t))^n 1 − ( 1 − F ( t ) ) n Q5 U ( 0 , 1 ) U(0,1) U ( 0 , 1 ) から n n n 個取った最大値の密度は?
A. n t n − 1 n t^{n-1} n t n − 1 B. n ( 1 − t ) n − 1 n(1-t)^{n-1} n ( 1 − t ) n − 1 C. t n t^n t n Q6 U ( 0 , 1 ) U(0,1) U ( 0 , 1 ) の n = 5 n=5 n = 5 個の最大値の期待値は?
A. 1 1 1 B. 5 / 6 ≈ 0.833 5/6\approx0.833 5/6 ≈ 0.833 C. 1 / 2 1/2 1/2 〜
(観測順、上段)を大きさの順に並べ替えると、最小値
から最大値
までの順序統計量になる(下段、両端を強調)。
標本を小さい順に並べた X ( 1 ) ≤ ⋯ ≤ X ( n ) X_{(1)}\le\dots\le X_{(n)} X ( 1 ) ≤ ⋯ ≤ X ( n ) 。最大値は「全部が t t t 以下」、最小値は「全部が t t t より大きい」と読み替えると、独立性から積で書ける。
で、丸括弧つきの添字が「小さい方から数えて何番目か」を表します(元の観測順を表す
X 1 , … , X n X_1,\dots,X_n X 1 , … , X n とは別物)。
最大値 X ( n ) X_{(n)} X ( n ) がある値x以下になるのは、n個全部がx以下のときだけなので
P ( X ( n ) ≤ x ) = P ( 全てのX i ≤ x ) = [ F ( x ) ] n P(X_{(n)}\le x)=P(\text{全てのX}_i\le x)=[F(x)]^n P ( X ( n ) ≤ x ) = P ( 全ての X i ≤ x ) = [ F ( x ) ] n
と、独立性から各 F ( x ) F(x) F ( x ) をn回掛けるだけで求まります。
最小値はこの逆を考えます。最小値がxより大きいのは、n個全部がxより大きいときだけなので
P ( X ( 1 ) > x ) = [ 1 − F ( x ) ] n P(X_{(1)}>x)=[1-F(x)]^n P ( X ( 1 ) > x ) = [ 1 − F ( x ) ] n 、よって
P ( X ( 1 ) ≤ x ) = 1 − [ 1 − F ( x ) ] n P(X_{(1)}\le x)=1-[1-F(x)]^n P ( X ( 1 ) ≤ x ) = 1 − [ 1 − F ( x ) ] n
です。「以下」と「より大きい」の向きを反転させる点が、最大値の式との一番の違いです。
最大値 X ( n ) = max ( X 1 , … , X n ) X_{(n)}=\max(X_1,\dots,X_n) X ( n ) = max ( X 1 , … , X n ) の分布は、順序統計量の一般論で F ( n ) ( t ) = [ F ( t ) ] n F_{(n)}(t)=[F(t)]^n F ( n ) ( t ) = [ F ( t ) ] n と与えられます。理由は明快で、『最大値が t t t 以下』 ⟺ \iff ⟺ 『全員が t t t 以下』であり、独立同分布なら各々が t t t 以下になる確率 F ( t ) F(t) F ( t ) の積で F ( t ) n F(t)^n F ( t ) n になるからです。ここではこの分布関数を出発点に、密度の導出と一様分布での具体計算まで踏み込みます。
密度は分布関数を t t t で微分して得ます。 F ( n ) ( t ) = F ( t ) n F_{(n)}(t)=F(t)^n F ( n ) ( t ) = F ( t ) n を合成関数の微分で t t t について微分すると、
f ( n ) ( t ) = d d t [ F ( t ) ] n = n [ F ( t ) ] n − 1 f ( t ) f_{(n)}(t)=\dfrac{d}{dt}[F(t)]^n=n[F(t)]^{n-1}f(t) f ( n ) ( t ) = d t d [ F ( t ) ] n = n [ F ( t ) ] n − 1 f ( t )
となります(外側の微分 n [ F ( t ) ] n − 1 n[F(t)]^{n-1} n [ F ( t ) ] n − 1 に、内側の F ′ ( t ) = f ( t ) F'(t)=f(t) F ′ ( t ) = f ( t ) が掛かる)。直感的には、最大値がちょうど t t t 付近にあるのは『誰か1人が t t t 付近にいて (f ( t ) f(t) f ( t ) ) 、残り n − 1 n-1 n − 1 人は全員 t t t 以下にいる ([ F ( t ) ] n − 1 [F(t)]^{n-1} [ F ( t ) ] n − 1 ) 』状況で、その人の選び方が n n n 通りある、と読めます。
一様分布 U ( 0 , 1 ) U(0,1) U ( 0 , 1 ) で具体化しましょう。 F ( t ) = t , f ( t ) = 1 F(t)=t,\ f(t)=1 F ( t ) = t , f ( t ) = 1 (0 < t < 1 0<t<1 0 < t < 1 ) なので、最大値の密度は f ( n ) ( t ) = n t n − 1 f_{(n)}(t)=n\,t^{n-1} f ( n ) ( t ) = n t n − 1 です。これは ベータ分布 B e t a ( n , 1 ) \mathrm{Beta}(n,1) Beta ( n , 1 ) 分布の密度に一致します。期待値は
E [ X ( n ) ] = ∫ 0 1 t ⋅ n t n − 1 d t = n ∫ 0 1 t n d t = n n + 1 E[X_{(n)}]=\int_0^1 t\cdot n t^{n-1}\,dt=n\int_0^1 t^{n}\,dt=\dfrac{n}{n+1} E [ X ( n ) ] = ∫ 0 1 t ⋅ n t n − 1 d t = n ∫ 0 1 t n d t = n + 1 n
となり、有名な結果 E = n / ( n + 1 ) E=n/(n+1) E = n / ( n + 1 ) が出ます。 n n n が大きいほど n / ( n + 1 ) n/(n+1) n / ( n + 1 ) は1に近づき、最大値が支持の上限1へ張り付いていくことが式からも読み取れます。保険の文脈では、 n n n 件の独立な保険金請求のうち最大の請求額の振る舞いがこれにあたり、契約件数 n n n が増えるほど『最大級の損害』が分布の上限近くまで伸びていく——再保険や巨大リスクの評価で効いてくる直感です。分散は V [ X ( n ) ] = n ( n + 1 ) 2 ( n + 2 ) V[X_{(n)}]=\dfrac{n}{(n+1)^2(n+2)} V [ X ( n ) ] = ( n + 1 ) 2 ( n + 2 ) n で、 n n n が増えると上限付近に密集して分散は小さくなります。
最小値の分布
F X ( 1 ) ( x ) = 1 − [ 1 − F ( x ) ] n F_{X_{(1)}}(x)=1-[1-F(x)]^n F X ( 1 ) ( x ) = 1 − [ 1 − F ( x ) ] n
一般のk番目の密度
f ( k ) ( x ) = n ! ( k − 1 ) ! ( n − k ) ! [ F ( x ) ] k − 1 [ 1 − F ( x ) ] n − k f ( x ) f_{(k)}(x)=\dfrac{n!}{(k-1)!(n-k)!}[F(x)]^{k-1}[1-F(x)]^{n-k}f(x) f ( k ) ( x ) = ( k − 1 )! ( n − k )! n ! [ F ( x ) ] k − 1 [ 1 − F ( x ) ] n − k f ( x )
指数分布の最小値の特例
独立な 指数分布 E x p ( λ ) \mathrm{Exp}(\lambda) Exp ( λ ) をn個とったときの最小値は 指数分布 E x p ( n λ ) \mathrm{Exp}(n\lambda) Exp ( nλ ) (指数分布の最小値の競合と同じ事実)。
範囲(レンジ)
R = X ( n ) − X ( 1 ) R=X_{(n)}-X_{(1)} R = X ( n ) − X ( 1 ) (最大値と最小値の差)。
ベータ分布との対応(一様分布)
一様分布 U ( 0 , 1 ) U(0,1) U ( 0 , 1 ) からのn個の標本のk番目の順序統計量 X ( k ) X_{(k)} X ( k ) は、ベータ分布 B e t a ( k , n − k + 1 ) \mathrm{Beta}(k,n-k+1) Beta ( k , n − k + 1 ) に従う。期待値 E [ X ( k ) ] = k n + 1 E[X_{(k)}]=\dfrac{k}{n+1} E [ X ( k ) ] = n + 1 k 、分散 V [ X ( k ) ] = k ( n − k + 1 ) ( n + 1 ) 2 ( n + 2 ) V[X_{(k)}]=\dfrac{k(n-k+1)}{(n+1)^2(n+2)} V [ X ( k ) ] = ( n + 1 ) 2 ( n + 2 ) k ( n − k + 1 ) 。
2つの順序統計量の共分散・相関(一様分布)
k<lのとき、 C o v ( X ( k ) , X ( l ) ) = k ( n − l + 1 ) ( n + 1 ) 2 ( n + 2 ) \mathrm{Cov}(X_{(k)},X_{(l)})=\dfrac{k(n-l+1)}{(n+1)^2(n+2)} Cov ( X ( k ) , X ( l ) ) = ( n + 1 ) 2 ( n + 2 ) k ( n − l + 1 ) 、相関係数 ρ ( X ( k ) , X ( l ) ) = k ( n − l + 1 ) l ( n − k + 1 ) \rho(X_{(k)},X_{(l)})=\sqrt{\dfrac{k(n-l+1)}{l(n-k+1)}} ρ ( X ( k ) , X ( l ) ) = l ( n − k + 1 ) k ( n − l + 1 ) 。2022年度は一様分布のこの共分散・相関係数を求める大問で出題された。
ベータ関数・ベータ積分(順序統計量の正規化定数)
B ( a , b ) = ∫ 0 1 t a − 1 ( 1 − t ) b − 1 d t = Γ ( a ) Γ ( b ) Γ ( a + b ) B(a,b) = \int_0^1 t^{a-1}(1-t)^{b-1}\,dt = \dfrac{\Gamma(a)\Gamma(b)}{\Gamma(a+b)} B ( a , b ) = ∫ 0 1 t a − 1 ( 1 − t ) b − 1 d t = Γ ( a + b ) Γ ( a ) Γ ( b ) 。一様分布 U ( 0 , 1 ) U(0,1) U ( 0 , 1 ) の第 k k k 順序統計量の密度 n ! ( k − 1 ) ! ( n − k ) ! t k − 1 ( 1 − t ) n − k \dfrac{n!}{(k-1)!\,(n-k)!}t^{k-1}(1-t)^{n-k} ( k − 1 )! ( n − k )! n ! t k − 1 ( 1 − t ) n − k は B e t a ( k , n − k + 1 ) \mathrm{Beta}(k,\, n-k+1) Beta ( k , n − k + 1 ) そのもので、正規化定数が 1 / B ( k , n − k + 1 ) 1/B(k,\, n-k+1) 1/ B ( k , n − k + 1 ) にあたる。 ∫ 0 1 t a − 1 ( 1 − t ) b − 1 d t \int_0^1 t^{a-1}(1-t)^{b-1}\,dt ∫ 0 1 t a − 1 ( 1 − t ) b − 1 d t の形の定積分は、部分積分を繰り返さずこの公式に当てはめて Γ \Gamma Γ の比に直すのが定石 (a , b a, b a , b が正の整数なら ( a − 1 ) ! ( b − 1 ) ! ( a + b − 1 ) ! \dfrac{(a-1)!\,(b-1)!}{(a+b-1)!} ( a + b − 1 )! ( a − 1 )! ( b − 1 )! )。
最大値:分布関数
F ( n ) ( t ) = [ F ( t ) ] n F_{(n)}(t)=[F(t)]^n F ( n ) ( t ) = [ F ( t ) ] n 。『最大値が t t t 以下』=『全員が t t t 以下』で、独立なら確率の積になる。
最大値:密度の導出
f ( n ) ( t ) = n [ F ( t ) ] n − 1 f ( t ) f_{(n)}(t)=n[F(t)]^{n-1}f(t) f ( n ) ( t ) = n [ F ( t ) ] n − 1 f ( t ) (F ( t ) n F(t)^n F ( t ) n を t t t で微分)。『1人が t t t 付近・残り n − 1 n-1 n − 1 人が t t t 以下』を n n n 通りで選ぶ形。
最大値:一様分布の特例
U ( 0 , 1 ) U(0,1) U ( 0 , 1 ) なら f ( n ) ( t ) = n t n − 1 f_{(n)}(t)=n t^{n-1} f ( n ) ( t ) = n t n − 1 で ベータ分布 B e t a ( n , 1 ) \mathrm{Beta}(n,1) Beta ( n , 1 ) に従う。期待値 E = n / ( n + 1 ) E=n/(n+1) E = n / ( n + 1 ) 、分散 = n / [ ( n + 1 ) 2 ( n + 2 ) ] =n/[(n+1)^2(n+2)] = n / [( n + 1 ) 2 ( n + 2 )] 。
最大値: n → ∞ n\to\infty n → ∞ の挙動
n / ( n + 1 ) → 1 n/(n+1)\to1 n / ( n + 1 ) → 1 。最大値は支持の上限(一様なら1)に張り付き、分散も縮む。最大級リスクの評価に効く。
最大値:順序統計量の端
最大値は k = n k=n k = n の順序統計量。一般の k k k 番目の密度はベータ型で、最大値・最小値はその両端の特例にあたる。
▶ 最大値: U ( 0 , 1 ) U(0,1) U ( 0 , 1 ) から n = 5 n=5 n = 5 個取った最大値。密度 f ( 5 ) ( t ) = 5 t 4 f_{(5)}(t)=5t^4 f ( 5 ) ( t ) = 5 t 4 (0 < t < 1 0<t<1 0 < t < 1 ) 。
E [ max ] = ∫ 0 1 t ⋅ 5 t 4 d t = 5 ∫ 0 1 t 5 d t = 5 / 6 ≈ 0.8333 E[\max]=\int_0^1 t\cdot5t^4\,dt=5\int_0^1 t^5\,dt=5/6\approx0.8333 E [ max ] = ∫ 0 1 t ⋅ 5 t 4 d t = 5 ∫ 0 1 t 5 d t = 5/6 ≈ 0.8333 (公式 n / ( n + 1 ) = 5 / 6 n/(n+1)=5/6 n / ( n + 1 ) = 5/6 と一致)。
最大値は ベータ分布 B e t a ( 5 , 1 ) \mathrm{Beta}(5,1) Beta ( 5 , 1 ) に従い、分散 = n ( n + 1 ) 2 ( n + 2 ) = 5 36 × 7 = 5 252 ≈ 0.0198 =\dfrac{n}{(n+1)^2(n+2)}=\dfrac{5}{36\times7}=\dfrac{5}{252}\approx0.0198 = ( n + 1 ) 2 ( n + 2 ) n = 36 × 7 5 = 252 5 ≈ 0.0198 。
X ( k ) X_{(k)} X ( k ) (k番目に小さい値)と元の X k X_k X k (k番目に観測した値)を混同する 独立性を使わずに P ( 全てがx以下 ) P(\text{全てがx以下}) P ( 全てが x 以下 ) を確率の積に分解できない、と思い込む 2つの順序統計量 X ( k ) , X ( l ) ( k < l ) X_{(k)},X_{(l)}\ (k<l) X ( k ) , X ( l ) ( k < l ) の同時密度を1変数の公式のまま求めようとしてしまう。実際は3区間に分割した f X ( k ) , X ( l ) ( x k , x l ) = n ! ( k − 1 ) ! ( l − k − 1 ) ! ( n − l ) ! [ F ( x k ) ] k − 1 [ F ( x l ) − F ( x k ) ] l − k − 1 [ 1 − F ( x l ) ] n − l f ( x k ) f ( x l ) f_{X_{(k)},X_{(l)}}(x_k,x_l)=\dfrac{n!}{(k-1)!(l-k-1)!(n-l)!}[F(x_k)]^{k-1}[F(x_l)-F(x_k)]^{l-k-1}[1-F(x_l)]^{n-l}f(x_k)f(x_l) f X ( k ) , X ( l ) ( x k , x l ) = ( k − 1 )! ( l − k − 1 )! ( n − l )! n ! [ F ( x k ) ] k − 1 [ F ( x l ) − F ( x k ) ] l − k − 1 [ 1 − F ( x l ) ] n − l f ( x k ) f ( x l ) が必要(2022年度で共分散・相関係数の導出に使用) 密度を出すとき指数を間違える (F ( t ) n F(t)^n F ( t ) n の微分は n [ F ( t ) ] n − 1 f ( t ) n[F(t)]^{n-1}f(t) n [ F ( t ) ] n − 1 f ( t ) なので、 f ( t ) f(t) f ( t ) を掛け忘れたり指数を n n n のままにしない) 期待値を n / ( n + 1 ) n/(n+1) n / ( n + 1 ) と覚えるだけで導出を忘れる (∫ 0 1 t ⋅ n t n − 1 d t = n / ( n + 1 ) \int_0^1 t\cdot n t^{n-1}dt=n/(n+1) ∫ 0 1 t ⋅ n t n − 1 d t = n / ( n + 1 ) 。 一様分布に限った値で、他分布では別計算) 最大値と最小値の式を混同する(最大は F ( t ) n F(t)^n F ( t ) n 、最小は 1 − [ 1 − F ( t ) ] n 1-[1-F(t)]^n 1 − [ 1 − F ( t ) ] n 。最大は『全員以下』、最小は『全員より大』の余事象) 離散分布 23
確率