Also, in a previous post, you suggest using 3dTstat to calculate a spikesum per voxel.
3dDespike -ssave spikiness -prefix despiked epi_r1+orig
3dcalc -a spikiness+orig. -expr 'step(a-2.5)' -prefix spike
3dTstat -sum -prefix spikesum spike+orig.
I'd also like to use the "spike+orig" dataset from 3dcalc to calculate the number of spikes across all voxels for each TR, because I'm noticing that certain TRs have a lot more spikes than others (which is to be expected). I've only ever done "voxel-wise" calculations and not "TR-wise" calculations on all brain voxels for a given TR. Can you tell me how to do this in AFNI.
Thanks!
Michael