Bukan jawaban yang lengkap, maaf, tetapi beberapa ide (untuk merindukan komentar). Perhatikan bahwa yang Anda miliki adalah produk dari iid variabel acak, di mana adalah variabel acak (rv) dengan distribusi poisson dengan parameter . Itu dapat digunakan untuk "pemeriksaan kewarasan" lain, simulasi (menggunakan eksponensial tingkat 1):K+1Kλ
set.seed(7*11*13)
N <- 1000000
prods <- rep(0, N)
ks <- rpois(N, 1)+1
for (i in 1:N) {
k <- ks[i]
prods[i] <- prod( rexp(k, 1))
}
qqplot( qexp(ppoints(N)), prods)
Hasilnya qqplot
(tidak diperlihatkan di sini) jauh dari garis lurus, jadi ini tidak terlihat sebagai eksponensial dari nilai 1. Rata-rata benar, varians ke besar, ada ekor kanan jauh lebih lama daripada eksponensial. Apa yang bisa dilakukan secara teoritis? Transformasi Mellin https://en.wikipedia.org/wiki/Mellin_transform disesuaikan dengan produk dari variabel acak independen. Saya akan menghitung hanya untuk eksponensial dengan laju 1. Transformasi Mellin dari kemudian adalah
sehingga transformasi Mellin dari produk eksponensial adalah
SejakV0
M1(s)=EVs0=∫∞0xse−xdx=Γ(s+1)
k+1Mk+1(s)=Γ(s+1)k+1
Kmemiliki distribusi poisson dengan parameter , transformasi Mellin dari produk acak dari bilangan acak faktor , adalah
tetapi saya tidak dapat menemukan kebalikan dari transformasi ini. Tetapi perhatikan bahwa jika adalah variabel acak nonnegatif dengan Mellin mengubah , maka mendefinisikan kita menemukan bahwa
sehingga transformasi Mellin dari adalah fungsi penghasil momen dari logaritma
λK+1M(s)=EMK+1(s)=EΓ(s+1)K+1=Γ(s+1)e−λ∑k=0∞λkk!Γ(s+1)k=e−λΓ(s+1)eλΓ(s+1)
XMX(t)Y=logXKY(t)=EetY=EetlogX=Eelog(Xt)=EXt=MX(y)
XY. Jadi, dengan menggunakan itu kita dapat memperkirakan distribusi dengan metode pendekatan saddlepoint,
Bagaimana cara kerja pendekatan saddlepoint? dan cari situs ini.
X