求助，如何在脚本中设计函数并求出极大值点？

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

才开始学R语言，因为是作业，写出的脚本总是报错，求高手指点。
问题是 如何在函数L中使用k函数，同时让L函数返回L0（上一步结果sumup乘上含h分式的最后表达式）
如何找到L的极大值点h0（L的极大值是上限不确定的正数）
多谢各位了。
set.seed(0)
x<-c(rnorm(1000))*100
k<-function(x){
  k0<-0.75*(1-x^2)
  return(k0)
}
p<-function(h){
  sum<-0
  pro<-1
  for(t in 1:1000){
    for(s in 1:1000){
      if(abs((x-x[t])/h)>1|s==t)
        next                                  //此处报错：“the condition has length > 1 and only the first element will be used”
        sum<-sum+k((x-x[t])/h)
    }
  pro<-pro*sum
  sum<-0
  }  
  return(pro)
}
L<-function(h){
  L0<-p(h)*((1/(999*h))^1000)
  return(L0)
}
curve(L,from=0,to=10000,n=101)

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：seed 高手指点 求高手 正整数 表达式