以下是引导程序的功能,但是如何使结果可再现?我尝试了set.seed(),但是那不起作用,因为每次lapply调用函数boot.lm.vector时,该函数仅生成一个模拟集并计算一次系数。 R中有什么功能可以像种子列表一样工作吗?或任何其他方法来使结果可重复?
boot.lm.vector <- function(index, inputData) {
d <- inputData[sample.int(nrow(inputData), replace = T),]
a <- ncol(inputData)-1
X <- d[, 1:a]
y <- d[, a+1]
solve(crossprod(X), crossprod(X,y))
}
rtest <- lapply(1:10000, fun = boot.lm.vector, inputData = boot_set)
rtestdf <- plyr::ldply(rtest)
答案 0 :(得分:0)
如果您在函数中使用{{ function }}
设置updated
,则应该能够重现它。下面的虚拟updated
功能:
seed