Scale the Data:

• • Scale the Data:
  • foreach run ( $runs )
  • 3dTstat -prefix rm.mean_r$run pb03.$subj.r$run.blur+orig
  • 3dcalc -a pb03.$subj.r$run.blur_orig -b rm.mean_r$run+orig \
  • -c full_mask.$subj+orig \
  • -expr 'c * min(200, a/b*100)' \
  • -prefix pb04.$subj.r$run.scale
  • end
  • dataset a : the blurred EPI time series (for a single run)
  • dataset b : a single sub-brick, where each voxel has the mean value for that run
  • dataset c : the full mask
  • -expr 'c * min(200, a/b*100)'
  • • compute a/b*100 (the EPI value 'a', as a percent of the mean 'b')
    • if that value is greater than 200, use 200
    • multiply by the mask value, which is 1 inside the mask, and 0 outside