[FieldTrip] Why are time-frequency results so much different from Beamformer source localization?

Tzvetan Popov tzvetan.popov at uni-konstanz.de
Fri Nov 13 12:31:41 CET 2020


Dear Patrick,

As you have already noticed the answer to your question is because you do the different things. I assume that what you did in the wavelet case was the following: you performed tf-decomposition on the epochs (correct and incorrect), applied a baseline correction and contrasted the baselined data. After subtraction of the two baselined conditions you found your effect in beta. If this is correct, in order to reproduce the same approach in source space I suggest you do the following.

1. Within each condition select the data into pre and post intervals of interest. For example, say your effect is in the latency from .5 to 1 sec after stimulus onset. Select for each condition a pre period (-1 to -.5) and post period (.5 to 1sec).
2. Append all the data using ft_appenddata and call ft_sourceanalysis as you did so far but with cfg.dics.keepfilter = ‘yes’;
3. Make 4 additional calls to ft_sourceanalysis for the pre and post intervals of the 2 conditions. Key here is to pass on the filters computed in the 2nd step by specifying cfg.sourcemodel.filter = outputfrom2ndstep.filter;
Why and how you should do this is further explained here https://www.fieldtriptoolbox.org/example/common_filters_in_beamforming/ <https://www.fieldtriptoolbox.org/example/common_filters_in_beamforming/>
4. Make a contrast baseline vs. Poststim activity per condition. This way you end up with two source maps, one for correct and one for incorrect. Use ft_math to compute a difference score and plot the result. 
The left-right flip has to do with the coordinates of the MR you are using. One easy fix is to follow the steps in the beamforming tutorial, i.e. ft_volumenormalise.

Good luck,
Tzvetan


> On 13 Nov 2020, at 11:20, Patrick Wiegel <patrick.wiegel at sport.uni-freiburg.de> wrote:
> 
> Dear Field-trip community,
> 
> I am working on EEG data from a motor learning task. in the first part of my analysis, I performed a time-frequency analysis (Morlet wavelet method) and cluster-based permutation test. I compared two different conditions (correct & incorrect movements) within subjects and found a significant difference in the beta band (20-35 Hz) when subjects were provided with feedback. I attached the topographical plot from one subject that illustrates the difference between the conditions at the selected time point and frequency (dif_morlet). As you can see, beta power was higher in one condition in lateral frontal electrodes.
> 
> Now, I would like to better spatially locate the effect for this single subject. For this purpose, I followed the Beamformer tutorial on the website. In the first step, I selected the raw data of interest (at the time where the effect was seen) and performed frequency analyses for the 20-35 Hz range (mtmfft, contrast between conditions shown in dif_mtmfft). There is an apparent difference between dif_morlet and dif_mtmfft and I am wondering why this is the case? One apparent reason is that the time-frequency data were baseline corrected and the frequency data from mtmfft not.
> 
> Using the power and CSD data from mtmfft, I performed the source analysis with a standard MRI and head model. The source was located in the left hemisphere (source plot), which is in stark contrast to the time-frequency data (dif_morlet). I am aware that the source analysis is based on a lot of assumptions (conductivity etc.) that influence the calculations but it is very difficult to make sense of such discrepancy. 
> Maybe it is necessary to also baseline correct the mtmfft data before inputting them to the source analysis? Or doing the contrasting source localization betwee a single condition and the corresponding baseline data before contrasting the 2 conditions? 
> 
> I am providing all necessary files and scripts (TF_analysis & BEAMFORMER) if anyone is interested to reproduce the results that I am showing and reporting for this single subject (https://github.com/PatrickWiegel/EEG-Beamformer-Source-localisation <https://github.com/PatrickWiegel/EEG-Beamformer-Source-localisation>). 
> 
> I would appreciate any help and discussion on that.
> 
> All the best,
> Patrick
> 
> 
> 
> 
> 
> Patrick Wiegel
> Department of Sport and Sport Science
> 
> University of Freiburg 
> Sandfangweg 4
> 79117 Freiburg i. Br . 
> 
> phone: +49 (0)761/ 203-4550
> email: patrick.wiegel at sport.uni-freiburg.de <mailto:patrick.wiegel at sport.uni-freiburg.de>
> web: www.sport.uni-freiburg.de
> 
> <dif_morlet.eps><dif_mtmfft.eps><sourceplot.eps>_______________________________________________
> fieldtrip mailing list
> https://mailman.science.ru.nl/mailman/listinfo/fieldtrip
> https://doi.org/10.1371/journal.pcbi.1002202

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.science.ru.nl/pipermail/fieldtrip/attachments/20201113/2a66423b/attachment.htm>


More information about the fieldtrip mailing list