2

我应用了该fitdistrplus软件包以适应经验分布。事实证明,最好的拟合是带参数的负二项分布:

size=0.6900788
mu=2.6522087

dnbinom(0:10, mu = 2.6522087, size =0.6900788)
[1] 0.33666338 0.18435650 0.12362301 0.08796440 0.06439416 0.04793144 0.03607044 0.02735574 0.02086667 0.01598815 0.01229390

我现在尝试在 EXCEL 上生成相同的数字,其中需要不同格式的参数:

NEGBINOMDIST(number_f,number_s,probability_s)

我该怎么做?有任何想法吗?非常感谢..

4

1 回答 1

3

根据 微软的文档,Excel 使用标准的“n 次失败前的绘制次数”定义;使用的参数化fitdistrplus是替代方法,?dnbinom称为:

另一种参数化(通常在生态学中使用)是 平均“mu”和“size”,即色散参数,其中“prob”=“size/(size+mu)”。在此参数化中,方差为“mu + mu^2/size”。

因此,如果您想从musize返回probsize(分别是 Excelprobability_snumber_s),您需要

number_s=size
probability_s=size/(size+mu)

muval <- 2.6522087
sizeval <- 0.6900788
(probval <- sizeval/(sizeval+muval))
## [1] 0.206469
all.equal(dnbinom(0:10,mu=muval,size=sizeval),
          dnbinom(0:10,prob=probval,size=sizeval))
## TRUE

但是,您还没有完成,因为(正如上面@James 评论的那样)Excel 只允许正整数number_s,并且上面的估计值为 0.69。您可能需要在与 Excel 相关的论坛上搜索/询问如何克服此限制......在最坏的情况下,由于 Excel 确实具有gamma 函数的实现,您可以使用给出的公式?dnbinom

Gamma(x+n)/(Gamma(n) x!) p^n (1-p)^x 

实现您自己的 NB 计算(此公式允许 的非整数值n)。最好使用GAMMLNExcel 中的函数来计算对数标度上的归一化常数的分子和分母......如果你幸运的话,那里有人会为你省去一些麻烦并已经实现了......

于 2013-05-28T13:40:28.143 回答