当前位置:   article > 正文

核密度图

核密度图左截断右拖尾
  • 核密度估计

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

  1. plot(density(x))
  2. #x是数值型向量

  • 核密度图
  1. #如下图
  2. d <- density(mtcars$mpg)
  3. plot(d)
  4. d <- density(mtcars$mpg)
  5. plot(d, main="Kernel Density of Miles Per Gallon") #添加标题
  6. polygon(d, col="red", border="blue") #ploygon()函数根据顶点的 x 与 y坐标(本例中由density()函数提供)绘制多边形
  7. rug(mtcars$mpg, col="brown") #添加棕色的轴须线

 

 
 
  • 可比较的核密度图
  1. #可比较的和密度图
  2. library(sm)
  3. attach(mtcars)
  4. cyl.f <-factor(cyl,levels = c(2,4,6),labels=c("4 cylinder","6 cylinder","8 cylinder"))
  5. sm.density.compare(mpg,cyl,xlab ="xxx") #绘制密度图
  6. title("yyy") #主标题
  7. colfill <- c(2:1+length(levels(cyl.f))) #得到值c(2,3,4)
  8. legend("topright",levels(cyl.f),fill=colfill) #locator(1)表示用鼠标点击想让图例出现的位置来交互式地放置图例
  9. #我修改为 topright
  10. #fill=colfill表示为每个水平指定了一种颜色

06113724_nMWD.png

从上图可以看出,核密度图的叠加也不失为一种在某个结果变量上跨组比较观测的强大方法,可以看到不同组所含值的分布状态,以及不同组之间的重叠程度

 

 

转载于:https://my.oschina.net/u/1785519/blog/1563079

声明:本文内容由网友自发贡献,不代表【wpsshop博客】立场,版权归原作者所有,本站不承担相应法律责任。如您发现有侵权的内容,请联系我们。转载请注明出处:https://www.wpsshop.cn/w/花生_TL007/article/detail/131998
推荐阅读
相关标签
  

闽ICP备14008679号