在R中拟合负二项分布时的警告

时间:2013-10-02 04:36:29

标签: r warnings binomial-coefficients

我正在尝试使用以下代码在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 

尽管有警告,我确实得到了系数估计值。我无法理解输出系数值是否值得信赖并且可以进一步使用。此外,这些警告背后的原因是什么,我怎么能取消它们?任何帮助将不胜感激。

感谢。

0 个答案:

没有答案