Hello.

I am not sure if this is the same problem, but i am trying to get iresp output in psc values.

I have tried scaling the data prior to 3dDeconvolve:

3dTstat -prefix ${ddir}/${RUN}mean_al \
${ddir}/${RUN}_al+orig
3dcalc -fscale -prefix ${ddir}/${RUN}_alpsc \
-a ${ddir}/${RUN}_al+orig -b ${ddir}/${RUN}mean_al+orig \
-expr "min(200, a/b*100)"

In the output, voxels within the brain have values aroun 100 and voxels outside the brain are around 0, but not null values (not sure if this could be the problem?)

The other way I've tried to scale is post-3dDeconvolve:
foreach TR (0 1 2 3 4 5 6 7 8 9 10)
3dcalc -fscale \
-a ${fmriID}${run}${cond}_${suf}.buc+orig'['0']' \
-b ${fmriID}${run}${cond}_${suf}.buc+orig'['2']' \
-c ${fmriID}${run}${cond}_${suf}iresp+orig'['${TR}']' \
-expr "astep(a,1000)*astep(b,1000)*((c*2)/(abs(a)+abs(b)))" \


where -a and -b (subbricks 0 and 2) are the baselines for 2 different functional runs, and -c (subbricks 0 to 10) are the TRs of the iresp file from 3dDeconvolve.

Both of these methods produce a dataset that is almost entirely full of zeros (null) values, except for a couple weird blobs at the perimeter of the frontal cortex.

My goal is simply to get this data into psc values for plotting. any help greatly appreciated.

thanks!
dharol