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目标搜寻

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