所使用的种子之间实际上存在一种单向关系set.seed()
以及中的信息.Random.seed
; as ?Random.seed
说,@MattTenenbaum 的回答显示,信息.Random.seed
可以保存和恢复。我很欣赏对一个函数的渴望,该函数可以从当前状态派生出一个简单的整数种子.Random.seed
,但在没有它的情况下,您必须保存并恢复完整信息......例如
set.seed(1001)
save(".Random.seed",file="random_state_seed1001.RData") ## save current state
runif(1)
## [1] 0.9856888
runif(1)
## [1] 0.4126285
runif(1)
## [1] 0.4295392
load("random_state_seed1001.RData") ## restore state just after set.seed()
runif(1)
## [1] 0.9856888
正如 @JoshuaUlrich 指出的,只有在保存和恢复之间不修改 RNG 的类型(以及正常偏差生成器的类型)时,这才有效/安全......
更完整的解决方案:
save_rng <- function(savefile=tempfile()) {
if (exists(".Random.seed")) {
oldseed <- get(".Random.seed", .GlobalEnv)
} else stop("don't know how to save before set.seed() or r*** call")
oldRNGkind <- RNGkind()
save("oldseed","oldRNGkind",file=savefile)
invisible(savefile)
}
restore_rng <- function(savefile) {
load(savefile)
do.call("RNGkind",as.list(oldRNGkind)) ## must be first!
assign(".Random.seed", oldseed, .GlobalEnv)
}
试试看:
set.seed(101)
RNGstore <- save_rng() ## save file name
runif(1)
## [1] 0.3721984
runif(10)
## [1] 0.04382482 0.70968402 0.65769040 0.24985572 0.30005483 0.58486663
## [7] 0.33346714 0.62201196 0.54582855 0.87979573
restore_rng(RNGstore)
runif(1)
## [1] 0.3721984
也可以看看:http://www.cookbook-r.com/Numbers/Saving_the_state_of_the_random_number_generator/ http://www.cookbook-r.com/Numbers/Saving_the_state_of_the_random_number_generator/