r - 优化求解器函数

Question

我正在尝试设置一个“求解器”函数，以将“gfc”的值优化为零，改变(并查找)下面等式中的变量“fc”。参数已给出。

f0 = 6
f1 = 1
k = 2
ft = 0.3

gfc = ft-((f0-fc)/k)+((f1/k)*ln((fc-f1)/(f0-f1)))

在 Excel 上求解这个函数，我发现 fc=5.504 的值。

Answer 1

您可以使用 uniroot 查找函数等于 0 的位置:

f0 = 6
f1 = 1
k = 2
ft = 0.3

gfc = function(fc) {
    ft - ((f0 - fc) / k) + ((f1 / k) * log((fc - f1) / (f0 - f1)))
}

uniroot(gfc, interval = c(f0, f1))
#> $root
#> [1] 5.504386
#> 
#> $f.root
#> [1] 6.72753e-09
#> 
#> $iter
#> [1] 5
#> 
#> $init.it
#> [1] NA
#> 
#> $estim.prec
#> [1] 6.103516e-05