<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
<meta name="Generator" content="Microsoft Word 15 (filtered medium)">
<style><!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0in;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri",sans-serif;}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:#0563C1;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:#954F72;
text-decoration:underline;}
p.MsoPlainText, li.MsoPlainText, div.MsoPlainText
{mso-style-priority:99;
mso-style-link:"Plain Text Char";
margin:0in;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri",sans-serif;}
span.PlainTextChar
{mso-style-name:"Plain Text Char";
mso-style-priority:99;
mso-style-link:"Plain Text";
font-family:"Calibri",sans-serif;}
.MsoChpDefault
{mso-style-type:export-only;
font-family:"Calibri",sans-serif;}
@page WordSection1
{size:8.5in 11.0in;
margin:1.0in 1.0in 1.0in 1.0in;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang="EN-US" link="#0563C1" vlink="#954F72">
<div class="WordSection1">
<p class="MsoPlainText"><a name="_MailEndCompose">, <o:p></o:p></a></p>
<p class="MsoPlainText"><span style="mso-bookmark:_MailEndCompose">1- I have done some research on the effect of ICA on brain connectivity in channel space two years ago. I have attached the results for your information. I found that applying ICA doesn’t
change the brain connectivity ( coherence and its imaginary part which is sensitive to phase). The analysis has been done using BrainStorm.<o:p></o:p></span></p>
<p class="MsoPlainText"><span style="mso-bookmark:_MailEndCompose">2- I would like to see more scientific evidences/references which are published in peer-reviewed journals (than conference papers) for statements such as “ However, superposition does not
hold when one uses a smaller number of ICA components to create a larger number of channels, for example, use 16 or 17 ICA components to create 19 channels”.<o:p></o:p></span></p>
<p class="MsoPlainText"><span style="mso-bookmark:_MailEndCompose">3- How do you know what is the “true” phase for instantaneous EEG considering the fact that it is a chaotic/stochastic signal. What is the meaning of “real phase” which we should preserve
it? You can do a very simple study. Generate some sinusoid signals with different phase delays. Add Gaussian noise to some of them. Apply ICA and remove the noisy component(s). Reconstruct the data and check out the phase.
<o:p></o:p></span></p>
<p class="MsoPlainText"><span style="mso-bookmark:_MailEndCompose">4- ICA has been used in audio processing for a long time. Suppose that there is a cocktail party and I am talking to you and I said “ Hello Rober” and you said, “ Oh Hi Iman” . Then in the case
of “huge” phase distortion after cleaning with ICA if someone was spying on us he would hear “Oh! Hi Iman … Hello, Robert”. I am exaggerating the case but this is the meaning of huge phase difference ( which in time domain means delay).
<o:p></o:p></span></p>
<p class="MsoPlainText"><span style="mso-bookmark:_MailEndCompose"> <o:p></o:p></span></p>
<p class="MsoPlainText"><span style="mso-bookmark:_MailEndCompose">Regards, <o:p>
</o:p></span></p>
<p class="MsoPlainText"><span style="mso-bookmark:_MailEndCompose">Iman<o:p></o:p></span></p>
<p class="MsoPlainText"><span style="mso-bookmark:_MailEndCompose"><o:p> </o:p></span></p>
<p class="MsoPlainText"><span style="mso-bookmark:_MailEndCompose">Iman Rezazadeh, Ph.D.<o:p></o:p></span></p>
<p class="MsoPlainText"><span style="mso-bookmark:_MailEndCompose">Research Scientist
<o:p></o:p></span></p>
<p class="MsoPlainText"><span style="mso-bookmark:_MailEndCompose">HRL Laboratories and UCLA David Geffen School of Medicine<o:p></o:p></span></p>
<p class="MsoPlainText"><span style="mso-bookmark:_MailEndCompose">Email: irezazadeh@mednet.ucla.edu<o:p></o:p></span></p>
<p class="MsoPlainText"><span style="mso-bookmark:_MailEndCompose"><o:p> </o:p></span></p>
<p class="MsoPlainText"><span style="mso-bookmark:_MailEndCompose"><o:p> </o:p></span></p>
<span style="mso-bookmark:_MailEndCompose"></span>
<p class="MsoPlainText">-----Original Message-----<br>
From: eeglablist-bounces@sccn.ucsd.edu [mailto:eeglablist-bounces@sccn.ucsd.edu] On Behalf Of Ramesh Srinivasan<br>
Sent: Saturday, June 17, 2017 1:28 PM<br>
To: Stefan Debener <stefan.debener@uni-oldenburg.de>; Arnaud Delorme <arno@UCSD.EDU>; Robert Thatcher <rwthatcher2@yahoo.com><br>
Cc: eeglablist <eeglablist@sccn.ucsd.edu>; pnunez@tulane.edu; r.srinivasan@uci.edu; Georges Otte <georges.otte@telenet.be><br>
Subject: Re: [Eeglablist] ICA Misinformation</p>
<p class="MsoPlainText"><o:p> </o:p></p>
<p class="MsoPlainText">Hi All -<o:p></o:p></p>
<p class="MsoPlainText"><o:p> </o:p></p>
<p class="MsoPlainText">Stefan's very first point is the heart of the matter, and he has expressed it much more clearly than me.<o:p></o:p></p>
<p class="MsoPlainText"><o:p> </o:p></p>
<p class="MsoPlainText">"We should not forget that many physiological artifacts are more or less continuous in nature. ECG artifacts for instance are usually not visible in the raw EEG, but they clearly exist - for as long as the subject your are recording
from is as alive. So, the important implication is that non-brain artifacts near continuoulsy contaminate the recording, even if you don't see them! Just because they are not visible does not mean that they don't exist."<o:p></o:p></p>
<p class="MsoPlainText"><o:p> </o:p></p>
<p class="MsoPlainText">There is no artifact-free EEG.<o:p></o:p></p>
<p class="MsoPlainText"><o:p> </o:p></p>
<p class="MsoPlainText">A natural consequence is that there is no actual gold standard test to compare artifact editing schemes that is meaningful. All of the approaches suggested as alternatives to ICA also have built-in assumptions, selection bias (for manual
editing), and almost certainly remove real brain data along with the artifact.<o:p></o:p></p>
<p class="MsoPlainText"><o:p> </o:p></p>
<p class="MsoPlainText">No EEG recording will ever be truly "clean", and it tells me that the focus in our work ought to be identifying robust effects. If you have a really worthwhile effect, then it should not be very sensitive to the artifact editing scheme.
Or for that matter data analysis choices, e.g., morlet wavelets versus FFT based spectrograms.<o:p></o:p></p>
<p class="MsoPlainText"><o:p> </o:p></p>
<p class="MsoPlainText">ramesh<o:p></o:p></p>
<p class="MsoPlainText"><o:p> </o:p></p>
<p class="MsoPlainText"><o:p> </o:p></p>
<p class="MsoPlainText">On 06/17/2017 12:36 PM, Stefan Debener wrote:<o:p></o:p></p>
<p class="MsoPlainText">> Hi there,<o:p></o:p></p>
<p class="MsoPlainText">><o:p> </o:p></p>
<p class="MsoPlainText">> I fully agree with Arno. Robert, you may consider the following thoughts:<o:p></o:p></p>
<p class="MsoPlainText">><o:p> </o:p></p>
<p class="MsoPlainText">> - We should not forget that many physiological artifacts are more or
<o:p></o:p></p>
<p class="MsoPlainText">> less continuous in nature. ECG artifacts for instance are usually not
<o:p></o:p></p>
<p class="MsoPlainText">> visible in the raw EEG, but they clearly exist - for as long as the
<o:p></o:p></p>
<p class="MsoPlainText">> subject your are recording from is as alive. So, the important
<o:p></o:p></p>
<p class="MsoPlainText">> implication is that non-brain artifacts near continuoulsy contaminate
<o:p></o:p></p>
<p class="MsoPlainText">> the recording, even if you don't see them! Just because they are not
<o:p></o:p></p>
<p class="MsoPlainText">> visible does not mean that they don't exist. This probably also holds
<o:p></o:p></p>
<p class="MsoPlainText">> for eye-related artifacts (among which many may contribute small
<o:p></o:p></p>
<p class="MsoPlainText">> voltage, but seem near always present, like microstates with eyes
<o:p></o:p></p>
<p class="MsoPlainText">> open, or eye ball rotations with eyes closed, etc...) and for
<o:p></o:p></p>
<p class="MsoPlainText">> muscle-related artifacts (e.g., neck muscles). They are not simply on
<o:p></o:p></p>
<p class="MsoPlainText">> or off, but more or less active, thus they contribute funny activation
<o:p></o:p></p>
<p class="MsoPlainText">> patterns to the EEG. It follows that an artifact-rejection approch
<o:p></o:p></p>
<p class="MsoPlainText">> alone may be misleading, regardless of whether it is implemented by
<o:p></o:p></p>
<p class="MsoPlainText">> visual inspection or using objective criteria). I personally don't
<o:p></o:p></p>
<p class="MsoPlainText">> believe that there is any such thing as artifact-free EEG recordings.<o:p></o:p></p>
<p class="MsoPlainText">> The artifact contributions are just more or less dominant...<o:p></o:p></p>
<p class="MsoPlainText">><o:p> </o:p></p>
<p class="MsoPlainText">> - The ICA model comes with a couple assumptions (like any statistical
<o:p></o:p></p>
<p class="MsoPlainText">> approach), and if the data do not adhere to the model assumptions then
<o:p></o:p></p>
<p class="MsoPlainText">> the resulting decomposition may be of very poor quality. In my view,
<o:p></o:p></p>
<p class="MsoPlainText">> it does not make much sense to praise or condemn a procedure without
<o:p></o:p></p>
<p class="MsoPlainText">> keeping in mind this fact. Keeping an eye on model assumptions could
<o:p></o:p></p>
<p class="MsoPlainText">> (and I think, should) guide the preprocessing steps. ICA decomposition
<o:p></o:p></p>
<p class="MsoPlainText">> quality depends clearly on how the data are preprocessed (e.g.;
<o:p></o:p></p>
<p class="MsoPlainText">> <a href="http://ieeexplore.ieee.org/document/7319296/?reload=true">
<span style="color:windowtext;text-decoration:none">http://ieeexplore.ieee.org/document/7319296/?reload=true</span></a>). You may
<o:p></o:p></p>
<p class="MsoPlainText">> also want to consider the SPR guidelines:<o:p></o:p></p>
<p class="MsoPlainText">> <a href="http://onlinelibrary.wiley.com/doi/10.1111/psyp.12147/abstract;jsessionid=ADC0B75ACEC7963AA2BCA8A10F317A3F.f02t04">
<span style="color:windowtext;text-decoration:none">http://onlinelibrary.wiley.com/doi/10.1111/psyp.12147/abstract;jsessionid=ADC0B75ACEC7963AA2BCA8A10F317A3F.f02t04</span></a>.<o:p></o:p></p>
<p class="MsoPlainText">><o:p> </o:p></p>
<p class="MsoPlainText">> - I remotely remember Robert or someones else stating in a previous
<o:p></o:p></p>
<p class="MsoPlainText">> mail that DTI, ADHD and other types of evidence would validate the
<o:p></o:p></p>
<p class="MsoPlainText">> functional connectivity patterns observed with inverse solutions of
<o:p></o:p></p>
<p class="MsoPlainText">> 19-channel raw EEG (I got a similar impression when exploring the
<o:p></o:p></p>
<p class="MsoPlainText">> appliedneuroscience.com website). Well, I find those associations to
<o:p></o:p></p>
<p class="MsoPlainText">> be not convincing. With or without ICA, the inverse problem should not
<o:p></o:p></p>
<p class="MsoPlainText">> be neglected, some care is required for the source level
<o:p></o:p></p>
<p class="MsoPlainText">> interpretation. Playing a little with parameter settings for source
<o:p></o:p></p>
<p class="MsoPlainText">> analysis is a good reminder that the resulting outcome could be
<o:p></o:p></p>
<p class="MsoPlainText">> dramatically different for different settings (same holds for
<o:p></o:p></p>
<p class="MsoPlainText">> channel-level analysis as well - outcome depends on your processing!).<o:p></o:p></p>
<p class="MsoPlainText">> In my opinion, source-level analysis can only confirm predictions
<o:p></o:p></p>
<p class="MsoPlainText">> (i.e., the result matching your priori expectations); I personally
<o:p></o:p></p>
<p class="MsoPlainText">> don't trust EEG source activations at unexpected locations, they
<o:p></o:p></p>
<p class="MsoPlainText">> appear, more often than not, spurious to me. In contrast, the use of
<o:p></o:p></p>
<p class="MsoPlainText">> ICA for eye blink correction seems much, much better validated to me!<o:p></o:p></p>
<p class="MsoPlainText">><o:p> </o:p></p>
<p class="MsoPlainText">> - I don't really get the concern of ICA messing up the phase of the
<o:p></o:p></p>
<p class="MsoPlainText">> continuous EEG signal. Of course, a spatial (or temporal) filter will
<o:p></o:p></p>
<p class="MsoPlainText">> modify the signal, that's it's purpose, and if the filter attenuates
<o:p></o:p></p>
<p class="MsoPlainText">> some portions of the signal, the residual signal may have a different
<o:p></o:p></p>
<p class="MsoPlainText">> phase, amplitude and/or topography. Because for real EEG recordings
<o:p></o:p></p>
<p class="MsoPlainText">> nobody knows the ground truth, the question should be whether a filter
<o:p></o:p></p>
<p class="MsoPlainText">> makes the data better or worse (for an excellent discussion on
<o:p></o:p></p>
<p class="MsoPlainText">> temporal filters making the data better/worse, see:<o:p></o:p></p>
<p class="MsoPlainText">> <a href="https://www.ncbi.nlm.nih.gov/pubmed/25128257">
<span style="color:windowtext;text-decoration:none">https://www.ncbi.nlm.nih.gov/pubmed/25128257</span></a>). Assuming that the raw
<o:p></o:p></p>
<p class="MsoPlainText">> data are a (hopefully linear) mixture of an unknown number of (brain
<o:p></o:p></p>
<p class="MsoPlainText">> and non-brain) sources, and, more likely than not, that the number of
<o:p></o:p></p>
<p class="MsoPlainText">> sources contributing are larger than the number of channels (at least
<o:p></o:p></p>
<p class="MsoPlainText">> for 19-ch this seems pretty obvious to me), NOT processing the data at
<o:p></o:p></p>
<p class="MsoPlainText">> all does not give a good reference for any comparison. In other words,
<o:p></o:p></p>
<p class="MsoPlainText">> the raw data do not qualify very well as the gold standard, because
<o:p></o:p></p>
<p class="MsoPlainText">> they may be messy (i.e., mixed!). More informative would be a
<o:p></o:p></p>
<p class="MsoPlainText">> comparison (of, say, two different filter approaches) with regard to a
<o:p></o:p></p>
<p class="MsoPlainText">> particular effect of interest (say, theta and working memory
<o:p></o:p></p>
<p class="MsoPlainText">> relationship, or any feature one has sufficient evidence to justify a
<o:p></o:p></p>
<p class="MsoPlainText">> clear prediction).<o:p></o:p></p>
<p class="MsoPlainText">><o:p> </o:p></p>
<p class="MsoPlainText">> Best,<o:p></o:p></p>
<p class="MsoPlainText">> Stefan<o:p></o:p></p>
<p class="MsoPlainText">><o:p> </o:p></p>
<p class="MsoPlainText">><o:p> </o:p></p>
<p class="MsoPlainText">> Am 17.06.17 um 19:42 schrieb Arnaud Delorme:<o:p></o:p></p>
<p class="MsoPlainText">>> Dear Robert,<o:p></o:p></p>
<p class="MsoPlainText">>><o:p> </o:p></p>
<p class="MsoPlainText">>> Thank you for your email and for your analysis. A few comments below:<o:p></o:p></p>
<p class="MsoPlainText">>><o:p> </o:p></p>
<p class="MsoPlainText">>> 1- I find it hard to believe all of the plots in your compressed
<o:p></o:p></p>
<p class="MsoPlainText">>> document compared artifact free regions of data. For example, you
<o:p></o:p></p>
<p class="MsoPlainText">>> observed 71% change over frontal channels in the delta band (see
<o:p></o:p></p>
<p class="MsoPlainText">>> attached screen capture in the PDF document from PCT DIFF PRE VS
<o:p></o:p></p>
<p class="MsoPlainText">>> POST-DELORME/PRE VS POST-DELORME_4.bmp). This is consistent with
<o:p></o:p></p>
<p class="MsoPlainText">>> selecting portions of data which contains eye-related artifacts as I
<o:p></o:p></p>
<p class="MsoPlainText">>> show in the rest of the attached PDF document.<o:p></o:p></p>
<p class="MsoPlainText">>><o:p> </o:p></p>
<p class="MsoPlainText">>> 2- The data you shared has a high density of blinks, there is only a
<o:p></o:p></p>
<p class="MsoPlainText">>> handful segments of data where the subject stops blinking for more
<o:p></o:p></p>
<p class="MsoPlainText">>> than 5 seconds. Even in these clean segments, in the frontal
<o:p></o:p></p>
<p class="MsoPlainText">>> channels, we can see some slow some activity that likely reflect eye
<o:p></o:p></p>
<p class="MsoPlainText">>> movement (or contamination from previous eye movements by the filter
<o:p></o:p></p>
<p class="MsoPlainText">>> you might have used). ICA is able to remove these as well. I am
<o:p></o:p></p>
<p class="MsoPlainText">>> attaching two examples in the PDF document. In between the two blinks
<o:p></o:p></p>
<p class="MsoPlainText">>> below (see document), you can see that the “clean” segment is not
<o:p></o:p></p>
<p class="MsoPlainText">>> that clean. It is still contaminated by eye artifacts at least on the
<o:p></o:p></p>
<p class="MsoPlainText">>> first 3 channels (in order from top to bottom FP1,FP2, F7, F3, Fz,
<o:p></o:p></p>
<p class="MsoPlainText">>> F4, F8). I do not think that anybody would argue that removing this
<o:p></o:p></p>
<p class="MsoPlainText">>> activity amounts to removing brain activity. So it is important to be
<o:p></o:p></p>
<p class="MsoPlainText">>> careful when selecting “clean” data segment.<o:p></o:p></p>
<p class="MsoPlainText">>><o:p> </o:p></p>
<p class="MsoPlainText">>> 3. I have removed all portions of data where eye activity was visible
<o:p></o:p></p>
<p class="MsoPlainText">>> preserving only 40 seconds out of 410 in the original file and I am
<o:p></o:p></p>
<p class="MsoPlainText">>> assessing below the data distortion. Spectrum was calculated on
<o:p></o:p></p>
<p class="MsoPlainText">>> non-overlapping 1 minute data segments (I am attaching the script). I
<o:p></o:p></p>
<p class="MsoPlainText">>> am attaching the scalp topography pre vs post in different frequency
<o:p></o:p></p>
<p class="MsoPlainText">>> bands.<o:p></o:p></p>
<p class="MsoPlainText">>><o:p> </o:p></p>
<p class="MsoPlainText">>> We can also see that the power at other channels is not affected.
<o:p></o:p></p>
<p class="MsoPlainText">>> Beyond 5 Hz, the power spectrum does not seem to be affected even in
<o:p></o:p></p>
<p class="MsoPlainText">>> frontal channels with very similar scalp topographies. In the
<o:p></o:p></p>
<p class="MsoPlainText">>> attached document, I show the scalp topography pre vs post in
<o:p></o:p></p>
<p class="MsoPlainText">>> different frequency bands. Note the difference in frontal channels at
<o:p></o:p></p>
<p class="MsoPlainText">>> 1Hz and 5Hz. This is due to what is mentioned in 2, that even clean
<o:p></o:p></p>
<p class="MsoPlainText">>> segments are not totally free of artifacts and I believe that these
<o:p></o:p></p>
<p class="MsoPlainText">>> represent true eye artifacts removed by ICA. We can also see that the
<o:p></o:p></p>
<p class="MsoPlainText">>> power at other channels is not affected. Beyond 5 Hz, the power
<o:p></o:p></p>
<p class="MsoPlainText">>> spectrum does not seem to be affected even in frontal channels with
<o:p></o:p></p>
<p class="MsoPlainText">>> very similar scalp topographies. There is no 71% difference in
<o:p></o:p></p>
<p class="MsoPlainText">>> absolute power as in your plot (in my case, it seems to be about 25%
<o:p></o:p></p>
<p class="MsoPlainText">>> at 1Hz over frontal channel and again, this correspond to removed
<o:p></o:p></p>
<p class="MsoPlainText">>> residual eye blinks).<o:p></o:p></p>
<p class="MsoPlainText">>><o:p> </o:p></p>
<p class="MsoPlainText">>> 4. P-value need to be corrected for multiple comparisons. You have 19
<o:p></o:p></p>
<p class="MsoPlainText">>> channels and are looking at 20 frequency bands. That’s 380
<o:p></o:p></p>
<p class="MsoPlainText">>> t-test/p-value. By chance at the 5% threshold, 19 of these are going
<o:p></o:p></p>
<p class="MsoPlainText">>> to be significant. This is an important limitation of the QEEG field
<o:p></o:p></p>
<p class="MsoPlainText">>> by the way and the reason why QEEG papers rarely make it to reputable
<o:p></o:p></p>
<p class="MsoPlainText">>> journals. And I have heard people say: well if I correct with
<o:p></o:p></p>
<p class="MsoPlainText">>> classical Bonferoni, nothing shows up significant anymore, but there
<o:p></o:p></p>
<p class="MsoPlainText">>> are less aggressive methods to correct for multiple comparison such
<o:p></o:p></p>
<p class="MsoPlainText">>> as False Discovery Rate, the cluster and max method. The QEEG field
<o:p></o:p></p>
<p class="MsoPlainText">>> should use these methods and they could easily be implemented in your
<o:p></o:p></p>
<p class="MsoPlainText">>> software.<o:p></o:p></p>
<p class="MsoPlainText">>><o:p> </o:p></p>
<p class="MsoPlainText">>> I am also attaching the script to reproduce my results from the raw
<o:p></o:p></p>
<p class="MsoPlainText">>> EDF data file you shared (it in the PDF).<o:p></o:p></p>
<p class="MsoPlainText">>> ICA is not a magical technique but it works pretty well for removing
<o:p></o:p></p>
<p class="MsoPlainText">>> eye artifacts with the least of distortion of the signal. See for
<o:p></o:p></p>
<p class="MsoPlainText">>> example this paper on real and simulated data (not from us<o:p></o:p></p>
<p class="MsoPlainText">>> <a href="http://journals.plos.org/plosone/article?id=10.1371/journal.pone.0003">
<span style="color:windowtext;text-decoration:none">http://journals.plos.org/plosone/article?id=10.1371/journal.pone.0003</span></a><o:p></o:p></p>
<p class="MsoPlainText">>> 004) concluding "ICA yielded almost perfect correction in all
<o:p></o:p></p>
<p class="MsoPlainText">>> conditions.”<o:p></o:p></p>
<p class="MsoPlainText">>> There are many other such papers). And I would be the first to agree
<o:p></o:p></p>
<p class="MsoPlainText">>> that there is a lot of subtleties of using ICA for artifact rejection
<o:p></o:p></p>
<p class="MsoPlainText">>> and lot of potential limitations (instabilities in some conditions,
<o:p></o:p></p>
<p class="MsoPlainText">>> issues with numerical precisions (double vs single), lack of exact
<o:p></o:p></p>
<p class="MsoPlainText">>> reproducibility in some other cases because of the random initial
<o:p></o:p></p>
<p class="MsoPlainText">>> conditions, multitude of algorithms). We have tried to address these
<o:p></o:p></p>
<p class="MsoPlainText">>> issues over the years. Even in your dataset, it seems that Infomax
<o:p></o:p></p>
<p class="MsoPlainText">>> ICA converge to at least 2 different solutions (I have run it 10
<o:p></o:p></p>
<p class="MsoPlainText">>> times on your data at least) yielding slightly different eye
<o:p></o:p></p>
<p class="MsoPlainText">>> artifacts (both seem valid but differ at very low frequencies < 3Hz
<o:p></o:p></p>
<p class="MsoPlainText">>> and very high frequencies above 60 Hz) and I will have to look at
<o:p></o:p></p>
<p class="MsoPlainText">>> that in more details (I suspect this has to do with the stopping rule
<o:p></o:p></p>
<p class="MsoPlainText">>> and the low number of channels in your dataset - I will run some
<o:p></o:p></p>
<p class="MsoPlainText">>> tests and might have to change the default stop threshold for
<o:p></o:p></p>
<p class="MsoPlainText">>> convergence). It remains that ICA is the best we have right now for
<o:p></o:p></p>
<p class="MsoPlainText">>> artifact removal, and I am convinced that a poll of scientist in our
<o:p></o:p></p>
<p class="MsoPlainText">>> field would show that more than 80% of EEG scientist agree with me.<o:p></o:p></p>
<p class="MsoPlainText">>><o:p> </o:p></p>
<p class="MsoPlainText">>> I think it is better to rely on published data than argue about a
<o:p></o:p></p>
<p class="MsoPlainText">>> particular dataset, although it is useful as well. At this point we
<o:p></o:p></p>
<p class="MsoPlainText">>> may agree to disagree and EEG users can compare and contrast
<o:p></o:p></p>
<p class="MsoPlainText">>> different approaches. They can download your data and use the script
<o:p></o:p></p>
<p class="MsoPlainText">>> I provide to reproduce my results.<o:p></o:p></p>
<p class="MsoPlainText">>><o:p> </o:p></p>
<p class="MsoPlainText">>> Best wishes,<o:p></o:p></p>
<p class="MsoPlainText">>><o:p> </o:p></p>
<p class="MsoPlainText">>> Arno<o:p></o:p></p>
<p class="MsoPlainText">>><o:p> </o:p></p>
<p class="MsoPlainText">>> PDF document<o:p></o:p></p>
<p class="MsoPlainText">>> <a href="http://sccn.ucsd.edu/~arno/download/icaphaserebutal.pdf">
<span style="color:windowtext;text-decoration:none">http://sccn.ucsd.edu/~arno/download/icaphaserebutal.pdf</span></a><o:p></o:p></p>
<p class="MsoPlainText">>> <<a href="http://sccn.ucsd.edu/%7Earno/download/icaphaserebutal.pdf"><span style="color:windowtext;text-decoration:none">http://sccn.ucsd.edu/%7Earno/download/icaphaserebutal.pdf</span></a>><o:p></o:p></p>
<p class="MsoPlainText">>><o:p> </o:p></p>
<p class="MsoPlainText">>><o:p> </o:p></p>
<p class="MsoPlainText">>>> On Jun 15, 2017, at 12:32 PM, Robert Thatcher <rwthatcher2@yahoo.com
<o:p></o:p></p>
<p class="MsoPlainText">>>> <<a href="mailto:rwthatcher2@yahoo.com"><span style="color:windowtext;text-decoration:none">mailto:rwthatcher2@yahoo.com</span></a>>> wrote:<o:p></o:p></p>
<p class="MsoPlainText">>>><o:p> </o:p></p>
<p class="MsoPlainText">>>> Dear Arnaud,<o:p></o:p></p>
<p class="MsoPlainText">>>><o:p> </o:p></p>
<p class="MsoPlainText">>>> I did statistical comparisons between 1 min & 40 seconds of
<o:p></o:p></p>
<p class="MsoPlainText">>>> artifact free EEG in the original EEG recording in with no eye
<o:p></o:p></p>
<p class="MsoPlainText">>>> movement artifact and the Win EEG ICA reconstruction and the ICA
<o:p></o:p></p>
<p class="MsoPlainText">>>> reconstruction that you did. I compared two different artifact
<o:p></o:p></p>
<p class="MsoPlainText">>>> rejection methods used on the original EEG: 1- manual selections of
<o:p></o:p></p>
<p class="MsoPlainText">>>> artifact free data and, 2- the automatic template method of artifact
<o:p></o:p></p>
<p class="MsoPlainText">>>> rejection where I hand selected a 10 second sample of artifact free
<o:p></o:p></p>
<p class="MsoPlainText">>>> EEG and then used an algorithm that matched the peak-to-peak
<o:p></o:p></p>
<p class="MsoPlainText">>>> amplitudes of the 10 second template to the remainder of the<o:p></o:p></p>
<p class="MsoPlainText">>>> record. There were no statistically significant differences
<o:p></o:p></p>
<p class="MsoPlainText">>>> between these two artifact rejection methods.<o:p></o:p></p>
<p class="MsoPlainText">>>><o:p> </o:p></p>
<p class="MsoPlainText">>>> Based on the time points of the artifact free data in the original
<o:p></o:p></p>
<p class="MsoPlainText">>>> EEG I selected the exact same time points in the Win EEG ICA
<o:p></o:p></p>
<p class="MsoPlainText">>>> reconstruction and in your ICA reconstruction. Therefore all three
<o:p></o:p></p>
<p class="MsoPlainText">>>> data files contained 1 minute & 40 seconds of the same time points.<o:p></o:p></p>
<p class="MsoPlainText">>>> I then computed percent differences as well as paired t-tests
<o:p></o:p></p>
<p class="MsoPlainText">>>> between the original EEG and the two ICA reconstructions. Here is
<o:p></o:p></p>
<p class="MsoPlainText">>>> a url to download the results:<o:p></o:p></p>
<p class="MsoPlainText">>>><o:p> </o:p></p>
<p class="MsoPlainText">>>> <a href="http://www.appliedneuroscience.com/STATISTICS">
<span style="color:windowtext;text-decoration:none">http://www.appliedneuroscience.com/STATISTICS</span></a> OF ARTIFACT FREE EEG
<o:p></o:p></p>
<p class="MsoPlainText">>>> VS POST ICA EEG.zip<o:p></o:p></p>
<p class="MsoPlainText">>>><o:p> </o:p></p>
<p class="MsoPlainText">>>> As you can see there were very large statistically significant
<o:p></o:p></p>
<p class="MsoPlainText">>>> differences between the artifact free EEG in the original recording<o:p></o:p></p>
<p class="MsoPlainText">>>> and the ICA reconstructions. Your reconstruction was less
<o:p></o:p></p>
<p class="MsoPlainText">>>> distorted than the Win EEG reconstruction but both were
<o:p></o:p></p>
<p class="MsoPlainText">>>> significantly different than the original artifact free EEG.<o:p></o:p></p>
<p class="MsoPlainText">>>><o:p> </o:p></p>
<p class="MsoPlainText">>>> I would be happy to send you the .edfs of the selected time points
<o:p></o:p></p>
<p class="MsoPlainText">>>> so that you can verify that the time points were identical and the
<o:p></o:p></p>
<p class="MsoPlainText">>>> original EEG did not have any eye movement artifact.<o:p></o:p></p>
<p class="MsoPlainText">>>><o:p> </o:p></p>
<p class="MsoPlainText">>>> These large magnitude of the differences between the original and
<o:p></o:p></p>
<p class="MsoPlainText">>>> unaltered data vs. the ICA altered data are similar to those that
<o:p></o:p></p>
<p class="MsoPlainText">>>> many WinEEG users find when they use the WinEEG ICA reconstruction<o:p></o:p></p>
<p class="MsoPlainText">>>> method. Therefore these large differences are not surprising and
<o:p></o:p></p>
<p class="MsoPlainText">>>> are commonly found especially when using the WinEEG ICA. For
<o:p></o:p></p>
<p class="MsoPlainText">>>> example, Georges or Robert Lawson and others.<o:p></o:p></p>
<p class="MsoPlainText">>>><o:p> </o:p></p>
<p class="MsoPlainText">>>> I also included screen captures of some of the waveforms showing
<o:p></o:p></p>
<p class="MsoPlainText">>>> visually detectable differences between the original and the ICA<o:p></o:p></p>
<p class="MsoPlainText">>>> reconstruction using the WinEEG ICA. The ICA that you used
<o:p></o:p></p>
<p class="MsoPlainText">>>> produced less visually obvious waveform changes but nonetheless
<o:p></o:p></p>
<p class="MsoPlainText">>>> there are some that are visually detectable. However, the best way
<o:p></o:p></p>
<p class="MsoPlainText">>>> to understand the alterations of the artifact free sections is by
<o:p></o:p></p>
<p class="MsoPlainText">>>> JTFA and/or FFT and statistics.<o:p></o:p></p>
<p class="MsoPlainText">>>><o:p> </o:p></p>
<p class="MsoPlainText">>>> Thank you for your patients in and dedication to exploring this
<o:p></o:p></p>
<p class="MsoPlainText">>>> important topic. It is an important topic because of the obvious
<o:p></o:p></p>
<p class="MsoPlainText">>>> discrepancies that will exist in the scientific literature between
<o:p></o:p></p>
<p class="MsoPlainText">>>> simple deletion of artifact vs ICA reconstruction going forward.<o:p></o:p></p>
<p class="MsoPlainText">>>> Also because the entire EEG record is modified the ability to
<o:p></o:p></p>
<p class="MsoPlainText">>>> replicate findings is reduced when using ICA reconstruction. Also,
<o:p></o:p></p>
<p class="MsoPlainText">>>> because there is some degree of decoupling between the underlying
<o:p></o:p></p>
<p class="MsoPlainText">>>> physiological origins of the EEG and a patient's brain then clinical
<o:p></o:p></p>
<p class="MsoPlainText">>>> correlation or effect size will be lower.<o:p></o:p></p>
<p class="MsoPlainText">>>><o:p> </o:p></p>
<p class="MsoPlainText">>>> Best regards,<o:p></o:p></p>
<p class="MsoPlainText">>>><o:p> </o:p></p>
<p class="MsoPlainText">>>> Robert<o:p></o:p></p>
<p class="MsoPlainText">>>><o:p> </o:p></p>
<p class="MsoPlainText">>>><o:p> </o:p></p>
<p class="MsoPlainText">>>><o:p> </o:p></p>
<p class="MsoPlainText">>>> On Thursday, June 15, 2017, 10:11:45 AM EDT, Robert Thatcher
<o:p></o:p></p>
<p class="MsoPlainText">>>> <<a href="mailto:rwthatcher2@yahoo.com%20%3cmailto:rwthatcher2@yahoo.com"><span style="color:windowtext;text-decoration:none">rwthatcher2@yahoo.com <mailto:rwthatcher2@yahoo.com</span></a>>> wrote:<o:p></o:p></p>
<p class="MsoPlainText">>>><o:p> </o:p></p>
<p class="MsoPlainText">>>><o:p> </o:p></p>
<p class="MsoPlainText">>>> Hi Arno,<o:p></o:p></p>
<p class="MsoPlainText">>>><o:p> </o:p></p>
<p class="MsoPlainText">>>> Thank you for your thoughtful post.<o:p></o:p></p>
<p class="MsoPlainText">>>><o:p> </o:p></p>
<p class="MsoPlainText">>>> “As far as phase distortion after removing ICA components (in my
<o:p></o:p></p>
<p class="MsoPlainText">>>> decomposition), I am not sure what you are referring to. Is it the
<o:p></o:p></p>
<p class="MsoPlainText">>>> minute shift when the red and black curve do not exactly superpose.”<o:p></o:p></p>
<p class="MsoPlainText">>>><o:p> </o:p></p>
<p class="MsoPlainText">>>> I am referring to the differences in phase between pairs of EEG
<o:p></o:p></p>
<p class="MsoPlainText">>>> channels. One can visually see differences in particular segments
<o:p></o:p></p>
<p class="MsoPlainText">>>> but it is best to use the Hilbert transform (cross-spectra) to
<o:p></o:p></p>
<p class="MsoPlainText">>>> compute instantaneous phase differences at any point in the record
<o:p></o:p></p>
<p class="MsoPlainText">>>> that one may want to average the absolute phase differences over
<o:p></o:p></p>
<p class="MsoPlainText">>>> some period of the record where there is no artifact and then<o:p></o:p></p>
<p class="MsoPlainText">>>> conduct t-tests to evaluate the large effect sizes. One can also
<o:p></o:p></p>
<p class="MsoPlainText">>>> compare the FFT spectra which is also an average, albeit more
<o:p></o:p></p>
<p class="MsoPlainText">>>> noisy. The alteration of phase differences are present no matter
<o:p></o:p></p>
<p class="MsoPlainText">>>> what measure one uses. The least reliable is a visual analysis
<o:p></o:p></p>
<p class="MsoPlainText">>>> although there are plenty of visual examples if one carefully
<o:p></o:p></p>
<p class="MsoPlainText">>>> reviews the traces.<o:p></o:p></p>
<p class="MsoPlainText">>>><o:p> </o:p></p>
<p class="MsoPlainText">>>> "I would argue that the data after removing ICA artifacts reflect
<o:p></o:p></p>
<p class="MsoPlainText">>>> more brain activity than before, and that the minute shift is due to
<o:p></o:p></p>
<p class="MsoPlainText">>>> removal of small eye movement activity. I agree that this would have
<o:p></o:p></p>
<p class="MsoPlainText">>>> to be demonstrated, and that you cannot take my word for it.”<o:p></o:p></p>
<p class="MsoPlainText">>>><o:p> </o:p></p>
<p class="MsoPlainText">>>> Myself and many others do not disagree that elimination of artifact
<o:p></o:p></p>
<p class="MsoPlainText">>>> is important what we disagree with is the ICA reconstruction method
<o:p></o:p></p>
<p class="MsoPlainText">>>> that adulterates the artifact free segments of the record. Why not
<o:p></o:p></p>
<p class="MsoPlainText">>>> simply delete the eye movement manually or like Neuroguide does with
<o:p></o:p></p>
<p class="MsoPlainText">>>> a signal detection algorithm that measures the voltage gradients
<o:p></o:p></p>
<p class="MsoPlainText">>>> produced by a blink or eye movement, etc? In this way all of the
<o:p></o:p></p>
<p class="MsoPlainText">>>> original digital data samples are unaltered.<o:p></o:p></p>
<p class="MsoPlainText">>>><o:p> </o:p></p>
<p class="MsoPlainText">>>><o:p> </o:p></p>
<p class="MsoPlainText">>>> “The EEG signal is extremely noisy.”<o:p></o:p></p>
<p class="MsoPlainText">>>><o:p> </o:p></p>
<p class="MsoPlainText">>>> The vast number of EEG experts would disagree with you that the “EEG
<o:p></o:p></p>
<p class="MsoPlainText">>>> is extremely noisy”. If this were true it would be obvious to
<o:p></o:p></p>
<p class="MsoPlainText">>>> every one with a total inability to replicate any EEG study and
<o:p></o:p></p>
<p class="MsoPlainText">>>> there would not be over 100,000 peer reviewed studies published in
<o:p></o:p></p>
<p class="MsoPlainText">>>> the National Library of Medicine. Simply visually examine the EEG
<o:p></o:p></p>
<p class="MsoPlainText">>>> traces showing well behaved and well organized alpha rhythms or
<o:p></o:p></p>
<p class="MsoPlainText">>>> theta rhythms or beta rhythms which reflect large synchronous LFPs.<o:p></o:p></p>
<p class="MsoPlainText">>>><o:p> </o:p></p>
<p class="MsoPlainText">>>><o:p> </o:p></p>
<p class="MsoPlainText">>>> “The phase of the signal at one electrode site and one given time is
<o:p></o:p></p>
<p class="MsoPlainText">>>> not representative of the underlying brain signal.”<o:p></o:p></p>
<p class="MsoPlainText">>>><o:p> </o:p></p>
<p class="MsoPlainText">>>> This also cannot be true because the phase difference between
<o:p></o:p></p>
<p class="MsoPlainText">>>> electrodes and/or sources are produced by the physiological
<o:p></o:p></p>
<p class="MsoPlainText">>>> foundations of the brain and networks and are due to differences in
<o:p></o:p></p>
<p class="MsoPlainText">>>> synaptic rise times, synaptic integration times, differences in
<o:p></o:p></p>
<p class="MsoPlainText">>>> conduction velocity, etc. This is the underlying brain signal and
<o:p></o:p></p>
<p class="MsoPlainText">>>> it is highly reproducible and clinically useful. If your belief were
<o:p></o:p></p>
<p class="MsoPlainText">>>> valid then there would be no clinical correlations to the EEG such
<o:p></o:p></p>
<p class="MsoPlainText">>>> as schizophrenia or ADHD or depression or epilepsy or drug effects,
<o:p></o:p></p>
<p class="MsoPlainText">>>> etc.<o:p></o:p></p>
<p class="MsoPlainText">>>><o:p> </o:p></p>
<p class="MsoPlainText">>>><o:p> </o:p></p>
<p class="MsoPlainText">>>> “, if you have a picture of a star, would you rather remove a visual
<o:p></o:p></p>
<p class="MsoPlainText">>>> artifact that is 10-fold the size of your original signal or
<o:p></o:p></p>
<p class="MsoPlainText">>>> continue to look at your original signal (not being able to see much
<o:p></o:p></p>
<p class="MsoPlainText">>>> because of the large artifact masking most of it).”<o:p></o:p></p>
<p class="MsoPlainText">>>><o:p> </o:p></p>
<p class="MsoPlainText">>>> This is an interesting take on my analogy and I agree that the
<o:p></o:p></p>
<p class="MsoPlainText">>>> 10-fold size artifact needs to be avoided or eliminated but not by
<o:p></o:p></p>
<p class="MsoPlainText">>>> using ICA reconstruction that effects the artifact free parts of the
<o:p></o:p></p>
<p class="MsoPlainText">>>> spectrum and thereby distorts the measurement not only of the one
<o:p></o:p></p>
<p class="MsoPlainText">>>> star that you are looking at but also all other stars and planets in
<o:p></o:p></p>
<p class="MsoPlainText">>>> the universe.<o:p></o:p></p>
<p class="MsoPlainText">>>><o:p> </o:p></p>
<p class="MsoPlainText">>>><o:p> </o:p></p>
<p class="MsoPlainText">>>> “Even if ICA was introducing minute distortion in phase”<o:p></o:p></p>
<p class="MsoPlainText">>>><o:p> </o:p></p>
<p class="MsoPlainText">>>> I wish that the distortion in phase difference was “minute” but the
<o:p></o:p></p>
<p class="MsoPlainText">>>> fact is that it is large and easily demonstrated as it has been by
<o:p></o:p></p>
<p class="MsoPlainText">>>> numerous scientists/clinicians over the last few years. For
<o:p></o:p></p>
<p class="MsoPlainText">>>> example, t-tests between the artifact free segments in the original
<o:p></o:p></p>
<p class="MsoPlainText">>>> EEG vs. the ICA reconstructed new time series are mostly significant
<o:p></o:p></p>
<p class="MsoPlainText">>>> at P < 0.00001. I will do some additional statistical comparisons so
<o:p></o:p></p>
<p class="MsoPlainText">>>> that you can better understand the large effect sizes of ICA phase
<o:p></o:p></p>
<p class="MsoPlainText">>>> difference distortion.<o:p></o:p></p>
<p class="MsoPlainText">>>><o:p> </o:p></p>
<p class="MsoPlainText">>>> This is an important dialog and I appreciate your dedication and
<o:p></o:p></p>
<p class="MsoPlainText">>>> willingness to consider these issues.<o:p></o:p></p>
<p class="MsoPlainText">>>><o:p> </o:p></p>
<p class="MsoPlainText">>>><o:p> </o:p></p>
<p class="MsoPlainText">>>> Best wishes,<o:p></o:p></p>
<p class="MsoPlainText">>>><o:p> </o:p></p>
<p class="MsoPlainText">>>><o:p> </o:p></p>
<p class="MsoPlainText">>>> Robert<o:p></o:p></p>
<p class="MsoPlainText">>>><o:p> </o:p></p>
<p class="MsoPlainText">>>><o:p> </o:p></p>
<p class="MsoPlainText">>>><o:p> </o:p></p>
<p class="MsoPlainText">>>> On Wednesday, June 14, 2017, 11:49:54 PM EDT, Arnaud Delorme
<o:p></o:p></p>
<p class="MsoPlainText">>>> <<a href="mailto:arno@ucsd.edu%20%3cmailto:arno@ucsd.edu"><span style="color:windowtext;text-decoration:none">arno@ucsd.edu <mailto:arno@ucsd.edu</span></a>>> wrote:<o:p></o:p></p>
<p class="MsoPlainText">>>><o:p> </o:p></p>
<p class="MsoPlainText">>>><o:p> </o:p></p>
<p class="MsoPlainText">>>> Dear Robert,<o:p></o:p></p>
<p class="MsoPlainText">>>><o:p> </o:p></p>
<p class="MsoPlainText">>>> There is no need to remove more components - except maybe for
<o:p></o:p></p>
<p class="MsoPlainText">>>> temporal muscles components (I would have to look again at your data
<o:p></o:p></p>
<p class="MsoPlainText">>>> to see if I can identify any). The procedure is to identify a
<o:p></o:p></p>
<p class="MsoPlainText">>>> handful of artifact components, remove them and then your data is
<o:p></o:p></p>
<p class="MsoPlainText">>>> cleaned of these artifacts. I personally rarely identify more than 4
<o:p></o:p></p>
<p class="MsoPlainText">>>> artifact component in a given subject (some other researchers have a
<o:p></o:p></p>
<p class="MsoPlainText">>>> more aggressive approach and remove more). I like to remove
<o:p></o:p></p>
<p class="MsoPlainText">>>> components I am sure of.<o:p></o:p></p>
<p class="MsoPlainText">>>><o:p> </o:p></p>
<p class="MsoPlainText">>>> ICA is a linear decomposition that isolate sources which are
<o:p></o:p></p>
<p class="MsoPlainText">>>> maximally independent. Blinks are mostly independent of brain
<o:p></o:p></p>
<p class="MsoPlainText">>>> activity (on first approximation) so ICA is able to isolate them.<o:p></o:p></p>
<p class="MsoPlainText">>>><o:p> </o:p></p>
<p class="MsoPlainText">>>> "You already showed that two ICA component removals results in more
<o:p></o:p></p>
<p class="MsoPlainText">>>> phase distortion than the removal of one ICA component.” In all of
<o:p></o:p></p>
<p class="MsoPlainText">>>> our exchanges I have always removed 2 components. I have never
<o:p></o:p></p>
<p class="MsoPlainText">>>> removed one ICA component.<o:p></o:p></p>
<p class="MsoPlainText">>>><o:p> </o:p></p>
<p class="MsoPlainText">>>> As far as phase distortion after removing ICA components (in my
<o:p></o:p></p>
<p class="MsoPlainText">>>> decomposition), I am not sure what you are referring to. Is it the
<o:p></o:p></p>
<p class="MsoPlainText">>>> minute shift when the red and black curve do not exactly superpose.<o:p></o:p></p>
<p class="MsoPlainText">>>> I have two comments on that.<o:p></o:p></p>
<p class="MsoPlainText">>>><o:p> </o:p></p>
<p class="MsoPlainText">>>> I would argue that the data after removing ICA artifacts reflect
<o:p></o:p></p>
<p class="MsoPlainText">>>> more brain activity than before, and that the minute shift is due to
<o:p></o:p></p>
<p class="MsoPlainText">>>> removal of small eye movement activity. I agree that this would have
<o:p></o:p></p>
<p class="MsoPlainText">>>> to be demonstrated, and that you cannot take my word for it.<o:p></o:p></p>
<p class="MsoPlainText">>>><o:p> </o:p></p>
<p class="MsoPlainText">>>> "Even a small amount of adulteration or distortion of EEG phase
<o:p></o:p></p>
<p class="MsoPlainText">>>> differences is not good and must be avoided at all costs. This is
<o:p></o:p></p>
<p class="MsoPlainText">>>> analogous to the use of telescopes that measure phase differences in
<o:p></o:p></p>
<p class="MsoPlainText">>>> the spectrum from stars moving in the universe."<o:p></o:p></p>
<p class="MsoPlainText">>>><o:p> </o:p></p>
<p class="MsoPlainText">>>> The data we are looking at on the scalp is a summation of millions
<o:p></o:p></p>
<p class="MsoPlainText">>>> of neuron activity and the phase we are observing a cumulative
<o:p></o:p></p>
<p class="MsoPlainText">>>> average of this signal (pondered by the geometry of the brain,
<o:p></o:p></p>
<p class="MsoPlainText">>>> difference in conductivity of different tissues etc...). The EEG
<o:p></o:p></p>
<p class="MsoPlainText">>>> signal is extremely noisy. The phase of the signal at one electrode
<o:p></o:p></p>
<p class="MsoPlainText">>>> site and one given time is not representative of the underlying
<o:p></o:p></p>
<p class="MsoPlainText">>>> brain signal. Even if ICA was introducing minute distortion of the
<o:p></o:p></p>
<p class="MsoPlainText">>>> "true phase" at given channels, properly removing artifacts (which
<o:p></o:p></p>
<p class="MsoPlainText">>>> are 10-fold the amplitude of brain EEG signal) like ICA does is more
<o:p></o:p></p>
<p class="MsoPlainText">>>> important than preserving the exact phase at a given time. In your
<o:p></o:p></p>
<p class="MsoPlainText">>>> analogy of looking at stars, if you have a picture of a star, would
<o:p></o:p></p>
<p class="MsoPlainText">>>> you rather remove a visual artifact that is 10-fold the size of your
<o:p></o:p></p>
<p class="MsoPlainText">>>> original signal or continue to look at your original signal (not
<o:p></o:p></p>
<p class="MsoPlainText">>>> being able to see much because of the large artifact masking most of
<o:p></o:p></p>
<p class="MsoPlainText">>>> it). Even if ICA was introducing minute distortion in phase (which I
<o:p></o:p></p>
<p class="MsoPlainText">>>> do not believe it does because it deals with instantaneous mixtures)
<o:p></o:p></p>
<p class="MsoPlainText">>>> , it is worth it given the advantage it provides.<o:p></o:p></p>
<p class="MsoPlainText">>>><o:p> </o:p></p>
<p class="MsoPlainText">>>> The exact phase at one electrode site is not informative in itself.
<o:p></o:p></p>
<p class="MsoPlainText">>>> Differences in phase between 2 electrode sites is not informative
<o:p></o:p></p>
<p class="MsoPlainText">>>> either because there may be dozens of possibility for activity
<o:p></o:p></p>
<p class="MsoPlainText">>>> within the brain to generate such phase difference. One must move to
<o:p></o:p></p>
<p class="MsoPlainText">>>> the source level, and this is what ICA is doing (although see also
<o:p></o:p></p>
<p class="MsoPlainText">>>> my previous message).<o:p></o:p></p>
<p class="MsoPlainText">>>><o:p> </o:p></p>
<p class="MsoPlainText">>>> Best wishes,<o:p></o:p></p>
<p class="MsoPlainText">>>><o:p> </o:p></p>
<p class="MsoPlainText">>>> Arno<o:p></o:p></p>
<p class="MsoPlainText">>>><o:p> </o:p></p>
<p class="MsoPlainText">>>><o:p> </o:p></p>
<p class="MsoPlainText">>>>> On Jun 14, 2017, at 6:50 PM, Robert Thatcher <rwthatcher2@yahoo.com
<o:p></o:p></p>
<p class="MsoPlainText">>>>> <<a href="mailto:rwthatcher2@yahoo.com"><span style="color:windowtext;text-decoration:none">mailto:rwthatcher2@yahoo.com</span></a>>> wrote:<o:p></o:p></p>
<p class="MsoPlainText">>>>><o:p> </o:p></p>
<p class="MsoPlainText">>>> Arnaud,<o:p></o:p></p>
<p class="MsoPlainText">>>> It is interesting to see of the amount of distortion of phase
<o:p></o:p></p>
<p class="MsoPlainText">>>> differences of the original artifact free segments of the EEG record
<o:p></o:p></p>
<p class="MsoPlainText">>>> by ICA is a linear or nonlinear function of the number of ICA
<o:p></o:p></p>
<p class="MsoPlainText">>>> components that are removed to do the reconstruction of a different
<o:p></o:p></p>
<p class="MsoPlainText">>>> time series. You already showed that two ICA component removals
<o:p></o:p></p>
<p class="MsoPlainText">>>> results in more phase distortion than the removal of one ICA
<o:p></o:p></p>
<p class="MsoPlainText">>>> component. If you were to remove three and then reconstruct and
<o:p></o:p></p>
<p class="MsoPlainText">>>> then four and then five, etc and reconstruct and then attach the
<o:p></o:p></p>
<p class="MsoPlainText">>>> .edf files and share them with the forum then we can plot the
<o:p></o:p></p>
<p class="MsoPlainText">>>> magnitude of phase distortion of the artifact free sections of the
<o:p></o:p></p>
<p class="MsoPlainText">>>> original record due to the ICA reconstructions. Based on embedding
<o:p></o:p></p>
<p class="MsoPlainText">>>> theory one would expect a linear relationship but there may be a
<o:p></o:p></p>
<p class="MsoPlainText">>>> nonlinear relationship with an asymptote at about two removals given
<o:p></o:p></p>
<p class="MsoPlainText">>>> there are only 19 channels.<o:p></o:p></p>
<p class="MsoPlainText">>>><o:p> </o:p></p>
<p class="MsoPlainText">>>> Please try this experiement with one or more EEG dataset, the one
<o:p></o:p></p>
<p class="MsoPlainText">>>> that was produced by ICA reconstruction in Australial is a good
<o:p></o:p></p>
<p class="MsoPlainText">>>> starting point but it will be good to do this experiment with two or
<o:p></o:p></p>
<p class="MsoPlainText">>>> three other EEG recordings.<o:p></o:p></p>
<p class="MsoPlainText">>>><o:p> </o:p></p>
<p class="MsoPlainText">>>> Thank you for honest interest in exploring the extent of phase
<o:p></o:p></p>
<p class="MsoPlainText">>>> difference distortion by ICA so that we can better understand it.<o:p></o:p></p>
<p class="MsoPlainText">>>><o:p> </o:p></p>
<p class="MsoPlainText">>>> Best regards,<o:p></o:p></p>
<p class="MsoPlainText">>>><o:p> </o:p></p>
<p class="MsoPlainText">>>> Robert<o:p></o:p></p>
<p class="MsoPlainText">>>><o:p> </o:p></p>
<p class="MsoPlainText">>>> On Wednesday, June 14, 2017, 9:00:58 PM EDT, Robert Thatcher
<o:p></o:p></p>
<p class="MsoPlainText">>>> <<a href="mailto:rwthatcher2@yahoo.com%20%3cmailto:rwthatcher2@yahoo.com"><span style="color:windowtext;text-decoration:none">rwthatcher2@yahoo.com <mailto:rwthatcher2@yahoo.com</span></a>>> wrote:<o:p></o:p></p>
<p class="MsoPlainText">>>><o:p> </o:p></p>
<p class="MsoPlainText">>>><o:p> </o:p></p>
<p class="MsoPlainText">>>> Arnaud,<o:p></o:p></p>
<p class="MsoPlainText">>>> It does not make any difference which components that the
<o:p></o:p></p>
<p class="MsoPlainText">>>> scientist/clinicians removed because your own analyses confirmed
<o:p></o:p></p>
<p class="MsoPlainText">>>> phase difference distortion by ICA when your removed your own
<o:p></o:p></p>
<p class="MsoPlainText">>>> components. Please try different ICA component removal and attach
<o:p></o:p></p>
<p class="MsoPlainText">>>> the edf files to see if you can create a reconstruction of the time
<o:p></o:p></p>
<p class="MsoPlainText">>>> series that DOES NOT distort or corrupt the phase differences
<o:p></o:p></p>
<p class="MsoPlainText">>>> between channels in the original EEG recording. Up to this point in
<o:p></o:p></p>
<p class="MsoPlainText">>>> time you have resoundeding proven that ICA reconstruction oes
<o:p></o:p></p>
<p class="MsoPlainText">>>> distort phase differences no matter what reconstruction is used.<o:p></o:p></p>
<p class="MsoPlainText">>>><o:p> </o:p></p>
<p class="MsoPlainText">>>> It is important to recognize and to pubically accept that phase or
<o:p></o:p></p>
<p class="MsoPlainText">>>> time differences between channels in the EEG is due to physiological
<o:p></o:p></p>
<p class="MsoPlainText">>>> processes like differences in synaptic rise times, differences in
<o:p></o:p></p>
<p class="MsoPlainText">>>> synaptic summation times and differences in conduction velocities,
<o:p></o:p></p>
<p class="MsoPlainText">>>> etc. Even a small amount of adulteration or distortion of EEG phase
<o:p></o:p></p>
<p class="MsoPlainText">>>> differences is not good and must be avoided at all costs.<o:p></o:p></p>
<p class="MsoPlainText">>>><o:p> </o:p></p>
<p class="MsoPlainText">>>> This is analogous to the use of telescopes that measure phase
<o:p></o:p></p>
<p class="MsoPlainText">>>> differences in the spectrum from stars moving in the universe. If
<o:p></o:p></p>
<p class="MsoPlainText">>>> ICA were used to distort the phase differences in the spectrum
<o:p></o:p></p>
<p class="MsoPlainText">>>> measured by telescopes because one believes that all telescopes have
<o:p></o:p></p>
<p class="MsoPlainText">>>> artifact then we would not know huge amounts about the nature and
<o:p></o:p></p>
<p class="MsoPlainText">>>> future of the universe. The same is true for the human EEG.<o:p></o:p></p>
<p class="MsoPlainText">>>><o:p> </o:p></p>
<p class="MsoPlainText">>>> Bob<o:p></o:p></p>
<p class="MsoPlainText">>>><o:p> </o:p></p>
<p class="MsoPlainText">>>><o:p> </o:p></p>
<p class="MsoPlainText">>>> On Wednesday, June 14, 2017, 7:23:43 PM EDT, Arnaud Delorme
<o:p></o:p></p>
<p class="MsoPlainText">>>> <<a href="mailto:arno@ucsd.edu%20%3cmailto:arno@ucsd.edu"><span style="color:windowtext;text-decoration:none">arno@ucsd.edu <mailto:arno@ucsd.edu</span></a>>> wrote:<o:p></o:p></p>
<p class="MsoPlainText">>>><o:p> </o:p></p>
<p class="MsoPlainText">>>><o:p> </o:p></p>
<p class="MsoPlainText">>>>> Thank you for attaching your ICA reconstructed edf file. It
<o:p></o:p></p>
<p class="MsoPlainText">>>>> involved removal of two ICA components and the magnitude of changes
<o:p></o:p></p>
<p class="MsoPlainText">>>>> in phase differences between channels is greater than the one
<o:p></o:p></p>
<p class="MsoPlainText">>>>> provided by the scientists/clinicians in Australia that deleted
<o:p></o:p></p>
<p class="MsoPlainText">>>>> only one ICA component. This is consistent with Taken's theorum
<o:p></o:p></p>
<p class="MsoPlainText">>>>> and also differential geometry theorums dealing with manifold
<o:p></o:p></p>
<p class="MsoPlainText">>>>> mapings and Lie groups etc. I know for certain that they used ICA
<o:p></o:p></p>
<p class="MsoPlainText">>>>> and not PCA.<o:p></o:p></p>
<p class="MsoPlainText">>>><o:p> </o:p></p>
<p class="MsoPlainText">>>> Yes, I meant that the data is usually preprocessed by PCA before
<o:p></o:p></p>
<p class="MsoPlainText">>>> doing ICA in commercial softwares, which could be the problem
<o:p></o:p></p>
<p class="MsoPlainText">>>> (although I do not think it was in that case). We would need to see
<o:p></o:p></p>
<p class="MsoPlainText">>>> which components were removed.<o:p></o:p></p>
<p class="MsoPlainText">>>> Best wishes,<o:p></o:p></p>
<p class="MsoPlainText">>>><o:p> </o:p></p>
<p class="MsoPlainText">>>> Arno<o:p></o:p></p>
<p class="MsoPlainText">>>><o:p> </o:p></p>
<p class="MsoPlainText">>>><o:p> </o:p></p>
<p class="MsoPlainText">>>>> ICA is excellent in feature detection and the brain operates by
<o:p></o:p></p>
<p class="MsoPlainText">>>>> highly efficient sub-clusters of neurons extracting features, e.g.,
<o:p></o:p></p>
<p class="MsoPlainText">>>>> face recognition by combining features like eye brows, head shape,
<o:p></o:p></p>
<p class="MsoPlainText">>>>> ears, chin, etc<o:p></o:p></p>
<p class="MsoPlainText">>>>><o:p> </o:p></p>
<p class="MsoPlainText">>>>> Here is a url to a recent study showing that only 206 neurons are
<o:p></o:p></p>
<p class="MsoPlainText">>>>> necessary to encode face recognition in monkeys:<o:p></o:p></p>
<p class="MsoPlainText">>>>><o:p> </o:p></p>
<p class="MsoPlainText">>>>> <a href="http://dx.doi.org/10.1016/j.cell.2017.05.011">
<span style="color:windowtext;text-decoration:none">http://dx.doi.org/10.1016/j.cell.2017.05.011</span></a><o:p></o:p></p>
<p class="MsoPlainText">>>>><o:p> </o:p></p>
<p class="MsoPlainText">>>>> The individual face components are like ICA face components for
<o:p></o:p></p>
<p class="MsoPlainText">>>>> face recognition. However, the anterior temporal lobes are just
<o:p></o:p></p>
<p class="MsoPlainText">>>>> one node among many nodes in a network so that the monkey can make
<o:p></o:p></p>
<p class="MsoPlainText">>>>> the correct adaptive decisions in very short periods of time by
<o:p></o:p></p>
<p class="MsoPlainText">>>>> network coherence and phase locking and phase shifting with other
<o:p></o:p></p>
<p class="MsoPlainText">>>>> nodes in networks.<o:p></o:p></p>
<p class="MsoPlainText">>>>><o:p> </o:p></p>
<p class="MsoPlainText">>>>> The problem with ICA is in its use in artifact rejection and then
<o:p></o:p></p>
<p class="MsoPlainText">>>>> reconstruction of a new time series that results in a new time
<o:p></o:p></p>
<p class="MsoPlainText">>>>> series that is disconnected from brain network connectivity
<o:p></o:p></p>
<p class="MsoPlainText">>>>> dynamics of phase shift and phase lock and coherence and
<o:p></o:p></p>
<p class="MsoPlainText">>>>> cross-frequency coupling and phase amplitude coupling, etc.<o:p></o:p></p>
<p class="MsoPlainText">>>>><o:p> </o:p></p>
<p class="MsoPlainText">>>>> Thank you again and lets continue to seek answers to how best to
<o:p></o:p></p>
<p class="MsoPlainText">>>>> use ICA for network dynamics without adulterating the original time
<o:p></o:p></p>
<p class="MsoPlainText">>>>> and phase relations between parts of the brain.<o:p></o:p></p>
<p class="MsoPlainText">>>>><o:p> </o:p></p>
<p class="MsoPlainText">>>>> Robert<o:p></o:p></p>
<p class="MsoPlainText">>>>><o:p> </o:p></p>
<p class="MsoPlainText">>>>><o:p> </o:p></p>
<p class="MsoPlainText">>>>> On Wednesday, June 14, 2017, 6:21:32 PM EDT, Arnaud Delorme
<o:p></o:p></p>
<p class="MsoPlainText">>>>> <<a href="mailto:arno@ucsd.edu%20%3cmailto:arno@ucsd.edu"><span style="color:windowtext;text-decoration:none">arno@ucsd.edu <mailto:arno@ucsd.edu</span></a>>> wrote:<o:p></o:p></p>
<p class="MsoPlainText">>>>><o:p> </o:p></p>
<p class="MsoPlainText">>>>><o:p> </o:p></p>
<p class="MsoPlainText">>>>> Hi Robert,<o:p></o:p></p>
<p class="MsoPlainText">>>>><o:p> </o:p></p>
<p class="MsoPlainText">>>>>> The Australian data was analyzed by two scientists/clinicians in
<o:p></o:p></p>
<p class="MsoPlainText">>>>>> the audience of a workshop that I was doing in 2014 and they are
<o:p></o:p></p>
<p class="MsoPlainText">>>>>> the ones that did the ICA component selection using commercial
<o:p></o:p></p>
<p class="MsoPlainText">>>>>> WinEEG software and not me.<o:p></o:p></p>
<p class="MsoPlainText">>>>><o:p> </o:p></p>
<p class="MsoPlainText">>>>> Most commercial EEG software preprocess the data using PCA to
<o:p></o:p></p>
<p class="MsoPlainText">>>>> reduce the dimensionality of the data. The idea behind this is that
<o:p></o:p></p>
<p class="MsoPlainText">>>>> users should not have to go through as many components as channels.<o:p></o:p></p>
<p class="MsoPlainText">>>>> It is easier to have them select components within 5 or 10
<o:p></o:p></p>
<p class="MsoPlainText">>>>> exemplars. However this PCA dimension reduction can bias the
<o:p></o:p></p>
<p class="MsoPlainText">>>>> reconstruction (we have data to back this up but it is not
<o:p></o:p></p>
<p class="MsoPlainText">>>>> published yet).<o:p></o:p></p>
<p class="MsoPlainText">>>>><o:p> </o:p></p>
<p class="MsoPlainText">>>>> However, I do not think PCA dimension reduction before running ICA
<o:p></o:p></p>
<p class="MsoPlainText">>>>> was responsible for what you observed (because your data is very
<o:p></o:p></p>
<p class="MsoPlainText">>>>> clean and even after PCA and the 2 artifact components have huge
<o:p></o:p></p>
<p class="MsoPlainText">>>>> contribution to the data variance, you would get very similar
<o:p></o:p></p>
<p class="MsoPlainText">>>>> components). I think the WinEEG users simply did not select the
<o:p></o:p></p>
<p class="MsoPlainText">>>>> correct artifact components, or maybe WinEEG failed to implement
<o:p></o:p></p>
<p class="MsoPlainText">>>>> ICA correctly.<o:p></o:p></p>
<p class="MsoPlainText">>>>><o:p> </o:p></p>
<p class="MsoPlainText">>>>>> You are welcome to download NeuroGuide and install and launch and
<o:p></o:p></p>
<p class="MsoPlainText">>>>>> then paste the key A into an email to me. I have posted a
<o:p></o:p></p>
<p class="MsoPlainText">>>>>> tutorial on our webpage but I can create a better tutorial to
<o:p></o:p></p>
<p class="MsoPlainText">>>>>> reduce the learning curve. Similarly when I am able to
<o:p></o:p></p>
<p class="MsoPlainText">>>>>> concentrate on EEGlab then you can tutor me to reduce my learning
<o:p></o:p></p>
<p class="MsoPlainText">>>>>> curve. Here is a url to the download webpage:<o:p></o:p></p>
<p class="MsoPlainText">>>>>> <a href="http://www.appliedneuroscience.com/Download_NeuroGuide.htm">
<span style="color:windowtext;text-decoration:none">http://www.appliedneuroscience.com/Download_NeuroGuide.htm</span></a><o:p></o:p></p>
<p class="MsoPlainText">>>>>><o:p> </o:p></p>
<p class="MsoPlainText">>>>>> At the end of the day together lets find ways to use the full
<o:p></o:p></p>
<p class="MsoPlainText">>>>>> power of ICA to explore network dynamics which is my favorite
<o:p></o:p></p>
<p class="MsoPlainText">>>>>> topic and also one that future science depends on.<o:p></o:p></p>
<p class="MsoPlainText">>>>><o:p> </o:p></p>
<p class="MsoPlainText">>>>> Yes, I agree on that view. Exploring network dynamics with ICA is
<o:p></o:p></p>
<p class="MsoPlainText">>>>> not an easy topic. The trend these days is not to use ICA for
<o:p></o:p></p>
<p class="MsoPlainText">>>>> connectivity analysis but instead define regions of interest and
<o:p></o:p></p>
<p class="MsoPlainText">>>>> compute pairwise connectivity between all brain regions as in this
<o:p></o:p></p>
<p class="MsoPlainText">>>>> recent paper <a href="https://www.ncbi.nlm.nih.gov/pubmed/28300640">
<span style="color:windowtext;text-decoration:none">https://www.ncbi.nlm.nih.gov/pubmed/28300640</span></a>. What can
<o:p></o:p></p>
<p class="MsoPlainText">>>>> be done is to use ICA components to define these regions and
<o:p></o:p></p>
<p class="MsoPlainText">>>>> compute activity in these regions. It is an open area of research.<o:p></o:p></p>
<p class="MsoPlainText">>>>><o:p> </o:p></p>
<p class="MsoPlainText">>>>> Best wishes,<o:p></o:p></p>
<p class="MsoPlainText">>>>><o:p> </o:p></p>
<p class="MsoPlainText">>>>> Arno<o:p></o:p></p>
<p class="MsoPlainText">>>>><o:p> </o:p></p>
<p class="MsoPlainText">>>>>><o:p> </o:p></p>
<p class="MsoPlainText">>>>>> Best regards,<o:p></o:p></p>
<p class="MsoPlainText">>>>>><o:p> </o:p></p>
<p class="MsoPlainText">>>>>> Robert<o:p></o:p></p>
<p class="MsoPlainText">>>>>><o:p> </o:p></p>
<p class="MsoPlainText">>>>>><o:p> </o:p></p>
<p class="MsoPlainText">>>>>> On Wednesday, June 14, 2017, 4:40:43 PM EDT, Arnaud Delorme
<o:p></o:p></p>
<p class="MsoPlainText">>>>>> <<a href="mailto:arno@ucsd.edu%20%3cmailto:arno@ucsd.edu"><span style="color:windowtext;text-decoration:none">arno@ucsd.edu <mailto:arno@ucsd.edu</span></a>>> wrote:<o:p></o:p></p>
<p class="MsoPlainText">>>>>><o:p> </o:p></p>
<p class="MsoPlainText">>>>>><o:p> </o:p></p>
<p class="MsoPlainText">>>>>> Dear Robert,<o:p></o:p></p>
<p class="MsoPlainText">>>>>><o:p> </o:p></p>
<p class="MsoPlainText">>>>>> There does seem to be a phase difference in your powerpoint.
<o:p></o:p></p>
<p class="MsoPlainText">>>>>> However, it is important to know which ICA component you removed
<o:p></o:p></p>
<p class="MsoPlainText">>>>>> to understand why this is the case. Are you sure these were
<o:p></o:p></p>
<p class="MsoPlainText">>>>>> artifactual components? Removing brain components may alter the
<o:p></o:p></p>
<p class="MsoPlainText">>>>>> phase of the signal recorded on the scalp (it would be as if you
<o:p></o:p></p>
<p class="MsoPlainText">>>>>> were removing from the scalp signal the contribution of a brain
<o:p></o:p></p>
<p class="MsoPlainText">>>>>> area). Without that information, it is not possible to figure out
<o:p></o:p></p>
<p class="MsoPlainText">>>>>> the origin of the phase difference.<o:p></o:p></p>
<p class="MsoPlainText">>>>>><o:p> </o:p></p>
<p class="MsoPlainText">>>>>> This seems to be the same data you shared yesterday. I have looked
<o:p></o:p></p>
<p class="MsoPlainText">>>>>> at it. Black is before ICA and red after removing the 2 eye
<o:p></o:p></p>
<p class="MsoPlainText">>>>>> components. You can see that there is no phase shift at 102.43
<o:p></o:p></p>
<p class="MsoPlainText">>>>>> second after I remove the two artifactual ICA components. I have
<o:p></o:p></p>
<p class="MsoPlainText">>>>>> provided the code in my email yesterday if you want to reproduce
<o:p></o:p></p>
<p class="MsoPlainText">>>>>> this result in EEGLAB.<o:p></o:p></p>
<p class="MsoPlainText">>>>>><o:p> </o:p></p>
<p class="MsoPlainText">>>>>> Best wishes,<o:p></o:p></p>
<p class="MsoPlainText">>>>>><o:p> </o:p></p>
<p class="MsoPlainText">>>>>> Arno<o:p></o:p></p>
<p class="MsoPlainText">>>>>><o:p> </o:p></p>
<p class="MsoPlainText">>>>>> <ICA_phase_example.png><o:p></o:p></p>
<p class="MsoPlainText">>>>>><o:p> </o:p></p>
<p class="MsoPlainText">>>>>>> On Jun 14, 2017, at 11:14 AM, Robert Thatcher <o:p>
</o:p></p>
<p class="MsoPlainText">>>>>>> <<a href="mailto:rwthatcher2@yahoo.com%20%3cmailto:rwthatcher2@yahoo.com"><span style="color:windowtext;text-decoration:none">rwthatcher2@yahoo.com <mailto:rwthatcher2@yahoo.com</span></a>>> wrote:<o:p></o:p></p>
<p class="MsoPlainText">>>>>>><o:p> </o:p></p>
<p class="MsoPlainText">>>>>>> <Example of Phase Differences at 1min & 46 seconds.pptx><o:p></o:p></p>
<p class="MsoPlainText">>>>>><o:p> </o:p></p>
<p class="MsoPlainText">>>>>> <ICA_phase_example.png><o:p></o:p></p>
<p class="MsoPlainText">>>>><o:p> </o:p></p>
<p class="MsoPlainText">>>><o:p> </o:p></p>
<p class="MsoPlainText">>>><o:p> </o:p></p>
<p class="MsoPlainText">>><o:p> </o:p></p>
<p class="MsoPlainText">>><o:p> </o:p></p>
<p class="MsoPlainText">>><o:p> </o:p></p>
<p class="MsoPlainText">>> _______________________________________________<o:p></o:p></p>
<p class="MsoPlainText">>> Eeglablist page: <a href="http://sccn.ucsd.edu/eeglab/eeglabmail.html">
<span style="color:windowtext;text-decoration:none">http://sccn.ucsd.edu/eeglab/eeglabmail.html</span></a><o:p></o:p></p>
<p class="MsoPlainText">>> To unsubscribe, send an empty email to <o:p></o:p></p>
<p class="MsoPlainText">>> <a href="mailto:eeglablist-unsubscribe@sccn.ucsd.edu">
<span style="color:windowtext;text-decoration:none">eeglablist-unsubscribe@sccn.ucsd.edu</span></a><o:p></o:p></p>
<p class="MsoPlainText">>> For digest mode, send an email with the subject "set digest mime" to
<o:p></o:p></p>
<p class="MsoPlainText">>> <a href="mailto:eeglablist-request@sccn.ucsd.edu"><span style="color:windowtext;text-decoration:none">eeglablist-request@sccn.ucsd.edu</span></a><o:p></o:p></p>
<p class="MsoPlainText">><o:p> </o:p></p>
<p class="MsoPlainText"><o:p> </o:p></p>
<p class="MsoPlainText">_______________________________________________<o:p></o:p></p>
<p class="MsoPlainText">Eeglablist page: <a href="http://sccn.ucsd.edu/eeglab/eeglabmail.html">
<span style="color:windowtext;text-decoration:none">http://sccn.ucsd.edu/eeglab/eeglabmail.html</span></a><o:p></o:p></p>
<p class="MsoPlainText">To unsubscribe, send an empty email to <a href="mailto:eeglablist-unsubscribe@sccn.ucsd.edu">
<span style="color:windowtext;text-decoration:none">eeglablist-unsubscribe@sccn.ucsd.edu</span></a><o:p></o:p></p>
<p class="MsoPlainText">For digest mode, send an email with the subject "set digest mime" to
<a href="mailto:eeglablist-request@sccn.ucsd.edu"><span style="color:windowtext;text-decoration:none">eeglablist-request@sccn.ucsd.edu</span></a><o:p></o:p></p>
</div>
</body>
</html>