Rでカーネル密度推定

基本はdensity関数を使えば
よい。

細かい調整が必要な場合
下のように実装できる。

矩形カーネル

#基本はこれでＯＫ．
plot(density(data))


#カーネル関数
k<-function(x){

	0.5*(abs(x)<=1)
}
#刻み幅
h=0.1
#密度関数
f<-function(x,h,data){
	X<-0
	for(i in 1:length(x)){	
		#X[i]<-1/h/length(data)*sum(dnorm(sd=1,mean=0,(x[i]-data)/h))	


		X[i]<-1/h/length(data)*sum(k((x[i]-data)/h))	
	}
	X
}


#curve( f(x));