2013-12-17 19 views
1

是否有将密度()函数返回的密度转换为计数的方法?我试着做这样的事情:从R密度图获取计数数据

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

但是,我不知道这个数学正确的,因为我并不完全了解的密度函数发生的事。

+2

你是指什么意思?你想要一个'hist'ogram吗? – James

回答

2

计数会总结到的输入点的数量,所以只是重新调整dens$y等于你输入向量的长度:

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

我觉得你的binwidth相乘的方法等同给我的除以和。而且可能更快,但是依赖于一个常规的bin大小。