Mengubah Statistik Pesanan

9

X1,...,XnY1,...,YnU(0,a)Zn=nlogmax(Y(n),X(n))min(Y(n),X(n))Exp(1)

Saya telah memulai masalah ini dengan mengatur Kemudian akan didistribusikan sebagai dan akan didistribusikan sebagai Kepadatan dapat ditemukan dengan mudah seperti dan{X1,...,Xn,Y1,...Yn}={Z1,...,Zn}max(Yn,Xn)=Z(2n)(za)2nmin(Yn,Xn)=Z(1)1(1za)2nfZ1(z)=(2n)(1za)2n11afZ(2n)(z)=(2n)(za)2n11a

Di sinilah saya mengalami kesulitan mengetahui ke mana harus pergi selanjutnya sekarang ini dihitung. Saya pikir itu harus melakukan sesuatu dengan transformasi, tapi saya tidak yakin ...

Susan
sumber
Tentunya Anda perlu berasumsi bahwa bukan hanya dan iid, tetapi juga tidak tergantung pada . Mengingat itu, apakah Anda berpikir untuk bekerja langsung dengan ? Y i X i Y j log ( Z i )XiYiXiYjlog(Zi)
whuber
@whuber pemikiran saya dari komentar Anda akan mengatur transformasi di mana saya memecahkan kepadatan n * log (Z )? i
Susan
Saya melakukan sedikit pemformatan ulang (terutama mengubah dan menjadi dan ) tetapi jika Anda tidak menyukainya, Anda dapat kembali ke versi sebelumnya (dengan mengklik tautan "diedit <x> lalu" di atas gravatar saya di bagian bawah posting Anda) dan kemudian mengklik tautan "putar kembali" di atas versi Anda sebelumnya. m i n log minlogminlogmin
Glen_b -Reinstate Monica
3
Susan, Anda tampaknya salah mengartikan / salah membaca pertanyaan. Pertanyaannya mencari rasio Penyebut mengacu pada : di mana adalah statistik urutan maksimum s, dan adalah statistik urutan maksimum s. Dengan kata lain, mencari min (maxX, maxY), BUKAN minimum semua dan , sehingga Anda tidak dapat menggunakan trik Z Anda untuk meratakan / menggabungkan semua nilai X dan Y. ....... min(Y(n),X(n))Y(n)YX(n)Xmin(Y(n),X(n))XY
max(Y(n),X(n))min(Y(n),X(n))
min(Y(n),X(n))Y(n)YX(n)Xmin(Y(n),X(n))XY
serigala
2
Dalam keadaan apa pun, dan sebagai masalah terpisah, tidak ada gunanya (seperti yang telah Anda lakukan) menghitung kepadatan , dan secara terpisah kepadatan , karena statistik urutan berbeda tidak umumnya mandiri. Untuk menemukan rasio , seseorang harus terlebih dahulu menemukan pdf gabungan dari , jika itu masalahnya di tangan (yang bukan). Z ( 2 n ) Z ( 2 n ) / Z ( 1 ) ( Z ( 1 ) , Z ( 2 n ) )Z(1)Z(2n)Z(2n)/Z(1)(Z(1),Z(2n))
serigala

Jawaban:

2

Masalah ini dapat diselesaikan dari definisi saja: satu-satunya perhitungan lanjutan adalah integral dari monomial.


Pengamatan awal

Kerja Mari kita dengan variabel dan seluruh: ini tidak mengubah tapi itu membuat IID dengan Uniform distribusi, menghilangkan semua penampilan mengganggu dari dalam perhitungan. Dengan demikian kita dapat mengasumsikan tanpa kehilangan keumuman.Y i / a Z n ( X 1 , ... , Y n ) ( 0 , 1 ) a a = 1Xi/aYi/aZn(X1,,Yn)(0,1)aa=1

Perhatikan bahwa independensi dan distribusi seragamnya menyiratkan bahwa untuk setiap angka yang , y 0 y 1Yiy0y1

Pr(yY(n))=Pr(yY1,,yYn)=Pr(yY1)Pr(yYn)=yn,

dengan hasil identik untuk . Untuk referensi di masa mendatang, ini memungkinkan kami untuk menghitungX(n)

E(2X(n)n)=012xnd(xn)=012nx2n1dx=1.

Larutan

Biarkan menjadi bilangan real positif. Untuk menemukan distribusi , gantilah definisi dan sederhanakan ketimpangan yang dihasilkan:Z ntZn

Pr(Zn>t)=Pr(Zn/n>t/n)=Pr(exp(Zn/n)>et/n)=Pr(max(X(n),Y(n))min(X(n),Y(n))>et/n)=Pr(et/nmax(X(n),Y(n))>min(X(n),Y(n))).

Peristiwa ini terbagi menjadi dua kasus yang dapat disetel, tergantung pada apakah atau adalah yang lebih kecil dari keduanya (dan persimpangan mereka, dengan probabilitas nol, dapat diabaikan). Jadi kita hanya perlu menghitung kemungkinan salah satu dari kasus ini (katakanlah di mana adalah yang lebih kecil) dan gandakan. Karena , , memungkinkan kita (setelah membiarkan untuk memainkan peran of ) untuk menerapkan perhitungan di bagian pendahuluan: Y ( n ) Y ( n ) t 0 0 e - t / n X ( n )1 e - t / n X ( n ) yX(n)Y(n)Y(n)t00et/nX(n)1et/nX(n)y

Pr(Zn>t)=2Pr(et/nX(n)>Y(n))=2E[(et/nX(n))n]=etE[2X(n)n]=et.

Itulah artinya bagi untuk memiliki distribusi Exp . ( 1 )Zn(1)

whuber
sumber
4

Saya akan membuat sketsa solusinya, di sini menggunakan sistem aljabar komputer untuk melakukan seluk-beluk ...

Larutan

Jika adalah sampel ukuran pada induk , maka pdf dari sampel maksimum adalah: dan juga untuk . n X Seragam ( 0 , a ) f n ( x ) = nX1,...,XnnXUniform(0,a)Y

fn(x)=nanxn1
Y

Pendekatan 1: Temukan pdf gabungan dari(X(n),Y(n))

Karena dan bersifat independen, pdf gabungan dari 2 sampel maksimum hanyalah produk dari 2 pdf, katakanlah :Y ( X ( n ) , Y ( n ) ) f ( n ) ( x , y )XY(X(n),Y(n))f(n)(x,y)

masukkan deskripsi gambar di sini

Diberikan . Kemudian, cdf dari adalah adalah: ZnP(Zn<z)Zn=nlogmax(Y(n),X(n))min(Y(n),X(n))ZnP(Zn<z)

masukkan deskripsi gambar di sini

di mana saya menggunakan Probfungsi dari paket mathStatica untuk Mathematica untuk mengotomatisasi. Membedakan cdf wrt menghasilkan pdf sebagai standar Eksponensial.Z nzZn


Pendekatan 2: Statistik pesanan

Kita dapat menggunakan statistik pesanan untuk 'by-pass' mekanisme karena harus berurusan dengan fungsi Max dan Min.

Sekali lagi: Jika adalah contoh ukuran pada induk , maka pdf dari sampel maksimum adalah, katakan, : n X Seragam ( 0 , a ) W = X ( n ) f n ( w )X1,...,XnnXUniform(0,a)W=X(n)fn(w)

masukkan deskripsi gambar di sini

Sampel maksimum dan hanyalah dua gambar independen dari distribusi ; yaitu statistik pesanan dan dari (dalam sampel ukuran 2) persis seperti yang kami cari:X(n)Y(n)W1st2ndW

  • W(1)=min(Y(n),X(n))

  • W(2)=max(Y(n),X(n))

Pdf gabungan dari , dalam sampel ukuran 2, misalkan , Adalah:(W(1),W(2))g(.,.)

masukkan deskripsi gambar di sini

Diberikan . Kemudian, cdf dari adalah adalah: ZnP(Zn<z)Zn=nlogmax(Y(n),X(n))min(Y(n),X(n))ZnP(Zn<z)

masukkan deskripsi gambar di sini

Keuntungan dari pendekatan ini adalah bahwa perhitungan probabilitas tidak lagi melibatkan fungsi max / min, yang dapat membuat derivasi (terutama dengan tangan) agak lebih mudah untuk diekspresikan.

Lain

Sesuai komentar saya di atas, tampaknya Anda salah menafsirkan pertanyaan ...

Kami diminta untuk menemukan:

Zn=nlogmax(Y(n),X(n))min(Y(n),X(n))

di mana penyebut adalah min (xmax, ymax), ... tidak minimal semua 's dan ' s.YXY

serigala
sumber
Mengikuti sketsa Anda, saya mengerti bagaimana saya salah mengartikan pertanyaan itu. Saya mengerti bagaimana cara menghitung pdf bersama dari dua sampel maksimum, tetapi saya masih tidak yakin bagaimana kita harus menafsirkan rasio maks / mnt.
Susan
Saya telah menambahkan derivasi alternatif menggunakan statistik pesanan, yang 'menghindari' max / min.
serigala
Jika Anda mulai dengan log data, Susan, maka Anda akan melihat perbedaan statistik urutan daripada rasio .
whuber
Saya tidak yakin menggunakan perhitungan formal komputer adalah cara terbaik untuk menjelaskan alasan mengapa rasio adalah variabel acak Exp (1).
Xi'an
1
Poin bagus ... kecuali OP tidak menanyakan alasannya ... tetapi untuk menunjukkan bahwa itu Exp [1]. Saya juga tidak yakin apakah ini adalah pekerjaan rumah (atau tugas) atau tidak ... dan itu sebenarnya satu keuntungan baik dari menggunakan komputer: seseorang memberikan langkah-langkah untuk diikuti, memverifikasi hasilnya, sehingga seseorang memiliki pendekatan yang tepat , tetapi mekanika masih diserahkan kepada OP. Akan menyenangkan bagi seseorang untuk mengeksplorasi saran @ whuber untuk mengambil log di awal.
serigala