<div dir="ltr">Thanks for the reference. In return, I also recommend reading Ciuparu and Mures an's recent rebuttal:<div><br></div><div>European Journal of Neuroscience, Vol. 43, pp. 861–869, 2016, <a href="http://dx.doi.org/10.1111/ejn.13179">doi:10.1111/ejn.13179</a><br></div><div><div>TECHNICAL SPOTLIGHT</div><div>Sources of bias in single-trial normalization procedures</div></div><div><br></div><div>They demonstrate that the positive bias Grandchamp and Delorme warned about with single-trial baseline normalization was, in fact, due to the correlated numerators and denominators in each of the baseline normalization procedures they tested, which was, in turn, caused by the skewed distributions of baseline power values. They demonstrate that this bias is reduced by using a much longer baseline, ideally incorporated into the experimental design, but when that's not possible, stitching together the baselines of many trials.</div><div><br></div><div>Neither article addresses my specific question of whether it would be even better to log-transform the raw power values before averaging, so I suppose I'll have to test it myself and write my own methods article! 🤓</div><div><br></div><div>I will go ahead and incorporate some of these alternative baseline normalization methods in my git fork of FieldTrip as I go along with my own analyses, so let me know if anyone else would find them useful, and I'll submit a pull request to contribute them to the master branch.</div><div><br></div><div>Thanks for the fruitful discussion, all!</div><div>~Teresa</div><div><br></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Dec 19, 2016 at 8:01 PM, Alik Widge <span dir="ltr"><<a href="mailto:alik.widge@gmail.com" target="_blank">alik.widge@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><p dir="ltr">Indeed, in a separate thread with Michael Cohen several months back he suggested precisely that paper. </p><div class="gmail-HOEnZb"><div class="gmail-h5">
<div class="gmail_extra"><br><div class="gmail_quote">On Dec 19, 2016 5:07 PM, "Nicholas A. Peatfield" <<a href="mailto:nick.peatfield@gmail.com" target="_blank">nick.peatfield@gmail.com</a>> wrote:<br type="attribution"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>I think this paper is relevant to this discussion. </div><div><br></div><div><span style="color:rgb(48,48,48);font-family:arial,helvetica,clean,sans-serif;font-size:13px">Grandchamp, R., & Delorme, A. (2011). Single-Trial Normalization for Event-Related Spectral Decomposition Reduces Sensitivity to Noisy Trials. </span><i style="color:rgb(48,48,48);font-family:arial,helvetica,clean,sans-serif;font-size:13px">Frontiers in Psychology</i><span style="color:rgb(48,48,48);font-family:arial,helvetica,clean,sans-serif;font-size:13px">, </span><i style="color:rgb(48,48,48);font-family:arial,helvetica,clean,sans-serif;font-size:13px">2</i><span style="color:rgb(48,48,48);font-family:arial,helvetica,clean,sans-serif;font-size:13px">, 236. <a href="http://doi.org/10.3389/fpsyg.2011.00236" target="_blank">http://doi.org/10.3389/fpsyg.2<wbr>011.00236</a></span><br></div><div><span style="color:rgb(48,48,48);font-family:arial,helvetica,clean,sans-serif;font-size:13px"><br></span></div><div><a href="https://www.ncbi.nlm.nih.gov/pmc/articles/PMC3183439/" target="_blank">https://www.ncbi.nlm.nih.gov/p<wbr>mc/articles/PMC3183439/</a><br></div><div><br></div><div><br></div><div><div class="gmail_extra"><br><div class="gmail_quote">On 19 December 2016 at 13:08, Teresa Madsen <span dir="ltr"><<a href="mailto:tmadsen@emory.edu" target="_blank">tmadsen@emory.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">I appreciate everyone's feedback, but I still wonder if something is being missed. I understand that the non-normally distributed power values may be less of an issue when performing non-parametric stats or even a paired-samples t-test that looks at difference values which may be normal even when the raw data isn't. However, my concern comes into play even before these statistical comparisons are made, whenever any averaging is done to freq-type data across times, frequencies, trials, electrodes, subjects, etc. That means any time any of these configuration options are used for any of these functions, and probably more:<div><br></div><div>ft_freqanalysis: cfg.keeptrials or cfg.keeptapers = 'no';</div><div>ft_freqgrandaverage: cfg.keepindividual = 'no';</div><div>ft_freqstatistics: cfg.avgoverchan, cfg.avgovertime, or cfg.avgoverfreq = 'yes';</div><div>ft_freqbaseline: cfg.baseline = anything but 'no'</div><div><br></div><div>In each case, if raw power values are averaged, the result will be positively skewed. Maybe it's not a huge problem if all of the data is treated identically, but the specific case that triggered my concern was in ft_freqbaseline, where the individual time-frequency bins are compared to the mean over time for the baseline period. For example, when using cfg.baselinetype = 'db', as Giuseppe Pellizzer suggested, the output freq data does indeed have a more normal distribution over time, but the mean over the baseline time period is performed *before* the log transform, when the distribution is still highly skewed:</div><div><br></div><div><div> meanVals = repmat(nanmean(data(:,:,baseli<wbr>neTimes), 3), [1 1 size(data, 3)]);</div></div><div><div> data = 10*log10(data ./ meanVals);</div></div><div><br></div><div>That's what I had originally done when analyzing data for my SfN poster, when I realized the background noise that shouldn't have changed much from baseline was mostly showing a decrease from baseline of about -3dB.</div><div><br></div><div>Now, I've realized I'm seeing this as more of a problem than others because of another tweak I made, which was to use a long, separate baseline recording to normalize my trial data, rather than a short pre-trial period as ft_freqbaseline is designed to do. Averaging a few hundred milliseconds for a baseline power estimate might be okay because overlapping time points in the original data are used to calculate those power values anyway, probably making them less skewed, but also (it seems to me) more arbitrary and prone to error. I already offered my custom function BLnorm.m to one person who was asking about this issue of normalizing to a separate baseline recording, and I would be happy to contribute it to FieldTrip if others would appreciate it.</div><div><br></div><div>Since a few people suggested using the median, and it is also suggested in <a href="https://mitpress.mit.edu/books/analyzing-neural-time-series-data" target="_blank">Cohen's textbook</a> as an alternative measure of the central tendency for skewed raw power values, I wonder if the simplest fix might be to add an option to select mean or median in each of the functions listed above. Another possibility would be adding an option to transform the power values upon output from ft_freqanalysis.</div><div><br></div><div>Would anyone else find such changes useful?</div><div><br></div><div>Thanks,</div><div>Teresa</div><div><br></div><div><div class="gmail_extra"><div><div class="gmail-m_-5284836939078097878m_5867596177854989507gmail-m_-6018713868620077266m_6572031705234916106m_-1732435530508580156h5"><br><div class="gmail_quote">On Wed, Dec 14, 2016 at 4:22 AM, Herring, J.D. (Jim) <span dir="ltr"><<a href="mailto:J.Herring@donders.ru.nl" target="_blank">J.Herring@donders.ru.nl</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<div lang="EN-US">
<div class="gmail-m_-5284836939078097878m_5867596177854989507gmail-m_-6018713868620077266m_6572031705234916106m_-1732435530508580156m_8485764955201434729gmail-m_-1930378958190960846WordSection1">
<p class="MsoNormal"><span style="font-size:11pt;font-family:calibri,sans-serif">In terms of statistics it is the distribution of values that you do the statistics on that matters. In case of a paired-samples t-test when comparing two conditions, it is
the distribution of difference values that has to be normally distributed. The distribution of difference values is often normal given two similarly non-normal distributions, offering no complications for a regular parametric test.
<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:calibri,sans-serif"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:calibri,sans-serif">The non-parametric tests offered in fieldtrip indeed do not assume normality, so you should have no problem there either.
<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:calibri,sans-serif"><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:11pt;font-family:calibri,sans-serif"><u></u> <u></u></span></p>
<p class="MsoNormal"><b><span style="font-size:11pt;font-family:calibri,sans-serif">From:</span></b><span style="font-size:11pt;font-family:calibri,sans-serif"> <a href="mailto:fieldtrip-bounces@science.ru.nl" target="_blank">fieldtrip-bounces@science.ru.n<wbr>l</a> [mailto:<a href="mailto:fieldtrip-bounces@science.ru.nl" target="_blank">fieldtrip-bounces@scie<wbr>nce.ru.nl</a>]
<b>On Behalf Of </b>Alik Widge<br>
<b>Sent:</b> Tuesday, December 13, 2016 3:10 PM<br>
<b>To:</b> FieldTrip discussion list <<a href="mailto:fieldtrip@science.ru.nl" target="_blank">fieldtrip@science.ru.nl</a>><br>
<b>Subject:</b> Re: [FieldTrip] impact of skewed power distributions on data analysis<u></u><u></u></span></p><div><div class="gmail-m_-5284836939078097878m_5867596177854989507gmail-m_-6018713868620077266m_6572031705234916106m_-1732435530508580156m_8485764955201434729gmail-h5">
<p class="MsoNormal"><u></u> <u></u></p>
<div>
<p class="MsoNormal">In this, Teresa is right and we have observed this in our own EEG data -- depending on one's level of noise and number of trials/patients, the mean can be a very poor estimator of central tendency. My students are still arguing about what
we really want to do with it, but at least one of them has shifted to using the median as a matter of course for baseline normalization.<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><br clear="all">
<u></u><u></u></p>
<div>
<div>
<div>
<p class="MsoNormal" style="margin-bottom:12pt">Alik Widge<br>
<a href="mailto:alik.widge@gmail.com" target="_blank">alik.widge@gmail.com</a><br>
<a href="tel:(206)%20866-5435" value="+12068665435" target="_blank">(206) 866-5435</a><u></u><u></u></p>
</div>
</div>
</div>
<p class="MsoNormal"><u></u> <u></u></p>
<div>
<p class="MsoNormal">On Mon, Dec 12, 2016 at 6:45 PM, Teresa Madsen <<a href="mailto:tmadsen@emory.edu" target="_blank">tmadsen@emory.edu</a>> wrote:<u></u><u></u></p>
<blockquote style="border-top:none;border-right:none;border-bottom:none;border-left:1pt solid rgb(204,204,204);padding:0in 0in 0in 6pt;margin-left:4.8pt;margin-right:0in">
<div>
<p class="MsoNormal">That may very well be true; to be honest, I haven't looked that deeply into the stats offerings yet. However, my plan is to express each electrode's experimental data in terms of change from their respective baseline recordings before attempting
any group averaging or statistical testing, and this problem shows up first in the baseline correction step, where FieldTrip averages raw power over time.
<br>
<span style="color:rgb(136,136,136)"><br>
<span class="gmail-m_-5284836939078097878m_5867596177854989507gmail-m_-6018713868620077266m_6572031705234916106m_-1732435530508580156m_8485764955201434729gmail-m_-1930378958190960846hoenzb">~Teresa </span></span><u></u><u></u></p>
</div>
<div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
<div>
<div>
<p class="MsoNormal">On Mon, Dec 12, 2016 at 4:56 PM Nicholas A. Peatfield <<a href="mailto:nick.peatfield@gmail.com" target="_blank">nick.peatfield@gmail.com</a>> wrote:<u></u><u></u></p>
</div>
<blockquote style="border-top:none;border-right:none;border-bottom:none;border-left:1pt solid rgb(204,204,204);padding:0in 0in 0in 6pt;margin-left:4.8pt;margin-right:0in">
<div>
<p class="MsoNormal">Correct me if I'm wrong, but, if you are using the non-parametric statistics implemented by fieldtrip, the data does not need to be normally distributed. <u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
<div>
<p class="MsoNormal">On 12 December 2016 at 13:39, Teresa Madsen <span class="gmail-m_-5284836939078097878m_5867596177854989507gmail-m_-6018713868620077266m_6572031705234916106m_-1732435530508580156m_8485764955201434729gmail-m_-1930378958190960846m3207961404418016072gmailmsg">
<<a href="mailto:tmadsen@emory.edu" target="_blank">tmadsen@emory.edu</a>></span> wrote:<u></u><u></u></p>
<blockquote style="border-top:none;border-right:none;border-bottom:none;border-left:1pt solid rgb(204,204,204);padding:0in 0in 0in 6pt;margin-left:4.8pt;margin-right:0in">
<div>
<p class="MsoNormal">No, sorry, that's not what I meant, but thanks for giving me the opportunity to clarify. Of course everyone is familiar with the 1/f pattern across frequencies, but the distribution across time (and according to the poster, also across
space), also has an extremely skewed, negative exponential distribution. I probably confused everyone by trying to show too much data in my figure, but each color represents the distribution of power values for a single frequency over time, using a histogram
and a line above with circles at the mean +/- one standard deviation. <br>
<br>
My main point was that the mean is not representative of the central tendency of such an asymmetrical distribution of power values over time. It's even more obvious which is more representative of their actual distributions when I plot e^mean(logpower) on the
raw plot and log(mean(rawpower)) on the log plot, but that made the figure even more busy and confusing.
<br>
<br>
I hope that helps,<br>
Teresa<u></u><u></u></p>
</div>
<div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
<div>
<div>
<p class="MsoNormal">On Mon, Dec 12, 2016 at 3:47 PM Nicholas A. Peatfield <<a href="mailto:nick.peatfield@gmail.com" target="_blank">nick.peatfield@gmail.com</a>> wrote:<u></u><u></u></p>
</div>
<blockquote style="border-top:none;border-right:none;border-bottom:none;border-left:1pt solid rgb(204,204,204);padding:0in 0in 0in 6pt;margin-left:4.8pt;margin-right:0in">
<div>
<div>
<p class="MsoNormal">Hi Teresa,<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<div>
<p class="MsoNormal">I think what you are discussing is the 1/f power scaling of the power spectrum. This is one of the reasons that comparisons are made within a band (i.e. alpha to alpha) and not between bands (i.e. alpha to gamma), as such the assumption
is that within bands there should be a relative change against baseline and this is what the statistics are performed on. That is, baseline correction is assumed to be the mean for a specific frequency and not a mean across frequencies.<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal"> And this leads to another point that when you are selecting a frequency range to do the non-parametric statistics on you should not do 1-64 Hz but break it up based on the bands.<u></u><u></u></p>
</div>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">Hope my interpretation of your point is correct. I sent in individually, as I wanted to ensure I followed your point.<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">Cheers,<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">Nick<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<div>
<p class="MsoNormal">On 12 December 2016 at 08:23, Teresa Madsen <span class="gmail-m_-5284836939078097878m_5867596177854989507gmail-m_-6018713868620077266m_6572031705234916106m_-1732435530508580156m_8485764955201434729gmail-m_-1930378958190960846m3207961404418016072m1027658596765936250m-2095170522474924584gmailmsg">
<<a href="mailto:tmadsen@emory.edu" target="_blank">tmadsen@emory.edu</a>></span> wrote:<u></u><u></u></p>
</div>
</div>
<div>
<div>
<blockquote style="border-top:none;border-right:none;border-bottom:none;border-left:1pt solid rgb(204,204,204);padding:0in 0in 0in 6pt;margin-left:4.8pt;margin-right:0in">
<div>
<div>
<div>
<div>
<p class="MsoNormal"><span class="gmail-m_-5284836939078097878m_5867596177854989507gmail-m_-6018713868620077266m_6572031705234916106m_-1732435530508580156m_8485764955201434729gmail-m_-1930378958190960846m3207961404418016072m1027658596765936250m-2095170522474924584gmailmsg"><span style="font-size:9.5pt">FieldTrippers,</span></span><u></u><u></u></p>
</div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:9.5pt"><u></u> <u></u></span></p>
</div>
<p class="MsoNormal"><span class="gmail-m_-5284836939078097878m_5867596177854989507gmail-m_-6018713868620077266m_6572031705234916106m_-1732435530508580156m_8485764955201434729gmail-m_-1930378958190960846m3207961404418016072m1027658596765936250m-2095170522474924584gmailmsg"><span style="font-size:9.5pt">While analyzing my data for the annual Society for Neuroscience meeting, I developed a concern that was quickly validated
by another poster (full abstract copied and linked below) focusing on the root of the problem: neural oscillatory power is not normally distributed across time, frequency, or space. The specific problem I had encountered was in baseline-correcting my experimental
data, where, regardless of cfg.baselinetype, ft_freqbaseline depends on the mean power over time. However, I found that the distribution of raw power over time is so skewed that the mean was not a reasonable approximation of the central tendency of the baseline
power, so it made most of my experimental data look like it had decreased power compared to baseline. The more I think about it, the more I realize that averaging is everywhere in the way we analyze neural oscillations (across time points, frequency bins,
electrodes, trials, subjects, etc.), and many of the standard statistics people use also rely on assumptions of normality. </span></span><u></u><u></u></p>
<div>
<p class="MsoNormal"><span style="font-size:9.5pt"><u></u> <u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span class="gmail-m_-5284836939078097878m_5867596177854989507gmail-m_-6018713868620077266m_6572031705234916106m_-1732435530508580156m_8485764955201434729gmail-m_-1930378958190960846m3207961404418016072m1027658596765936250m-2095170522474924584gmailmsg"><span style="font-size:9.5pt">The most obvious solution for me was to log transform the data first, as it appears to be fairly log normal, and I always
use log-scale visualizations anyway. Erik Peterson, middle author on the poster, agreed that this would at least "restore (some) symmetry to the error distribution." I used a natural log transform, sort of arbitrarily to differentiate from the standard decibel
transform included in FieldTrip as cfg.baselinetype = 'db'. The following figures compare the 2 distributions across several frequency bands (using power values from a wavelet spectrogram obtained from a baseline LFP recorded in rat prelimbic cortex). The
lines at the top represent the mean +/- one standard deviation for each frequency band, and you can see how those descriptive stats are much more representative of the actual distributions in the log scale.</span></span><span style="font-size:9.5pt"><u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:9.5pt"><u></u> <u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:9.5pt"><img border="0" width="822" height="725" style="width: 8.5625in; height: 7.552in;" id="gmail-m_-5284836939078097878m_5867596177854989507gmail-m_-6018713868620077266m_6572031705234916106m_-1732435530508580156m_8485764955201434729gmail-m_-1930378958190960846_x0000_i1025" src="cid:image001.png@01D255F3.787B5C10"><br>
<span class="gmail-m_-5284836939078097878m_5867596177854989507gmail-m_-6018713868620077266m_6572031705234916106m_-1732435530508580156m_8485764955201434729gmail-m_-1930378958190960846m3207961404418016072m1027658596765936250m-2095170522474924584gmailmsg"></span><u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span class="gmail-m_-5284836939078097878m_5867596177854989507gmail-m_-6018713868620077266m_6572031705234916106m_-1732435530508580156m_8485764955201434729gmail-m_-1930378958190960846m3207961404418016072m1027658596765936250m-2095170522474924584gmailmsg"><span style="font-size:9.5pt">For my analysis, I also calculated a z-score on the log transformed power to assess how my experimental data compared to the
variability of the noise in a long baseline recording from before conditioning, rather than a short pre-trial baseline period, since I find that more informative than any of FieldTrip's built-in baseline types. I'm happy to share the custom functions I wrote
for this if people think it would be a useful addition to FieldTrip. I can also share more about my analysis and/or a copy of the poster, if anyone wants more detail - I just didn't want to make this email too big.</span></span><span style="font-size:9.5pt"><u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:9.5pt"><u></u> <u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:9.5pt">Mostly, I'm just hoping to start some discussion here as to how to address this. I searched the <a href="http://www.fieldtriptoolbox.org/development/zscores" target="_blank">wiki</a>, <a href="https://mailman.science.ru.nl/pipermail/fieldtrip/2006-December/000773.html" target="_blank">listserv</a> <a href="https://mailman.science.ru.nl/pipermail/fieldtrip/2010-March/002718.html" target="_blank">archives</a>,
and <a href="http://bugzilla.fieldtriptoolbox.org/show_bug.cgi?id=1574" target="_blank">bugzilla</a> for anything related and came up with a few topics surrounding normalization and baseline correction, but only skirting this issue. It seems important, so
I want to find out whether others agree with my approach or already have other ways of avoiding the problem, and whether FieldTrip's code needs to be changed or just documentation added, or what?<u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:9.5pt"><u></u> <u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:9.5pt">Thanks for any insights,<u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:9.5pt">Teresa<u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:9.5pt"><u></u> <u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:9.5pt"><u></u> <u></u></span></p>
</div>
<div>
<h1 style="margin:0in 0in 0.0001pt;line-height:21pt"><span style="font-size:18pt;font-family:arial,sans-serif;color:rgb(45,163,189)"><a href="http://www.abstractsonline.com/pp8/#!/4071/presentation/24150" target="_blank">271.03 / LLL17 - Neural oscillatory
power is not Gaussian distributed across time</a><u></u><u></u></span></h1>
<div>
<p class="MsoNormal" style="line-height:15pt"><b><span style="font-size:10.5pt;font-family:arial,sans-serif;color:black">Authors<u></u><u></u></span></b></p>
<p class="MsoNormal" style="margin-right:0in;margin-bottom:3.75pt;margin-left:0.5in;line-height:15pt">
<b><span style="font-size:10.5pt;font-family:arial,sans-serif;color:black">*L. IZHIKEVICH</span></b><span style="font-size:10.5pt;font-family:arial,sans-serif;color:black">, E. PETERSON, B. VOYTEK; <br>
Cognitive Sci., UCSD, San Diego, CA<u></u><u></u></span></p>
<p class="MsoNormal" style="margin-left:7.5pt;line-height:15pt"><b><span style="font-size:10.5pt;font-family:arial,sans-serif;color:black">Disclosures<u></u><u></u></span></b></p>
<p class="MsoNormal" style="margin-right:0in;margin-bottom:3.75pt;margin-left:0.5in;line-height:15pt">
<span style="font-size:10.5pt;font-family:arial,sans-serif;color:black"> <b>L. Izhikevich:</b> None. <b>E. Peterson:</b> None. <b>B. Voytek:</b> None.<u></u><u></u></span></p>
<p class="MsoNormal" style="margin-left:15pt;line-height:15pt"><b><span style="font-size:10.5pt;font-family:arial,sans-serif;color:black">Abstract<u></u><u></u></span></b></p>
<p class="MsoNormal" style="margin-right:0in;margin-bottom:3.75pt;margin-left:0.5in;line-height:15pt">
<span style="font-size:10.5pt;font-family:arial,sans-serif;color:black">Neural oscillations are important in organizing activity across the human brain in healthy cognition, while oscillatory disruptions are linked to numerous disease states. Oscillations
are known to vary by frequency and amplitude across time and between different brain regions; however, this variability has never been well characterized. We examined human and animal EEG, LFP, MEG, and ECoG data from over 100 subjects to analyze the distribution
of power and frequency across time, space and species. We report that between data types, subjects, frequencies, electrodes, and time, an inverse power law, or negative exponential distribution, is present in all recordings. This is contrary to, and not compatible
with, the Gaussian noise assumption made in many digital signal processing techniques. The statistical assumptions underlying common algorithms for power spectral estimation, such as Welch's method, are being violated resulting in non-trivial misestimates
of oscillatory power. Different statistical approaches are warranted.<u></u><u></u></span></p>
</div>
</div>
</div>
<div>
<p class="MsoNormal" style="margin-left:22.5pt"><span style="color:rgb(136,136,136)"><u></u> <u></u></span></p>
</div>
<p class="MsoNormal" style="margin-left:22.5pt"><span class="gmail-m_-5284836939078097878m_5867596177854989507gmail-m_-6018713868620077266m_6572031705234916106m_-1732435530508580156m_8485764955201434729gmail-m_-1930378958190960846m3207961404418016072m1027658596765936250m-2095170522474924584m-2898402136910438041m4827568291528697442gmail-hoenzb"><span style="color:rgb(136,136,136)">--
</span><u></u><u></u></span></p>
<div>
<div>
<div>
<div>
<div>
<div>
<p class="MsoNormal" style="margin-left:22.5pt"><span style="font-size:13.5pt;font-family:garamond,serif;color:blue">Teresa E. Madsen, PhD</span><span style="font-size:13.5pt;color:blue"><br>
</span><span style="font-family:garamond,serif;color:blue">Research Technical Specialist: <i>in vivo
</i>electrophysiology & data analysis</span><u></u><u></u></p>
</div>
<div>
<p class="MsoNormal" style="margin-left:22.5pt"><span style="font-family:garamond,serif;color:blue">Division of Behavioral Neuroscience and Psychiatric Disorders<br>
Yerkes National Primate Research Center</span><span style="color:rgb(136,136,136)"><u></u><u></u></span></p>
<div>
<p class="MsoNormal" style="margin-left:22.5pt"><span style="font-family:garamond,serif;color:blue">Emory University</span><span style="color:rgb(136,136,136)"><u></u><u></u></span></p>
<div>
<p class="MsoNormal" style="margin-left:22.5pt"><span style="font-family:garamond,serif;color:blue">Rainnie Lab, NSB 5233<br>
954 Gatewood Rd. NE<br>
Atlanta, GA 30329</span><span style="color:rgb(136,136,136)"><u></u><u></u></span></p>
</div>
<div>
<div>
<p class="MsoNormal" style="margin-left:22.5pt"><span style="font-family:garamond,serif;color:blue"><a href="tel:(770)%20296-9119" target="_blank">(770) 296-9119</a></span><span style="color:rgb(136,136,136)"><u></u><u></u></span></p>
</div>
</div>
</div>
<div>
<p class="MsoNormal" style="margin-left:22.5pt"><span style="font-family:garamond,serif;color:blue"><a href="mailto:braingirl@gmail.com" target="_blank">braingirl@gmail.com</a></span><span style="color:rgb(136,136,136)"><u></u><u></u></span></p>
</div>
<div>
<div>
<p class="MsoNormal" style="margin-left:22.5pt"><span style="font-family:garamond,serif;color:blue"><a href="https://www.linkedin.com/in/temadsen" target="_blank">https://www.linkedin.com/in/te<wbr>madsen</a><u></u><u></u></span></p>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
<p class="MsoNormal" style="margin-left:22.5pt"><u></u> <u></u></p>
</div>
<p class="MsoNormal" style="margin-left:22.5pt"><u></u> <u></u></p>
</blockquote>
</div>
</div>
<div>
<div>
<blockquote style="border-top:none;border-right:none;border-bottom:none;border-left:1pt solid rgb(204,204,204);padding:0in 0in 0in 6pt;margin-left:4.8pt;margin-right:0in">
<p class="MsoNormal" style="margin-left:22.5pt">______________________________<wbr>_________________<br>
fieldtrip mailing list<br>
<a href="mailto:fieldtrip@donders.ru.nl" target="_blank">fieldtrip@donders.ru.nl</a><br>
<a href="https://mailman.science.ru.nl/mailman/listinfo/fieldtrip" target="_blank">https://mailman.science.ru.nl/<wbr>mailman/listinfo/fieldtrip</a><u></u><u></u></p>
</blockquote>
</div>
</div>
<div>
<p class="MsoNormal" style="margin-left:22.5pt"><br>
<br clear="all">
<u></u><u></u></p>
<div>
<p class="MsoNormal" style="margin-left:22.5pt"><u></u> <u></u></p>
</div>
<p class="MsoNormal" style="margin-left:22.5pt">-- <u></u><u></u></p>
<div>
<div>
<div>
<div>
<div>
<p class="MsoNormal" style="margin-left:22.5pt"><span style="font-size:10pt;font-family:arial,sans-serif">Nicholas Peatfield, PhD</span><u></u><u></u></p>
</div>
<div>
<p class="MsoNormal" style="margin-left:22.5pt"><u></u> <u></u></p>
</div>
</div>
</div>
</div>
</div>
</div>
</blockquote>
</div>
</div>
</div>
</blockquote>
</div>
<p class="MsoNormal" style="margin-left:22.5pt"><br>
<br clear="all">
<u></u><u></u></p>
<div>
<p class="MsoNormal" style="margin-left:22.5pt"><u></u> <u></u></p>
</div>
</div>
<div>
<p class="MsoNormal" style="margin-left:22.5pt">-- <u></u><u></u></p>
<div>
<div>
<div>
<div>
<div>
<p class="MsoNormal" style="margin-left:22.5pt"><span style="font-size:10pt;font-family:arial,sans-serif">Nicholas Peatfield, PhD</span><u></u><u></u></p>
</div>
<div>
<p class="MsoNormal" style="margin-left:22.5pt"><u></u> <u></u></p>
</div>
</div>
</div>
</div>
</div>
</div>
</blockquote>
</div>
</div>
</div>
<p class="MsoNormal" style="margin-left:22.5pt"><br>
______________________________<wbr>_________________<br>
fieldtrip mailing list<br>
<a href="mailto:fieldtrip@donders.ru.nl" target="_blank">fieldtrip@donders.ru.nl</a><br>
<a href="https://mailman.science.ru.nl/mailman/listinfo/fieldtrip" target="_blank">https://mailman.science.ru.nl/<wbr>mailman/listinfo/fieldtrip</a><u></u><u></u></p>
</blockquote>
</div>
<p class="MsoNormal" style="margin-left:22.5pt"><u></u> <u></u></p>
</div>
</div></div></div>
</div>
<br>______________________________<wbr>_________________<br>
fieldtrip mailing list<br>
<a href="mailto:fieldtrip@donders.ru.nl" target="_blank">fieldtrip@donders.ru.nl</a><br>
<a href="https://mailman.science.ru.nl/mailman/listinfo/fieldtrip" rel="noreferrer" target="_blank">https://mailman.science.ru.nl/<wbr>mailman/listinfo/fieldtrip</a><br></blockquote></div><br><br clear="all"><div><br></div>-- <br></div></div><div class="gmail-m_-5284836939078097878m_5867596177854989507gmail-m_-6018713868620077266m_6572031705234916106m_-1732435530508580156m_8485764955201434729gmail_signature"><div dir="ltr"><font color="#0000ff" face="garamond, serif"><font size="4">Teresa E. Madsen, PhD<br></font><span>Division of Behavioral Neuroscience and Psychiatric Disorders<br>Yerkes National Primate Research Center</span></font><span><div><font color="#0000ff" face="garamond, serif">Emory University<br></font><div><font color="#0000ff" face="garamond, serif">Rainnie Lab, NSB 5233<br>954 Gatewood Rd. NE<br>Atlanta, GA 30329</font></div><div><font color="#0000ff" face="garamond, serif"><a href="tel:(770)%20296-9119" value="+17702969119" target="_blank">(770) 296-9119</a></font></div></div></span></div></div>
</div></div></div>
<br>______________________________<wbr>_________________<br>
fieldtrip mailing list<br>
<a href="mailto:fieldtrip@donders.ru.nl" target="_blank">fieldtrip@donders.ru.nl</a><br>
<a href="https://mailman.science.ru.nl/mailman/listinfo/fieldtrip" rel="noreferrer" target="_blank">https://mailman.science.ru.nl/<wbr>mailman/listinfo/fieldtrip</a><br></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail-m_-5284836939078097878m_5867596177854989507gmail-m_-6018713868620077266m_6572031705234916106m_-1732435530508580156gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><font face="arial, helvetica, sans-serif" size="2">Nicholas Peatfield, PhD</font></div><div><br></div></div></div></div></div>
</div></div></div>
<br>______________________________<wbr>_________________<br>
fieldtrip mailing list<br>
<a href="mailto:fieldtrip@donders.ru.nl" target="_blank">fieldtrip@donders.ru.nl</a><br>
<a href="https://mailman.science.ru.nl/mailman/listinfo/fieldtrip" rel="noreferrer" target="_blank">https://mailman.science.ru.nl/<wbr>mailman/listinfo/fieldtrip</a><br></blockquote></div></div>
</div></div><br>______________________________<wbr>_________________<br>
fieldtrip mailing list<br>
<a href="mailto:fieldtrip@donders.ru.nl">fieldtrip@donders.ru.nl</a><br>
<a href="https://mailman.science.ru.nl/mailman/listinfo/fieldtrip" rel="noreferrer" target="_blank">https://mailman.science.ru.nl/<wbr>mailman/listinfo/fieldtrip</a><br></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature"><div dir="ltr"><font color="#0000ff" face="garamond, serif"><font size="4">Teresa E. Madsen, PhD<br></font>Division of Behavioral Neuroscience and Psychiatric Disorders<br>Yerkes National Primate Research Center</font><div><font color="#0000ff" face="garamond, serif">Emory University<br></font><div><font color="#0000ff" face="garamond, serif">Rainnie Lab, NSB 5233<br>954 Gatewood Rd. NE<br>Atlanta, GA 30329</font></div><div><font color="#0000ff" face="garamond, serif">(770) 296-9119</font></div></div></div></div>
</div></div>