我使用excel目标搜索功能计算了所需的安全库存。下面是图像。 但是现在我想用R来做同样的事情。
当我进入SafetyStock&amp ;;时,以下功能给了我相同的excel结果。 SD。现在我需要进行反向计算(每当我提供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)
}
答案 0 :(得分:0)
为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