2013-12-17 50 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大小。