关于set.seed的疑问

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

给位坛友，我看set.seed的help，还是不太明白。所以向大家请教一下。
我的问题是，我在一个程序中要做3次bootstrap的检验，为了使结果可重复。是不是只要在第一个test前设定种子（set.seed）还是每一个检验都得重新设定。还有set.seed（）的整数值可以选的范围是？不同的数有什么不同的？
谢谢！

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏1 回帖

关键词：seed see Set Bootstrap Bootstra 疑问

就在最开始设个set.seed就好了

qoiqpwqr 发表于 2013-5-30 21:15
就在最开始设个set.seed就好了

接版大的话，关于set.seed这个function一直觉得是一个学术造假的利器，这个种子事实上将会决定你未来生成的一批随机数（因为计算机生成随机数都是依赖于算法的）。所以你如果不是要重新开始模拟的话正如版大所说设一次就够了（当然重新开始模拟就要重新设了）
关于种子的取值，不知道R3.0.0之后的有没有改进，2.15.x的都是正负2^31-1，如果你设的种子带小数的话R会自动取整，也就是你设3和3.4是一样的
不同的数当然结果就不同咯
之前还在网上看到过一个反向获取种子的函数，lz可以玩一玩(只能是最简单的情况)

1. find.seed<-function(x,FUN = rnorm,seed = 0:10000, ...){
   
2.   res = NULL
   
3.    for (i in seed) {
   
4.      set.seed(i)
   
5.       rx = FUN(length(x), ...)
   
6.         if (all(x == rx)) {
   
7.             res = i
   
8.             break
   
9.         }
   
10.     }
    
11.     res
    
12. }
    
13. > set.seed(3.5);
    
14. > a<-rnorm(1);
    
15. > find.seed(a);
    
16. [1] 3

复制代码

qoiqpwqr 发表于 2013-5-30 21:15
就在最开始设个set.seed就好了

非常感谢！

求证1加1 发表于 2013-5-31 10:07
接版大的话，关于set.seed这个function一直觉得是一个学术造假的利器，这个种子事实上将会决定你未来生成 ...

谢谢！我还有点疑惑，如果只在第一次test之前设定了种子。那后面两个test生成的还是随机数么

kerrydu 发表于 2013-5-31 11:18
谢谢！我还有点疑惑，如果只在第一次test之前设定了种子。那后面两个test生成的还是随机数么

这本身只是一列伪随机数，意识就是说在一个特定的种子下，你之后产生的每一个随机数都已经被确定了
我的意思是你如果还想从头再来一遍的话就要重新再设了