古い S-plus スクリプトがあり、R で結果を再現したいと考えています。私が抱えている唯一の問題は、ランダム シードです。疑似乱数の生成に異なるアルゴリズムを使用していることは知っています。S-plus ファイルでは、次を使用してシードが設定されました。
set.seed(337)
私が見つけた関連情報は
このドキュメントから、S-plus は疑似乱数の生成に「Super-duper」アルゴリズムを使用したようです。Rにはこのオプションがありますが、2つの整数を受け取ると書かれていますが、S-plusは0から1000までの整数を1つしか必要としません。さらに、Rのドキュメントは言う
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 で同じ結果を得ようとしています。