楼主: 疏楼梦妖
879 1

[问答] 有大佬知道这个程序错在哪里么 [推广有奖]

  • 0关注
  • 0粉丝

高中生

55%

还不是VIP/贵宾

-

威望
0
论坛币
661 个
通用积分
0.4500
学术水平
0 点
热心指数
0 点
信用等级
0 点
经验
2010 点
帖子
7
精华
0
在线时间
39 小时
注册时间
2017-3-31
最后登录
2024-1-12

+2 论坛币
k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群
赵安豆老师微信:zhaoandou666

经管之家联合CDA

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

感谢您参与论坛问题回答

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

+2 论坛币
> bw <- function(para,lb,ub){
+ set.seed(10000)
+ mydata <<- rnorm(1000)
+ fit = nlminb(start = para, objective= llh, lower=lb, upper=ub)
+ fit$par
+ }
> epkernel <- function(z){(3/4)*(1-z^2)*(abs(z)<=1)}
> llh <- function(h){
+ tt <- length(mydata)
+ fhat <- sapply(mydata,function(x){(1/((tt-1)*h))*(sum(epkernel((mydata-x)/h))-epkernel(0))})-sum(log(fhat[fhat!=0]))
+ }
> optbw = bw(1,1,1)
Error in objective(.par, ...) : 找不到对象'fhat'
>
一路都没问题,最后说找不到对象fhat,但是这个fhat不是就在它上面一点点么
楼主是正在学R的小白,程序是按照老师给的范本修改出来的,哪里写错了大佬们轻喷

二维码

扫码加我 拉你入群

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

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


找不到“对象”系列[titter

使用道具

您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注cda
拉您进交流群

京ICP备16021002-2号 京B2-20170662号 京公网安备 11010802022788号 论坛法律顾问:王进律师 知识产权保护声明   免责及隐私声明

GMT+8, 2024-4-23 14:33