# source band induced power

**URL:** <https://mne.discourse.group/t/source-band-induced-power/1501>\
**Category:** Mailing List Archive (read-only)\
**Tags:** list-archive\
**Created:** [April 30, 2018, 8:32pm UTC](https://mne.discourse.group/t/source-band-induced-power/1501 "2018-04-30T20:32:43Z")\
**Posts on this page:** 2\
**Page:** 1

<div class="post-metadata">

**Author:** ![system](https://global.discourse-cdn.com/free1/uploads/mne/original/1X/85cc6bd2b69cb698a166dc6d880fb550510d0144.jpeg) [@system](https://mne.discourse.group/u/system)\
**Post date:** [April 30, 2018, 8:32pm UTC](https://mne.discourse.group/t/source-band-induced-power/1501/1 "2018-04-30T20:32:43Z")

</div>

Dear group,

I'm trying to estimate the mean power over time in different bands for different labels.

I've tried the two following approaches:

1) Running mne.minimum\_norm.source\_band\_induced\_power on the epochs, (which is eating almost all the memory of our RAM monster)

2) First calculating the labels time course from the stc, then calculating the power using scipy.signal.welch

I've noticed, that for example for one of the labels and one of the bands, Welch gives me 1.886, wherein the first approach, the stc.data.mean(0) is between 45,260 and 145,023.

What am I'm missing?

Thanks!

Noam  
-------------- next part --------------  
An HTML attachment was scrubbed...  
URL: [http://mail.nmr.mgh.harvard.edu/pipermail/mne\_analysis/attachments/20180430/2bd4c1d5/attachment.html](http://mail.nmr.mgh.harvard.edu/pipermail/mne_analysis/attachments/20180430/2bd4c1d5/attachment.html)

---

<div class="post-metadata">

**Author:** ![system](https://global.discourse-cdn.com/free1/uploads/mne/original/1X/85cc6bd2b69cb698a166dc6d880fb550510d0144.jpeg) [@system](https://mne.discourse.group/u/system)\
**Post date:** [May 3, 2018, 7:48pm UTC](https://mne.discourse.group/t/source-band-induced-power/1501/2 "2018-05-03T19:48:01Z")

</div>

External Email - Use Caution

hi Noam,

> I'm trying to estimate the mean power over time in different bands for different labels.
> 
> I've tried the two following approaches:
> 
> 1) Running mne.minimum\_norm.source\_band\_induced\_power on the epochs, (which is eating almost all the memory of our RAM monster)
> 
> 2) First calculating the labels time course from the stc, then calculating the power using scipy.signal.welch
> 
> I've noticed, that for example for one of the labels and one of the bands, Welch gives me 1.886, wherein the first approach, the stc.data.mean(0) is between 45,260 and 145,023.
> 
> What am I'm missing?

I am not sure it's the reason but source\_band\_induced\_power using  
Morlet wavelets while welch  
is using an averaged periodogram with fixed window size.

Note that source\_band\_induced\_power has not much been improved since  
it was written  
7 years ago...

Alex
