我需要适合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中是否存在可用于获取此内容的内置代码吗？