[Mne_analysis] Compute source power spectral density (PSD) in a label

[Alexandre Gramfort]

        External Email - Use Caution        

yes this should be correct

Alex

On Fri, Dec 27, 2019 at 11:38 PM Vivek Sharma
<vivek.sharma1510 at gmail.com> wrote:
>
>         External Email - Use Caution
>
> Hi,
> In continuation with the above query.
> If during the compute of psd the option db is set as false then how it can be plotted in db.
> stc = compute_source_psd(raw, inverse_operator, lambda2=1. / 9., method="dSPM",
>                          tmin=tmin, tmax=tmax, fmin=fmin, fmax=fmax,
>                          pick_ori="normal", n_fft=n_fft, label=label,
>                          dB=False)
>
> To plot power in db I did the following change in code:
> plt.plot(stc.times, 10*np.log10(stc.data.T))
> I need validation for the changes I have made.
>
> On Thu, Dec 26, 2019 at 2:08 AM Vivek Sharma <vivek.sharma1510 at gmail.com> wrote:
>>
>> Hi,
>> In the following example
>> (https://mne.tools/stable/auto_examples/time_frequency/plot_source_power_spectrum.html)
>> why are we multiplying stc.times and stc.data.T by 1e3?
>>
>> Secondly, I would like to bring it to notice that in the same example the plot that is being obtained shows frequency as 10000, 20000, etc. and the unit is represented as Hz which I believe is wrong as it should be millihertz.
>> I try to keep the frequency in Hz by making the following change in the code.
>> plt.plot(stc.times, (1e3 * stc.data.T))
>> is the change valid?
>>
>> --
>> Vivek Sharma
>
>
>
> --
> Vivek Sharma
> _______________________________________________
> Mne_analysis mailing list
> Mne_analysis at nmr.mgh.harvard.edu
> https://mail.nmr.mgh.harvard.edu/mailman/listinfo/mne_analysis