根据 微软的文档,Excel 使用标准的“n 次失败前的绘制次数”定义;使用的参数化fitdistrplus
是替代方法,?dnbinom
称为:
另一种参数化(通常在生态学中使用)是
平均“mu”和“size”,即色散参数,其中“prob”=“size/(size+mu)”。在此参数化中,方差为“mu + mu^2/size”。
因此,如果您想从mu
和size
返回prob
和size
(分别是 Excelprobability_s
和number_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
)。最好使用GAMMLN
Excel 中的函数来计算对数标度上的归一化常数的分子和分母......如果你幸运的话,那里有人会为你省去一些麻烦并已经实现了......