Density plot

도수 table을 이산으로 표시하는 것은 앞서 해 봤습니다. 이번에는 density plot(도수의 전체 합을 1로 했을 경우 각각의 예상 분포)을 그려봅시다. 그린 다음에 아랫부분을 붉게 채워보겠습니다.

> x1 <- rep(c(0,1,2,3,4,5,6),c(7,14,13,8,4,2,2))
> plot(table(x1),xlab='x1',ylab='frequency',main='test',type='h',lwd=10)
> plot(density(x1),xlab='x1',ylab="density",col="red")
> polygon(density(x1),col='red',border='red')

결과는

density()로 밀도를 구하고 이를 plot으로 그립니다. 이후 polygon을 써서 density의 모든 영역을 촘촘히 색칠합니다.
col : 색상
border : 경계 색상