Turunkan distribusi Poisson bivariat


13

Saya baru-baru ini menemukan distribusi Poisson bivariat, tapi saya agak bingung bagaimana itu bisa diturunkan.

Distribusi diberikan oleh:

P(X=x,Y=y)=e(θ1+θ2+θ0)θ1xx!θ2yy!i=0min(x,y)(xi)(yi)i!(θ0θ1θ2)i

Dari apa yang saya dapat kumpulkan, istilah θ0 adalah ukuran korelasi antara X dan Y ; karenanya, ketika X dan Y bersifat independen, θ0=0 dan distribusinya hanya menjadi produk dari dua distribusi Poisson univariat.

Mengingat hal ini, kebingungan saya didasarkan pada jangka penjumlahan - Aku menduga istilah ini menjelaskan korelasi antara X dan Y .

Tampak bagi saya bahwa ringkasan merupakan semacam produk fungsi distribusi kumulatif binomial di mana probabilitas "sukses" diberikan oleh (θ0θ1θ2) dan probabilitas "kegagalan" diberikan oleh i!1min(x,y)i , karena (i!1min(x,y)i!)(min(x,y)i)=i!, tapi aku bisa pergi dengan ini.

Bisakah seseorang memberikan bantuan tentang bagaimana distribusi ini dapat diturunkan? Juga, jika itu dapat dimasukkan dalam jawaban apa pun bagaimana model ini dapat diperluas ke skenario multivarian (katakan tiga atau lebih variabel acak), itu akan bagus!

(Akhirnya, saya telah mencatat bahwa ada pertanyaan serupa yang diposting sebelumnya ( Memahami distribusi bivariat Poisson ), tetapi derivasi itu sebenarnya tidak dieksplorasi.)


2
Bukankah istilah pertama dengan eksponen adalah alih-alih ? e θ 1 + θ 2 + θ 0e(θ1+θ2+θ0)eθ1+θ2+θ0
Gilles

1
@ Giles Maaf, saya salah membaca komentar Anda awalnya - ya, Anda benar; istilah tersebut harus dibaca . Terima kasih telah menangkap itu! e(θ1+θ2+θ0)
user9171

3
Secara umum, ini bukan "the" untuk versi multivariate dari distribusi univariat, dengan beberapa pengecualian konvensional ("the" multivariate normal misalnya). Ada banyak cara untuk mendapatkan ekstensi multivarian, tergantung pada fitur mana yang paling penting untuk dimiliki. Penulis yang berbeda mungkin memiliki versi multivariat yang berbeda untuk distribusi univariat yang umum. Jadi secara umum, orang mungkin mengatakan sesuatu seperti " sebuah Poisson multivariat", atau 'Jadi-dan-jadi itu bivariat Poisson" Ini adalah salah satu yang cukup alami, tapi bukan satu-satunya..
Glen_b -Reinstate Monica

2
(ctd) ... mis. beberapa penulis mencari distribusi multivariat yang mampu ketergantungan negatif, kemampuan yang tidak dimiliki yang satu ini.
Glen_b -Reinstate Monica

Jawaban:


17

Dalam presentasi slide , Karlis dan Ntzoufras mendefinisikan Poisson bivariat sebagai distribusi mana secara independen memiliki distribusi Poisson . Ingatlah bahwa memiliki sarana distribusi seperti ituX i θ i(X,Y)=(X1+X0,X2+X0)Xiθi

Pr(Xi=k)=eθiθikk!

untukk=0,1,2,.

Peristiwa adalah penyatuan yang terpisah dari peristiwa tersebut(X,Y)=(x,y)

(X0,X1,X2)=(i,xi,yi)

untuk semua yang membuat ketiga komponen bilangan bulat non-negatif, dari mana kita dapat menyimpulkan bahwa . Karena independen probabilitasnya berlipat ganda, dari mana0 i min ( x , y ) X ii0imin(x,y)Xi

F(θ0,θ1,θ2)(x,y)=Pr((X,Y)=(x,y))=i=0min(x,y)Pr(X0=i)Pr(X1=xi)Pr(X2=yi).

Ini adalah formula; kita sudah selesai. Tetapi untuk melihat bahwa itu setara dengan rumus dalam pertanyaan, gunakan definisi distribusi Poisson untuk menulis probabilitas ini dalam hal parameter dan (dengan asumsi tidak ada dari adalah nol) secara aljabar agar terlihat sebanyak mungkin seperti produk :θ 1 , θ 2 Pr ( X 1 = x ) Pr ( X 2 = y )θiθ1,θ2Pr(X1=x)Pr(X2=y)

F(θ0,θ1,θ2)(x,y)=i=0min(x,y)(eθ0θ0ii!)(eθ1θ1xi(xi)!)(eθ2θ2yi(yi)!)=e(θ1+θ2)θ1xx!θ2yy!(eθ0i=0min(x,y)θ0ii!x!θ1i(xi)!y!θ2i(yi)!).

Jika Anda benar-benar ingin - ini agak sugestif - Anda dapat mengekspresikan kembali istilah dalam jumlah menggunakan koefisien binomial Dan , menghasilkan(xi)=x!/((xi)!i!)(yi)

F(θ0,θ1,θ2)(x,y)=e(θ0+θ1+θ2)θ1xx!θ2yy!i=0min(x,y)i!(xi)(yi)(θ0θ1θ2)i,

persis seperti dalam pertanyaan.


Skenario generalisasi ke multivarian dapat dilanjutkan dalam beberapa cara, tergantung pada fleksibilitas yang dibutuhkan. Yang paling sederhana akan merenungkan distribusi

(X1+X0,X2+X0,,Xd+X0)

untuk terdistribusi Poisson independen . Untuk lebih fleksibel, variabel tambahan dapat diperkenalkan. Misalnya, gunakan variabel Poisson independen dan pertimbangkan distribusi multivariat dari ,X0,X1,,XdηiY1,,YdXi+(Yi+Yi+1++Yd)i=1,2,,d.


1
pujian! Btw, bukankah seharusnya dalam tanda kurung besar sebelum langkah terakhir adalah ? eθ0eθ2
Gilles

1
@Gilles Terima kasih telah menangkap kesalahan ketik - Saya memperbaikinya. Eksponen awal harus ; yang dalam kurung adalah benar. θ0+θ1θ1+θ2eθ0
whuber

@whuber Terima kasih banyak! Itu jawaban yang sempurna!
user9171

@whuber Jawaban yang bagus! Saya masih tidak melihat mengapa acara harus menjadi persatuan acara yang terpisah . Saya kira ini hanya berlaku untuk . Mungkin yang Anda maksud (komponen-bijaksana)? Tetapi apakah itu cukup untuk menandai fungsi distribusi? (X,Y)=(x,y)(X0,X1,X2)=(i,xi,yi)i=0(X,Y)(x,y)
vanguard2k

@ vanguard2k Saya tidak mengerti komentar Anda. Apakah Anda menegaskan peristiwa itu tidak terputus - putus? (Namun mereka harus, karena mereka memiliki nilai .) Atau apakah Anda menyatakan mereka tidak lengkap? (Jika demikian, berapa nilai ( menurut Anda belum termasuk?)X0(X,Y)
whuber

4

Berikut adalah cara untuk menurunkan distribusi bivariat poisson.

Biarkan menjadi variabel acak poisson independen dengan parameter . Kemudian kita mendefinisikan . Variabel , umum untuk dan , menyebabkan pasanganX0,X1,X2θ0,θ1,θ2Y1=X0+X1,Y2=X0+X2X0Y1Y2(Y1,Y2) berkorelasi. Maka kita harus menghitung probabilitas massa funtion:

P(Y1=y1,Y2=y2)=P(X0+X1=y1,X0+X2=y2)=x0=0min(y1,y2)P(X0=x0)P(X1=y1x0)P(X2=y2y0)=x0=0min(y1,y2)eθ0θ0x0x0!eθ1θ1y1x0(y1x0)!eθ2θ2y2x0(y2x0)!=eθ0θ1θ2θ1y1θ2y2x0=0min(y1,y2)(θ0θ1θ2)x0x0!(y1x0)(y2x0)

Semoga ini bisa membantu!

1
Hai Kjetil - Saya memperbaiki masalah dengan format (tetapi, ingin mengubah sesedikit mungkin, meninggalkan beberapa kesalahan ketik utuh). Saya tidak mengerti mengapa Anda memposting replika derivasi pada jawaban saya sebelumnya, terutama ketika Anda kehilangan beberapa faktor penting di sepanjang jalan yang menyebabkan hasil akhir salah. Apakah ada poin tertentu yang ingin Anda sampaikan? TEX
whuber

1
whuber: Saya mulai menulis jawaban saya sebelum jawaban Anda diposting! kalau tidak, saya tidak akan menulisnya.
kjetil b halvorsen
Dengan menggunakan situs kami, Anda mengakui telah membaca dan memahami Kebijakan Cookie dan Kebijakan Privasi kami.
Licensed under cc by-sa 3.0 with attribution required.