從核密度生成隨機樣本時的不同結果

Question

library(ks) 

x<-rnorm(1000) 
hist(x, col="red") 

y <- rkde(kde(x), n=1000) 
hist(y, col="green") 

y <- rkde(density(x), n=1000) 
hist(y, col="blue") 

最後的直方圖是錯誤的。我之前使用過density，我發現它對於更復雜的分佈是準確的。爲什麼在這種情況下表現如此糟糕？謝謝

Answer 1

因爲你使用的功能錯了。 rkde需要類kde的一個對象。 density不會返回kde對象，並且結構不同。

這就像告訴某人拍攝他們的手槍並遞給他們獵槍彈，然後想知道爲什麼當他們開槍時並沒有真正做任何事情。