如何摆脱在R中找不到函数do()的错误

Question

我正在做bootstrap采样。我使用了以下代码：

RNGkind(sample.kind = "Rejection")
set.seed(7377)
B = do(1000)*mean(resample(c(rep(1,258),rep(0,352-258)),352))

但得到以下错误:错误在do(1000)：找不到函数"do“

Answer 1

您可能需要replicate.

但是，对于gdata::resample的默认设置replace=FALSE，您的样本将始终具有相同的均值。我建议您在resample命令中添加replace=TRUE：

library(gdata)
set.seed(7377)

B <- replicate(1000, mean(resample(x=c(rep(1,258), rep(0,352-258)), size=352, replace=TRUE)))
B

#[1] 0.6988636 0.7386364 0.7187500 0.7698864 0.7329545 0.7244318