任何用于加权引导程序的 R 包？

Question

在这个领域很新，因此欢迎任何建议。 假设我有一个关于 n=200 个样本的数据。我必须给每个样本赋予权重。

假设，一开始我生成数据y=rnorm(200)，权重为w=rexp(200)。此权重需要应用于我的原始样本 (y)。我必须这样做 1000 次并重新生成 1000 个数据（引导程序）。我怎么能用简单的 R 函数来做呢？有包裹吗？或者如何手动完成？ 我知道我可以在 R 中使用带有“replace=TRUE”的“sample”执行非参数引导程序。但我不知道如何进行加权引导程序。

Answer 1

这样的事情应该有效：

y <- rnorm(200)
w <- rexp(200)
w <- w/sum(w)
y_samp <- sample(y, length(y), replace=TRUE, prob = w)

prob 参数允许您对采样值应用权重。