1.构造的NW估计(写成函数形式);
> X<-runif(100)
> epsilon<-rnorm(100,0,1)
> Y<-sin(2*pi*X)+epsilon
>h=0.08
>kh<- function(x){
suma <- 0
for(i in 1:100){
if((abs(x-X)/h)<=1
suma <- suma+1/2/h
}
return(suma)
}
>fh<- function(x){
sumb <- 0
for(i in 1:100){
if(abs((x-X)/h)<=1)
sumb <- sumb+1/2*Y/h
}
return(sumb)}
> mh<-function(x){fh(x)/kh(x)}
>f=function(x) mh(x)
>curve(f,0,1)
程序可以输出结果,但是无法画图,提示条件长度大于1,只能用其第一个元素,不知道应该怎么改