小编use*_*rLx的帖子

我正在使用optim()(和optimize()) 尝试查找二项式分布的分位数，但是对于 N ~ 2000 (N = 2135)，函数不会给出正确的值。

optim(21, function(x) abs(1 - pbinom(x, 2135, 21/2135) - 0.1), 
      method = "Brent", lower = 1, upper = 2135) 

optimize(function(x) abs(1 - pbinom(x, 2135, 21/2135) - 0.1), c(1,2135))

PS：我也尝试将min参数设置为等于概率，但我仍然得到错误的答案。