<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><font size="2" face="Courier New"><font size="2" face="Courier New">
<p></p></font></font><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><font face="Courier New"><font face="Courier New">
<p><font size="1"></font></p></font></font><font face="Courier New"><font face="Courier New">
<p><font size="1">cfg.output = </font></p></font></font><font color="#a020f0" face="Courier New"><font color="#a020f0" face="Courier New"><font color="#a020f0" size="1" face="Courier New">'powandcsd'</font></font></font><font face="Courier New"><font face="Courier New"><font size="1">;</font> 
<p><font size="1">cfg.foi = 75;</font></p>
<p><font size="1">cfg.foilim = [75 75];</font></p>
<p><font size="1">cfg.tapsmofrq = 15; <font color="#228b22">% 75 +- 15 = 60 - 90 Hz</font></font></p>
<p><font size="1">cfg.pad = </font></p></font></font><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><font face="Courier New"><font face="Courier New"> 
<p><font size="1">cfg.method = </font></p></font></font><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><font face="Courier New"><font face="Courier New"> 
<p><font size="1">cfg.width = 5; </font></p></font></font><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><font color="#228b22" size="1" face="Courier New"></font> </p><font face="Courier New"><font size="1" face="Courier New"><font color="#228b22" size="2" face="Courier New"><font color="#228b22" size="2" face="Courier New"><font color="#228b22" size="2" face="Courier New">
<p>2) % Perform the actual beamforming</p>
<p>% load grids of MRI volume</p></font></font></font><font size="2" face="Courier New"><font size="2" face="Courier New">
<p>fullnameGrid = [ CurSubject.mripath CurSubject.MRIModel];</p>
<p>cfg = load(fullnameGrid); </p>
<p>cfg.vol = cfg.hdm;</p>
<p></p></font></font><font color="#228b22" size="2" face="Courier New"><font color="#228b22" size="2" face="Courier New"><font color="#228b22" size="2" face="Courier New">
<p></p></font></font></font><font size="2" face="Courier New"><font size="2" face="Courier New">
<p>cfg.frequency = 75;</p>
<p>cfg.method = </p></font></font><font color="#a020f0" size="2" face="Courier New"><font color="#a020f0" size="2" face="Courier New"><font color="#a020f0" size="2" face="Courier New">'dics'</font></font></font><font size="2" face="Courier New"><font size="2" face="Courier New">;</font></font><font size="2" face="Courier New"><font size="2" face="Courier New"> 
<p>cfg.keepfilter = </p></font></font><font color="#a020f0" size="2" face="Courier New"><font color="#a020f0" size="2" face="Courier New"><font color="#a020f0" size="2" face="Courier New">'yes'</font></font></font><font size="2" face="Courier New"><font size="2" face="Courier New">; </font></font><font color="#228b22" size="2" face="Courier New"><font color="#228b22" size="2" face="Courier New"><font color="#228b22" size="2" face="Courier New">% keep the beamformer filters (virtual electrodes) for later use</font></font></font><font size="2" face="Courier New"><font size="2" face="Courier New"> 
<p>cfg.keeptrial = </p></font></font><font color="#a020f0" size="2" face="Courier New"><font color="#a020f0" size="2" face="Courier New"><font color="#a020f0" size="2" face="Courier New">'no'</font></font></font><font size="2" face="Courier New"><font size="2" face="Courier New">; 
<p>cfg.keepmom = </p></font></font><font color="#a020f0" size="2" face="Courier New"><font color="#a020f0" size="2" face="Courier New"><font color="#a020f0" size="2" face="Courier New">'yes'</font></font></font><font size="2" face="Courier New"><font size="2" face="Courier New">; </font></font><font color="#228b22" size="2" face="Courier New"><font color="#228b22" size="2" face="Courier New"><font color="#228b22" size="2" face="Courier New">% useful for connectivity calculations. Dipole moments..</font></font></font><font size="2" face="Courier New"><font size="2" face="Courier New"> 
<p>cfg.projectmom = </p></font></font><font color="#a020f0" size="2" face="Courier New"><font color="#a020f0" size="2" face="Courier New"><font color="#a020f0" size="2" face="Courier New">'yes'</font></font></font><font size="2" face="Courier New"><font size="2" face="Courier New">; 
<p>cfg.projectnoise = </p></font></font><font color="#a020f0" size="2" face="Courier New"><font color="#a020f0" size="2" face="Courier New"><font color="#a020f0" size="2" face="Courier New">'yes'</font></font></font><font size="2" face="Courier New"><font size="2" face="Courier New">; 
<p>cfg.latency = 1;</p></font></font></font></font></div></div><br>