Hi, I have a question about the option -tdiff in 3dTstat. According the instruction message,

-tdiff: Means to take the first difference of each time

series before further processing.

It is not clear what kind of further processing? I checked the relevant source code, it looks like this:

/* RWC: first difference here [25 May 2011] */

if( do_tdiff ){

float tsm , tss ;

ts_dif = (float*)calloc(npts, sizeof(float)) ;

for( ii=1 ; ii < npts ; ii++ )

ts_dif[ii-1] = ts[ii] - ts[ii-1] ;

get_linear_trend( npts-1 , ts_dif , &tsm , &tss ) ;

ts_mean = (double)tsm ;

ts_slope = (double)tss ;

npts-- ;

ts = ts_dif ;

}

It looks like the outputs are ts_mean and ts_slope. But I only get one output (the mean?). If I have time series of sin(t), would I get mean=0 and slope=0? Such results would not provide any information of the changes of the time series. Could any one please tell me a mathematical expression of the -tdiff output?

Thank you.