Question about the output of 'psd_multitaper'

✉️ Mailing List Archive (read-only)
1

External Email - Use Caution

Hi all,
       When I used the command *mne.time_frequency.psd_multitaper, *my
input is raw data (n_channel,n_timepoints), However, my output was not
(n_channel,n_freqs) but was (n_channel,n_freqs*time(s)). Why did this
output come?
       Thank you for your attention
       Tao
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mail.nmr.mgh.harvard.edu/pipermail/mne_analysis/attachments/20190315/20511711/attachment.html

2

External Email - Use Caution

hi,

I am not sure to understand the problem. This seems to work:

import os.path as op
import mne

data_path = op.join(mne.datasets.sample.data_path(), 'MEG', 'sample')
raw = mne.io.read_raw_fif(op.join(data_path, 'sample_audvis_raw.fif'),
preload=True)
raw.crop(0, 15)

psd, freqs = mne.time_frequency.psd_multitaper(raw)

print(psd.shape)
# (364, 4506)

print(freqs.shape)
# (4506,)