我有一个旧的 S-plus 脚本,我想在 R 中重现结果。我遇到的唯一问题是随机种子。我知道他们使用不同的算法来生成伪随机数。在 S-plus 文件中,使用以下方式设置种子:
set.seed(337)
我能找到的相关信息是
从这个文档看来,S-plus 使用了“Super-duper”算法来生成伪随机数。R有这个选项,但它说它需要2个整数,而S-plus只需要0到1000之间的1个整数。此外,R doc说
The two seeds are the Tausworthe and congruence long integers, respectively.
A one-to-one mapping to S's .Random.seed[1:12] is possible but we will not
publish one, not least as this generator is not exactly the same as that in
recent versions of S-PLUS.
我不太确定这意味着什么。那么有人知道是否可以复制结果吗?
R 邮件列表上的旧帖子试图在 S/R/Splus 中获得相同的结果。