r - 如何在 R 中最小化此功能？

Question

我正在尝试编写一个公式来确定 a 的值，该值最小化函数输出 myfun(即 a-fptotal ). MWE:

c <- as.matrix(c(.25,.5,.25))
d <- as.matrix(c(10000,12500,15000))
e <- 700
f <- 1.1
tr <- .30

myfun <- function(a) {
   b <- max(a-e,0)
   df <- data.frame(u1=c(c*b*.40),u2=c(c*b*.60))
   df$year <- 1:nrow(df)
   df$factor <- 1/(f)^df$year
   df$d <- d
   df$t <- (df$d-df$u1-df$u2)*tr
   df$f <- df$d-df$t
   df$fp <- df$f*df$factor
   fptotal <- sum(df$fp)
   return(a-fptotal) # what is a good approach to minimize this?
}

我试过这样的事情但没有成功:

o2 <- optimize(myfun, lower = 0, upper = 30000)

正确答案应该在 28355 左右:

myfun(28355)  
## [1] -1.04151

任何想法表示赞赏。谢谢!

Answer 1

我想你想最小化 a-fptotal 的 square ...

ff <- function(x) myfun(x)^2
> optimize(ff,lower=0,upper=30000)
$minimum
[1] 28356.39

$objective
[1] 1.323489e-23

或者找到根(即 myfun(x)==0):

uniroot(myfun,interval=c(0,30000))
$root
[1] 28356.39

$f.root
[1] 1.482476e-08

$iter
[1] 4

$init.it
[1] NA

$estim.prec
[1] 6.103517e-05