我想在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只在我的函数中设置种子?
答案 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