Mne_analysis Digest, Vol 101, Issue 15

✉️ Mailing List Archive (read-only)
1

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
        Mne_analysis Info Page
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
> Mne_analysis Info Page
>
>
> 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
> MyComplianceReport.com: Compliance and Ethics Reporting . 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
Mne_analysis Info Page

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