给定 R 中所需的答案，如何找到函数参数的最佳值？

Question

我想知道如何找到 t 的最精确值，因为我希望下面整个函数调用的答案在 R 中变为 3？

if(!require(BayesFactor)){install.packages('BayesFactor')} ; library(BayesFactor)

# 3 = function call below:
unname(ttest.tstat( t, 20, 20, rscale = sqrt(2)/2, simple = TRUE)[[1]])

Answer 1

这是一个使用 optim 的解决方案。 par 是使您最接近 3 的 t 值。其他方法可能会使您更接近，但那是在我的驾驶室之外。

library(BayesFactor)

f1 <- function(t){
  abs(ttest.tstat( t, 20, 20, rscale = sqrt(2)/2, simple = TRUE)[[1]]-3)
}

f2 <- function(t){
  ttest.tstat( t, 20, 20, rscale = sqrt(2)/2, simple = TRUE)[[1]]
}

optim(par=-1, fn = f1)

#$par
#[1] -2.440637
#
#$value
#[1] 1.140988e-07
#
#$counts
#function gradient 
#      56       NA 
#
#$convergence
#[1] 0
# 
#$message
#NULL

编辑：

考虑到 Nelder-Mead 对于一维优化不可靠的警告，优化 t 值的其他两种方法。警告只是意味着有/有更好的方法来优化单变量函数。

optim(par=-1, fn = f1, method="Brent", lower = -100, upper=100)
optimize(f1, interval = c(-100,100))