从R密度图获得计数数据

Question

是否有方法将density（）函数返回的密度转换为计数？我尝试过这样的事情：

dens = density(x)
plot(dens$x,dens$y*dens$n*(dens$x[2]-dens$x[1]))

但是，我不确定这个数学是否正确，因为我完全不知道密度函数下发生了什么。

Answer 1

计数将与输入点的数量相加，因此只需将dens$y重新缩放为等于输入向量的长度：

> x=rnorm(1000,5)
> 
> dens=density(x)
> plot(dens$x,dens$n*dens$y/sum(dens$y))

我认为你乘以binwidth的方法相当于我的除以总和。并且可能更快，但确实依赖于常规的bin大小。