Welch PSD in source space

✉️ Mailing List Archive (read-only)
1

Dear MNE-users,

I would like to compute PSD in source space using Welch's method. I did not
find how to do it with the currently implemented functions:
- compute_source_psd is doing it but from raw file, yet I'ld like to
perform the analysis only on some parts of epochs,
- compute_source_psd_epochs is using multitaper,
- psd_welch had epochs/evoked instances as argument and not source data

How could I do it?

Thanks a lot,
Best,
Laetitia G.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mail.nmr.mgh.harvard.edu/pipermail/mne_analysis/attachments/20161107/e1c0932b/attachment.html

2

Hi Laetitia I would simply source localize epochs and compute the Welch on
the output. I happen to write some code for this right now, I'm happy to
share it with you. We can then make a PR to MNE. Currently the API there
does not seam to support this yet.
Denis

3

For now you can use `mne.time_frequency.psd._psd_welch` and pass sfreq from
your info and data from the stc.

4

Yes, I'ld be happy to use your code. I was also thinking about source
reconstructing PSD computed on the epochs, do you think it is valid?

2016-11-07 16:03 GMT+01:00 Denis-Alexander Engemann <
denis.engemann at gmail.com>:

Hi Laetitia I would simply source localize epochs and compute the Welch on
the output. I happen to write some code for this right now, I'm happy to
share it with you. We can then make a PR to MNE. Currently the API there
does not seam to support this yet.
Denis

Dear MNE-users,

I would like to compute PSD in source space using Welch's method. I did
not find how to do it with the currently implemented functions:
- compute_source_psd is doing it but from raw file, yet I'ld like to
perform the analysis only on some parts of epochs,
- compute_source_psd_epochs is using multitaper,
- psd_welch had epochs/evoked instances as argument and not source data

How could I do it?

Thanks a lot,
Best,
Laetitia G.
_______________________________________________
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

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.

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mail.nmr.mgh.harvard.edu/pipermail/mne_analysis/attachments/20161107/db54f3d9/attachment-0001.html

5

No only for the Fourier coefficients is valid once you compute the power
it's non-linear.

Use apply_inverse_epochs and do a welch on the output with the private
function and accumulate the result.

That's anyways what I will do.

6

Great, thanks!

2016-11-07 16:11 GMT+01:00 Denis-Alexander Engemann <
denis.engemann at gmail.com>:

No only for the Fourier coefficients is valid once you compute the power
it's non-linear.

Use apply_inverse_epochs and do a welch on the output with the private
function and accumulate the result.

That's anyways what I will do.

Yes, I'ld be happy to use your code. I was also thinking about source
reconstructing PSD computed on the epochs, do you think it is valid?

2016-11-07 16:03 GMT+01:00 Denis-Alexander Engemann <
denis.engemann at gmail.com>:

Hi Laetitia I would simply source localize epochs and compute the Welch
on the output. I happen to write some code for this right now, I'm happy to
share it with you. We can then make a PR to MNE. Currently the API there
does not seam to support this yet.
Denis

Dear MNE-users,

I would like to compute PSD in source space using Welch's method. I did
not find how to do it with the currently implemented functions:
- compute_source_psd is doing it but from raw file, yet I'ld like to
perform the analysis only on some parts of epochs,
- compute_source_psd_epochs is using multitaper,
- psd_welch had epochs/evoked instances as argument and not source data

How could I do it?

Thanks a lot,
Best,
Laetitia G.
_______________________________________________
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

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

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

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.

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mail.nmr.mgh.harvard.edu/pipermail/mne_analysis/attachments/20161107/f9dab915/attachment-0001.html