Ниже приведены функции для начальной загрузки, но как я могу сделать воспроизводимый результат? Я пробовал 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)