如何在納入檢查時在R中擬合負二項分佈

Question

我需要擬合Y_ij ~ NegBin(m_ij,k)，因此對於計數爲負二項分佈。然而，我觀察到的數據被審查，我知道y_ij的價值，但它可能超過這個價值。書面方式向下對數似然這個問題會是：

ll = \sum_{i=1}^n w_i (c_i log(P(Y_ij=y_ij|X_ij)) + (1- c_i) log(1- \sum_{k=1}^32 P(Y_ij = k|X_ij))) 

凡X_ij表示設計矩陣（與感興趣的協變量），w_i爲每個觀測重量，y_ij是響應變量和P(Y_ij=y_ij|Xij)是負二項式分佈，其中m_ij=exp(X_ij \beta)和\alpha是超分散參數。

有人知道在R中是否存在可用於獲取該內置代碼的內置代碼？

Answer 1

檢查本文出來：Regression Models for Count Data in R