Есть ли способ найти вероятность p с R для pbinom

Question

Я думаю, что это довольно распространенная проблема, но я не смог найти решения.

Я хочу решить следующее уравнение: pbinom (18,25, p) = 0,05.

Есть ли способ найти неизвестный p с помощью программы R?

Любая помощь приветствуется.

Answer 1

Root вывод:

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

Этот код оценивает pdf для различных значений p, пока вы не «достаточно близко» к 0,05. Здесь "достаточно близко" означает диапазон 0,001.

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