w*******9
发帖数: 1433 | 1 Hi all, I am using the mclapply function in the multicore package in R. I am
doing a simulation study which requires 10000 replicates. How should I set
the seed so that the other people can reproduce my result? I tried the
following code but it gives me different output every time I run it.
set.seed(123)
mclapply(1:3, rnorm, mc.set.seed=T) |
|
t*****w
发帖数: 254 | 2 set.seed(123)
x<-rnorm(10000)
mclapply(1:3, x, mc.set.seed=T)
am
set |
|
A*****n
发帖数: 243 | 3 MCMC没法并行化，但是你的20k个subject可以并行化，而不是使用loop
最简单的一个方法就是用parallel包的mclapply函数，如果
机器的cpu比较多的话，时间就可以省下啦了。
apply的速度和for loop没有太多的区别，google一下可以发现有很多的
apply vs. for速度比较的blog。 |
|
