如何在合并检查的同时拟合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

查阅本文： R中计数数据的回归模型