# SNR estimate for single trials

**URL:** <https://mne.discourse.group/t/snr-estimate-for-single-trials/1199>\
**Category:** Mailing List Archive (read-only)\
**Tags:** list-archive\
**Created:** [March 30, 2017, 10:47am UTC](https://mne.discourse.group/t/snr-estimate-for-single-trials/1199 "2017-03-30T10:47:32Z")\
**Posts on this page:** 6\
**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:** [March 30, 2017, 10:47am UTC](https://mne.discourse.group/t/snr-estimate-for-single-trials/1199/1 "2017-03-30T10:47:32Z")

</div>

Hello everyone,

I want to ask what would be the recommended SNR amplitude estimate when  
computing single trial source estimates with  
mne-python(apply\_inverse/apply\_inverse\_epochs)?

I looked at this tutorial  
[http://martinos.org/mne/stable/auto\_examples/inverse/plot\_compute\_mne\_inverse\_epochs\_in\_label.html](http://martinos.org/mne/stable/auto_examples/inverse/plot_compute_mne_inverse_epochs_in_label.html)

and for single trial SNR is kept the same as for the evoked data -  
default 3. Would it not make sense to reduce it to less (e.g. 1, I have  
been using previously with mne-c) if we know that single trials are much  
noisier?

Thanks for your help,  
Ana

---

<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:** [March 30, 2017, 11:09am UTC](https://mne.discourse.group/t/snr-estimate-for-single-trials/1199/2 "2017-03-30T11:09:24Z")

</div>

hi Ana,

yes 1 is typically what is used if you estimate things on single trials.

I think we should update this example to clarify / change this.

See [https://github.com/mne-tools/mne-python/issues/4131](https://github.com/mne-tools/mne-python/issues/4131)

Alex

---

<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:** [March 30, 2017, 12:44pm UTC](https://mne.discourse.group/t/snr-estimate-for-single-trials/1199/3 "2017-03-30T12:44:57Z")

</div>

Hi Alex,

Thanks for the clarification.

Ana

---

<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 2, 2017, 2:19pm UTC](https://mne.discourse.group/t/snr-estimate-for-single-trials/1199/4 "2017-04-02T14:19:13Z")

</div>

Hi Ana and Alex,

I am not so sure about specifying a low SNR for single trials... Yes, it's true of course that single trials are noisier, but often single trial processing in the source space is followed by combining across trials.. So in many (but not all) instances, it is useful to have a consistent inverse operator (specifying a fixed SNR... and for dSPM scaling the noise cov by a fixed number of trials).

For instance, if you did inverse first on single trials and then averaging in source space vs. average across trials in sensor space first and then do inverse, you would get different answers if you specified different SNRs for the two cases.

This issue came up before, and the 3.0 was used to make them consistent:

[https://mail.nmr.mgh.harvard.edu/pipermail//mne\_analysis/2014-April/002066.html](https://mail.nmr.mgh.harvard.edu/pipermail//mne_analysis/2014-April/002066.html)

[https://github.com/mne-tools/mne-python/pull/1237](https://github.com/mne-tools/mne-python/pull/1237)

Just something to keep in mind depending on how you are using single trial inverses.

- Hari

---

<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 2, 2017, 4:34pm UTC](https://mne.discourse.group/t/snr-estimate-for-single-trials/1199/5 "2017-04-02T16:34:48Z")

</div>

Hi All,

I second Hari, I always used 3.0 in my old Matlab stream for processing  
single trial data for event related connectivity.

Thanks

Sheraz

---

<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 3, 2017, 2:07am UTC](https://mne.discourse.group/t/snr-estimate-for-single-trials/1199/6 "2017-04-03T02:07:42Z")

</div>

ok fair enough. Can we just add some comments in the examples then?

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