来自离散命名向量的R-概率密度

Question

我有一个命名向量modVect，它包含每个名称的概率值（-.01，-.02 ....这些是离散的错误值）。

head(modVect)
-.01         -.02         -.03         -.04         -.05         -.06 
0.0006927649 0.0033267561 0.0080302663 0.0065196742 0.0018305703 0.0003929703 

我想根据这些概率模拟错误值。如何使用概率密度函数或核密度估计来做到这一点？到目前为止，我已经完成了这个

xseq<-as.numeric(names(modVect))
plot(xseq, modVect)
y<-kde(modVect)

但这是不正确的，因为它不会将modVect值解释为概率。我希望有一个函数可以根据概率创建一个pdf。

由于

Answer 1

要使用probs进行抽样，请将sample与probs一起使用。

N <- 1000
sample(xseq, N, replace=TRUE, probs=modVect)