[Mne_analysis] Cortical labels power spectrum different approaches

Peled, Noam NPELED at mgh.harvard.edu
Wed Aug 29 14:32:39 EDT 2018
Search archives:

Hey all,

I'm calculating MEG cortical labels power spectrum (for resting state data) in two different ways. Can you help me understand the differences? The power spectrums are quite different (see attached).


1) Go through the source space time-series:

stcs = mne.minimum_norm.apply_inverse_epochs(epochs, ...)
labels_ts = mne.extract_label_time_course(stcs, labels, ...)
for label_ts:
  psds, freqs = mne.time_frequency.psd_array_welch(label_ts, ...)
  psds = 10 * np.log10(psds)



2) Compute the PSD from the epochs:


for label_ind, label in enumerate(labels):
  stcs = mne.minimum_norm.compute_source_psd_epochs(epochs, ...)
  for stc in stcs:
    psds = np.mean(stc.data, axis=0)



Thanks!

Noam

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mail.nmr.mgh.harvard.edu/pipermail/mne_analysis/attachments/20180829/f20f7f3d/attachment-0001.html 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: via_epochs.png
Type: image/png
Size: 108924 bytes
Desc: via_epochs.png
Url : http://mail.nmr.mgh.harvard.edu/pipermail/mne_analysis/attachments/20180829/f20f7f3d/attachment-0002.png 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: via_ts.png
Type: image/png
Size: 74593 bytes
Desc: via_ts.png
Url : http://mail.nmr.mgh.harvard.edu/pipermail/mne_analysis/attachments/20180829/f20f7f3d/attachment-0003.png 


More information about the Mne_analysis mailing list