#======================================================================
# re-scale each
run's baseline to 100.
# If baseline is
100, and result of 3dcalc on one voxel is 106, then
# we can say that
at that voxel shows a 6% increase is signal activity
# relative to
baseline.
# Use full_mask
to remove non-brain
#======================================================================
foreach run (
`count -digits 1 1 10` )
3dTstat -prefix
mean_r{$run} {$subj}_r{$run}_vr_bl+orig
3dcalc -a
{$subj}_r{$run}_vr_bl+orig \
-b
mean_r{$run}+orig \
-c
full_mask+orig
\
-expr
"(a/b * 100) * c" \
-prefix scaled_r{$run}
/bin/rm
mean_r{$run}+orig*
end
¥ Divide
each voxel value (ÔaÕ) by its temporal mean (ÔbÕ) and scale by 100
¥ Result
will have temporal mean of 100
¥ Voxels
not in the mask will be set to 0 (by ÔcÕ)