求助R语言NW核回归，CV选带宽，程序出错！-经管之家官网！

source("d:/R脚本/cvbandNW.R")

K=function(x)

0.75*(1-x^2)*(abs(x)<=1)

m=function(X)

0.5*exp(X^2)

k=20

grid=rep(0,k)

a=-1;b=1

for(i in 1:k)

grid=a+(i-1)*(b-a)/(k-1)

Mestimated.NWE1=rep(0,k)

n=100

M=500

rmse1=rep(0,M)

for(j in 1:M){

X=runif(n,min=a,max=b)

e=rep(0,n)

for(i in 1:n)

e=rnorm(1,0,2*abs(X))

Y=m(X)+e

h1=cvbandNW(X,Y)$h

estimated.NWE1=rep(0,k)

for(i in 1:k){

X1=X[abs(X-grid)<=h1]

Y1=Y[abs(X-grid)<=h1]

Kvec=K((X1-grid))/h1)

lengh(Y1)=length(Kvec)

S=sum(Kvec)

if(S==0) S=S+0.0001

estimated.NWE1=sum(Kvec*Y1)/S

}

Mestimated.NWE1= Mestimated.NWE1+ estimated.NWE1

rmse1[j]=sqrt(mean((estimated.NWE1-m(grid))^2))

Estimated.NWE1=Mestimated.NWE1/M

op=par(mfrow=c(2,2))

plot(grid,m(grid),type=”l”,xlab=”x”,ylab=”m(x)”),xlim=(-1,1),ylim=(.4,1.4))

lines(grid,estimated.NWE1,lty=2,col=2)

title(sub=”(a)”)

boxplot(remse1,xlab=””,ylab=”RMSEs”,range=3.5,names=”expression(REMSE[1])”)

title(sub=”(b)”)

扫码或添加微信号：坛友素质互助

「经管之家」APP：经管人学习、答疑、交友，就上经管之家！
免流量费下载资料----在经管之家app可以下载论坛上的所有资源，并且不额外收取下载高峰期的论坛币。
涵盖所有经管领域的优秀内容----覆盖经济、管理、金融投资、计量统计、数据分析、国贸、财会等专业的学习宝库，各类资料应有尽有。
来自五湖四海的经管达人----已经有上千万的经管人来到这里，你可以找到任何学科方向、有共同话题的朋友。
经管之家（原人大经济论坛），跨越高校的围墙，带你走进经管知识的新世界。
扫描下方二维码下载并注册APP