我正在尝试使用以下代码在R中使用负二项式ditribution:
negfit <- fitdistr(Probhist$DemandQuantity,"negative binomial")
coef(negfit)
数据如下所示:
Probhist <- data.frame(
DemandQuantity=rep(
c(0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 13, 14, 17),
c(8632, 2242, 707, 272, 132, 65, 32, 22, 12, 10, 1, 2, 2, 2, 1)
)
)
当我在上面的数据上运行fitdistr时,我收到以下警告:
> negfit <- fitdistr(Probhist$DemandQuantity,"negative binomial")
Warning messages:
1: In densfun(x, parm[1], parm[2], ...) : NaNs produced
2: In densfun(x, parm[1], parm[2], ...) : NaNs produced
3: In densfun(x, parm[1], parm[2], ...) : NaNs produced
4: In densfun(x, parm[1], parm[2], ...) : NaNs produced
5: In densfun(x, parm[1], parm[2], ...) : NaNs produced
> coef(negfit)
size mu
0.4692958 0.4911798
尽管有警告,我确实得到了系数估计值。我无法理解输出系数值是否值得信赖并且可以进一步使用。此外,这些警告背后的原因是什么,我怎么能取消它们?任何帮助将不胜感激。
感谢。