Maaf untuk jawaban yang terlambat, tetapi ini menyadap saya dan saya menemukan jawabannya. Distribusi memang Dirichlet-Multinomial dan neg individu. distribusi binomial bahkan tidak perlu identik, selama faktor Fano mereka (rasio varian terhadap rata-rata) identik.
Jawaban panjang:
Jika Anda parameterkan NB sebagai:
p(X=x|λ,θ)=NB(x|λ,θ)=(θ−1λ+x−1x)(11+θ−1)x(θ−11+θ−1)θ−1λ
Kemudian dan danE(X)=λVar(X)=λ(1+θ)
∀i:Xi∼NB(λi,θ) menyiratkan
∑Xi∼NB(∑λi,θ)
Kemudian mengambil probabilitas dengan memberi jumlah:
∏NB(xi|λi,θ)NB(∑xi|∑λi,θ)=(11+θ−1)∑xi(θ−11+θ−1)θ−1∑λi∏(θ−1λi+xi−1xi)(11+θ−1)∑xi(θ−11+θ−1)θ−1∑λi(θ−1∑λi+∑xi−1∑xi)==Γ(∑xi+1)Γ(θ−1∑λi)Γ(θ−1∑λi+∑xi)∏Γ(θ−1λi+xi)Γ(xi+1)Γ(θ−1λi)=DM(x1,...,xn|θ−1λ1,...,θ−1λn)
di mana adalah kemungkinan Dirichlet-Multinomial. Ini hanya dihasilkan dari fakta bahwa kecuali untuk koefisien multinomial, banyak istilah dalam fraksi di sisi kiri dibatalkan, sehingga Anda hanya dengan istilah fungsi gamma yang kebetulan sama dengan kemungkinan DM.DM
Perhatikan juga bahwa parameter model ini tidak dapat diidentifikasi sebagai peningkatan dengan penurunan simultan di semua menghasilkan kemungkinan yang persis sama.θλi
Referensi terbaik yang saya miliki untuk bagian ini adalah bagian 2 hingga 3.1 dari Guimarães & Lindrooth (2007): Mengontrol penyebaran berlebih dalam model logit bersyarat bersyarat: Aplikasi sederhana untuk regresi Dirichlet-multinomial - sayangnya paywall, tetapi saya tidak dapat menemukan referensi yang tidak dibayar di dinding.