- 核密度估计
用于估计随机变量概率密函数的一种非参数方法
- plot(density(x))
- #x是数值型向量
例
- 核密度图
- #如下图
- 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") #ploygon()函数根据顶点的 x 与 y坐标(本例中由density()函数提供)绘制多边形
- rug(mtcars$mpg, col="brown") #添加棕色的轴须线
- 可比较的核密度图
- #可比较的和密度图
- library(sm)
- attach(mtcars)
- cyl.f <-factor(cyl,levels = c(2,4,6),labels=c("4 cylinder","6 cylinder","8 cylinder"))
- sm.density.compare(mpg,cyl,xlab ="xxx") #绘制密度图
- title("yyy") #主标题
-
- colfill <- c(2:1+length(levels(cyl.f))) #得到值c(2,3,4)
- legend("topright",levels(cyl.f),fill=colfill) #locator(1)表示用鼠标点击想让图例出现的位置来交互式地放置图例
- #我修改为 topright
- #fill=colfill表示为每个水平指定了一种颜色
从上图可以看出,核密度图的叠加也不失为一种在某个结果变量上跨组比较观测的强大方法,可以看到不同组所含值的分布状态,以及不同组之间的重叠程度