Bir Poisson dağılımını takip etmeyi bekleyebileceğim bir veri setim var, ancak yaklaşık 3 kat fazla dağılmış durumda. Şu anda, bu overdispersion R aşağıdaki kod gibi bir şey kullanarak modelleme.
## assuming a median value of 1500
med = 1500
rawdist = rpois(1000000,med)
oDdist = rawDist + ((rawDist-med)*3)
Görme, bu ampirik verilerime çok iyi uyuyor gibi görünüyor. Eğer uyumdan memnunsam , burada açıklandığı gibi negatif bir binom dağılımı kullanmak gibi daha karmaşık bir şey yapmamın bir nedeni var mı? (Öyleyse, bunu yapan herhangi bir işaretçi veya bağlantı çok takdir edilecektir).
Oh, ve bunun biraz pürüzlü bir dağılım yarattığının farkındayım (üçe çarpma nedeniyle), ancak bu benim uygulama için önemli olmamalı.
Güncelleme: Bu soruyu araştıran ve bulan başka herkes uğruna, negatif bir binom dağılımı kullanarak aşırı dağılmış bir poisson modellemek için basit bir R işlevi. D'yi istenen ortalama / varyans oranına ayarlayın:
rpois.od<-function (n, lambda,d=1) {
if (d==1)
rpois(n, lambda)
else
rnbinom(n, size=(lambda/(d-1)), mu=lambda)
}
(R posta listesinden: https://stat.ethz.ch/pipermail/r-help/2002-Haziran/022425.html )