[Mne_analysis] DICS source power
Alexandre Gramfort
alexandre.gramfort at telecom-paristech.fr
Sun Jun 26 18:27:15 EDT 2016
hi,
> 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.
>
> Could you please help me to understand why the matrix filter (‘Wk’) is
> divided by ‘Ck’ in the case of fixed source orientation (_dics.py , line
> 398).
with fixed ori Ck is a float so
Wk[:] = np.dot(linalg.pinv(Ck, 0.1), Wk)
is replaced by
Wk /= Ck
ok?
if you think there is a but please with us on github a gist so we can
replicate the pb.
Alex
> Ck = Gk.T.dot(Cm_inv).dot(Gk)
>
> where G is the forward solution and Cm_inv is the pseudoinverse of the CSD.
>
>
> Any help or suggestion would be appreciated!
>
>
> Many thanks.
>
>
> Alexandre F
>
>
> _______________________________________________
> Mne_analysis mailing list
> Mne_analysis at nmr.mgh.harvard.edu
> https://mail.nmr.mgh.harvard.edu/mailman/listinfo/mne_analysis
>
>
> The information in this e-mail is intended only for the person to whom it is
> addressed. If you believe this e-mail was sent to you in error and the
> e-mail
> contains patient information, please contact the Partners Compliance
> HelpLine at
> http://www.partners.org/complianceline . If the e-mail was sent to you in
> error
> but does not contain patient information, please contact the sender and
> properly
> dispose of the e-mail.
>
More information about the Mne_analysis
mailing list