Hi Lisa,
For creating unit vectors of the time series, you might consider
what @compute_gcor does:
3dTnorm -polort 0 -prefix OUTPUT INPUT
The default scaling is the L2 norm, setting the sum of squares
equal to 1. If instead you want the standard deviation, then there
is a factor of sqrt(NT-1) to scale in. Maybe I should consider
adding such an option...
- rick