Ð63Ð
Script for
Deconvolution - 4
#======================================================================
# 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
Mean
of the runth dataset,
through
time: run=1..10
¥ 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Õ)