我有一些计数数据,有点像泊松分布,但分布过大。我在python中使用统计模型拟合了一个负二项式glm模型,并为alpha选择了一个值(即分散参数)。然后,我将该模型应用于测试集,以进行一些预测。我现在想计算累积分布函数,这样我就可以计算随机变量X<=9给出了我的输出预测μ(即7.8)和阿尔法的预定值(即0.2)。
scipy文档( https://docs.scipy.org/doc/scipy/reference/generated/scipy.stats.nbinom.html )建议我可以用nbinom.cdf(k,n,p)在python中计算这一点,但当我只有mu和alpha时,我如何获得n和p的值?
p和n与alpha和mu有关,您提供的文档中给出了这种关系:
sigma_squared = mu + alpha * mu**2 p = mu / sigma_squared n = mu**2 / (sigma_squared - mu)
其相当于:
p = 1 / (1 + alpha * mu) n = 1 / alpha