adhd200preproc
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
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
| Title | Author | Date |
|---|---|---|
| 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 | |
| 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 | |
