合并截尾时如何拟合R中的负二项分布_R

合并截尾时如何拟合R中的负二项分布

合并截尾时如何拟合R中的负二项分布,r,R,我需要拟合Y_ij~NegBinm_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是响应变量，PY_ij=y_ij | Xij是负二项

我需要拟合Y_ij~NegBinm_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是响应变量，PY_ij=y_ij | Xij是负二项分布，其中m_ij=expX_ij\beta和\alpha是过度分散参数

有人知道R中是否存在可用于获取此信息的内置代码吗？

请查看这篇文章：

您的问题是相当统计性的，因此您可以尝试在新的统计信息交换网站上提问