0

我有一个命名向量 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。

谢谢

4

1 回答 1

0

要使用 probs 进行采样,请使用samplewith probs

N <- 1000
sample(xseq, N, replace=TRUE, probs=modVect)
于 2013-07-30T04:31:26.023 回答