[Mne_analysis] Mne_analysis Digest, Vol 101, Issue 15

[Alexandre FABRE]

Hi

Yes, it's okay

Thank you so much

2016-06-27 18:00 GMT+02:00 <mne_analysis-request at nmr.mgh.harvard.edu>:

> Send Mne_analysis mailing list submissions to
>         mne_analysis at nmr.mgh.harvard.edu
>
> To subscribe or unsubscribe via the World Wide Web, visit
>         https://mail.nmr.mgh.harvard.edu/mailman/listinfo/mne_analysis
> or, via email, send a message with subject or body 'help' to
>         mne_analysis-request at nmr.mgh.harvard.edu
>
> You can reach the person managing the list at
>         mne_analysis-owner at nmr.mgh.harvard.edu
>
> When replying, please edit your Subject line so it is more specific
> than "Re: Contents of Mne_analysis digest..."
>
>
> Today's Topics:
>
>    1. Re: DICS source power (Alexandre Gramfort)
>
>
> ----------------------------------------------------------------------
>
> Message: 1
> Date: Mon, 27 Jun 2016 00:27:15 +0200
> From: Alexandre Gramfort <alexandre.gramfort at telecom-paristech.fr>
> Subject: Re: [Mne_analysis] DICS source power
> To: Discussion and support forum for the users of MNE Software
>         <mne_analysis at nmr.mgh.harvard.edu>
> Message-ID:
>         <
> CADeotZoN5b38KPjG1Rfm6o1ak3J-dPHU5zJD8Dz5gOvvY+ONpw at mail.gmail.com>
> Content-Type: text/plain; charset=UTF-8
>
> 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.
> >
>
>
>
> ------------------------------
>
> _______________________________________________
> Mne_analysis mailing list
> Mne_analysis at nmr.mgh.harvard.edu
> https://mail.nmr.mgh.harvard.edu/mailman/listinfo/mne_analysis
>
> End of Mne_analysis Digest, Vol 101, Issue 15
> *********************************************
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mail.nmr.mgh.harvard.edu/pipermail/mne_analysis/attachments/20160628/9a51f4a1/attachment.html