R语言学习系列(画向量的密度直方图)

最新推荐文章于 2024-07-20 09:40:11 发布

广州接入

最新推荐文章于 2024-07-20 09:40:11 发布

阅读量1.3w

点赞数 1

分类专栏： R语言 4、学习笔记文章标签：语言 plot function 图形 c 解密

本文链接：https://blog.csdn.net/hawksoft/article/details/7581524

版权

本文介绍如何利用R语言的plot和lines函数创建密度直方图，重点在于理解核密度函数及其背后的概率计算。通过公式1/(n*bw)解释了密度值的计算，并指出当样本数量足够多、分段足够细时，自画的密度直方图能逼近真实概率密度函数。

摘要由CSDN通过智能技术生成

利用R语言画密度直方图比较方便，但为了理解密度函数的意义和如何计算密度值，于是用plot,lines两个画图函数来自己实现密度直

方图的画图程序脚本如下：

DrawDensity = function(x,bw=5)
{
   if(any(bw<=0))
   {
     bw <- 5
    
   }
   #print(bw)
   n <- length(x)
   if(any(n>0))
   {
     mn <- min(x)
     mx <- max(x)
     rmn <- (mn %/% bw)
     begx <-  rmn * bw
     rmn <- (mx %/% bw)
     endx <- (rmn) * bw
     if ((mx-rmn * bw) > 0)
     {
        endx <- (rmn+1) * bw
     }
     
     segs <- (endx-begx)  %/% bw
     segcounts <- numeric(segs)
     validsegs <- 0
     print(segs)
     for(i in 1:segs)
     {
        slw <- begx+(i-1) * bw
        sup <- begx+(i) * bw
        tmp<- x[x > slw & x<=sup];
        segcounts[i]<- length(tmp)
        if(length(tmp)>0)
        {
           validsegs <- validsegs+1
        }
     }
     #计算基本密度: