
<html><head></head><body bgcolor='#FFFFFF' style='font-size:12px;background-color:#FFFFFF;font-family:Verdana, Arial, sans-serif;'>Hi Tobias,<br/><br/>it could be that the temporal activity at the sensor level is of muscular origin. Given all runs well in beamforming you shouldn't see too much of that insode the head - have you considered this possibility?<br/><br/>Michael<br/><br/><br/><br/><blockquote style="border-left: 2px solid blue; margin-left: 5px; padding-left: 5px; padding-top: 5px;"><hr/><b>Von:</b> "Tobias Navarro Schroeder" <tobiasnavarroschroeder@googlemail.com><br/><b>Gesendet:</b> May 30, 2011 5:58:32 PM<br/><b>An:</b> fieldtrip@donders.ru.nl<br/><b>Betreff:</b> [FieldTrip] Beamforming on resting state data<br/><br/><div class="gmail_quote"><p>Hi Fieldtrippers,</p><div> </div><div>I performed DICS beamforming on resting-state data (i.e. no task, no baseline, eyes closed) of a clinical population and controls.</div><div>According to the sensor data the patients have more gamma-band (60-90 Hz) activity over temporal areas (somewhat spread-out).</div><div>Curiously, after beamforming and group comparisons with 'indepsamplesT' permutation statistics the difference map shows more gamma-band (60-90 Hz) activity for controls in the left, middle-temporal gyrus.</div><div>Hence, the source-level results are opposite to the sensor-level results (although more spatially confined). This is <u>not</u> a problem of the design matrix, or confusing the groups.</div><div> </div><div>Question 1: Does anybody have experience with opposite results on source and sensor level?</div><div> </div><div>The individual raw beamforming results (interpolated on the subjects' MRIs) look suspicios as well. Due to the 'unit gain constraint' there should be an overestimation of power in the center of the head, but only a minority of subjects show that. Often I see temporal or parietal activity.</div><div> </div><div>Question 2:</div><div> </div><div>Would a lack of a blob in the center of the head in individual, raw, beamforming data indicate an erroneous source analysis?</div><div> </div><div>Thanks a lot for help and suggestions.</div><div> </div><div>Best, Tobias Navarro Schröder</div><div> </div><div>ps: The configurations for calculating the covariance and the beamforming is added below:</div><div> </div><div><p> </p><font color="#228b22" face="Courier New"><font color="#228b22" face="Courier New"><font color="#228b22" face="Courier New"><font size="1">1) % Compute the covariance matrices </font></font></font></font><p> </p><p><font face="Courier New"><font face="Courier New"><font size="1">cfg.output = </font></font></font></p><font color="#a020f0" face="Courier New"><font color="#a020f0" face="Courier New"><font color="#a020f0" face="Courier New" size="1">'powandcsd'</font></font></font><font face="Courier New"><font face="Courier New"><font size="1">;</font></font></font><p><font face="Courier New"><font face="Courier New"><font size="1">cfg.foi = 75;</font></font></font></p><p><font face="Courier New"><font face="Courier New"><font size="1">cfg.foilim = [75 75];</font></font></font></p><p><font face="Courier New"><font face="Courier New"><font size="1">cfg.tapsmofrq = 15; <font color="#228b22">% 75 +- 15 = 60 - 90 Hz</font></font></font></font></p><p><font face="Courier New"><font face="Courier New"><font size="1">cfg.pad = </font></font></font></p><font size="1"><font color="#a020f0" face="Courier New"><font color="#a020f0" face="Courier New"><font color="#a020f0" face="Courier New">'maxperlen'</font></font></font><font face="Courier New"><font face="Courier New">; </font></font><font color="#228b22" face="Courier New"><font color="#228b22" face="Courier New"><font color="#228b22" face="Courier New">% default</font></font></font></font><p><font face="Courier New"><font face="Courier New"><font size="1">cfg.method = </font></font></font></p><font size="1"><font color="#a020f0" face="Courier New"><font color="#a020f0" face="Courier New"><font color="#a020f0" face="Courier New">'wltconvol'</font></font></font><font face="Courier New"><font face="Courier New">; </font></font><font color="#228b22" face="Courier New"><font color="#228b22" face="Courier New"><font color="#228b22" face="Courier New">% wavelet analysis</font></font></font></font><p><font face="Courier New"><font face="Courier New"><font size="1">cfg.width = 5; </font></font></font></p><font color="#228b22" face="Courier New"><font color="#228b22" face="Courier New"><font color="#228b22" face="Courier New"><font size="1">% Width of the wavelet. Determines the temporal and spectral resolution of the analysis constant.</font></font></font></font><p> </p><p><font face="Courier New"><font face="Courier New" size="1"><font color="#228b22" face="Courier New" size="2"><font color="#228b22" face="Courier New" size="2"><font color="#228b22" face="Courier New" size="2">2) % Perform the actual beamforming</font></font></font></font></font></p><p><font face="Courier New"><font face="Courier New" size="1"><font color="#228b22" face="Courier New" size="2"><font color="#228b22" face="Courier New" size="2"><font color="#228b22" face="Courier New" size="2">% load grids of MRI volume</font></font></font></font></font></p><p><font face="Courier New"><font face="Courier New" size="1"><font face="Courier New" size="2"><font face="Courier New" size="2">fullnameGrid = [ CurSubject.mripath CurSubject.MRIModel];</font></font></font></font></p><p><font face="Courier New"><font face="Courier New" size="1"><font face="Courier New" size="2"><font face="Courier New" size="2">cfg = load(fullnameGrid);</font></font></font></font></p><p><font face="Courier New"><font face="Courier New" size="1"><font face="Courier New" size="2"><font face="Courier New" size="2">cfg.vol = cfg.hdm;</font></font></font></font></p><p> </p><p> </p><p><font face="Courier New"><font face="Courier New" size="1"><font face="Courier New" size="2"><font face="Courier New" size="2">cfg.frequency = 75;</font></font></font></font></p><p><font face="Courier New"><font face="Courier New" size="1"><font face="Courier New" size="2"><font face="Courier New" size="2">cfg.method =</font></font></font></font></p><font face="Courier New"><font face="Courier New" size="1"><font color="#a020f0" face="Courier New" size="2"><font color="#a020f0" face="Courier New" size="2"><font color="#a020f0" face="Courier New" size="2">'dics'</font></font></font><font face="Courier New" size="2"><font face="Courier New" size="2">;</font></font></font></font><p><font face="Courier New"><font face="Courier New" size="1"><font face="Courier New" size="2"><font face="Courier New" size="2">cfg.keepfilter =</font></font></font></font></p><font face="Courier New"><font face="Courier New" size="1"><font color="#a020f0" face="Courier New" size="2"><font color="#a020f0" face="Courier New" size="2"><font color="#a020f0" face="Courier New" size="2">'yes'</font></font></font><font face="Courier New" size="2"><font face="Courier New" size="2">; </font></font><font color="#228b22" face="Courier New" size="2"><font color="#228b22" face="Courier New" size="2"><font color="#228b22" face="Courier New" size="2">% keep the beamformer filters (virtual electrodes) for later use</font></font></font></font></font><p><font face="Courier New"><font face="Courier New" size="1"><font face="Courier New" size="2"><font face="Courier New" size="2">cfg.keeptrial =</font></font></font></font></p><font face="Courier New"><font face="Courier New" size="1"><font color="#a020f0" face="Courier New" size="2"><font color="#a020f0" face="Courier New" size="2"><font color="#a020f0" face="Courier New" size="2">'no'</font></font></font><font face="Courier New" size="2"><font face="Courier New" size="2">;</font></font></font></font><p><font face="Courier New"><font face="Courier New" size="1"><font face="Courier New" size="2"><font face="Courier New" size="2">cfg.keepmom =</font></font></font></font></p><font face="Courier New"><font face="Courier New" size="1"><font color="#a020f0" face="Courier New" size="2"><font color="#a020f0" face="Courier New" size="2"><font color="#a020f0" face="Courier New" size="2">'yes'</font></font></font><font face="Courier New" size="2"><font face="Courier New" size="2">; </font></font><font color="#228b22" face="Courier New" size="2"><font color="#228b22" face="Courier New" size="2"><font color="#228b22" face="Courier New" size="2">% useful for connectivity calculations. Dipole moments..</font></font></font></font></font><p><font face="Courier New"><font face="Courier New" size="1"><font face="Courier New" size="2"><font face="Courier New" size="2">cfg.projectmom =</font></font></font></font></p><font face="Courier New"><font face="Courier New" size="1"><font color="#a020f0" face="Courier New" size="2"><font color="#a020f0" face="Courier New" size="2"><font color="#a020f0" face="Courier New" size="2">'yes'</font></font></font><font face="Courier New" size="2"><font face="Courier New" size="2">;</font></font></font></font><p><font face="Courier New"><font face="Courier New" size="1"><font face="Courier New" size="2"><font face="Courier New" size="2">cfg.projectnoise =</font></font></font></font></p><font face="Courier New"><font face="Courier New" size="1"><font color="#a020f0" face="Courier New" size="2"><font color="#a020f0" face="Courier New" size="2"><font color="#a020f0" face="Courier New" size="2">'yes'</font></font></font><font face="Courier New" size="2"><font face="Courier New" size="2">;</font></font></font></font><p><font face="Courier New"><font face="Courier New" size="1"><font face="Courier New" size="2"><font face="Courier New" size="2">cfg.latency = 1;</font></font></font></font></p></div></div></blockquote></body></html>