adhd200preproc > RE: Question about liner trend and filter
Jun 2, 2011  06:06 AM | Cameron Craddock
RE: Question about liner trend and filter
Alright,

I have figured out the problem.  

I looked at sfnwmrda1779922_session_1_rest_1.nii.gz and its corresponding AAL timecourses.  The first of which is for ROI #2001. In the attached plot, the black line corresponds to the mean timecourse for this ROI extracted from snwmrda1779922_session_1_rest_1.nii.gz, which is correctly detrended.  The red line corresponds to the averaged band passed filtered time course from sfnwmrda1779922_session_1_rest_1_AAL_TCs.1D.  The blue line is the correctly bandpassed filtered time course.

The problem is with how I am calling 3dFourier to do the bandpass filtering, I used the -retrend option, which is described as:
"Any mean and linear trend are removed before filtering. This will restore the trend after filtering.".  I used this option as a carryover from the days when I bandpassed filtered (haven't used that in a while) and also wanted to preserve the image detail.  This is what is causing the odd behavior.

The odd thing is that the file that is the input to 3dFourier is nwmrda1779922_session_1_rest_1.nii.gz, which has already been detrended.  I verified this by re-detrending the black line in the plot, and got the black line again.  So somehow 3dFourier is finding a bogus trend in the data.  My guess is that what it is doing is fitting a straight line through the first point, which is -1 and the last point, which is 2 and ignoring everything in between.  

Interestingly, there is a new discussion on the AFNI forum about 3dFourier weirdness.

Thank you very much Cheng for finding this problem.  I will fix the data and put it on NITRC soon.

Cheers,
Cameron

Threaded View

TitleAuthorDate
cheng wei Jun 1, 2011
cheng wei Jul 12, 2011
cheng wei Jul 11, 2011
Cameron Craddock Jul 12, 2011
cheng wei Jul 14, 2011
Cameron Craddock Jul 11, 2011
cheng wei Jul 11, 2011
Cameron Craddock Jul 11, 2011
Cameron Craddock Jun 22, 2011
Cameron Craddock Jun 1, 2011
cheng wei Jun 2, 2011
Cameron Craddock Jun 2, 2011
Cameron Craddock Jun 2, 2011
Cameron Craddock Jun 2, 2011
Pierre Bellec Jun 2, 2011
Cameron Craddock Jun 2, 2011
Cameron Craddock Jun 2, 2011
RE: Question about liner trend and filter
Cameron Craddock Jun 2, 2011
cheng wei Jun 6, 2011
Carlton Chu Jun 7, 2011
Carlton Chu Jul 27, 2011
cheng wei Jun 8, 2011
Cameron Craddock Jun 8, 2011
cheng wei Jun 22, 2011
cheng wei Jun 16, 2011