R中的Excel目标搜索

Question

我已经使用 excel 目标搜索功能计算了所需的安全库存。下面是图片。 但现在我想用 R 做同样的事情。

当我输入 SafetyStock & SD 时，下面的函数给出了相同的 Excel 结果。现在我需要进行反向计算（每当我提供 x 和 SD 时，我都需要 SS）。有人可以帮助我吗？

我尝试了 Optix 和其他类似的 R 包，但没有成功。

opt<-function(SS,SD){

  x=-SS*(1-pnorm(SS/SD)) + SD * dnorm(SS/SD,mean=0,sd =1,0)
  print(x)

}

Excel Goal seek

Answer 1

为 x 求解 f(x)=c 与求解 f(x)-c=0 相同。您可以使用uniroot 查找根目录：

f <- function(SS, SD, ESC) {
  -SS*(1-pnorm(SS/SD)) + SD * dnorm(SS/SD,mean=0,sd =1,0) - ESC
}

zero <- uniroot(f,c(0,1000),SD=600,ESC=39.3)
zero$root

第二个参数是要搜索的区间：0 到 1000 之间。返回

674.0586 

zero 结构有更多有趣的信息：

$root
[1] 674.0586

$f.root
[1] 1.933248e-08

$iter
[1] 8

$init.it
[1] NA

$estim.prec
[1] 6.103516e-05