<html><head></head><body><div style="color:#000; background-color:#fff; font-family:HelveticaNeue, Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif;font-size:14px"><div id="yui_3_16_0_1_1448896819639_3038" dir="ltr">Dear <span class="" id="yui_3_16_0_1_1448896819639_3033">Frédéri,</span></div><div id="yui_3_16_0_1_1448896819639_3039" dir="ltr"><br><span class="" id="yui_3_16_0_1_1448896819639_3033"></span></div><div id="yui_3_16_0_1_1448896819639_3204" dir="ltr"><span class="" id="yui_3_16_0_1_1448896819639_3033">maybe you have the solution already.</span></div><div id="yui_3_16_0_1_1448896819639_3267" dir="ltr"><br><span class="" id="yui_3_16_0_1_1448896819639_3033"></span></div><div id="yui_3_16_0_1_1448896819639_3268" dir="ltr"><span class="" id="yui_3_16_0_1_1448896819639_3033">In case not, could you explain how did you get the last two plots? Are spectra computed after averaging across all virtual channel spectra or you are representing the spectrum of the source with maximal power?</span></div><div id="yui_3_16_0_1_1448896819639_3559" dir="ltr"><br><span class="" id="yui_3_16_0_1_1448896819639_3033"></span></div><div id="yui_3_16_0_1_1448896819639_3560" dir="ltr"><span class="" id="yui_3_16_0_1_1448896819639_3033">Best,</span></div><div id="yui_3_16_0_1_1448896819639_3563" dir="ltr"><span class="" id="yui_3_16_0_1_1448896819639_3033">Maite<br></span></div><br><div class="qtdSeparateBR"><br><br></div><div style="display: block;" class="yahoo_quoted"> <div style="font-family: HelveticaNeue, Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif; font-size: 14px;"> <div style="font-family: HelveticaNeue, Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif; font-size: 16px;"> <div dir="ltr"><font size="2" face="Arial"> El Jueves 26 de noviembre de 2015 13:55, Frédéric Roux <f.roux@bcbl.eu> escribió:<br></font></div> <br><br> <div class="y_msg_container">Dear all,<br><br>I am observing a counter-intuitive result after applying the spatial filters obtained from DICS beamforming to my MEG data.<br><br>The plot in the attachment summarizes the issue: The raw data shows a clear peak of activity at 10 Hz, however, after applying the spatial filter to the MEG data the spectrum of the virtual channels shows that power in this frequency band is reduced. The same also happens when I apply a spatial filter in the gamma range for power at 70 Hz to the raw MEG signals.<br><br>Is this something that falls out of the beamformer math due to the fact how the algorithm suppresses power at all locations by minimizing the noise level, or is there something wrong with my code? I (maybe naively) assumed that the results should come out the opposite. On the other hand, when plotting the noise normalized maps (NAI) of both the alpha and gamma beamformers the sources look quite accurate. Should I normalize with the level of noise first before plotting the spectrum to see the expected results (ie an upregulation of power instead of a downregulation of power)?<br> <br>The code that I use to generate the figures is:<br><br>% 1) compute CSD matrix<br>cfg = [];<br>cfg.method = 'mtmfft';<br>cfg.output = 'powandcsd';<br>cfg.foilim = [10 10];<br>cfg.tapsmofrq = 2;<br>cfg.pad = 'maxperlen';<br>cfg.taper = 'dpss';<br><br>[csd] = ft_freqanalysis(cfg,meg_data);<br><br><br>% 2) compute DICS filter<br>cfg = [];<br>cfg.method = 'dics';<br>cfg.grad = meg_data.grad;<br>cfg.headmodel = vol;<br>cfg.frequency = csd.freq;<br><br>cfg.dics.realfiter = 'yes';<br>cfg.dics.fixedori = 'yes';<br><br>[alpha_filter]= ft_sourceanalysis(cfg,csd);<br><br>% 3) compute the virtual channel data<br><br>VC = meg_data;<br>VC.trial = cell(1,length(meg_data.trial));<br>VC.label = cell(1,length(alpha_filter.avg.pow));<br><br>for i = 1:length(VC.trial)<br> for j = 1:length(VC.label)<br><br> VC.trial{i}(j,:) = alpha_filter{j}*meg_data.trial{i};<br> VC.label(j) = {['virtual_channel',num2str(j)]};<br> end;<br>end;<br><br>% 4) compute the spectrum of the virtual channel data<br>cfg = [];<br>cfg.method = 'mtmfft';<br>cfg.pad = 'maxperlen';<br>cfg.taper = 'dpss';<br>cfg.tapsmofrq = 1;<br>cfg.foi = 0.1:100;<br><br>[pow] = ft_freqanalysis(cfg,VC);<br><br><br>Any help or suggestions would be greatly appreciated.<br><br>Fred<br><br>--<br>Frédéric Roux<br>Postdoctoral Scientist, Marie-Curie fellow<br>BCBL. Basque Center on Cognition, Brain & Language.<br><br><a ymailto="mailto:f.roux@bcbl.eu" href="mailto:f.roux@bcbl.eu">f.roux@bcbl.eu</a><br>Tel: +34 943 309 300 Ext 211<br>Fax: +34 943 309 052<br><br>Legal disclaimer/Aviso legal/Lege-oharra: www.bcbl.eu/legal-disclaimer<br>---------------------------------------------------------------------------<br><br>“The probability of success is difficult to estimate; but if we never search the chance of success is zero.”<br>_______________________________________________<br>fieldtrip mailing list<br><a ymailto="mailto:fieldtrip@donders.ru.nl" href="mailto:fieldtrip@donders.ru.nl">fieldtrip@donders.ru.nl</a><br><a href="http://mailman.science.ru.nl/mailman/listinfo/fieldtrip" target="_blank">http://mailman.science.ru.nl/mailman/listinfo/fieldtrip</a><br><br></div> </div> </div> </div></div></body></html>