[R시각화] 조건부 커널밀도추정함수(KDE) 그래프 그리는 방법 (densityplot)

lattice 패키지의 densityplot 함수를 이용하면 조건부 커널밀도추정함수를 그릴 수 있습니다. 커널밀도추정은 커널함수를 이용하여 밀도추정을 하는 것인데요. 커널함수와 밀도추정이 무엇인지 간단히 알아봅시다.

 

커널함수 : 원점을 중심, 대칭, 적분값이 1, non-negative인 함수

 

밀도추정 : 데이터들의 분포를 이용하여 변수의 분포를 추정하는 것. 변수의 확률밀도함수를 구하는 것.

 

iris 데이터에 densityplot을 적용해봅시다. 꽃의 '종(species)' 별로 꽃받침 길이(Sepal.Length)의 커널밀도함수를 그려주는 것입니다.

 

library("lattice")
densityplot( ~ Sepal.Width |Species, data=iris)

 

 

설정할 수 있는 옵션들이 많이 있는데 몇가지 주요한 옵션을 알아봅시다. 

 

layout=c(1,3) : 출력레이아웃

bw="nrd0" : bandwidth인데 rull of thumb 에 의해 정해진 값이다. 마치 p값이 0.05인 것처럼. 

 

옵션을 설정해봅시다. 레이아웃 옵션만 설정하겠습니다. 

 

library("lattice")
densityplot( ~ Sepal.Width |Species, data=iris,
             layout=c(3,1))