Dear R HPC People:

Here is a function that I have put together:

library(foreach)
library(doMPI)
cl <- startMPIcluster()
registerDoMPI(cl)
zz <- numeric(length=10000)
foreach(i=1:10000,.combine=cbind) %dopar% {
 x <- arima.sim(list(order=c(1,0,0),ar=-0.9),n=120)
 z[i] <- arima(x,order=c(1,0,0))$model$phi
 }
assign("zz",z,envir=.GlobalEnv)
     # and shut it down
     closeCluster(cl)

but it doesn't save the zz variable.

Does anything look out of whack, please?

Thanks,
Erin



Erin M. Hodgess, PhD
Associate Professor
Department of Computer and Mathematical Sciences
University of Houston - Downtown
mailto: hodgesse@uhd.edu


	[[alternative HTML version deleted]]

