> 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的小白,程序是按照老师给的范本修改出来的,哪里写错了大佬们轻喷