Question

我是R新手。

让我说我有代码行

set.seed(123456)

然后我想得到值123456所以我可以打印出来用于文档目的，如果需要在将来的某个时间重新输入值。那我怎么得到那粒子呢？

请注意，我可能会将上面的代码行注释掉，因此我实际上不会知道种子在123456处是什么。因此我需要r将当前种子的位置打印为单个整数，而不是626个整数的列表。

Answer 1

set.seed()中使用的种子与.Random.seed中的信息之间存在单向关系;如?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/

Answer 2

调用set.seed后，您应该可以引用.Random.seed（有关详细信息，请参阅http://stat.ethz.ch/R-manual/R-patched/library/base/html/Random.html）。

一个简单的例子：

set.seed(123)
tmp <- .Random.seed
> runif(10)
 [1] 0.2875775 0.7883051 0.4089769 0.8830174 0.9404673 0.0455565 0.5281055 0.8924190 0.5514350 0.4566147
.Random.seed <- tmp
> runif(10)
 [1] 0.2875775 0.7883051 0.4089769 0.8830174 0.9404673 0.0455565 0.5281055 0.8924190 0.5514350 0.4566147

Answer 3

这些在实践中应该可以正常工作（但不应该用于加密）：

new.seed <- as.integer(runif(1)*2e9)
cat("Random seed: ", new.seed, "\n")
set.seed(new.seed)

你需要实际设置种子以使其工作，所以它不是原始问题的答案，但可能是你正在寻找的，因为它可以让你有一个容易指定的随机种子，而无需手动设置带有预先指定的整数。

打印当前随机种子，以便稍后可以使用set.seed（）输入它

3 个答案: