如何在R中找到多变量函数的最大值

Question

我有一个共同的可能性，我需要最大化，这由大约25个不同的变量决定。我希望有一种方法可以找到最大化此功能的变量值组合，而不需要求助于25＆＃39;循环遍历每个可能的值。

这是一个可能性小得多的例子，使用＆＃39;来最大化。循环方法，其中＆＃39; temp＆＃39;计算似然的值，test是一个记录所用变量值的向量：

h1=23
h2=31

test=c(NA,NA,NA,NA,0)
for (N in seq(60,150,10)) {
  for (p1 in seq(0.01,1,.01)) {
    for (p2 in seq(0.01,1,.01)) {
      for (S1 in seq(0.5,1,.005)) {
        temp=factorial(N)/(factorial(h1)*factorial(h2)*factorial(N-h1-h2))*(p1)^h1*((1-p1)*S1*p2)^h2*(1-(p1+(1-p1)*S1*p2))^(N-h1-h2)
        if (temp>test[5]) test=c(p1,p2,S1,N,temp)
      }
    }
  }
}
test

提前致谢， 塞吉