Dear Pierre,
thank you for your detailed response.
The default cutoff makes perfect sense now, thank you. Additional
information about this issue is also present in the thread referenced by
Donald.
Concerning the filtering approach, point (2) makes sense, although
should be of no concern when inference is at the group level. Concerning
point (1), I wonder whether a possible approach to minimizing edge
effects with non-DCT filters could be to simply acquire an additional N
baseline scans at the beginning and end of the fMRI time series of
interest, with N = order of FIR filter, and then discard those 2*N scans
after forward-reverse filtering.
I am also curious about the comparisons of different filtering
approaches reported in this page:
http://nipy.sourceforge.net/nitime/examples/filtering_fmri.html
Different approaches appear to lead to different levels of distortion of
the spectral band of interest. It would be interesting to include the
DCT approach in the comparison. If this has been done before, I would
appreciate if someone could point me to the relevant literature.
Best,
Bruno
On 05/07/2013 7:22 PM, Pierre Bellec wrote:
> Dear Bruno,
>
> I am aware of at least two reasons. First edges effects. Standard FIR
> filters (e.g. Butterworth) with off-the-shelf parameters are very poorly
> adapted to fMRI time series. This is because the TR in fMRI (~1s) is
> very long compared to typical applications in signal processing. If you
> filter time series over the full brain, normalize the time series to
> zero mean and unit variance and then compute standard deviation across
> the brain, you'll find larger deviations at the edges (or at least
> that's what I remember from a few experiments I ran years ago). The DCT
> solves that problem, because it relies on a periodization of the signal
> that does not introduce discontinuities at the edges. It may be possible
> to tweak a FIR filter better for fMRI time series, but I don't have
> experience with that. (2) the regression of discrete cosines fits in the
> linear model framework and it thus properly handled in the statistical
> inference.
>
> Otherwise, regarding the amount of noise in different frequency band, I
> guess that at least part of the reasoning is the 1/f spectral power of
> fMRI time series (http://www.citeulike.org/user/pbellec/article/956093).
> There is lot of variance below 0.01 Hz , and that variance is mostly
> noise. That can be seen by observing the spatial distribution of the
> relative variance in ultra-low frequencies, which follows notably the
> edges of the brain. We also showed a couple years ago that simulation of
> motion generates slow time drifts (below 0.01 Hz) in an amount somewhat
> similar to what is actually observed in real data
> (http://www.ncbi.nlm.nih.gov/pubmed/19570641, Fig. 6a). The simulation
> model was very crude, but it does suggest that this frequency band is
> very noisy.
>
> If you find some good reference down the road about this issue, please
> share, I'd be interested. I hope this helps, best regards,
>
> Pierre Bellec, PhD
> Researcher / Chercheur, Research Centre of the Montreal Geriatric Institute
> Professeur adjoint sous octroi, Department of Computer Science and
> Operations Research
> University of Montreal, Québec, Canada
> 1 514 340 3540 #3367
> http://simexp-lab.org/brainwiki/doku.php?id=pierrebellec
>
>
> 2013/7/4 Bruno L. Giordano <[log in to unmask] <mailto:[log in to unmask]>>
>
> Hello,
>
> I am interested in the implementation of the temporal filtering of
> fMRI time series. In particular, I am trying to understand the
> advantages associated with the method based on the discrete cosine
> transform (DCT) basis set: are they simply computational, or there
> is something else that makes it preferable to a more classical
> forward-reverse filtering approach based on, e.g., FIR filters?
>
> Additionally, could someone please point me towards studies of the
> amount of noise in different frequency bands of fMRI time series?
> What's the reasoning behind the default of a 128-s period high pass
> in SPM?
>
> Thank you,
>
> Bruno
>
> ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~__~~
> Bruno L. Giordano, PhD
> Institute of Neuroscience and Psychology
> 58 Hillhead Street, University of Glasgow
> Glasgow, G12 8QB, Scotland
> T +44 (0) 141 330 5484 <tel:%2B44%20%280%29%20141%20330%205484>
> Www: http://www.brunolgiordano.net
>
>
--
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
Bruno L. Giordano, PhD
Institute of Neuroscience and Psychology
58 Hillhead Street, University of Glasgow
Glasgow, G12 8QB, Scotland
T +44 (0) 141 330 5484
Www: http://www.brunolgiordano.net
|