1 R 패키지 설치하고 불러오기

if (! ("devtools" %in% rownames(installed.packages()))) 
  { install.packages("devtools") }
base::require("devtools")
if (! ("roxygen2" %in% rownames(installed.packages()))) 
  { install.packages("roxygen2") }
base::require("roxygen2")

2 github에서 패키지 설치하고 부르기

github에서 설치하는 kde 패키지는 커널 함수를 추정하고 그래프를 그리는 \[ \hat{f}(x;h) = n^{-1} \displaystyle\sum_{i=1}^{n}K_h(x-X_i) \] 이다. 여기서 \(K_h(u) = h^{-1}K(u/h)\) 이고 \(K\) 커널 함수는 \(\int K(x) dx=1\)을 만족한다. 또한 \(h\)는 양수이며 광대역(bandwidth) 또는 창 너비(window width)로 부른다.

install_github(repo="hmkang98/kde")
library(kde) 

3 예제

example(kde)
## 
## kde> Xi <- c(3,4.5,5,8,9)
## 
## kde> x <- seq(min(Xi)-3, max(Xi)+3,by=0.1)
## 
## kde> sx <- seq(-3,3,by=0.1)
## 
## kde> haty <- kde(x,Xi)
## 
## kde> plot(x,haty,type="l",ylim=c(0,max(haty)))

## 
## kde> for(i in 1:length(Xi)){
## kde+    lines(sx+Xi[i], dnorm(sx)/length(Xi),lty=2) 
## kde+ }