slicedata <- read.fmrisample(slice=3)
fsave <- "/tmp/simultest.sav"
answ <- readline("Run MCMC simulation first ? ")
run <- FALSE
if (substr(answ, 1, 1) == "y") { run <- TRUE }
if(run) {
out <- cudaMultireg.slice(slicedata, fsave=fsave )
} else {
load(file=fsave)
cat("loaded",fsave,"\n")
}
pmeans <- pmeans.hcoef(out$betadraw)
Run the code above in your browser using DataLab