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

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

1 个答案: