AFNI program: 3dTstat

Usage: 3dTstat [options] dataset
Computes one or more voxel-wise statistics for a 3D+time dataset
and stores them in a bucket dataset.

Options:
 -mean   = compute mean of input voxels [DEFAULT]
 -slope  = compute mean slope of input voxels vs. time
 -stdev  = compute standard deviation of input voxels
             [N.B.: this is computed after    ]
             [      the slope has been removed]
 -cvar   = compute coefficient of variation of input
             voxels = stdev/fabs(mean)
   **N.B.: You can add NOD to the end of the above 2
           options to turn off detrending, as in
             -stdevNOD or -cvarNOD

 -MAD    = compute MAD (median absolute deviation) of
             input voxels = median(|voxel-median(voxel)|)
             [N.B.: the trend is NOT removed for this]
 -DW    = compute Durbin-Watson Statistic of
             input voxels
             [N.B.: the trend is removed for this]
 -median = compute median of input voxels  [undetrended]
 -min    = compute minimum of input voxels [undetrended]
 -max    = compute maximum of input voxels [undetrended]

 -prefix p = use string 'p' for the prefix of the
               output dataset [DEFAULT = 'stat']
 -datum d  = use data type 'd' for the type of storage
               of the output, where 'd' is one of
               'byte', 'short', or 'float' [DEFAULT=float]
 -autocorr n = compute autocorrelation function and return
               first n coefficients
 -autoreg n = compute autoregression coefficients and return
               first n coefficients
    [N.B.: -autocorr 0 and/or -autoreg 0 will return coefficients
           equal to the length of the input data

The output is a bucket dataset.  The input dataset
may use a sub-brick selection list, as in program 3dcalc.

INPUT DATASET NAMES
-------------------
This program accepts datasets that are modified on input according to the
following schemes:
  'r1+orig[3..5]'                                    {sub-brick selector}
  'r1+orig<100.200>'                                 {sub-range selector}
  'r1+orig[3..5]<100..200>'                          {both selectors}
  '3dcalc( -a r1+orig -b r2+orig -expr 0.5*(a+b) )'  {calculation}
For the gruesome details, see the output of 'afni -help'.