Perkiraan statistik pesanan untuk variabel acak normal

Adakah formula terkenal untuk statistik urutan distribusi acak tertentu? Khususnya statistik urutan pertama dan terakhir dari variabel acak normal, tetapi jawaban yang lebih umum juga akan dihargai.

Sunting: Untuk memperjelas, saya mencari rumus perkiraan yang dapat lebih atau kurang dievaluasi secara eksplisit, bukan ekspresi integral yang tepat.

Sebagai contoh, saya telah melihat dua perkiraan berikut untuk statistik orde pertama (yaitu minimum) dari rv normal:

$e_{1:n} \geq \mu - \frac{n-1}{\sqrt{2n-1}}\sigma$

dan

$e_{1:n} \approx \mu + \Phi^{-1} \left( \frac{1}{n+1} \right)\sigma$

Yang pertama, untuk $n=200$ , memberikan sekitar $e_{1:200} \geq \mu - 10\sigma$ yang tampak seperti ikatan yang sangat longgar.

Yang kedua memberi $e_{1:200} \approx \mu - 2.58\sigma$ sedangkan Monte Carlo yang cepat memberi , jadi itu bukan perkiraan yang buruk tapi tidak bagus juga, dan yang lebih penting aku tidak ' Saya tidak memiliki intuisi tentang dari mana asalnya.$e_{1:200} \approx \mu - 2.75\sigma$

Ada bantuan?

Referensi klasik adalah Royston (1982) [1] yang memiliki algoritma melampaui rumus eksplisit. Ini juga mengutip formula terkenal oleh Blom (1958): $E(r:n) \approx \mu + \Phi^{-1}(\frac{r-\alpha}{n-2\alpha+1})\sigma$dengan$\alpha=0.375$. Formula ini memberikan pengali -2,73 untuk$n=200, r=1$.

[1]: Algoritma AS 177: Statistik Orde Normal yang Diharapkan (Tepat dan Perkiraan) JP Royston. Jurnal Masyarakat Statistik Kerajaan. Seri C (Statistik Terapan) Vol. 31, No. 2 (1982), hlm. 161-165

$$\newcommand{\Pr}{\mathrm{Pr}}\newcommand{\Beta}{\mathrm{Beta}}\newcommand{\Var}{\mathrm{Var}}$$ Distribusi statistik urutan ke-i dari setiap variabel acak kontinu dengan PDF diberikan oleh distribusi senyawa "beta-F". Cara intuitif untuk berpikir tentang distribusi ini, adalah untuk mempertimbangkan urutan engan statistik dalam sampel . Sekarang agar nilai statistik urutan ke-i dari variabel acak X menjadi sama dengan x kita membutuhkan 3 kondisi:$N$$X$$x$

1. bawah x , ini memiliki probabilitas F X ( x ) untuk setiap pengamatan, di mana F X ( x ) = P r ( X < x ) adalah CDF dari variabel acak X.$i-1$$x$$F_{X}(x)$$F_X(x)=\Pr(X<x)$
2. atas x , ini memiliki probabilitas 1 - F X ( x $N-i$$x$$1-F_{X}(x)$
3. 1 nilai di dalam interval sangat kecil yang mengandung , ini memiliki probabilitas f X ( x ) d x di mana f X ( x ) d x = d F X ( x ) = P r ( x < X < x + d x ) adalah PDF dari variabel acak $x$$f_{X}(x)dx$$f_{X}(x)dx=dF_{X}(x)=\Pr(x<X<x+dx)$$X$

Ada cara untuk membuat pilihan ini, jadi kami memiliki:${N \choose 1}{N-1 \choose i-1}$

$$f_{i}(x_{i})=\frac{N!}{(i-1)!(N-i)!}f_{X}(x_{i})\left[1-F_{X}(x_{i})\right]^{N-i}\left[F_{X}(x_{i})\right]^{i-1}dx$$

Sunting di posting asli saya, saya melakukan upaya yang sangat buruk untuk melangkah lebih jauh dari titik ini, dan komentar di bawah mencerminkan ini. Saya telah berusaha untuk memperbaikinya di bawah ini

Jika kita mengambil nilai rata-rata dari pdf ini, kita mendapatkan:

$$E(X_{i})=\int_{-\infty}^{\infty} x_{i}f_{i}(x_{i})dx_{i}$$

Dan dalam integral ini, kami membuat perubahan berikut variabel (mengambil petunjuk @ henry), dan integral menjadi:$p_{i}=F_{X}(x_{i})$

$$E(X_{i})=\int_{0}^{1} F_{X}^{-1}(p_{i})\Beta(p_{i}|i,N-i+1)dp_{i}=E_{\Beta(p_{i}|i,N-i+1)}\left[F_{X}^{-1}(p_{i})\right]$$

Jadi ini adalah nilai yang diharapkan dari CDF terbalik, yang dapat didekati dengan baik menggunakan metode delta untuk memberikan:

$$E_{\Beta(p_{i}|i,N-i+1)}\left[F_{X}^{-1}(p_{i})\right]\approx F_{X}^{-1}\left[E_{\Beta(p_{i}|i,N-i+1)}\right]=F_{X}^{-1}\left[\frac{i}{N+1}\right]$$

To make a better approximation, we can expand to 2nd order (prime denoting differentiation), and noting that the second derivative of an inverse is:

$$\frac{\partial^{2}}{\partial a^{2}}F_{X}^{-1}(a)=-\frac{F_{X}^{''}(F_{X}^{-1}(a))}{\left[F_{X}^{'}(F_{X}^{-1}(a))\right]^{3}}=-\frac{f_{X}^{'}(F_{X}^{-1}(a))}{\left[f_{X}(F_{X}^{-1}(a))\right]^{3}}$$

Let $\nu_{i}=F_{X}^{-1}\left[\frac{i}{N+1}\right]$. Then We have:

$$E_{\Beta(p_{i}|i,N-i+1)}\left[F_{X}^{-1}(p_{i})\right]\approx F_{X}^{-1}\left[\nu_{i}\right]-\frac{\Var_{\Beta(p_{i}|i,N-i+1)}\left[p_{i}\right]}{2}\frac{f_{X}^{'}(\nu_{i})}{\left[f_{X}(\nu_{i})\right]^{3}}$$ $$=\nu_{i}-\frac{\left(\frac{i}{N+1}\right)\left(1-\frac{i}{N+1}\right)}{2(N+2)}\frac{f_{X}^{'}(\nu_{i})}{\left[f_{X}(\nu_{i})\right]^{3}}$$

Now, specialising to normal case we have

$$f_{X}(x)=\frac{1}{\sigma}\phi(\frac{x-\mu}{\sigma})\rightarrow f_{X}^{'}(x)=-\frac{x-\mu}{\sigma^{3}}\phi(\frac{x-\mu}{\sigma})=-\frac{x-\mu}{\sigma^{2}}f_{X}(x)$$ $$F_{X}(x)=\Phi(\frac{x-\mu}{\sigma})\implies F_{X}^{-1}(x)=\mu+\sigma\Phi^{-1}(x)$$

Note that $f_{X}(\nu_{i})=\frac{1}{\sigma}\phi\left[\Phi^{-1}\left(\frac{i}{N+1}\right)\right]$ And the expectation approximately becomes:

$$E[x_{i}]\approx \mu+\sigma\Phi^{-1}\left(\frac{i}{N+1}\right)+\frac{\left(\frac{i}{N+1}\right)\left(1-\frac{i}{N+1}\right)}{2(N+2)}\frac{\sigma\Phi^{-1}\left(\frac{i}{N+1}\right)}{\left[\phi\left[\Phi^{-1}\left(\frac{i}{N+1}\right)\right]\right]^{2}}$$

And finally:

$$E[x_{i}]\approx \mu+\sigma\Phi^{-1}\left(\frac{i}{N+1}\right)\left[1+\frac{\left(\frac{i}{N+1}\right)\left(1-\frac{i}{N+1}\right)}{2(N+2)\left[\phi\left[\Phi^{-1}\left(\frac{i}{N+1}\right)\right]\right]^{2}}\right]$$

Although as @whuber has noted, this will not be accurate in the tails. In fact I think it may be worse, because of the skewness of a beta with different parameters

Jawaban Aniko bergantung pada formula terkenal Blom yang melibatkan pilihan $\alpha = 3/8$. Ternyata rumus ini sendiri merupakan perkiraan belaka dari jawaban yang tepat karena G. Elfving (1947), Distribusi asimtot dari rentang sampel dari populasi normal , Biometrika, Vol. 34, hlm. 111-119. Rumus Elfving ditujukan pada minimum dan maksimum sampel, untuk pilihan alpha yang tepat$\pi/8$. Formula hasil Blom ketika kami memperkirakan$\pi$ oleh $3$.

Dengan menggunakan rumus Elfving daripada perkiraan Blom, kita mendapatkan pengganda -2,744165. Angka ini lebih dekat dengan jawaban pasti Erik P. (-2.746) dan perkiraan Monte Carlo (-2.75) daripada perkiraan Blom (-2.73), sementara lebih mudah diterapkan daripada rumus yang tepat.

Depending on what you want to do, this answer may or may not help - I got the following exact formula from Maple's Statistics package.

with(Statistics):
X := OrderStatistic(Normal(0, 1), 1, n):
m := Mean(X):
m;

$$\int _{-\infty }^{\infty }\!1/2\,{\frac {{\it \_t0}\,n!\,\sqrt {2}{ {\rm e}^{-1/2\,{{\it \_t0}}^{2}}} \left( 1/2-1/2\, {{\rm erf}\left(1/2\,{\it \_t0}\,\sqrt {2}\right)} \right) ^{-1+n}}{ \left( -1+n \right) !\,\sqrt {\pi }}}{d{\it \_t0}}$$

By itself this isn't very useful (and it could probably be derived fairly easily by hand, since it's the minimum of $n$ random variables), but it does allow for quick and very accurate approximation for given values of $n$ - much more accurate than Monte Carlo:

evalf(eval(m, n = 200));
evalf[25](eval(m, n = 200));

gives -2.746042447 and -2.746042447451154492412344, respectively.

(Full disclosure - I maintain this package.)