r-如何在循环中调用种子向量,每次循环运行时将结果并排保存



我有两个问题:

  1. 我正在进行模拟,我编写了必要的函数,并在循环中生成了我需要的东西,其中包含以下代码:
sim_seeds<- as.vector(sample(1:30000, 5, replace = FALSE), mode = "numeric")
save(sim_seeds,file = "~/Desktop/untitled folder/sim_seeds.Rda")
load(file = "~/Desktop/untitled folder/sim_seeds.Rda")
for (i in 1:5) {
load(file = "~/Desktop/untitled folder/sim_seeds.Rda")
seeds=set.seed(sim_seeds[i])
#Generating data using functions had been wrote before
data<-generate_data(seed =seeds )
Y_1<- mean(data$Y)
#estimation
weights<-generate_weights(T1=S1~Year+growth, T2=R1~Age+Sex+HIV, data=data)
w<-weights$w
g<-g_est(data=data)
p1<-g$p1
Q<-Q_est(data=data,w , p1)
mu1_Q<-Q$mu1

#Results
results <- rbind(seeds,Y_1,mu1_Q)
results
}

我的问题是种子部分!我想做的是生成5个不同的数据集,但每次都需要一个单独的种子;对于";运行,所以我想创建一个种子向量,然后在每次循环运行时调用第I个值,但当我想在循环内部调用它时,它会给出一个NULL值!

  1. 其他问题是,我希望将最终结果保存并并排打印,以便进行比较。所以简单地说,我生成了5个不同的数据集,所以我在";结果">

有人知道吗?

1.:在开始时设置一个种子就足以在每次循环迭代中获得不同的结果。编辑后用下面的示例进行说明

2.:例如:

set.seed(123)
seeds = 1:5
results = c()
for (i in 1:5){
set.seed(1) #this is wrong, produces the same set of values at each loop iteration
# set.seed(seeds[i]) # works, but is unecessary, setting one seed at the beginning is ok
# commenting the two lines above is fine, and will generate 5 different vectors of random numbers. 
x = rnorm(5,0,1)
results = rbind(results, x)
}
results

最新更新