从核密度生成随机样本时的结果不同

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对象，结构不同。

这就像告诉别人射击他们的手枪并递给他们的猎枪弹，然后想知道为什么当他们开枪时它并没有真正做任何事情。