对于set.seed()
语句,如果我没有显式地设置种子,我可以在运行一些代码之后获得它吗?
我重新运行了一些代码(交互式/在控制台上)包含随机化输入数据样本的函数(该函数是kohonen
包的一部分).用它玩了一段时间后,看看输出的各种(这是一个“不稳定”的问题),我注意到一个结果相当有趣。我当然没有使用set.seed()
,但是我想知道在运行代码后是否可以得到种子来重现结果?
在?set.seed
中我看到
.Random.seed保存均匀随机数生成器的种子集
但我不知道这有什么用。
5条答案
按热度按时间y3bcpkx11#
如果你没有保留种子,那么在你观察到随机抽取之后,没有通用的方法可以将随机数生成器“回滚”到之前的状态。接下来,你可能想做的是保存
.Random.seed
的值沿着你的计算结果,类似于这样。然后,您可以按如下方式重置PRNG:
result2
应与result
相同。6tr1vspr2#
Hong的回答是可靠的。对于快速和肮脏的解决方案,我只是重新执行整个脚本,直到我得到有趣的行为,我随机选择一个整数,打印出来,然后使用它作为种子。如果我的特定运行有有趣的行为,我注意到种子:
pengsaosao3#
为了补充mpettis给出的答案,如果您不想手动重新执行脚本--每次迭代都生成新的随机种子--您可以执行以下操作:
其中变量“runs”是一个正整数,指示要运行代码的次数。
通过这种方式,您可以快速生成大量输出,并为每个迭代提供单独的种子以实现可重复性。
bis0qfac4#
zynd9foi5#
下面是修复
setSeed
函数没有反getSeed
函数的问题的尝试。大约十二小时前,我在Using R, how to get.seed()?上发布了一个类似的问题,该问题被关闭,因为它被归类为“重复”。我“破解”了一个带有种子内存的解决方案,它需要一个全局变量
.random.seed.memory
。计时非常重要,因为我必须使用
set.seed
“生成种子函数
setSeed
的行为通常与set.seed
类似,但传递给set.seed
的任何超出整数(种类、正常种类、样本种类)的自定义参数都需要在args.set
中列出,因为setSeed
的省略号...
用于将参数传递给initSeed(...)
,initSeed(...)
是一个内部函数,使setSeed
和getSeed
能够工作。我还编写了一个C标准的
rand()
函数,它传入min,max,n方法,这就是我如何生成一个“整数”以提供给setSeed
并存储在内存中。我使用Sys.time()
作为默认种子生成的最小值/最大值(min = -1*as.integer(Sys.time())
和max = as.integer(Sys.time())
)。sample
是个坏主意因为它必须在该范围内创建一个向量来计算单个值,但它是rand()
的一个方法选项,它馈送initSeed
。我发现默认的“高-低”比“地板”稍快。用法
默认情况下,它将种子值存储到全局列表中名为“last”的元素中...这使您能够根据正在运行的进程跟踪不同的内存种子。在下面的示例中,我专门访问“last”和“nsim”...存储在内存中的第二个种子...
当然,它可能有错误,我将感谢任何建议或错误发现。
--太平洋标准时间2021年2月19日凌晨5点左右--
当然,传递固定种子的能力也是可能的。