R语言绘制 核密度图

核密度估计是用于估计随机变
量概率密度函数的一种非参数方法。虽然其数学细节已经超出了本书的范畴，但从总体上讲，核
密度图不失为一种用来观察连续型变量分布的有效方法。绘制密度图的方法（不叠加到另一幅图
上方）为：
plot(density(x))
其中的x是一个数值型向量。由于plot()函数会创建一幅新的图形，所以要向一幅已经存在的图
形上叠加一条密度曲线，可以使用lines()函数

par(mfrow=c(2,1)) 
 d <- density(mtcars$mpg) 
 plot(d) 
 d <- density(mtcars$mpg)
 plot(d, main="Kernel Density of Miles Per Gallon")
 polygon(d, col="red", border="blue")
 rug(mtcars$mpg, col="brown")

polygon()函数根据顶点的x和y坐标（本例中由density()函数提供）绘制了多边形。
核密度图可用于比较组间差异。可能是由于普遍缺乏方便好用的软件，这种方法其实完全没