Here are a few simple things you can do with 3dDespike's spikiness. 3dcalc is uses to identify a voxel at a point in time as spike using the 2.5 threshold. Then the number of spikes are summed across time with 3dTstat. 3dBrickStat gives the maximum number of spikes and the total number of spikes within a range of values (from 2 to 32 spikes at that 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.
3dBrickStat -max spikesum+orig
32
3dBrickStat -count -non-zero spikesum+orig'<2..32>'
33839
You may also want to check 3dTqual and 3dToutcount for other statistical quality measures.