Hi Caroline,

The TENT betas should be contiguous, ordered by mean response, then modulators.

So in this case, the 8 mean response betas for "hour_target" should be #0..#7, [50..64(2)], then come the 8 for the first modulator (RT), then the 8 for confidence.

The gltsym, when not provided with a condition index list, will use the sum of all betas. I think the output of the 3dDeconvolve execution should indicate that, if you would like to verify. For example, the all_targets GLT is probably 0.25 times the sum of 96 betas(!), 24 for each condition.

To get just the sum of the mean TENT response betas use ' .25*hour_targets[0..7] +...'

I am not sure about the follow-up. Do you mean you would like to look at the 8 betas of each modulation term, or something like an F-stat for it? Would you like to extract each of the 8 betas into a new dataset? If so, using 3dbucket with a selector like the above, [50..64(2)], should work. That means to extract volume index 50 though 64, with a step of 2 (i.e. 50, 52, 54, ..., 62, 64).

Alternatively, you could run 3dDeconvolve with the -iresp option, but that would give you a TR-grid time series of only the mean response.

Do you want the betas, stats or some GLT?

- rick