我需要拟合，即负二项分布来计算数据。但是，我观察到的数据是经过审查的——我知道的值，但它可能超过该值。对数似然是$Y_{ij} \sim NegBin(m_{ij},k)$$y_{ij}$

$$\begin{equation} 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}))) \end{equation}$$

其中表示设计矩阵（带有感兴趣的协变量），是每个观察的权重，是响应变量，是负二项分布，其中和是过度分散参数。$X_{ij}$$w_i$$y_{ij}$$P(Y_{ij}=y_{ij}|X_{ij})$$m_{ij}=exp(X_{ij} \beta)$$\alpha$

有谁知道解决这个问题的 R 包？