カーネル密度推定
# [-h,h]の範囲の一様分布をカーネル関数とする my.rect <- function(x,h){ ret <- rep(0,length(x)) ret[which(abs(x)<=1)] <- 1 ret } hs <- c(0.1,0.25,0.5,1) y <- seq(from=-10,to=10,length=1000) x.k <- matrix(0,length(y),length(hs)) for(i in 1:l…
# [-h,h]の範囲の一様分布をカーネル関数とする my.rect <- function(x,h){ ret <- rep(0,length(x)) ret[which(abs(x)<=1)] <- 1 ret } hs <- c(0.1,0.25,0.5,1) y <- seq(from=-10,to=10,length=1000) x.k <- matrix(0,length(y),length(hs)) for(i in 1:l…