## load simulation
fsave <- "/tmp/simultest1.sav"
load(file=fsave)
cat("loaded",fsave,"\n")
vreg <- 2
pmeans <- pmeans.hcoef(out$betadraw)
px <- regpostsim(pmeans, vreg=vreg)
spma <- px$spma # active voxels
spmn <- px$spmn # non-active voxels
plot(out$betadraw,spmname="activated",spm=spma)
plot(out$betadraw,spmname="non-activated",spm=spmn)
Run the code above in your browser using DataLab