I was denoising fMRI data using GLM and had a question regarding the “data” of fMRI in general.
(For reference you can look at the “GLM” part of https://fmriprep.org/en/stable/outputs.html or https://www.sciencedirect.com/science/article/pii/S1053811917310972)
Below is a time series of the before/after GLM for a given voxel:
As you can see, the raw values (actually we went through a preprocessing step prior to GLM denoising) had a value of around $600\pm 50$, while the residual (i.e. denoised data) had a value around $0\pm 50$. In other words, it seems that via GLM, the raw data was somewhat “centered”. At first I thought this was OK, as I thought that only the fluctuations of the data matter (i.e. I had this assumption that mean centering doesn’t remove any information in the data). However, I realized that when we did this, the spatial distribution of signal values before and after GLM changes, as the pictures below show (the picture below is the distribution of fMRI signal for a given time point before and after GLM denoising) :
Is this OK?
In other words,
- does only the “fluctuations” (across time) of the data matter in fMRI (i.e. is adding a constant to an fMRI time series data irrelevant?).
- if so, then since adding a constant to an voxel’s time series data is irrelevant, it stands to reason that if I add different constant values to different voxels (for example as part of a detrending step), the spatial distribution of values before and after would change. Is this OK? (In other words, does the distribution of fMRI signals for a given time point meaningless?)
Thank you :)