如何使样本和lapply的组合可重现?

时间:2018-10-31 00:36:57

标签: r lapply

以下是引导程序的功能,但是如何使结果可再现?我尝试了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)

1 个答案:

答案 0 :(得分:0)

如果您在函数中使用{{ function }}设置updated,则应该能够重现它。下面的虚拟updated 功能

seed