核密度图-CSDN博客

2019独角兽企业重金招聘Python工程师标准>>>

核密度估计

用于估计随机变量概率密函数的一种非参数方法

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表示为每个水平指定了一种颜色