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