Apakah Poisson nol terpotong dan Poisson dasar bersarang atau tidak bersarang?

9

Saya telah melihat banyak yang membahas apakah regresi Poisson dasar adalah versi bersarang dari regresi Poisson nol-meningkat. Misalnya situs ini berpendapat bahwa itu adalah, karena yang terakhir mencakup parameter tambahan untuk memodelkan nol tambahan, tetapi sebaliknya menyertakan parameter regresi Poisson yang sama dengan yang sebelumnya, meskipun halaman tersebut tidak menyertakan referensi yang tidak setuju.

Apa yang saya tidak dapat menemukan informasi tentang apakah Poisson nol terpotong dan Poisson dasar bersarang. Jika Poisson nol terpotong hanya Poisson dengan ketentuan tambahan bahwa probabilitas nol dihitung adalah nol, maka saya kira kedengarannya seperti itu, tapi saya berharap untuk jawaban yang lebih pasti.

Alasan saya bertanya-tanya adalah bahwa itu akan mempengaruhi apakah saya harus menggunakan tes Vuong (untuk model non-bersarang), atau tes chi-square yang lebih mendasar berdasarkan perbedaan dalam kemungkinan loglikel (untuk model bersarang).

Wilson (2015) berbicara tentang apakah tes Vuong sesuai untuk membandingkan regresi nol-inflasi dengan yang dasar, tetapi saya tidak dapat menemukan sumber yang membahas data nol-terpotong.

Justin
sumber

Jawaban:

4

Temukan ini sekarang. Untuk menghindari kebingungan, saya adalah Wilson dari Wilson (2015) yang dirujuk dalam pertanyaan asli, yang menanyakan apakah model Poisson dan Poisson terpotong bersarang, tidak bersarang, dll. Menyederhanakan, model yang lebih kecil bersarang dalam model yang lebih besar jika lebih besar model berkurang ke yang lebih kecil jika subset parameternya ditetapkan pada nilai yang dinyatakan; dua model tumpang tindih jika keduanya mengurangi ke model yang sama ketika himpunan bagian dari parameter masing-masing ditetapkan ke nilai-nilai tertentu, mereka tidak bersarang jika tidak peduli bagaimana parameter diperbaiki satu tidak dapat mengurangi yang lain. Menurut definisi ini Poisson terpotong dan Poisson standar tidak bersarang. NAMUN, dan ini adalah poin yang tampaknya telah diabaikan oleh banyak orang, teori distribusi Vuong mengacu pada STRICTLY nested, STRICTLY non-nested, dan SANGAT tumpang tindih. "STRICTLY" mengacu pada penambahan enam pembatasan pada definisi dasar bertumpuk dll. Pembatasan ini tidak persis sederhana, tetapi mereka, antara lain, berarti bahwa hasil Vuong tentang distribusi rasio kemungkinan log tidak berlaku dalam kasus di mana model / distribusi bersarang pada batas ruang parameter (seperti halnya Poisson / zero inflated Poisson dengan tautan identitas untuk parameter nol-inflasi) atau ketika satu model cenderung ke yang lain ketika parameter cenderung tak terhingga, seperti adalah kasus dengan Poisson / zero-inflated Poisson ketika tautan logit digunakan untuk memodelkan parameter nol-inflasi. Vuong tidak mengajukan teori tentang distribusi rasio kemungkinan log dalam keadaan ini. Sayangnya di sini,

Kode R berikut akan mensimulasikan distribusi poisson dan rasio loglikelihood Poisson terpotong. Itu membutuhkan VGAMpaket.

n<-30   
lambda1<-1
H<-rep(999,10000)
for(i in 1:10000){
  print(i)
  y<-rpospois(n, lambda1)
  fit1 <- vglm(y ~ 1, pospoisson)
  fit2<-glm(y~1, family=poisson(link="log"))
  H[i]<-logLik(fit1)-logLik(fit2)
}

hist(H,col="lemonchiffon")
Pauljw11
sumber
4

Poisson dasar dapat dianggap bersarang di dalam bentuk yang lebih umum:

p(x)=(1p)eλλxx!+p1(x=0)

Ketika , kita memiliki Poisson dasar. Ketika , kita memiliki Poisson tanpa terpotong. Ketika , kami memiliki Poisson yang dikurangi nol. Ketika , kami memiliki Poisson nol-inflasi, dan kami memiliki distribusi yang merosot pada .p=0p=exp{λ}/(1exp{λ})exp{λ}/(1exp{λ})<p<00<p<1p=1

Jadi menurut saya versi uji Vuong yang bersarang, atau chi-square seperti yang Anda sarankan, akan sesuai untuk kasus Anda. Perhatikan, bagaimanapun, bahwa chi-square dapat memiliki masalah karena probabilitas kecil dari pengamatan "besar" (relatif terhadap ). Anda mungkin ingin menggunakan bootstrap untuk mendapatkan nilai p untuk statistik chi-square alih-alih mengandalkan asimptotik kecuali Anda memiliki banyak data.λ

Jbowman
sumber
Terima kasih @jbowman - itu semacam jawaban yang lebih ketat yang saya harapkan. Namun saya tidak jelas: Saya pikir inti dari tes Vuong adalah untuk model yang tidak bersarang, jadi meskipun melampaui postingan asli saya, dapatkah Anda memberikan sedikit informasi lebih lanjut tentang "versi tes Vuong yang bersarang". Untuk menjadi jelas tentang sumber kebingungan saya: sampai saat ini saya hanya mengetahui vuongfungsi dalam paket pscldi R yang mengatakan itu untuk model non-bersarang. Saya baru saja googled dan menemukan fungsi vuongtestdalam paket nonnest2yang termasuk argumen 'bersarang'. Itu saja?
Justin
Iya itu. Sebenarnya, halaman Wikipedia en.wikipedia.org/wiki/Vuong%27s_closeness_test pada tes Vuong sedikit membantu (seringkali tidak terlalu banyak) dalam menjelaskan perbedaannya.
jbowman
1
NB Kedua Poisson & Poisson nol-dipotong adalah kasus khusus dari distribusi yang telah Anda tetapkan. Satu tidak bersarang di yang lain. Jadi Anda tidak dapat menggunakan teorema Wilks untuk mendapatkan distribusi chi-squared asimptotik untuk dua kali rasio kemungkinan log, mana yang Anda anggap sebagai hipotesis nol. (Saya pikir ada beberapa kondisi keteraturan untuk tes Vuong juga.)
Scortchi - Reinstate Monica
3
@Scortchi Saya ingin tahu tentang definisi "bersarang" yang Anda terapkan. Walaupun saya tidak setuju dengan kesimpulan Anda, saya sampai pada sudut pandang yang sedikit berbeda: ya, Poisson bersarang di dalam keluarga ini (karena muncul dengan membatasi ) tetapi berbagai kesimpulan tentang distribusi asimptotik MLE estimasi parameter untuk tidak berlaku karena nilai terletak pada batas keluarga. Apakah saya kehilangan beberapa perbedaan penting? p=0pp
whuber
2
@whuber, saya akan berkomentar / memberikan jawaban tentang hal yang sama. The Link direferensikan melakukan catatan: "... meskipun distribusi chi-kuadrat mungkin perlu beberapa penyesuaian karena pembatasan pada batas ruang parameter"
Ben Bolker