slicedata <- read.fmrislice(fbase="fmri", slice=3, swap=TRUE)
ymaskdata <- premask(slicedata)
fsave1 <- "/tmp/simultest1.sav"
out <- cudaMultireg.slice(slicedata, ymaskdata, R=2000, keep=5, nu.e=nu1, fsave=fsave1,
zprior=FALSE )
vreg <- 2
post.shrinkage.minmax(out, slicedata$X, vreg=vreg)Run the code above in your browser using DataLab