## Not run:
# slicedata <- read.fmrislice(fbase="fmri", slice=3, swap=FALSE)
# ymaskdata <- premask(slicedata)
# fsave <- paste(tempdir(),"/simultest1",fileext = ".sav", sep="")
# answ <- readline("Run MCMC simulation first ? ")
# run <- FALSE
# if (substr(answ, 1, 1) == "y") { run <- TRUE }
# if(run) {
# out <- cudaMultireg.slice(slicedata, ymaskdata, R=2000, keep=5, nu.e=3,
# fsave=fsave, zprior=FALSE)
# } else {
# load(file=fsave)
# cat("loaded",fsave,"\n")
# }
# pmeans <- pmeans.hcoef(out$betadraw)
# ## End(Not run)
Run the code above in your browser using DataLab