求助neldermead包里neldermead函数的用法

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

大家好，我在看一段代码的时间

1. # The Nelder-Mead simplex method
   
2. f <- function(x, y) ((x - y) ^ 2 + (x - 2) ^ 2 + (y - 3) ^ 4) / 10
   
3. x <- seq(0, 5, len=20)
   
4. y <- seq(0, 5, len=20)
   
5. z <- outer(x, y, f)
   
6. contour(x, y, z)
   
7. 
   
8. library(neldermead)
   
9. x <- matrix(c(0, 0, 2, 0, 2, 0), 3, 2)
   
10. polygon(x)
    
11. for (i in 1:9) {
    
12.   cat(i,":")
    
13.   x <- neldermead(x,f)
    
14.   polygon(x)
    
15.   text(rbind(apply(x, 2, mean)), labels=i)
    
16. }
    
17. x

复制代码

运行不出来，显示错误如下：
Error: unknownValueForOption: Unknown value 0 for optbase optionunknownValueForOption: Unknown value 0 for optbase optionunknownValueForOption: Unknown value 2 for optbase optionunknownValueForOption: Unknown value 0 for optbase optionunknownValueForOption: Unknown value 2 for optbase optionunknownValueForOption: Unknown value 0 for optbase option
In addition: Warning message:
In if (is(value) != "optimbase") { :
  the condition has length > 1 and only the first element will be used

没有用过这种算法，不知道怎么回事，麻烦高人解释下.
这段代码来自于A First Course in Statistical Programming with R （2nd Edition）的7.3节最后一部分。
谢谢！

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：Elder Mead der EAD ERM