Hi Mary,

To be sure, are you running 3dDeconvolve/3dREMLfit using -stim_times_AM2, where each event has an attached modulator in {0,1,2,3} (or {-1.5, -.5, .5, 1.5})?

If so, that assumes a linear relationship between the modulator and the amplitude of the BOLD response, where the output (at each voxel) is a beta for the mean response, and a beta for the incremental effect of the modulator.

Since the output is a mean and slope for the modulated conditions, there isn't really so much to plot.

What you might want to do is model all 4 levels separately, and then run a simple linear fit to that. That could be compared with the outputs from 3dDeconvolve (they should be similar), and it could be viewed as a 4-point time series, or even as scatter plots over ROIs/clusters, say.

Does this seem reasonable?

- rick