如何在纳入检查时在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