非参数核密度估计R代码

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

###########################################################


n=1000
h=0.45


#生成data###############################################


xi=rnorm(n, mean=0, sd=1)

##h=1.06*sd(xi)*n^(-1/5)

###############ker#############################################

ker=function(u)
{
kern=0.75*(1-u^2)*(abs(u)<=1)
return(kern)
}


#############密度估计####################################

hatfx=function(x)
{
hf=sum(ker((xi-x)/h))/(n*h)

return(hf)
}
#################plot######################################


ira=30

a=range(xi)[1]

b=range(xi)[2]

x0=seq(a,b,length=ira)

hatf=c()

for(i in 1:ira)
hatf[i]=hatfx(x0[i])


plot(x0,hatf,type="l",col=1,xlim=c(a,b),ylim=c(0,0.5),xlab="x",lwd=2,
ylab="hatf(x)",cex.lab=1.5,cex.axis=1.5)

lines(x0,(1/sqrt(2*pi))*exp(-x0^2/2),lty=3,col=1,lwd=2)

扫码加我 拉你入群

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

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

分享0 收藏2 回帖

关键词：非参数核密度估计 核密度估计 非参数 R代码 核密度