<div dir="ltr"><span style="font-size:12.8px">Thanks for the additional feedback, Mike and Matt.</span><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">I agree that the skewness of the power distribution is "trivial" in that it is skewed <i>by definition</i> as a squared value, but you're right that my concern is whether it has a nontrivial impact on the analysis of real data. I do intend to do further simulations to test this and various potential solutions, but I need to get another paper submitted first. I will report my findings back to this listserv when possible.</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">Thanks again,</div><div style="font-size:12.8px">Teresa</div><div><br></div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Jan 13, 2017 at 9:51 AM, Matt Craddock <span dir="ltr"><<a href="mailto:m.p.craddock@leeds.ac.uk" target="_blank">m.p.craddock@leeds.ac.uk</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">Hi all,<br>
<br>
I'm fiddling around with simulations at the moment but not got a lot of free time, and every time I do fiddle with them the more complicated I seem to make them...<br>
<br>
correct me if i'm wrong, but there are a couple of different issues here.<br>
<br>
About the skewness of the distribution - you're not usually comparing two skewed distributions statistically (which even then would probably be ok if they're similary skewed), usually it's a bunch of means from those distributions, and the distribution of the means will often be normally distributed (as will their differences), and thus ok for parametric stats. So the main problem is really whether the mean is a good summary statistic for capturing differences between the underlying distributions or not (similar debates about reaction times) - it's certainly not robust to outliers. Median might be better in some ways, but it's a biased estimator of the population median and might have consequences on stat power - that's one for the simulations at some point.<br>
<br>
Taking the log first won't necessarily help. As the code Michael Cohen provides shows, taking the log can still result in skewed distributions - power is right-skewed but its log is left-skewed. Taking the square root would obviously be better in this case but that's because it just returns it to the original unsquared normal distribution - whether that's true of real data needs checking out. In addition if you take the log you're no longer testing the same thing when you compare the means - you're testing the geometric mean instead of the arithmetic mean. Differences in the arithmetic means across conditions will not necessarily translate into differences in the geometric means and vice versa.<br>
<br>
If you use a subtractive baseline, using the median or mean of power will only change what constant you subtract from the data. If you baseline correct using the average of all trials rather than condition specific baselines, it should make no difference to the stats if you use the mean or the median of baseline power (or indeed if you don't use baseline correction at all). If you test differences across electrodes as well as condition, and use electrode but not condition specific baselines - which seems reasonable, since you might expect bigger baseline differences across than within electrodes, it'll only change the main electrode effect, not the condition effect or the interaction between condition and electrode. If you use condition and electrode specific baselines, then it influences *everything*, and you also no longer know if the differences between conditions/electrodes are due to baseline or post-stimulus differences.<br>
<br>
Using the median over the mean will change the stats if you use divisive baselines, but whether that's a good or bad difference is another question for simulations.<br>
<br>
Cheers,<br>
Matt<div><div class="gmail-h5"><br>
<br>
On 12/01/2017 13:37, Mike X Cohen wrote:<br>
</div></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div><div class="gmail-h5">
Interesting discussion here. I think we should take a step back and<br>
distinguish between trivial and nontrivial causes and consequences for<br>
the skewed distribution. To some extent, the non-normal distribution is<br>
simply the result of defining power as a squared distance -- distances<br>
are always positive and squaring them means big values become really<br>
big. Consider the following:<br>
<br>
d = randn(10000,1); % random numbers<br>
subplot(311), hist(d,500); % their distribution<br>
subplot(312), hist(d.^2,500); % "power" distribution, also try a<br>
log-scaled y-axis<br>
subplot(313), hist(log(d.^2),500); % log-power distribution<br>
<br>
The fact that power values have a power-law-like distribution is<br>
therefore trivial.<br>
<br>
But this leads to two non-trivial questions:<br>
1) Is this distribution meaningful for brain function (beyond simply the<br>
result of taking squared values)? People who study "the log-brain" and<br>
fractal-like (or scale-free) organization of brain function would argue<br>
that these distributions reveal meaningful insights into brain function,<br>
and therefore it is not really an artifact for analyses. In other words,<br>
large values are large for a reason; they might not be outliers that we<br>
should attempt to compress.<br>
<br>
2) Does it matter for real data analysis? I think this is Teresa's<br>
initial concern. I'm inclined to think that it doesn't really matter,<br>
but that's just based on the idea (hope!) that if it did really matter<br>
and the way we do it is wrong, the field would have discovered this a<br>
long time ago. On the other hand, this wouldn't be the first time that<br>
people have gotten things wrong for decades.<br>
<br>
I think the best way to investigate #2 would be with simulated data,<br>
showing that a "true" effect is missed when not first computing<br>
log-power, presumably because the effect was overshadowed by<br>
large-amplitude "noise" (but somehow this would have to be physiological<br>
noise that wouldn't get rejected during data cleaning). In empirical<br>
data, all you'd be able to do is show a difference without knowing the<br>
right answer.<br>
<br>
btw, make sure to be careful with baselining log-power -- any divisions<br>
(e.g., dB or percent change) will be unstable/run off to infinity when<br>
baseline power is close to zero, i.e., raw power is close to 1. The sign<br>
might also get weird. Probably best to use a baseline subtraction.<br>
<br>
Mike<br>
<br>
<br>
<br>
<br>
<br>
------------------------------<wbr>--------------><--------------<wbr>------------------------------<wbr>----<br>
<br>
Thanks for the reference. In return, I also recommend reading Ciuparu and<br>
Mures an's recent rebuttal:<br>
<br>
European Journal of Neuroscience, Vol. 43, pp. 861–869, 2016,<br>
doi:10.1111/ejn.13179 <<a href="http://dx.doi.org/10.1111/ejn.13179" rel="noreferrer" target="_blank">http://dx.doi.org/10.1111/ejn<wbr>.13179</a>><br>
TECHNICAL SPOTLIGHT<br>
Sources of bias in single-trial normalization procedures<br>
<br>
They demonstrate that the positive bias Grandchamp and Delorme warned about<br>
with single-trial baseline normalization was, in fact, due to the<br>
correlated numerators and denominators in each of the baseline<br>
normalization procedures they tested, which was, in turn, caused by the<br>
skewed distributions of baseline power values. They demonstrate that this<br>
bias is reduced by using a much longer baseline, ideally incorporated into<br>
the experimental design, but when that's not possible, stitching together<br>
the baselines of many trials.<br>
<br>
Neither article addresses my specific question of whether it would be even<br>
better to log-transform the raw power values before averaging, so I suppose<br>
I'll have to test it myself and write my own methods article! 🤓<br>
<br>
I will go ahead and incorporate some of these alternative baseline<br>
normalization methods in my git fork of FieldTrip as I go along with my own<br>
analyses, so let me know if anyone else would find them useful, and I'll<br>
submit a pull request to contribute them to the master branch.<br>
<br>
Thanks for the fruitful discussion, all!<br>
~Teresa<br>
<br>
<br>
On Mon, Dec 19, 2016 at 8:01 PM, Alik Widge <alik.widge at <a href="http://gmail.com" rel="noreferrer" target="_blank">gmail.com</a><br></div></div><span class="gmail-">
<<a href="http://gmail.com" rel="noreferrer" target="_blank">http://gmail.com</a>>> wrote:<br>
<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
Indeed, in a separate thread with Michael Cohen several months back he<br>
suggested precisely that paper.<br>
<br>
On Dec 19, 2016 5:07 PM, "Nicholas A. Peatfield" <nick.peatfield at<br>
</blockquote>
</span><a href="http://gmail.com" rel="noreferrer" target="_blank">gmail.com</a> <<a href="http://gmail.com" rel="noreferrer" target="_blank">http://gmail.com</a>>><span class="gmail-"><br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
wrote:<br>
<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
I think this paper is relevant to this discussion.<br>
<br>
Grandchamp, R., & Delorme, A. (2011). Single-Trial Normalization for<br>
Event-Related Spectral Decomposition Reduces Sensitivity to Noisy<br>
</blockquote></blockquote>
Trials. *Frontiers<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
in Psychology*, *2*, 236. <a href="http://doi.org/10.3389/fpsyg.2011.00236" rel="noreferrer" target="_blank">http://doi.org/10.3389/fpsyg.2<wbr>011.00236</a><br>
<br>
<a href="https://www.ncbi.nlm.nih.gov/pmc/articles/PMC3183439/" rel="noreferrer" target="_blank">https://www.ncbi.nlm.nih.gov/p<wbr>mc/articles/PMC3183439/</a><br>
<br>
<br>
<br>
On 19 December 2016 at 13:08, Teresa Madsen <tmadsen at <a href="http://emory.edu" rel="noreferrer" target="_blank">emory.edu</a><br>
</blockquote></blockquote></span><div><div class="gmail-h5">
<<a href="http://emory.edu" rel="noreferrer" target="_blank">http://emory.edu</a>>> wrote:<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
I appreciate everyone's feedback, but I still wonder if something is<br>
being missed. I understand that the non-normally distributed power<br>
</blockquote></blockquote></blockquote>
values<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
may be less of an issue when performing non-parametric stats or even a<br>
paired-samples t-test that looks at difference values which may be<br>
</blockquote></blockquote></blockquote>
normal<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
even when the raw data isn't. However, my concern comes into play even<br>
before these statistical comparisons are made, whenever any averaging is<br>
done to freq-type data across times, frequencies, trials, electrodes,<br>
subjects, etc. That means any time any of these configuration<br>
</blockquote></blockquote></blockquote>
options are<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
used for any of these functions, and probably more:<br>
<br>
ft_freqanalysis: cfg.keeptrials or cfg.keeptapers = 'no';<br>
ft_freqgrandaverage: cfg.keepindividual = 'no';<br>
ft_freqstatistics: cfg.avgoverchan, cfg.avgovertime, or<br>
cfg.avgoverfreq = 'yes';<br>
ft_freqbaseline: cfg.baseline = anything but 'no'<br>
<br>
In each case, if raw power values are averaged, the result will be<br>
positively skewed. Maybe it's not a huge problem if all of the data is<br>
treated identically, but the specific case that triggered my concern<br>
</blockquote></blockquote></blockquote>
was in<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
ft_freqbaseline, where the individual time-frequency bins are<br>
</blockquote></blockquote></blockquote>
compared to<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
the mean over time for the baseline period. For example, when using<br>
cfg.baselinetype = 'db', as Giuseppe Pellizzer suggested, the output<br>
</blockquote></blockquote></blockquote>
freq<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
data does indeed have a more normal distribution over time, but the mean<br>
over the baseline time period is performed *before* the log<br>
</blockquote></blockquote></blockquote>
transform, when<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
the distribution is still highly skewed:<br>
<br>
meanVals = repmat(nanmean(data(:,:,baseli<wbr>neTimes), 3), [1 1<br>
size(data, 3)]);<br>
data = 10*log10(data ./ meanVals);<br>
<br>
That's what I had originally done when analyzing data for my SfN poster,<br>
when I realized the background noise that shouldn't have changed<br>
</blockquote></blockquote></blockquote>
much from<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
baseline was mostly showing a decrease from baseline of about -3dB.<br>
<br>
Now, I've realized I'm seeing this as more of a problem than others<br>
because of another tweak I made, which was to use a long, separate<br>
</blockquote></blockquote></blockquote>
baseline<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
recording to normalize my trial data, rather than a short pre-trial<br>
</blockquote></blockquote></blockquote>
period<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
as ft_freqbaseline is designed to do. Averaging a few hundred<br>
</blockquote></blockquote></blockquote>
milliseconds<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
for a baseline power estimate might be okay because overlapping time<br>
</blockquote></blockquote></blockquote>
points<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
in the original data are used to calculate those power values anyway,<br>
probably making them less skewed, but also (it seems to me) more<br>
</blockquote></blockquote></blockquote>
arbitrary<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
and prone to error. I already offered my custom function BLnorm.m<br>
</blockquote></blockquote></blockquote>
to one<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
person who was asking about this issue of normalizing to a separate<br>
baseline recording, and I would be happy to contribute it to<br>
</blockquote></blockquote></blockquote>
FieldTrip if<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
others would appreciate it.<br>
<br>
Since a few people suggested using the median, and it is also suggested<br>
in Cohen's textbook<br>
<<a href="https://mitpress.mit.edu/books/analyzing-neural-time-series-data" rel="noreferrer" target="_blank">https://mitpress.mit.edu/book<wbr>s/analyzing-neural-time-series<wbr>-data</a>> as<br>
an alternative measure of the central tendency for skewed raw power<br>
</blockquote></blockquote></blockquote>
values,<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
I wonder if the simplest fix might be to add an option to select mean or<br>
median in each of the functions listed above. Another possibility<br>
</blockquote></blockquote></blockquote>
would be<br>
</div></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div><div class="gmail-h5">
adding an option to transform the power values upon output from<br>
ft_freqanalysis.<br>
<br>
Would anyone else find such changes useful?<br>
<br>
Thanks,<br>
Teresa<br>
<br>
<br>
On Wed, Dec 14, 2016 at 4:22 AM, Herring, J.D. (Jim) <<br></div></div><span class="gmail-">
J.Herring at <a href="http://donders.ru.nl" rel="noreferrer" target="_blank">donders.ru.nl</a> <<a href="http://donders.ru.nl" rel="noreferrer" target="_blank">http://donders.ru.nl</a>>> wrote:<br>
<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
In terms of statistics it is the distribution of values that you do the<br>
statistics on that matters. In case of a paired-samples t-test when<br>
comparing two conditions, it is the distribution of difference<br>
</blockquote></span></blockquote></blockquote></blockquote><span class="gmail-">
values that<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
has to be normally distributed. The distribution of difference<br>
</blockquote></blockquote></blockquote></blockquote>
values is<br>
</span><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><span class="gmail-">
often normal given two similarly non-normal distributions, offering no<br>
complications for a regular parametric test.<br>
<br>
<br>
<br>
The non-parametric tests offered in fieldtrip indeed do not assume<br>
normality, so you should have no problem there either.<br>
<br>
<br>
<br>
<br>
<br></span>
*From:* fieldtrip-bounces at <a href="http://science.ru.nl" rel="noreferrer" target="_blank">science.ru.nl</a> <<a href="http://science.ru.nl" rel="noreferrer" target="_blank">http://science.ru.nl</a>><br>
</blockquote></blockquote></blockquote></blockquote>
[mailto:<a href="mailto:fieldtrip-bounces" target="_blank">fieldtrip-bounces</a> <mailto:<a href="mailto:fieldtrip-bounces" target="_blank">fieldtrip-bounces</a>> at scie<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<a href="http://nce.ru.nl" rel="noreferrer" target="_blank">nce.ru.nl</a> <<a href="http://nce.ru.nl" rel="noreferrer" target="_blank">http://nce.ru.nl</a>>] *On Behalf Of *Alik Widge<span class="gmail-"><br>
*Sent:* Tuesday, December 13, 2016 3:10 PM<br>
*To:* FieldTrip discussion list <fieldtrip at <a href="http://science.ru.nl" rel="noreferrer" target="_blank">science.ru.nl</a><br>
</span></blockquote></blockquote></blockquote></blockquote>
<<a href="http://science.ru.nl" rel="noreferrer" target="_blank">http://science.ru.nl</a>>><span class="gmail-"><br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
*Subject:* Re: [FieldTrip] impact of skewed power distributions on<br>
data analysis<br>
<br>
<br>
<br>
In this, Teresa is right and we have observed this in our own EEG data<br>
-- depending on one's level of noise and number of trials/patients, the<br>
mean can be a very poor estimator of central tendency. My students are<br>
still arguing about what we really want to do with it, but at least<br>
</blockquote></blockquote></blockquote></blockquote>
one of<br>
</span><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><span class="gmail-">
them has shifted to using the median as a matter of course for baseline<br>
normalization.<br>
<br>
<br>
Alik Widge<br></span>
alik.widge at <a href="http://gmail.com" rel="noreferrer" target="_blank">gmail.com</a> <<a href="http://gmail.com" rel="noreferrer" target="_blank">http://gmail.com</a>><span class="gmail-"><br>
<a href="tel:%28206%29%20866-5435" value="+12068665435" target="_blank">(206) 866-5435</a><br>
<br>
<br>
<br>
On Mon, Dec 12, 2016 at 6:45 PM, Teresa Madsen <tmadsen at<br>
</span></blockquote></blockquote></blockquote></blockquote>
<a href="http://emory.edu" rel="noreferrer" target="_blank">emory.edu</a> <<a href="http://emory.edu" rel="noreferrer" target="_blank">http://emory.edu</a>>><span class="gmail-"><br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
wrote:<br>
<br>
That may very well be true; to be honest, I haven't looked that deeply<br>
into the stats offerings yet. However, my plan is to express each<br>
electrode's experimental data in terms of change from their respective<br>
baseline recordings before attempting any group averaging or<br>
</blockquote></blockquote></blockquote></blockquote>
statistical<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
testing, and this problem shows up first in the baseline correction<br>
</blockquote></blockquote></blockquote></blockquote>
step,<br>
</span><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><span class="gmail-">
where FieldTrip averages raw power over time.<br>
<br>
~Teresa<br>
<br>
<br>
<br>
On Mon, Dec 12, 2016 at 4:56 PM Nicholas A. Peatfield <<br></span><span class="gmail-">
nick.peatfield at <a href="http://gmail.com" rel="noreferrer" target="_blank">gmail.com</a> <<a href="http://gmail.com" rel="noreferrer" target="_blank">http://gmail.com</a>>> wrote:<br>
<br>
Correct me if I'm wrong, but, if you are using the non-parametric<br>
statistics implemented by fieldtrip, the data does not need to be<br>
</span></blockquote></blockquote></blockquote></blockquote><span class="gmail-">
normally<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
distributed.<br>
<br>
<br>
<br>
On 12 December 2016 at 13:39, Teresa Madsen <tmadsen at <a href="http://emory.edu" rel="noreferrer" target="_blank">emory.edu</a><br>
</blockquote></blockquote></blockquote></blockquote></span><span class="gmail-">
<<a href="http://emory.edu" rel="noreferrer" target="_blank">http://emory.edu</a>>> wrote:<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<br>
No, sorry, that's not what I meant, but thanks for giving me the<br>
opportunity to clarify. Of course everyone is familiar with the 1/f<br>
</blockquote></blockquote></blockquote></blockquote>
pattern<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
across frequencies, but the distribution across time (and according<br>
</blockquote></blockquote></blockquote></blockquote>
to the<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
poster, also across space), also has an extremely skewed, negative<br>
exponential distribution. I probably confused everyone by trying to<br>
</blockquote></blockquote></blockquote></blockquote>
show<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
too much data in my figure, but each color represents the<br>
</blockquote></blockquote></blockquote></blockquote>
distribution of<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
power values for a single frequency over time, using a histogram<br>
</blockquote></blockquote></blockquote></blockquote>
and a line<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
above with circles at the mean +/- one standard deviation.<br>
<br>
My main point was that the mean is not representative of the central<br>
tendency of such an asymmetrical distribution of power values over<br>
</blockquote></blockquote></blockquote></blockquote>
time.<br>
</span><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><span class="gmail-">
It's even more obvious which is more representative of their actual<br>
distributions when I plot e^mean(logpower) on the raw plot and<br>
log(mean(rawpower)) on the log plot, but that made the figure even more<br>
busy and confusing.<br>
<br>
I hope that helps,<br>
Teresa<br>
<br>
<br>
<br>
On Mon, Dec 12, 2016 at 3:47 PM Nicholas A. Peatfield <<br></span><span class="gmail-">
nick.peatfield at <a href="http://gmail.com" rel="noreferrer" target="_blank">gmail.com</a> <<a href="http://gmail.com" rel="noreferrer" target="_blank">http://gmail.com</a>>> wrote:<br>
<br>
Hi Teresa,<br>
<br>
<br>
<br>
I think what you are discussing is the 1/f power scaling of the power<br>
spectrum. This is one of the reasons that comparisons are made within<br>
a band (i.e. alpha to alpha) and not between bands (i.e. alpha to<br>
</span></blockquote></blockquote></blockquote></blockquote><span class="gmail-">
gamma),<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
as such the assumption is that within bands there should be a relative<br>
change against baseline and this is what the statistics are<br>
</blockquote></blockquote></blockquote></blockquote>
performed on.<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
That is, baseline correction is assumed to be the mean for a specific<br>
frequency and not a mean across frequencies.<br>
<br>
<br>
<br>
And this leads to another point that when you are selecting a<br>
frequency range to do the non-parametric statistics on you should<br>
</blockquote></blockquote></blockquote></blockquote>
not do<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
1-64 Hz but break it up based on the bands.<br>
<br>
<br>
<br>
Hope my interpretation of your point is correct. I sent in<br>
individually, as I wanted to ensure I followed your point.<br>
<br>
<br>
<br>
Cheers,<br>
<br>
<br>
<br>
Nick<br>
<br>
<br>
<br>
<br>
<br>
On 12 December 2016 at 08:23, Teresa Madsen <tmadsen at <a href="http://emory.edu" rel="noreferrer" target="_blank">emory.edu</a><br>
</blockquote></blockquote></blockquote></blockquote></span><div><div class="gmail-h5">
<<a href="http://emory.edu" rel="noreferrer" target="_blank">http://emory.edu</a>>> wrote:<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<br>
FieldTrippers,<br>
<br>
<br>
<br>
While analyzing my data for the annual Society for Neuroscience<br>
meeting, I developed a concern that was quickly validated by<br>
</blockquote></blockquote></blockquote></blockquote>
another poster<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
(full abstract copied and linked below) focusing on the root of the<br>
problem: neural oscillatory power is not normally distributed<br>
</blockquote></blockquote></blockquote></blockquote>
across time,<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
frequency, or space. The specific problem I had encountered was in<br>
baseline-correcting my experimental data, where, regardless of<br>
cfg.baselinetype, ft_freqbaseline depends on the mean power over time.<br>
However, I found that the distribution of raw power over time is so<br>
</blockquote></blockquote></blockquote></blockquote>
skewed<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
that the mean was not a reasonable approximation of the central<br>
</blockquote></blockquote></blockquote></blockquote>
tendency of<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
the baseline power, so it made most of my experimental data look<br>
</blockquote></blockquote></blockquote></blockquote>
like it<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
had decreased power compared to baseline. The more I think about<br>
</blockquote></blockquote></blockquote></blockquote>
it, the<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
more I realize that averaging is everywhere in the way we analyze<br>
</blockquote></blockquote></blockquote></blockquote>
neural<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
oscillations (across time points, frequency bins, electrodes, trials,<br>
subjects, etc.), and many of the standard statistics people use<br>
</blockquote></blockquote></blockquote></blockquote>
also rely<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
on assumptions of normality.<br>
<br>
<br>
<br>
The most obvious solution for me was to log transform the data first,<br>
as it appears to be fairly log normal, and I always use log-scale<br>
visualizations anyway. Erik Peterson, middle author on the poster,<br>
</blockquote></blockquote></blockquote></blockquote>
agreed<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
that this would at least "restore (some) symmetry to the error<br>
distribution." I used a natural log transform, sort of arbitrarily to<br>
differentiate from the standard decibel transform included in<br>
</blockquote></blockquote></blockquote></blockquote>
FieldTrip as<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
cfg.baselinetype = 'db'. The following figures compare the 2<br>
</blockquote></blockquote></blockquote></blockquote>
distributions<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
across several frequency bands (using power values from a wavelet<br>
spectrogram obtained from a baseline LFP recorded in rat prelimbic<br>
cortex). The lines at the top represent the mean +/- one standard<br>
deviation for each frequency band, and you can see how those<br>
</blockquote></blockquote></blockquote></blockquote>
descriptive<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
stats are much more representative of the actual distributions in<br>
</blockquote></blockquote></blockquote></blockquote>
the log<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
scale.<br>
<br>
<br>
<br>
<br>
<br>
<br>
For my analysis, I also calculated a z-score on the log transformed<br>
power to assess how my experimental data compared to the<br>
</blockquote></blockquote></blockquote></blockquote>
variability of the<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
noise in a long baseline recording from before conditioning, rather<br>
</blockquote></blockquote></blockquote></blockquote>
than a<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
short pre-trial baseline period, since I find that more informative<br>
</blockquote></blockquote></blockquote></blockquote>
than<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
any of FieldTrip's built-in baseline types. I'm happy to share the<br>
</blockquote></blockquote></blockquote></blockquote>
custom<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
functions I wrote for this if people think it would be a useful<br>
</blockquote></blockquote></blockquote></blockquote>
addition to<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
FieldTrip. I can also share more about my analysis and/or a copy<br>
</blockquote></blockquote></blockquote></blockquote>
of the<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
poster, if anyone wants more detail - I just didn't want to make<br>
</blockquote></blockquote></blockquote></blockquote>
this email<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
too big.<br>
<br>
<br>
<br>
Mostly, I'm just hoping to start some discussion here as to how to<br>
address this. I searched the wiki<br>
<<a href="http://www.fieldtriptoolbox.org/development/zscores" rel="noreferrer" target="_blank">http://www.fieldtriptoolbox.o<wbr>rg/development/zscores</a>>, listserv<br>
<br>
</blockquote></blockquote></blockquote></blockquote>
<<a href="https://mailman.science.ru.nl/pipermail/fieldtrip/2006-December/000773.html" rel="noreferrer" target="_blank">https://mailman.science.ru.nl<wbr>/pipermail/fieldtrip/2006-Dece<wbr>mber/000773.html</a>><br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
archives<br>
<br>
</blockquote></blockquote></blockquote></blockquote>
<<a href="https://mailman.science.ru.nl/pipermail/fieldtrip/2010-March/002718.html" rel="noreferrer" target="_blank">https://mailman.science.ru.nl<wbr>/pipermail/fieldtrip/2010-Marc<wbr>h/002718.html</a>>,<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
and bugzilla<br>
<<a href="http://bugzilla.fieldtriptoolbox.org/show_bug.cgi?id=1574" rel="noreferrer" target="_blank">http://bugzilla.fieldtriptool<wbr>box.org/show_bug.cgi?id=1574</a>> for<br>
anything related and came up with a few topics surrounding<br>
</blockquote></blockquote></blockquote></blockquote>
normalization<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
and baseline correction, but only skirting this issue. It seems<br>
</blockquote></blockquote></blockquote></blockquote>
important,<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
so I want to find out whether others agree with my approach or<br>
</blockquote></blockquote></blockquote></blockquote>
already have<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
other ways of avoiding the problem, and whether FieldTrip's code<br>
</blockquote></blockquote></blockquote></blockquote>
needs to<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
be changed or just documentation added, or what?<br>
<br>
<br>
<br>
Thanks for any insights,<br>
<br>
Teresa<br>
<br>
<br>
<br>
<br>
271.03 / LLL17 - Neural oscillatory power is not Gaussian distributed<br>
across time<br>
<<a href="http://www.abstractsonline.com/pp8/#!/4071/presentation/24150" rel="noreferrer" target="_blank">http://www.abstractsonline.co<wbr>m/pp8/#!/4071/presentation/241<wbr>50</a>><br>
<br>
*Authors*<br>
<br>
**L. IZHIKEVICH*, E. PETERSON, B. VOYTEK;<br>
Cognitive Sci., UCSD, San Diego, CA<br>
<br>
*Disclosures*<br>
<br>
*L. Izhikevich:* None. *E. Peterson:* None. *B. Voytek:* None.<br>
<br>
*Abstract*<br>
<br>
Neural oscillations are important in organizing activity across the<br>
human brain in healthy cognition, while oscillatory disruptions are<br>
</blockquote></blockquote></blockquote></blockquote>
linked<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
to numerous disease states. Oscillations are known to vary by<br>
</blockquote></blockquote></blockquote></blockquote>
frequency and<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
amplitude across time and between different brain regions; however,<br>
</blockquote></blockquote></blockquote></blockquote>
this<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
variability has never been well characterized. We examined human<br>
</blockquote></blockquote></blockquote></blockquote>
and animal<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
EEG, LFP, MEG, and ECoG data from over 100 subjects to analyze the<br>
distribution of power and frequency across time, space and species. We<br>
report that between data types, subjects, frequencies, electrodes, and<br>
time, an inverse power law, or negative exponential distribution, is<br>
present in all recordings. This is contrary to, and not compatible<br>
</blockquote></blockquote></blockquote></blockquote>
with,<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
the Gaussian noise assumption made in many digital signal processing<br>
techniques. The statistical assumptions underlying common<br>
</blockquote></blockquote></blockquote></blockquote>
algorithms for<br>
</div></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div><div class="gmail-h5">
power spectral estimation, such as Welch's method, are being violated<br>
resulting in non-trivial misestimates of oscillatory power. Different<br>
statistical approaches are warranted.<br>
<br>
<br>
<br>
--<br>
<br>
Teresa E. Madsen, PhD<br>
Research Technical Specialist: *in vivo *electrophysiology & data<br>
analysis<br>
<br>
Division of Behavioral Neuroscience and Psychiatric Disorders<br>
Yerkes National Primate Research Center<br>
<br>
Emory University<br>
<br>
Rainnie Lab, NSB 5233<br>
954 Gatewood Rd. NE<br>
Atlanta, GA 30329<br>
<br>
<a href="tel:%28770%29%20296-9119" value="+17702969119" target="_blank">(770) 296-9119</a><br>
<br></div></div>
braingirl at <a href="http://gmail.com" rel="noreferrer" target="_blank">gmail.com</a> <<a href="http://gmail.com" rel="noreferrer" target="_blank">http://gmail.com</a>><span class="gmail-"><br>
<br>
<a href="https://www.linkedin.com/in/temadsen" rel="noreferrer" target="_blank">https://www.linkedin.com/in/te<wbr>madsen</a><br>
<br>
<br>
<br>
<br>
<br>
______________________________<wbr>_________________<br>
fieldtrip mailing list<br></span>
fieldtrip at <a href="http://donders.ru.nl" rel="noreferrer" target="_blank">donders.ru.nl</a> <<a href="http://donders.ru.nl" rel="noreferrer" target="_blank">http://donders.ru.nl</a>><span class="gmail-"><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>
<br>
<br>
<br>
<br>
<br>
--<br>
<br>
Nicholas Peatfield, PhD<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
--<br>
<br>
Nicholas Peatfield, PhD<br>
<br>
<br>
<br>
<br>
______________________________<wbr>_________________<br>
fieldtrip mailing list<br></span>
fieldtrip at <a href="http://donders.ru.nl" rel="noreferrer" target="_blank">donders.ru.nl</a> <<a href="http://donders.ru.nl" rel="noreferrer" target="_blank">http://donders.ru.nl</a>><span class="gmail-"><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>
<br>
<br>
<br>
______________________________<wbr>_________________<br>
fieldtrip mailing list<br></span>
fieldtrip at <a href="http://donders.ru.nl" rel="noreferrer" target="_blank">donders.ru.nl</a> <<a href="http://donders.ru.nl" rel="noreferrer" target="_blank">http://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>
<br>
</blockquote><span class="gmail-">
<br>
<br>
<br>
--<br>
Teresa E. Madsen, PhD<br>
Division of Behavioral Neuroscience and Psychiatric Disorders<br>
Yerkes National Primate Research Center<br>
Emory University<br>
Rainnie Lab, NSB 5233<br>
954 Gatewood Rd. NE<br>
Atlanta, GA 30329<br>
<a href="tel:%28770%29%20296-9119" value="+17702969119" target="_blank">(770) 296-9119</a><br>
<br>
______________________________<wbr>_________________<br>
fieldtrip mailing list<br></span>
fieldtrip at <a href="http://donders.ru.nl" rel="noreferrer" target="_blank">donders.ru.nl</a> <<a href="http://donders.ru.nl" rel="noreferrer" target="_blank">http://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>
<br>
</blockquote><span class="gmail-">
<br>
<br>
<br>
--<br>
Nicholas Peatfield, PhD<br>
<br>
<br>
______________________________<wbr>_________________<br>
fieldtrip mailing list<br></span>
fieldtrip at <a href="http://donders.ru.nl" rel="noreferrer" target="_blank">donders.ru.nl</a> <<a href="http://donders.ru.nl" rel="noreferrer" target="_blank">http://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>
<br>
</blockquote>
<br>
______________________________<wbr>_________________<br>
fieldtrip mailing list<br>
fieldtrip at <a href="http://donders.ru.nl" rel="noreferrer" target="_blank">donders.ru.nl</a> <<a href="http://donders.ru.nl" rel="noreferrer" target="_blank">http://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>
<br>
</blockquote><span class="gmail-">
<br>
<br>
<br>
--<br>
Teresa E. Madsen, PhD<br>
Division of Behavioral Neuroscience and Psychiatric Disorders<br>
Yerkes National Primate Research Center<br>
Emory University<br>
Rainnie Lab, NSB 5233<br>
954 Gatewood Rd. NE<br>
Atlanta, GA 30329<br>
<a href="tel:%28770%29%20296-9119" value="+17702969119" target="_blank">(770) 296-9119</a><br>
<br>
--<br>
Mike X Cohen, PhD<br>
</span><a href="http://mikexcohen.com" rel="noreferrer" target="_blank">mikexcohen.com</a> <<a href="http://mikexcohen.com" rel="noreferrer" target="_blank">http://mikexcohen.com</a>><span class="gmail-"><br>
<br>
<br>
______________________________<wbr>_________________<br>
fieldtrip mailing list<br>
<a href="mailto:fieldtrip@donders.ru.nl" target="_blank">fieldtrip@donders.ru.nl</a><br>
</span><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>
<br>
</blockquote>
______________________________<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></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>