# \donttest{
oldpar <- par(no.readonly = TRUE)
library(cmR)
data(cmrsim)
mask=array(NA,c(30,30))
space.mbf=space.ci=array(NA,c(30,30,3))
for (i in 1:3){
mask=array(NA,c(30,30))
mask[cmrdata_sim[,,i,1]!=0]=1
temp=cmr.space(cmrdata_sim[,,i,], mask, input_sim)
space.mbf[,,i]=t(as.matrix(temp$mbf))
space.ci[,,i]=t(as.matrix(temp$ci))
}
par(mfrow=c(2,1))
imageMBF(maxresp_sim, zlim=c(0,5))
imageMBF(space.mbf, zlim=c(0,5))
imageMBF(space.ci, zlim=c(0,0.8))
par(oldpar)
# }
Run the code above in your browser using DataLab