<div dir="ltr">

<p class="MsoNormal">Hello !</p><p class="MsoNormal"><br></p>

<p class="MsoNormal">I’m trying to modify the dics_source_power function (mne/beamformer/_dics.py)
and the compute_epochs_csd function (mne/time_frequency/csd.py) to apply DICS on
each source and each epoch individually. I’m trying to understand how these
functions work but I have a problem with DICS.</p><br>

<p class="MsoNormal"><span style="mso-spacerun:yes"> </span>Could you please help
me to understand<span style="mso-spacerun:yes">  </span>why the matrix filter
(‘Wk’) is divided by ‘Ck’ in the case of fixed source orientation (_dics.py ,
line 398).</p><p class="MsoNormal"><br></p>

<p class="MsoNormal">Ck = Gk.T.dot(Cm_inv).dot(Gk) <br></p><p class="MsoNormal">where G is the forward solution and Cm_inv<span style="mso-spacerun:yes">  </span>is the pseudoinverse of the CSD.</p><p class="MsoNormal"><br><span style="mso-ansi-language:EN-US" lang="EN-US"></span>

</p><p class="MsoNormal" style="line-height:normal">Any help or suggestion would be appreciated! <br></p><p class="MsoNormal" style="line-height:normal"><br></p>

<p class="MsoNormal" style="line-height:normal">Many thanks.</p><p class="MsoNormal" style="line-height:normal"><br><span style="mso-ansi-language:EN-US" lang="EN-US"></span></p>

<p class="MsoNormal" style="line-height:normal">Alexandre F</p>

</div>