在R中本地(非全局)设置种子

时间:2013-01-14 18:08:52

标签: r global local seeding

我想在R中仅在本地(内部函数)设置种子,但似乎R不仅在本地设置种子,而且在全局设置种子。这是我正在尝试(不)做的一个简单示例。

myfunction <- function () {
  set.seed(2)
}

# now, whenever I run the two commands below I'll get the same answer
myfunction()
runif(1)

所以,我的问题是:为什么R在全局设置种子,而不仅仅是在我的函数中?我怎样才能让R只在我的函数中设置种子?

2 个答案:

答案 0 :(得分:36)

这样的事情对我有用:

myfunction <- function () {
  old <- .Random.seed
  set.seed(2)
  res <- runif(1)
  .Random.seed <<- old
  res
}

或许更优雅:

myfunction <- function () {
  old <- .Random.seed
  on.exit( { .Random.seed <<- old } )
  set.seed(2)
  runif(1)
}

例如:

> myfunction()
[1] 0.1848823
> runif(1)
[1] 0.3472722
> myfunction()
[1] 0.1848823
> runif(1)
[1] 0.4887732

答案 1 :(得分:0)

使用@Romain Francois的答案,将其概括为函数:

withRandom <- function(expr, seed = 1) {
    old <- .Random.seed
    on.exit({.Random.seed <<- old})
    set.seed(seed)
    expr
}

用法:

runif(2)
withRandom(seed = 2, {
    runif(1)
    runif(1)
})
runif(2)
withRandom(seed = 2, runif(2))
runif(2)

输出:

> runif(2)
[1] 0.5776099 0.6309793
> withRandom(seed = 2, {
+     runif(1)
+     runif(1)
+ })
[1] 0.702374
> runif(2)
[1] 0.5120159 0.5050239
> withRandom(seed = 2, runif(2))
[1] 0.1848823 0.7023740
> runif(2)
[1] 0.5340354 0.5572494

相关问题