如何从估计的密度函数上随机取样？

随便给个向量，如下，用density()方程求这个向量元素的概率密度函数，即根据已知有限几个数据估计一个适用于这个规律的分布。想从这个估计的分布上随机取值怎么取？
如下面程序，A即是估计的密度函数，假如想从这个密度函数上随机取10000个数，如何操作？

Data = c(2,3,8)
A=density(Data, kernel = "gaussian", n=512)

如果A是个数据，那就直接 sample(A, 10000, replace=FALSE), 但是A是函数不是数据集所以行不通啊。

1. Data = c(2,3,8)
   
2. A=density(Data, kernel = "gaussian", n=512)
   
3. n <- 10000
   
4. sam <- sample(Data, n, replace=TRUE) + rnorm(n, 0, A$bw)
   
5. hist(sam, prob=T)
   
6. lines(A)

复制代码

刚才试了下，可以啊

qoiqpwqr 发表于 2014-10-9 09:38

qoiqpwqr 发表于 2014-10-8 22:19

多谢！

难道 "rnorm(n, 0, A$bw)" 就是对应KDE的 kernel 部分？

sty1882 发表于 2014-10-9 19:55
多谢！

难道 "rnorm(n, 0, A$bw)" 就是对应KDE的 kernel 部分？

对，因为你用的是gaussian kernel

qoiqpwqr 发表于 2014-10-9 19:58
对，因为你用的是gaussian kernel

多谢！