核密度图也是一种用来观察连续型变量分布的有效方法。
考虑mtcars数据集的mpg变量
> plot(density(mtcars$mpg))
从图中我们可以更详细地观察mpg变量的分布情况。
核密度图可用于比较组间的差异。如果想要比较cyl=4,6,8这三种不同组的汽车mpg变量分布的差异。可以使用sm包。
> library(sm)
> cyl.f <- factor(mtcars$cyl,levels=c(4,6,8),labels=c("4 cyl","6 cyl","8 cyl"))
> sm.density.compare(mtcars$mpg,mtcars$cyl)
> legend(locator(1),levels(cyl.f),fill=c(2,3,4))
从图中,可以明显看出,根据cyl数量分成的三组汽车的mpg变量分布差异。