有没有办法用 R 找到 pbinom 的概率 p

Question

我认为这是一个相当普遍的问题，但我找不到解决方案。

我想解下面的方程： pbinom(18,25,p)=0.05。

有没有办法用程序R找到未知的p？

感谢您的帮助。

Answer 1

根发现：

print(
  res <- uniroot(function(p) pbinom(18,25,p) - 0.05, c(0, 1), 
          tol = .Machine$double.eps)
)

pbinom(18,25,res$root)
#[1] 0.05

Answer 2

蛮力：

p = 0.0001 # starting point
while (abs(pbinom(18,25,p) -  0.05) > 0.001) p <- p + 0.001

此代码针对不同的 p 值评估 pdf，直到您 "close enough" 达到 0.05。这里 "close enough" 表示在 0.001 范围内。

> p
[1] 0.8601
> pbinom(18,25,0.8601)
[1] 0.05070763

Is there a way to find the probability p with R for pbinom