Hi FSL experts,
I'm trying to understand how FSL does prewhitening and am writing MATLAB implementations of some of the prewhitening procedures in the article "Temporal Autocorrelation in Univariate Linear Modeling of FMRI Data". The paper states that for an autocovariance matrix V=KK^T :
"Given an estimate of V, V^{-1} and hence K^{-1} can be obtained by inverting in the spectral domain (some autocorrelation models, e.g., au toregressive, have simple parameterized forms for K^{-1}, and hence inversion in the spectral domain is not necessary)."
Using some of the non-parametric methods for autocorrelation estimation, I've found that I can't invert in the time domain. Do you have a reference for how inverting in the spectral domain is done?
Cheers,
-Garren
|