Excel中的求解器

我一直在试图find类似于Excel中的excel求解器。

例如，我有以下程序：

 L=2.4 UCL=L; LCL=-L ARL=1/(pnorm(LCL, mean = 0, sd = 1, lower.tail = TRUE)+pnorm(UCL, mean = 0, sd = 1, lower.tail = FALSE)) 

我想要search一个“L”的值，使得我的结果“ARL”接近例如200。

有没有办法在r？

Answer 1

 f <- function(L) { UCL=L; LCL=-L ARL=1/(pnorm(LCL, mean = 0, sd = 1, lower.tail = TRUE)+pnorm(UCL, mean = 0, sd = 1, lower.tail = FALSE)) ARL-200 } uniroot(f,c(1,10)) 

它会返回

 $root [1] 2.807 $f.root [1] -0.0003725 $iter [1] 13 $init.it [1] NA $estim.prec [1] 6.532e-05 

所以解决scheme是L = 2.807，它需要13次迭代。