对单个变量的fminsearch

Question

将R在fminsearch上的帮助页面示例作为起点：

# Rosenbrock function
rosena <- function(x, a) 100*(x[2]-x[1]^2)^2 + (a-x[1])^2  # min: (a, a^2)

fminsearch(rosena, c(-1.2, 1), a = sqrt(2))
# x = (1.414214 2.000010) , fval = 1.239435e-11

我想评估类似的情况，但是只有一个变量，例如：

 rosena <- function(x, a) 100*(x[1]-x[1]^2)^2 + (a-x[1])^2 

但当我跑

fminsearch(rosena, c(1), a = sqrt(2))

它给出了X2:d1中的错误:错误的维数

fminsearch似乎想要一个长度大于或等于2的向量，但是对于这个例子来说，向量需要长度1。

注意: fminsearch是在"pracma“软件包中。

Answer 1

它看起来像是pracma包中的一个bug。

anms函数正在向下标的相关摘录中删除维度：

d <- length(x0) # i.e. 1
d1 <- d + 1 # i.e. 2
...
X <- matrix(0, nrow = d1, ncol = d)
...
X <- X[o, ] # could put drop = FALSE here

我认为你应该向软件包的作者发布一个bug。