<html><head></head><body><div style="font-family:Helvetica Neue, Helvetica, Arial, sans-serif;font-size:small;"><div><div>Hi Joe,</div></div><div> As Ramesh mentioned and also Paul Nunez "Electrical Fields of the Brain" 1981 eplain that the physics of the electrical potential between the retina and cornea decreases with distance from its source but is inextricably superimposed upon the EEG waves especially in the delta and theta bands and less in the alpha2 to beta to gamma bands. No regression method can remove the superimposition and it is best to simply use improved recording hygiene methods and delete and sections of the EEG record that have artifact no matter what the artifact is.</div><div><br></div><div>I am not sure what methods that you are using ("<span style="color: rgb(38, 40, 42); font-size: 13px;">Multiple Algorithm Artifact Correction or MAAC"?) </span>but it is important to recognize that Mathematical Regression methods based on adjacent time points or smoothing or decomposition and reconstruction with refducefd components or with less information to then alter the entire EEG record incuding the physics of the artifact free parts of the recoring is bad science in my opinion. Most scientists insist that one must preserve the physics of the brain at all costs and not conduct analyses based on an altered time series. There is no equivalence of the physics ground truth by an artificial replacement base on an hypothesis or belief that all EEG is noise or artifact.</div><div><br></div><div><br></div><div>Best wishes.</div><div><br></div><div>Robert</div><div><br></div><div id="ydpe431f14dyahoo_quoted_9723898929" class="ydpe431f14dyahoo_quoted"><div style="font-family:'Helvetica Neue', Helvetica, Arial, sans-serif;font-size:13px;color:#26282a;"><div>On Sunday, July 9, 2017, 7:34:33 PM EDT, Joseph Dien <jdien07@mac.com> wrote:</div><div><br></div><div><br></div><div><div id="ydpe431f14dyiv0019821171"><div><div class="ydpe431f14dyiv0019821171">I’ve been working on an artifact correction manuscript with an EEG plus eye-tracker dataset. The corneo-retinal dipole (CRD) artifact reflects eye position not eye movement so it’ll be constantly present even if the eyes maintain fixation. My take is that each artifact is best corrected with a method tailored to its unique characteristics. I’ve got a working algorithm (Multiple Algorithm Artifact Correction or MAAC) implemented in my EP Toolkit (<a shape="rect" class="ydpe431f14dyiv0019821171" href="https://sourceforge.net/projects/erppcatoolkit/?source=navbar" rel="nofollow" target="_blank">https://sourceforge.net/projects/erppcatoolkit/?source=navbar</a>) but I need to finish testing and validation. Ideally it also needs to be tested out on different systems and montages to see how well it generalizes.</div><div class="ydpe431f14dyiv0019821171"><br clear="none" class="ydpe431f14dyiv0019821171"></div><div class="ydpe431f14dyiv0019821171">Cheers!</div><div class="ydpe431f14dyiv0019821171"><br clear="none" class="ydpe431f14dyiv0019821171"></div><div class="ydpe431f14dyiv0019821171">Joe</div><br clear="none" class="ydpe431f14dyiv0019821171"><div class="ydpe431f14dyiv0019821171yqt0940008496" id="ydpe431f14dyiv0019821171yqt90600"><div><blockquote class="ydpe431f14dyiv0019821171" type="cite"><div class="ydpe431f14dyiv0019821171">On Jul 5, 2017, at 07:37, Clayton Hickey <<a shape="rect" class="ydpe431f14dyiv0019821171" href="mailto:cmhickey@gmail.com" rel="nofollow" target="_blank">cmhickey@gmail.com</a>> wrote:</div><br clear="none" class="ydpe431f14dyiv0019821171Apple-interchange-newline"><div class="ydpe431f14dyiv0019821171"></div></blockquote></div></div></div><div class="ydpe431f14dyiv0019821171yqt0940008496" id="ydpe431f14dyiv0019821171yqt34220"><div><div class="ydpe431f14dyiv0019821171" style="word-wrap:break-word;"><div class="ydpe431f14dyiv0019821171">Hi Marius, Hi All, </div><div class="ydpe431f14dyiv0019821171"><br clear="none" class="ydpe431f14dyiv0019821171"></div><div class="ydpe431f14dyiv0019821171">In reading your message, Marius, it struck me that I might make a small contribution to this discussion. </div><div class="ydpe431f14dyiv0019821171"><br clear="none" class="ydpe431f14dyiv0019821171"></div><div class="ydpe431f14dyiv0019821171">The idea that EEG can be free of eye movement artifacts is something that can be empirically studied… what is needed is an independent measure of eye movements. Many of us have this in concurrent eye-tracking data. </div><div class="ydpe431f14dyiv0019821171"><br clear="none" class="ydpe431f14dyiv0019821171"></div><div class="ydpe431f14dyiv0019821171">In my group we use an eye-link 1000 and the wonderful EYE-EEG toolbox to co-register tracker data with concurrently recorded 64-channel EEG (Dimigen et al. & Kliegl, 2011, JEP:GEN). As per the pipeline suggested by Dimigen et al., we derive ICA components from the EEG data and identify artifactual components as those that strongly correlate with the tracker data. This removes some of the ambiguity inherent to identification of artifacts from EOG, topography etc, because the IR-camera-based tracker data is independent of the encephalic voltage data. </div><div class="ydpe431f14dyiv0019821171"><br clear="none" class="ydpe431f14dyiv0019821171"></div><div class="ydpe431f14dyiv0019821171">For the fun of it (procrastination, get thee behind me) I just had a quick look at a few datasets, extracting intervals where participants maintained fixation. I ran infomax ICA on this data and identified ‘artifactual’ components as those that correlated with the eye-tracking data. I consistently got at least a couple of components that strongly correlated with the tracker signal. This in spite of the fact that the data was putatively ‘artifact-free’ and collected during consistent maintenance of fixation. </div><div class="ydpe431f14dyiv0019821171"><br clear="none" class="ydpe431f14dyiv0019821171"></div><div class="ydpe431f14dyiv0019821171">Some of this variance may come from the brain. Dilation and microsaccades are likely to have related neural activity that could be confounded with muscle / rotation-of-eyeball activity by ICA. But the brain activity seems to be a very small source of variance. The components load heavily on frontal sensors close to the eyes, sharply tapering off along the anterior-to-posterior axis with very low or zero weights at posterior electrodes. I would expect microsaccades and dilation to elicit brain activity in subcortex (diffuse topography) or over early visual areas (posterior topography). </div><div class="ydpe431f14dyiv0019821171"><br clear="none" class="ydpe431f14dyiv0019821171"></div><div class="ydpe431f14dyiv0019821171">So... the eyes are doing something to the EEG signal during the maintenance of fixation, and ICA can pick this up. This seems to be primarily artifact from muscle / eyeball rotation. </div><div class="ydpe431f14dyiv0019821171"><br clear="none" class="ydpe431f14dyiv0019821171"></div><div class="ydpe431f14dyiv0019821171">best, clayton</div><div class="ydpe431f14dyiv0019821171"><br clear="none" class="ydpe431f14dyiv0019821171"></div><div class="ydpe431f14dyiv0019821171"><br clear="none" class="ydpe431f14dyiv0019821171"></div><div class="ydpe431f14dyiv0019821171"><br clear="none" class="ydpe431f14dyiv0019821171"></div><div class="ydpe431f14dyiv0019821171"><br clear="none" class="ydpe431f14dyiv0019821171"></div><br clear="none" class="ydpe431f14dyiv0019821171"><div class="ydpe431f14dyiv0019821171"><blockquote class="ydpe431f14dyiv0019821171" type="cite"><div class="ydpe431f14dyiv0019821171"> 8. Re: ICA Misinformation (Marius Klug)<br clear="none" class="ydpe431f14dyiv0019821171"><br clear="none" class="ydpe431f14dyiv0019821171"><div class="ydpe431f14dyiv0019821171" style="margin-top:0px;margin-right:0px;margin-bottom:0px;margin-left:0px;"><span class="ydpe431f14dyiv0019821171" style="font-family:-webkit-system-font, Helvetica Neue, Helvetica, sans-serif;"><b class="ydpe431f14dyiv0019821171">From: </b></span><span class="ydpe431f14dyiv0019821171" style="font-family:-webkit-system-font, Helvetica Neue, Helvetica, sans-serif;">Marius Klug <<a shape="rect" class="ydpe431f14dyiv0019821171" href="mailto:marius.s.klug@gmail.com" rel="nofollow" target="_blank">marius.s.klug@gmail.com</a>><br clear="none" class="ydpe431f14dyiv0019821171"></span></div><div class="ydpe431f14dyiv0019821171" style="margin-top:0px;margin-right:0px;margin-bottom:0px;margin-left:0px;"><span class="ydpe431f14dyiv0019821171" style="font-family:-webkit-system-font, Helvetica Neue, Helvetica, sans-serif;"><b class="ydpe431f14dyiv0019821171">Subject: </b></span><span class="ydpe431f14dyiv0019821171" style="font-family:-webkit-system-font, Helvetica Neue, Helvetica, sans-serif;"><b class="ydpe431f14dyiv0019821171">Re: [Eeglablist] ICA Misinformation</b><br clear="none" class="ydpe431f14dyiv0019821171"></span></div><div class="ydpe431f14dyiv0019821171" style="margin-top:0px;margin-right:0px;margin-bottom:0px;margin-left:0px;"><span class="ydpe431f14dyiv0019821171" style="font-family:-webkit-system-font, Helvetica Neue, Helvetica, sans-serif;"><b class="ydpe431f14dyiv0019821171">Date: </b></span><span class="ydpe431f14dyiv0019821171" style="font-family:-webkit-system-font, Helvetica Neue, Helvetica, sans-serif;">June 30, 2017 at 3:02:52 PM GMT+2<br clear="none" class="ydpe431f14dyiv0019821171"></span></div><div class="ydpe431f14dyiv0019821171" style="margin-top:0px;margin-right:0px;margin-bottom:0px;margin-left:0px;"><span class="ydpe431f14dyiv0019821171" style="font-family:-webkit-system-font, Helvetica Neue, Helvetica, sans-serif;"><b class="ydpe431f14dyiv0019821171">To: </b></span><span class="ydpe431f14dyiv0019821171" style="font-family:-webkit-system-font, Helvetica Neue, Helvetica, sans-serif;">otte georges <<a shape="rect" class="ydpe431f14dyiv0019821171" href="mailto:georges.otte@pandora.be" rel="nofollow" target="_blank">georges.otte@pandora.be</a>><br clear="none" class="ydpe431f14dyiv0019821171"></span></div><div class="ydpe431f14dyiv0019821171" style="margin-top:0px;margin-right:0px;margin-bottom:0px;margin-left:0px;"><span class="ydpe431f14dyiv0019821171" style="font-family:-webkit-system-font, Helvetica Neue, Helvetica, sans-serif;"><b class="ydpe431f14dyiv0019821171">Cc: </b></span><span class="ydpe431f14dyiv0019821171" style="font-family:-webkit-system-font, Helvetica Neue, Helvetica, sans-serif;">eeglablist <<a shape="rect" class="ydpe431f14dyiv0019821171" href="mailto:eeglablist@sccn.ucsd.edu" rel="nofollow" target="_blank">eeglablist@sccn.ucsd.edu</a>>, <a shape="rect" class="ydpe431f14dyiv0019821171" href="mailto:pnunez@tulane.edu" rel="nofollow" target="_blank">pnunez@tulane.edu</a>, <a shape="rect" class="ydpe431f14dyiv0019821171" href="mailto:r.srinivasan@uci.edu" rel="nofollow" target="_blank">r.srinivasan@uci.edu</a>, Georges Otte <<a shape="rect" class="ydpe431f14dyiv0019821171" href="mailto:georges.otte@telenet.be" rel="nofollow" target="_blank">georges.otte@telenet.be</a>><br clear="none" class="ydpe431f14dyiv0019821171"></span></div><br clear="none" class="ydpe431f14dyiv0019821171"><br clear="none" class="ydpe431f14dyiv0019821171"><div class="ydpe431f14dyiv0019821171" dir="ltr"><div class="ydpe431f14dyiv0019821171"><div class="ydpe431f14dyiv0019821171">Dear Georges, all,<br clear="none" class="ydpe431f14dyiv0019821171"><br clear="none" class="ydpe431f14dyiv0019821171"></div>I am not sure if I missed something, this is the last email about the subject to my knowledge, so I'm answering to it now since I feel it should get one. My answers are in the text below, I think this will also be my first and last email on this topic. Everything I write is my own personal opinion.<br clear="none" class="ydpe431f14dyiv0019821171"></div><div class="ydpe431f14dyiv0019821171"><div class="ydpe431f14dyiv0019821171"><div class="ydpe431f14dyiv0019821171gmail_extra"><br clear="none" class="ydpe431f14dyiv0019821171"><div class="ydpe431f14dyiv0019821171gmail_quote">2017-06-25 22:18 GMT+02:00 otte georges <span class="ydpe431f14dyiv0019821171" dir="ltr"><<a shape="rect" class="ydpe431f14dyiv0019821171" href="mailto:georges.otte@pandora.be" rel="nofollow" target="_blank">georges.otte@pandora.be</a>></span>:<br clear="none" class="ydpe431f14dyiv0019821171"><blockquote class="ydpe431f14dyiv0019821171gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex;">Dear all<br clear="none" class="ydpe431f14dyiv0019821171">
<br clear="none" class="ydpe431f14dyiv0019821171">
<br clear="none" class="ydpe431f14dyiv0019821171">
Although I find the subject very interesting, the "cynicism" embedded in some responses feels a bit vexatious and this is something that we should try to avoid in important scientific discussions not unlike this on concerning the possible side effects of ICA signal reconstruction and the model that ICA approach of EEG is based upon<br clear="none" class="ydpe431f14dyiv0019821171"></blockquote><div class="ydpe431f14dyiv0019821171"> </div><div class="ydpe431f14dyiv0019821171">I find it strange that you attribute non-scientific way of arguing mostly to the ICA experts' side, having followed the complete discussion. I found Robert to be rather bold in his statements and kind of immune to well-formulated and thought through counter-arguments to his statements, even if they are backed up by data and plots. No wonder the tone in the discussion has turned a bit darker after a while, if you ask me. We all have a lot of things on our plate and must allocate our resources, some of the writers have allocated a lot of theirs to the topic and their time and reasoning has been met with denial disguised as interest, and statements have been repeated even though they were addressed already.<br clear="none" class="ydpe431f14dyiv0019821171"><br clear="none" class="ydpe431f14dyiv0019821171"></div><blockquote class="ydpe431f14dyiv0019821171gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex;">
<br clear="none" class="ydpe431f14dyiv0019821171">
An example: In Prof Debners reply to Dr. Thatcher he invokes the question of the possibility of invisible ECG artefact contamination of the EEG recordings and then continues by saying " because you don't see heart-electrical activity in your "original" 19-ch 10/20 recordings, can you seriously claim that the heart of your participant was not beating during recording?" Ha.. Ha ... ha<br clear="none" class="ydpe431f14dyiv0019821171"></blockquote><div class="ydpe431f14dyiv0019821171"> </div><div class="ydpe431f14dyiv0019821171">I suspect that not even to be specifically humorous, but an actual argument. Robert and other in his side of argument several times spoke of "artifact-free data segments" which would be distorted by taking out an ICA component, even if it was artifact-free. The argument here is, and I can't believe that I have to write this again, since it's been written so often already, that no such data exists, end of story. EEG is just an electrical recording. Electrical signals are abundant and generated constantly during the recording by a wide number of sources, not only the brain. One example here being heart beats. Now obviously those signals do have an amplitude and a phase, and since they are superpositioning the brain signals and thus part of the EEG recording, they will - no way around this! - have an impact on the phase (and amplitude) of the data set. The data set is thus _continuously_ distorted by the artifact. So by saying that a data set has artifact-free time periods, Robert implies that the heart has stopped beating, which - as Stefan pointed out - is rather unlikely (now here is a bit of humor on my side, I admit ;-) ). <br clear="none" class="ydpe431f14dyiv0019821171"><br clear="none" class="ydpe431f14dyiv0019821171">The thing is: This applies not only to heart beats but to each and every single artifact source that can be recorded by the EEG! So, even if you have no super strong eye blink artifacts, the eyeballs and eye muscles will still continuously contaminate the data - a bit at least. Since the IC for eyes in a 19 channel data set does likely contain both blinks and movements, taking it out will also take out the continuous eye movement artifacts produced by all kinds of drifts, saccades, and microsaccades, so it will take out parts of the data in the complete set! Now, since the eye signal again has a phase and an amplitude, the resulting signal will be distorted - but to the better, not to the worse! The new signal does not contain the artifacts of the eyes, AT ALL TIME POINTS, at least to the degree that ICA was able to separate them. So yes, the phase will be distorted, AND THAT IS A GOOD THING! The IC topography should, however, not extend to the parietal regions, so there should not be any distortions there. You can also check the continuous IC activity (meaning eye activity) and at times where it is very low, the distortion in the reconstructed sensor channels should be equally low.<br clear="none" class="ydpe431f14dyiv0019821171"><br clear="none" class="ydpe431f14dyiv0019821171"></div><blockquote class="ydpe431f14dyiv0019821171gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex;">
<br clear="none" class="ydpe431f14dyiv0019821171">
The question of possible ECG contamination that Stefan raioses is indeed very important but instead of just joking about it, perhaps it would be even better to (also) propose a possible way to solve this question. </blockquote><div class="ydpe431f14dyiv0019821171"><br clear="none" class="ydpe431f14dyiv0019821171"></div><div class="ydpe431f14dyiv0019821171">The proposed way is to take out the artifact via ICA. That's doable with high-density EEG, I find ECG ICs in my subjects without trouble usually. We use a neck band with electrodes in addition, that does help a lot, I suppose...<br clear="none" class="ydpe431f14dyiv0019821171"></div><div class="ydpe431f14dyiv0019821171"> </div><blockquote class="ydpe431f14dyiv0019821171gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex;">I am not in the same scientific league that Prof Debner is in but should that question pop up I would try to search for a test or setup that could clarify it or lead to a solution. Is that not the core of science? I can imagine that at Stefan's university, open hart operations are being performed and it would perhaps be feasible to record EEG before and during the period that the heart function is stopped and a heart-lung machine takes over. I agree that this is not a trivial task due to -among others- the electrical and other potential interferences in the operational theatre and it would require a good research capacity but I am confident that Stefan is up to that task. </blockquote><div class="ydpe431f14dyiv0019821171"><br clear="none" class="ydpe431f14dyiv0019821171"></div><div class="ydpe431f14dyiv0019821171">This is indeed a non-trivial task and you must be joking to casually assume that Stefan will just lay down other projects at hand to have an experiment to prove small subleties of heart beat artifact attenuation via ICA - that are obvious to the skilled user and will likely not convince Robert anyway - in a setup that is likely to contain artifacts from all the machines around that are so strong that not much else can be recorded as a start. I can see the heart beat in the neck electrodes even visually WHILE RECORDING without any preprocessing and ICA, there is no way to deny this.<br clear="none" class="ydpe431f14dyiv0019821171"> <br clear="none" class="ydpe431f14dyiv0019821171"></div><blockquote class="ydpe431f14dyiv0019821171gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex;">I think the results would seriously help the community to come closer to a more definite answer. The same could be done for eye movement artefacts. In ophthalmologic clinics all kind of eye diseased are treated. What about people with no corpus vitreum: do they have other eye movement artefacts before and after operation? It is often repeated that eye movement potentials are caused by the cornea-retina dipole. Is that so? Have You people tested this or do You just assumed this is so because it sounds so very plausible. I heard at least one ICA specialist state that this is a mémé and that the corpus vitreous is the originator of those eye movement potentials. I know that all handbooks of EEG claim otherwise so why not test this? We should be very critical with plausible models and try to test the stuff that we all take for granted.<br clear="none" class="ydpe431f14dyiv0019821171"></blockquote><div class="ydpe431f14dyiv0019821171"><br clear="none" class="ydpe431f14dyiv0019821171"></div><div class="ydpe431f14dyiv0019821171">Inhowfar does it matter if the eye ball movement artifact is generated in to corpus vitreum? I don't see this to be a counter to the notion of a dipolar eye ball, it's just a specification. It doesn't matter. The eye moves -> signals detectable in the EEG which should be taken out. You're giving a straw man argument here.<br clear="none" class="ydpe431f14dyiv0019821171"></div><div class="ydpe431f14dyiv0019821171"> </div><blockquote class="ydpe431f14dyiv0019821171gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex;">
<br clear="none" class="ydpe431f14dyiv0019821171">
Scientists with access to large university hospital facilities are empowered to solve these kinds of question. They s are the people we clinicians need to get questions solved and personal beliefs or opinions are less relevant in those matters. I can of course understand humor even with a sniff of cynicism but it should not stop there. Humor should not serve to avoid tricky questions as a rapier in a duel but should be followed by an indepth endeavor to prove one point of view in a way that everybody can agree as solid evidence.. I do not want to start a flame war nor do I have revendicate intentions but as a clinician I feel that we deserve a more serious solution from the side of You scientists.<br clear="none" class="ydpe431f14dyiv0019821171">
The challenge is open and until now unmet.<br clear="none" class="ydpe431f14dyiv0019821171"></blockquote><div class="ydpe431f14dyiv0019821171"><br clear="none" class="ydpe431f14dyiv0019821171"></div><div class="ydpe431f14dyiv0019821171">Rest assured that the majority of "us scientists" are giving our best to investigate our own methods and continuously try to find better ways to analyze our data, since we are usually our own most harsh critic and overthink everything (at least I am and do)... You however should not have stopped with your assessment of cynical humor above but thought further why this assumed humorous argument was indeed spot on.<br clear="none" class="ydpe431f14dyiv0019821171"><br clear="none" class="ydpe431f14dyiv0019821171"></div><div class="ydpe431f14dyiv0019821171">Our serious solution is right there: High-density EEG with spatial filtering of some kind (ICA, spatiospectral decomposition, joint decorrelation, ...), source localization, source-level analysis or back-projection of sources of interest (without artifacts) into the sensor-level. Investigating all ICs and seeing the obvious: There are noise and artifact signals in the EEG data all the time. Of course it would be nice to have and provide you with a really fancy way to have perfect brain-only data with 19 electrodes, but at least for the time being, that is not possible. Also, I am wondering if you would then still say that it is different from the raw data and thus bad, ignoring the arguments given... Which brings me to my last point:<br clear="none" class="ydpe431f14dyiv0019821171"></div><div class="ydpe431f14dyiv0019821171"> </div><blockquote class="ydpe431f14dyiv0019821171gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex;">
<br clear="none" class="ydpe431f14dyiv0019821171">
When Dr Thatcher launches a potential serious problem, it did not help me very much to hear his data being considered flawed, to hear preemptive suppositions that his choice of ICA algorithm was wrong, that his selection was maybe erroneous etc. One cannot solve a question or a problem by attacking the man who asks the question. It feels like shooting at the pianist because one dislikes the melody. Please bring forward more sturdy arguments and proof that ICA reconstruction does not alter phase.<br clear="none" class="ydpe431f14dyiv0019821171"></blockquote><div class="ydpe431f14dyiv0019821171"><br clear="none" class="ydpe431f14dyiv0019821171"></div><div class="ydpe431f14dyiv0019821171">I don't remember anyone attacking Robert as a person, except for the harsh tone and bold statements. The "preemptive suppositions" were indeed important and correct: He did not provide any information about the ICA decomposition and the ICs that were taken out. Since eye movements should not extend to the parietal parts of the EEG, but the phase has been altered in those channels after back projection, there must have been other ICs taken out (possibly containing brain signals) or the decomposition was seriously bad (mixing eye and brain signals), which led to doubts of the decomposition algorithm, the choice of ICs, and to the general statement that 19 channels are rather few for a good ICA decomposition in general. In fact, the discussion about his data began when Arno could clearly show that taking out eye movement artifacts did NOT distort phase to a relevant degree in times other than eye artifacts occurring in one of his earlier emails. How Robert can interpret those figures to suit his own argument is a mystery to me. Arno was, so to say, not able to replicate your bug and then the list proceeded to search for other reasons, but we were stuck by the fact, that there was no movement on Roberts side, no careful examination of the facts and arguments that have been laid out.<br clear="none" class="ydpe431f14dyiv0019821171"><br clear="none" class="ydpe431f14dyiv0019821171"></div><div class="ydpe431f14dyiv0019821171">And, again, I can't believe that I have to write this again: YES, TAKING OUT ICs WILL ALTER PHASE! This is not even an argument! The thing is a) that this is something good for the time points where the artifacts occur, because in fact the original data has been distorted by a the artifacts, eye in this case, and taking out the artifacts will restore the more correct brain signal phases and b) the extent and spatial distribution of the distortion presented by Robert led to the conclusion that something must have been done incorrectly or at least not with a lot of scrutiny. Who knows... workshop participants might make mistakes, they learn, that's what they are there for. But this was not addressed by Robert at all.<br clear="none" class="ydpe431f14dyiv0019821171"><br clear="none" class="ydpe431f14dyiv0019821171"></div><div class="ydpe431f14dyiv0019821171">The arguments and proof have been as sturdy as it can get and the fact that you and Robert don't see this just proves that your heads are sturdier, sorry to be so clear in my expression. I have nothing personal against either you or Robert, I don't know you at all and I don't judge you as persons, to be clear, we might get along well over a beer in a bar. I just judge the facts and your reactions to them in the emails I read.<br clear="none" class="ydpe431f14dyiv0019821171"></div><div class="ydpe431f14dyiv0019821171"> </div><blockquote class="ydpe431f14dyiv0019821171gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex;">
<br clear="none" class="ydpe431f14dyiv0019821171">
Considering the EEG as a chaotic multivariate time series continuously contaminated by visible and invisible artefacts is a model that if true justifies the use of ICA but as such a model that needs rock solid scientific confirmation and not just authority based replication of opinions.<br clear="none" class="ydpe431f14dyiv0019821171"></blockquote><div class="ydpe431f14dyiv0019821171"><br clear="none" class="ydpe431f14dyiv0019821171"></div><div class="ydpe431f14dyiv0019821171">This is the thing. That's what Stefan meant with his "cynical humor argument": You cannot seriously argue that physiological signals all of a sudden cease to exist every now and then and then come back to life again later. Let me be clear: <br clear="none" class="ydpe431f14dyiv0019821171"><br clear="none" class="ydpe431f14dyiv0019821171">As long as your subjects have eyes that move, muscles that work, sweat glands on their skin, and a beating heart, there are biological artifacts in your EEG that continuously contaminate your data by superpositioning with the brain signals. And as long as you have electricity in the house that you record in, especially if it is close to the electrodes, and no faraday cage, you have other noise and artifacts that continuously contaminate your signals. This is as true as water is wet and Evolution is a thing. <br clear="none" class="ydpe431f14dyiv0019821171"><br clear="none" class="ydpe431f14dyiv0019821171">And the best thing is: In times where this is not the case, the continuous activation of the IC components would be zero (as long as the decomposition is perfect - otherwise at least very low) and there would be no data alteration at those times if you subtract that IC! I feel like there should be a bit more in-depth understanding ofbiophysics of EEG, ICA, spatial filters as a method in general, linear mixture models, etc., before such a strong critique against this methods is put forward. Robert may or may not have this, but at least to me, from his arguments, it doesn't look like it. I don't say this as an attack but rather with the last spark of hope that spatial filtering as a useful method gets more understood and not bashed into oblivion in the clinics and diagnostics and in the end patients can benefit from the advances neuroscience has made in the last 30 years.<br clear="none" class="ydpe431f14dyiv0019821171"><br clear="none" class="ydpe431f14dyiv0019821171"></div><div class="ydpe431f14dyiv0019821171">Frankly, I am kinda sad about this discussion: It could have been really intersting and fruitful, and brought the investigation of ICA as a method forward (which would be good, since it's far from perfect and many things still need to be tested and learned!) but as someone stated earlier, it's mostly interesting from a sociological point of view...<br clear="none" class="ydpe431f14dyiv0019821171"></div><div class="ydpe431f14dyiv0019821171"> <br clear="none" class="ydpe431f14dyiv0019821171"></div><div class="ydpe431f14dyiv0019821171">I wish you all the best,<br clear="none" class="ydpe431f14dyiv0019821171"></div><div class="ydpe431f14dyiv0019821171">Marius<br clear="none" class="ydpe431f14dyiv0019821171"></div><div class="ydpe431f14dyiv0019821171"><br clear="none" class="ydpe431f14dyiv0019821171"></div><blockquote class="ydpe431f14dyiv0019821171gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex;">
<br clear="none" class="ydpe431f14dyiv0019821171">
<br clear="none" class="ydpe431f14dyiv0019821171">
Sincerely<br clear="none" class="ydpe431f14dyiv0019821171">
<span class="ydpe431f14dyiv0019821171gmail-m_-1304450643944854106HOEnZb"><font class="ydpe431f14dyiv0019821171" color="#888888"><br clear="none" class="ydpe431f14dyiv0019821171">
Georges Otte<br clear="none" class="ydpe431f14dyiv0019821171">
</font></span><span class="ydpe431f14dyiv0019821171gmail-m_-1304450643944854106im ydpe431f14dyiv0019821171gmail-m_-1304450643944854106HOEnZb"><br clear="none" class="ydpe431f14dyiv0019821171">
<br clear="none" class="ydpe431f14dyiv0019821171">
<br clear="none" class="ydpe431f14dyiv0019821171">
<br clear="none" class="ydpe431f14dyiv0019821171">
<br clear="none" class="ydpe431f14dyiv0019821171">
-----Original Message-----<br clear="none" class="ydpe431f14dyiv0019821171">
From: eeglablist [mailto:<a shape="rect" class="ydpe431f14dyiv0019821171" href="mailto:eeglablist-bounces@sccn.ucsd.edu" rel="nofollow" target="_blank">eeglablist-bounces@scc n.ucsd.edu</a>] On Behalf Of Stefan Debener<br clear="none" class="ydpe431f14dyiv0019821171">
Sent: Saturday, June 24, 2017 12:53 PM<br clear="none" class="ydpe431f14dyiv0019821171">
To: Robert Thatcher <<a shape="rect" class="ydpe431f14dyiv0019821171" href="mailto:rwthatcher2@yahoo.com" rel="nofollow" target="_blank">rwthatcher2@yahoo.com</a>>; Arnaud Delorme <<a shape="rect" class="ydpe431f14dyiv0019821171" href="mailto:arno@ucsd.edu" rel="nofollow" target="_blank">arno@ucsd.edu</a>><br clear="none" class="ydpe431f14dyiv0019821171">
Cc: eeglablist <<a shape="rect" class="ydpe431f14dyiv0019821171" href="mailto:eeglablist@sccn.ucsd.edu" rel="nofollow" target="_blank">eeglablist@sccn.ucsd.edu</a>>; <a shape="rect" class="ydpe431f14dyiv0019821171" href="mailto:pnunez@tulane.edu" rel="nofollow" target="_blank">pnunez@tulane.edu</a>; <a shape="rect" class="ydpe431f14dyiv0019821171" href="mailto:r.srinivasan@uci.edu" rel="nofollow" target="_blank">r.srinivasan@uci.edu</a>; Georges Otte <<a shape="rect" class="ydpe431f14dyiv0019821171" href="mailto:georges.otte@telenet.be" rel="nofollow" target="_blank">georges.otte@telenet.be</a>><br clear="none" class="ydpe431f14dyiv0019821171">
</span><span class="ydpe431f14dyiv0019821171gmail-m_-1304450643944854106im ydpe431f14dyiv0019821171gmail-m_-1304450643944854106HOEnZb">Subject: Re: [Eeglablist] ICA Misinformation<br clear="none" class="ydpe431f14dyiv0019821171">
<br clear="none" class="ydpe431f14dyiv0019821171">
</span><div class="ydpe431f14dyiv0019821171gmail-m_-1304450643944854106HOEnZb"><div class="ydpe431f14dyiv0019821171gmail-m_-1304450643944854106h5">Dear Robert,<br clear="none" class="ydpe431f14dyiv0019821171">
<br clear="none" class="ydpe431f14dyiv0019821171">
Ok, I guess I have to give up on you. Of course you ascribe the "original" times series some magic, and this is one of several flaws in your reasoning. To cite your own published paper:<br clear="none" class="ydpe431f14dyiv0019821171">
<br clear="none" class="ydpe431f14dyiv0019821171">
" PROBLEMS WITH RE-MONTAGING AND DISTORTIONS OF THE ORIGINAL TIME SERIES If the original EEG/event-related potential (ERP) time series is transformed into a second time series by using the average reference then the original phase differences from three electrode locations may be scrambled and lost. For example, with an average reference the entire surface of the brain is not measured, thus the averaging does not create a true zero potential at each instant of time." (citation taken from <a shape="rect" class="ydpe431f14dyiv0019821171" href="http://www.appliedneuroscience.com/Coh_phasediff&phase_resetinEEG-ERP.pdf" rel="nofollow" target="_blank">http://www.appliedneuroscience .com/Coh_phasediff&phase_ resetinEEG-ERP.pdf</a>).<br clear="none" class="ydpe431f14dyiv0019821171">
<br clear="none" class="ydpe431f14dyiv0019821171">
Your phrase "distortion of the original time series" does make sense only if you believe that the original time series represents somehow a gold standard, something special (something "magic", forgive my poor use of English) that is magically close to the contributions from brain generators to the surface-recorded signal. In my view this is not justified, for instance because it disregards the fact that recording settings will influence how the data will be recorded, that is, they determine phase and amplitude! I argued before that attributing something special to the "original" time series is highly misleading, the "original" time series is not closer to the brain signal, in contrast it may be pretty far away from it, not only because phase of mixed brain generators does not make much sense, but also becasue of all the artifactual influences not accounted for. Now if you really believe that one particular reference scheme is appropriate in bringing the original data close to the brain while all others are rubbish, then I wish you good luck in trying to convince the community. I predict that nobody will take you serious. The reference discussion has been going on for ages, and there are good reasons to change the reference (online or<br clear="none" class="ydpe431f14dyiv0019821171">
offline) depending on what the purpose of the study/analysis is.<br clear="none" class="ydpe431f14dyiv0019821171">
<br clear="none" class="ydpe431f14dyiv0019821171">
I argue that any post-recording signal processing that changes the morphology of a time series will change the phase values as well, NOT just the average reference, and NOT just ICA. In your case, the average reference is far from zero because 19 10-20 channels, which are located on the top half of the head sphere, can hardly sum up to zero. Only a full equidistant spatial coverage of the head sphere would make the spatial average approximate zero potential. Check this paper for in-depth discussion on the average reference and the bias it introduces:<br clear="none" class="ydpe431f14dyiv0019821171">
<a shape="rect" class="ydpe431f14dyiv0019821171" href="https://www.ncbi.nlm.nih.gov/pubmed/10402104" rel="nofollow" target="_blank">https://www.ncbi.nlm.nih.gov/p ubmed/10402104</a>). If you don't like the average reference, fair enough, then simply take any other reference, but you will observe the same "distortion" on phase values. By the way, non-invasive EEG does not measure signals from the brain, there is a skull in between and a couple other layers, all with different conductivies...just another reason to be aware of the inverse problem and trat signals originally recorded from the scalp with great care.<br clear="none" class="ydpe431f14dyiv0019821171">
<br clear="none" class="ydpe431f14dyiv0019821171">
The second flaw in your reasoning is that you believe there are "artifact-free" intervals. The report I included showed (far from perfect, but clearly evident!) heart-electrical activity. Recording parameters strongly determine whether such "EKG" contributions show up in ICA decompositions or not (such as sub-10/20 spatial coverage). Now, only because you don't see heart-electrical activity in your "original"<br clear="none" class="ydpe431f14dyiv0019821171">
19-ch 10/20 recordings, can you seriously claim that the heart of your participant was not beating during recording? Of course not! All it says is that the influcene may be stronger or weaker represented in your recordings (depending on individual differences, and, again, recording parameters). With your philosophy, you gonna miss it if it does not jump into your eyes.<br clear="none" class="ydpe431f14dyiv0019821171">
<br clear="none" class="ydpe431f14dyiv0019821171">
Anyway, I hope others join the discussion, I am giving up.<br clear="none" class="ydpe431f14dyiv0019821171">
<br clear="none" class="ydpe431f14dyiv0019821171">
Best,<br clear="none" class="ydpe431f14dyiv0019821171">
Stefan<br clear="none" class="ydpe431f14dyiv0019821171">
<br clear="none" class="ydpe431f14dyiv0019821171">
<br clear="none" class="ydpe431f14dyiv0019821171">
Am 23.06.17 um 19:16 schrieb Robert Thatcher:<br clear="none" class="ydpe431f14dyiv0019821171">
><br clear="none" class="ydpe431f14dyiv0019821171">
> Dear Stefan,<br clear="none" class="ydpe431f14dyiv0019821171">
><br clear="none" class="ydpe431f14dyiv0019821171">
> Thank you for your dilligence and dedication to this important<br clear="none" class="ydpe431f14dyiv0019821171">
> issue. I am pleased that you are in agreement with myself and<br clear="none" class="ydpe431f14dyiv0019821171">
> Georges and Gert and scienific publications that "a spatial filter<br clear="none" class="ydpe431f14dyiv0019821171">
> operation such as ICA or other, the phase differences may indeed be<br clear="none" class="ydpe431f14dyiv0019821171">
> different." I would add that we have not yet found an instance where<br clear="none" class="ydpe431f14dyiv0019821171">
> ICA reconstruction did not alter phase differences between channel and<br clear="none" class="ydpe431f14dyiv0019821171">
> there I would change the word "may" and phrase to "are indeed different".<br clear="none" class="ydpe431f14dyiv0019821171">
><br clear="none" class="ydpe431f14dyiv0019821171">
> You sated: "As a toy example I include the common average reference. "<br clear="none" class="ydpe431f14dyiv0019821171">
><br clear="none" class="ydpe431f14dyiv0019821171">
> In Neuroguide we do not allow one to compute phase diffferences or<br clear="none" class="ydpe431f14dyiv0019821171">
> coherence when using a common reference. We found with simulation and<br clear="none" class="ydpe431f14dyiv0019821171">
> real data that the common reference mixes the phases differences<br clear="none" class="ydpe431f14dyiv0019821171">
> between all channels and itself distorts phase differences and often<br clear="none" class="ydpe431f14dyiv0019821171">
> in strange ways, for example, if one or two channels has a suddent<br clear="none" class="ydpe431f14dyiv0019821171">
> large amplitude alpha or theta rhythm then the phase differences<br clear="none" class="ydpe431f14dyiv0019821171">
> between channels that do not have a alpha or beta or theta rhythm are<br clear="none" class="ydpe431f14dyiv0019821171">
> altered. If one uses a single common reference then if an alpha<br clear="none" class="ydpe431f14dyiv0019821171">
> rhythm appeas for example in O1/2 then there is not change in phase<br clear="none" class="ydpe431f14dyiv0019821171">
> differences in channels where there is no alpha. We also compared two<br clear="none" class="ydpe431f14dyiv0019821171">
> sine waves at different phase differences, e.g., 30 deg, 60 deg, 90<br clear="none" class="ydpe431f14dyiv0019821171">
> deg, etc and mixed different amounts of white noise in one of the<br clear="none" class="ydpe431f14dyiv0019821171">
> channels we found linear reductions in coherence (i.e., the phase<br clear="none" class="ydpe431f14dyiv0019821171">
> stability over time) as a function of the SNR and the mean phase<br clear="none" class="ydpe431f14dyiv0019821171">
> differences were stable until the noise was too high and when<br clear="none" class="ydpe431f14dyiv0019821171">
> coherencer was near zero. In contrast, uses the averge refence and<br clear="none" class="ydpe431f14dyiv0019821171">
> repeats the same signal and noise test with different phase<br clear="none" class="ydpe431f14dyiv0019821171">
> differences then the mean phase difference is quickly lost and there<br clear="none" class="ydpe431f14dyiv0019821171">
> is no valid measure of coherence. Here is a url to a publication<br clear="none" class="ydpe431f14dyiv0019821171">
> that discusses this topic:<br clear="none" class="ydpe431f14dyiv0019821171">
> <a shape="rect" class="ydpe431f14dyiv0019821171" href="http://www.appliedneuroscience.com/Coh_phasediff&phase_resetinEEG-ERP" rel="nofollow" target="_blank">http://www.appliedneuroscience .com/Coh_phasediff&phase_ resetinEEG-ERP</a>.<br clear="none" class="ydpe431f14dyiv0019821171">
> pdf<br clear="none" class="ydpe431f14dyiv0019821171">
><br clear="none" class="ydpe431f14dyiv0019821171">
> You stated: "Your claim that ICA has somehow corrupted the data such<br clear="none" class="ydpe431f14dyiv0019821171">
> that previously super reliable clinical effects all over a sudden<br clear="none" class="ydpe431f14dyiv0019821171">
> vanished is not convincing either."<br clear="none" class="ydpe431f14dyiv0019821171">
><br clear="none" class="ydpe431f14dyiv0019821171">
> I never said this. There is no ICA reconstruction in NeuroGuide and<br clear="none" class="ydpe431f14dyiv0019821171">
> the over 3,000 users of Neuroguide have not complained about coherence<br clear="none" class="ydpe431f14dyiv0019821171">
> or phase and they always obtain repeatable measures when the retest<br clear="none" class="ydpe431f14dyiv0019821171">
> patients. It was only the WinEEG users that are worried about ICA<br clear="none" class="ydpe431f14dyiv0019821171">
> and they never said that they got "super reliable clinical effects".<br clear="none" class="ydpe431f14dyiv0019821171">
> They just noticed that coherence using the WinEEG after ICA<br clear="none" class="ydpe431f14dyiv0019821171">
> reconstruction was totally different than when they use NxLink or SKIL<br clear="none" class="ydpe431f14dyiv0019821171">
> or Neurorep or Neuroguide or Brindx or other software, etc.<br clear="none" class="ydpe431f14dyiv0019821171">
><br clear="none" class="ydpe431f14dyiv0019821171">
> You stated: "Now which phase values are valid, those obtained by one<br clear="none" class="ydpe431f14dyiv0019821171">
> particular reference scheme or those by another? In my view they are<br clear="none" class="ydpe431f14dyiv0019821171">
> both arbitraty"<br clear="none" class="ydpe431f14dyiv0019821171">
><br clear="none" class="ydpe431f14dyiv0019821171">
> See my earlier reply regarding average references and the Laplacian<br clear="none" class="ydpe431f14dyiv0019821171">
> reference in regard to the accurate and reproducible measures of phase<br clear="none" class="ydpe431f14dyiv0019821171">
> differences as opposed to using a single common reference. I am<br clear="none" class="ydpe431f14dyiv0019821171">
> recopying the link here:<br clear="none" class="ydpe431f14dyiv0019821171">
> <a shape="rect" class="ydpe431f14dyiv0019821171" href="http://www.appliedneuroscience.com/Coh_phasediff&phase_resetinEEG-ERP" rel="nofollow" target="_blank">http://www.appliedneuroscience .com/Coh_phasediff&phase_ resetinEEG-ERP</a>.<br clear="none" class="ydpe431f14dyiv0019821171">
> pdf here is another review on this topic:<br clear="none" class="ydpe431f14dyiv0019821171">
> <a shape="rect" class="ydpe431f14dyiv0019821171" href="http://www.appliedneuroscience.com/Brain%20Connectivity-A%20Tutorial.p" rel="nofollow" target="_blank">http://www.appliedneuroscience .com/Brain%20Connectivity-A% 20Tutorial.p</a><br clear="none" class="ydpe431f14dyiv0019821171">
> df<br clear="none" class="ydpe431f14dyiv0019821171">
><br clear="none" class="ydpe431f14dyiv0019821171">
> You stated: “there is no such magically clean raw brainsignal<br clear="none" class="ydpe431f14dyiv0019821171">
> available in the first place!”<br clear="none" class="ydpe431f14dyiv0019821171">
><br clear="none" class="ydpe431f14dyiv0019821171">
> No one says that “magic” is involved in the EEG. However, the physics<br clear="none" class="ydpe431f14dyiv0019821171">
> of EEG is involved and during a 5 minute to 20 minute EEG recording<br clear="none" class="ydpe431f14dyiv0019821171">
> there is plenty of artifact free data. There are over 100,000 QEEG<br clear="none" class="ydpe431f14dyiv0019821171">
> publications in the National Library of Medicine with high effect<br clear="none" class="ydpe431f14dyiv0019821171">
> sizes and high test retest reliability and highly reproducible<br clear="none" class="ydpe431f14dyiv0019821171">
> findings. My concern and the concern of many others is that ICA<br clear="none" class="ydpe431f14dyiv0019821171">
> reconstruction alters phase differences in an entire EEG recording<br clear="none" class="ydpe431f14dyiv0019821171">
> even if there are only a few instances of artifact. The reason that<br clear="none" class="ydpe431f14dyiv0019821171">
> QEEG has been so successful and so widely used since the late 1950s is<br clear="none" class="ydpe431f14dyiv0019821171">
> because people have been successful in avoiding or deleting artifact<br clear="none" class="ydpe431f14dyiv0019821171">
> and selecting multiple artifact free parts of the record and achieving<br clear="none" class="ydpe431f14dyiv0019821171">
> 0.95 or higher test retest reliability. Over reaction to the<br clear="none" class="ydpe431f14dyiv0019821171">
> presence of small amounts of artifact is not a justification for<br clear="none" class="ydpe431f14dyiv0019821171">
> altering the electricity of the brain including network dynamics such<br clear="none" class="ydpe431f14dyiv0019821171">
> as average synaptic rise times and conduction velocities and couplings<br clear="none" class="ydpe431f14dyiv0019821171">
> between groups of neurons.<br clear="none" class="ydpe431f14dyiv0019821171">
><br clear="none" class="ydpe431f14dyiv0019821171">
> You stated: “Artifacts not accounted for adulterate EEG phase values”<br clear="none" class="ydpe431f14dyiv0019821171">
><br clear="none" class="ydpe431f14dyiv0019821171">
> We agree on this but this is not what the discussion is about. The<br clear="none" class="ydpe431f14dyiv0019821171">
> concerns of many is that the phase differences in the artifact free<br clear="none" class="ydpe431f14dyiv0019821171">
> sections are altered. Rarely if ever do clinicians/scientists bother<br clear="none" class="ydpe431f14dyiv0019821171">
> computing the phase differences during an eye movement artifact.<br clear="none" class="ydpe431f14dyiv0019821171">
> Usually phase differences are zero and this physics fact plus the<br clear="none" class="ydpe431f14dyiv0019821171">
> electrical gradients from the eyes is how many people detect eye<br clear="none" class="ydpe431f14dyiv0019821171">
> movement artifact and then omit the artifact from analyses without<br clear="none" class="ydpe431f14dyiv0019821171">
> using ICA. The main focus in QEEG is the parts of the recording<br clear="none" class="ydpe431f14dyiv0019821171">
> where there is no artifact and the electrical potentials are generated<br clear="none" class="ydpe431f14dyiv0019821171">
> by the brain inside the skull.<br clear="none" class="ydpe431f14dyiv0019821171">
><br clear="none" class="ydpe431f14dyiv0019821171">
> Also, thank you for your use of the Hilbert transform, this is only of<br clear="none" class="ydpe431f14dyiv0019821171">
> the tools that we use and it allows one to evaluate phase differences<br clear="none" class="ydpe431f14dyiv0019821171">
> in every individual time sample in the artifact free sections and<br clear="none" class="ydpe431f14dyiv0019821171">
> prove that each and every one of the phase differences for all<br clear="none" class="ydpe431f14dyiv0019821171">
> channel combinations is altered by ICA reconstruction. You may be<br clear="none" class="ydpe431f14dyiv0019821171">
> interested that we use the Hilbert transform to measure phase shift<br clear="none" class="ydpe431f14dyiv0019821171">
> and phase lock duration that have a high correlation with Autism and<br clear="none" class="ydpe431f14dyiv0019821171">
> intelligence in short distance connections and with use the Hilbert<br clear="none" class="ydpe431f14dyiv0019821171">
> transform to measure the magnitude of information flow (phase slope<br clear="none" class="ydpe431f14dyiv0019821171">
> index) and intelligence. Here are some hyperlinks to these studies:<br clear="none" class="ydpe431f14dyiv0019821171">
> <a shape="rect" class="ydpe431f14dyiv0019821171" href="http://www.appliedneuroscience.com/Intelligence-phase_reset_Nature.pdf" rel="nofollow" target="_blank">http://www.appliedneuroscience .com/Intelligence-phase_reset_ Nature.pdf</a><br clear="none" class="ydpe431f14dyiv0019821171">
><br clear="none" class="ydpe431f14dyiv0019821171">
> <a shape="rect" class="ydpe431f14dyiv0019821171" href="http://www.appliedneuroscience.com/Autism%20Thatcher%20et%20al.pdf" rel="nofollow" target="_blank">http://www.appliedneuroscience .com/Autism%20Thatcher%20et% 20al.pdf</a><br clear="none" class="ydpe431f14dyiv0019821171">
><br clear="none" class="ydpe431f14dyiv0019821171">
> <a shape="rect" class="ydpe431f14dyiv0019821171" href="http://www.appliedneuroscience.com/Default_Network_LORETA_Phase_Reset-" rel="nofollow" target="_blank">http://www.appliedneuroscience .com/Default_Network_LORETA_ Phase_Reset-</a><br clear="none" class="ydpe431f14dyiv0019821171">
> Thatcher_et_al.pdf<br clear="none" class="ydpe431f14dyiv0019821171">
><br clear="none" class="ydpe431f14dyiv0019821171">
> <a shape="rect" class="ydpe431f14dyiv0019821171" href="http://www.appliedneuroscience.com/Intelligence%20&%20information%20fl" rel="nofollow" target="_blank">http://www.appliedneuroscience .com/Intelligence%20&%20inform ation%20fl</a><br clear="none" class="ydpe431f14dyiv0019821171">
> ow-Thatcher%20et%20al%202016.p df<br clear="none" class="ydpe431f14dyiv0019821171">
><br clear="none" class="ydpe431f14dyiv0019821171">
> We are also using the Hilbert transforms for cross-frequency network<br clear="none" class="ydpe431f14dyiv0019821171">
> dynamics including phase-amplitude coupling. I do not believe that<br clear="none" class="ydpe431f14dyiv0019821171">
> we would have discovered these important network correlations if we<br clear="none" class="ydpe431f14dyiv0019821171">
> had used ICA reconstruction.<br clear="none" class="ydpe431f14dyiv0019821171">
><br clear="none" class="ydpe431f14dyiv0019821171">
> Robert<br clear="none" class="ydpe431f14dyiv0019821171">
><br clear="none" class="ydpe431f14dyiv0019821171">
><br clear="none" class="ydpe431f14dyiv0019821171">
> On Friday, June 23, 2017, 9:53:01 AM EDT, Stefan Debener<br clear="none" class="ydpe431f14dyiv0019821171">
> <<a shape="rect" class="ydpe431f14dyiv0019821171" href="mailto:stefan.debener@uni-oldenburg.de" rel="nofollow" target="_blank">stefan.debener@uni-oldenburg. de</a>> wrote:<br clear="none" class="ydpe431f14dyiv0019821171">
><br clear="none" class="ydpe431f14dyiv0019821171">
><br clear="none" class="ydpe431f14dyiv0019821171">
> Dear Robert,<br clear="none" class="ydpe431f14dyiv0019821171">
><br clear="none" class="ydpe431f14dyiv0019821171">
> I have expanded my illustration and now consider the phase differences<br clear="none" class="ydpe431f14dyiv0019821171">
> between two channels, slides 13 to 16 of the updated pdf:<br clear="none" class="ydpe431f14dyiv0019821171">
> <a shape="rect" class="ydpe431f14dyiv0019821171" href="https://www.dropbox.com/s/e70qhf91dgc5anu/Thatcher_summary_2.pdf?dl=0" rel="nofollow" target="_blank">https://www.dropbox.com/s/e70q hf91dgc5anu/Thatcher_summary_ 2.pdf?dl=0</a><br clear="none" class="ydpe431f14dyiv0019821171">
><br clear="none" class="ydpe431f14dyiv0019821171">
> Note that phase values were derived by the Hilbert transform of the<br clear="none" class="ydpe431f14dyiv0019821171">
> bandpass filtered signal, as explained by W Freeman here:<br clear="none" class="ydpe431f14dyiv0019821171">
> <a shape="rect" class="ydpe431f14dyiv0019821171" href="http://www.scholarpedia.org/article/Hilbert_transform_for_brain_waves" rel="nofollow" target="_blank">http://www.scholarpedia.org/ar ticle/Hilbert_transform_for_br ain_waves</a><br clear="none" class="ydpe431f14dyiv0019821171">
><br clear="none" class="ydpe431f14dyiv0019821171">
> More details on the particular implementation I used are here:<br clear="none" class="ydpe431f14dyiv0019821171">
> <a shape="rect" class="ydpe431f14dyiv0019821171" href="https://de.mathworks.com/help/signal/ref/hilbert.html" rel="nofollow" target="_blank">https://de.mathworks.com/help/ signal/ref/hilbert.html</a><br clear="none" class="ydpe431f14dyiv0019821171">
><br clear="none" class="ydpe431f14dyiv0019821171">
> If you measure phase differences between two channels, consider the<br clear="none" class="ydpe431f14dyiv0019821171">
> result as your gold standard, and then apply a spatial filter<br clear="none" class="ydpe431f14dyiv0019821171">
> operation such as ICA or other, the phase differences may indeed be<br clear="none" class="ydpe431f14dyiv0019821171">
> different. I assume any spatial filter (that effectively spatially<br clear="none" class="ydpe431f14dyiv0019821171">
> filters the data) changes phase values and phase difference values. As<br clear="none" class="ydpe431f14dyiv0019821171">
> a toy example I include the common average reference. If you apply a<br clear="none" class="ydpe431f14dyiv0019821171">
> common average reference to the raw data, then bandpass filter as<br clear="none" class="ydpe431f14dyiv0019821171">
> before, and compare the phase difference values to your "gold<br clear="none" class="ydpe431f14dyiv0019821171">
> standard", then the phase differences will change as well. Now which<br clear="none" class="ydpe431f14dyiv0019821171">
> phase values are valid, those obtained by one particular reference<br clear="none" class="ydpe431f14dyiv0019821171">
> scheme or those by another? In my view they are both arbitraty, since<br clear="none" class="ydpe431f14dyiv0019821171">
> recording settings as well as preprocessing steps may have a strong<br clear="none" class="ydpe431f14dyiv0019821171">
> impact on the actually measured phase. There is no reason to assume<br clear="none" class="ydpe431f14dyiv0019821171">
> that a change in phase, or in phase differences, "adulterates" a<br clear="none" class="ydpe431f14dyiv0019821171">
> magically clean phase signal obtained from the raw data - simply<br clear="none" class="ydpe431f14dyiv0019821171">
> because there is no such magically clean raw brain signal available in the first place!<br clear="none" class="ydpe431f14dyiv0019821171">
><br clear="none" class="ydpe431f14dyiv0019821171">
> Your claim that ICA has somehow corrupted the data such that<br clear="none" class="ydpe431f14dyiv0019821171">
> previously super reliable clinical effects all over a sudden vanished<br clear="none" class="ydpe431f14dyiv0019821171">
> is not convincing either. Artifacts such as eye blinks and lateral eye<br clear="none" class="ydpe431f14dyiv0019821171">
> movements are very common, I hope you can agree at least here. Now,<br clear="none" class="ydpe431f14dyiv0019821171">
> keep in mind that they contribute fixed spatial patterns - as long as<br clear="none" class="ydpe431f14dyiv0019821171">
> the electrodes cap does not shift during acquisition the projections<br clear="none" class="ydpe431f14dyiv0019821171">
> of the sources of those artifacts do not change. My illustrations<br clear="none" class="ydpe431f14dyiv0019821171">
> above show very clearly how artifacs indeed adulterate phase values,<br clear="none" class="ydpe431f14dyiv0019821171">
> just as Arnos illustrations do! Now, if you disregard artifactual<br clear="none" class="ydpe431f14dyiv0019821171">
> influences you may end up with highly reliable connectivity effects -<br clear="none" class="ydpe431f14dyiv0019821171">
> but they tell you very little about brain function! Even more<br clear="none" class="ydpe431f14dyiv0019821171">
> troubling, if you compare two individuals EEGs (say, one "healthy",<br clear="none" class="ydpe431f14dyiv0019821171">
> one "abnormal"), then a different amount of artifacts in the data, if<br clear="none" class="ydpe431f14dyiv0019821171">
> not carefully taken care of during preprocessing, will produce<br clear="none" class="ydpe431f14dyiv0019821171">
> spurious results that are falsely attributed to differences in brain<br clear="none" class="ydpe431f14dyiv0019821171">
> function. Actually, given that many artifacts often contribute much<br clear="none" class="ydpe431f14dyiv0019821171">
> more variance to that raw signals than (reasonably well validated)<br clear="none" class="ydpe431f14dyiv0019821171">
> brain signals, such as fronto-midline theta, this is actually very likely! So, what we learn is that:<br clear="none" class="ydpe431f14dyiv0019821171">
><br clear="none" class="ydpe431f14dyiv0019821171">
> Artifacts not accounted for adulterate EEG phase values<br clear="none" class="ydpe431f14dyiv0019821171">
><br clear="none" class="ydpe431f14dyiv0019821171">
> Best,<br clear="none" class="ydpe431f14dyiv0019821171">
><br clear="none" class="ydpe431f14dyiv0019821171">
> Stefan<br clear="none" class="ydpe431f14dyiv0019821171">
><br clear="none" class="ydpe431f14dyiv0019821171">
><br clear="none" class="ydpe431f14dyiv0019821171">
><br clear="none" class="ydpe431f14dyiv0019821171">
> Am 22.06.17 um 20:30 schrieb Robert Thatcher:<br clear="none" class="ydpe431f14dyiv0019821171">
> > Dear Stefan,<br clear="none" class="ydpe431f14dyiv0019821171">
> > The attachment did not contain any measures of phase differences<br clear="none" class="ydpe431f14dyiv0019821171">
> > between channels. It is very difficult to visually see differences<br clear="none" class="ydpe431f14dyiv0019821171">
> > in phase differences. One must use the cross-spectrum to calculate<br clear="none" class="ydpe431f14dyiv0019821171">
> > phase differences and compare phase differences in degrees. Phase<br clear="none" class="ydpe431f14dyiv0019821171">
> > difference varies from -180 to 180 degrees and one must look at the<br clear="none" class="ydpe431f14dyiv0019821171">
> > numbers. Below is a url to the two power points that also show<br clear="none" class="ydpe431f14dyiv0019821171">
> > visually similar EEG tracings but also computed the instantaneous<br clear="none" class="ydpe431f14dyiv0019821171">
> > phase differences using the Hilbert transform (complex demodulation).<br clear="none" class="ydpe431f14dyiv0019821171">
> > Four identical time points were selected and they demonstrated<br clear="none" class="ydpe431f14dyiv0019821171">
> > totally different phase differences with respect to the O1 channel<br clear="none" class="ydpe431f14dyiv0019821171">
> > and the other 18 channels. No matter what reference channel one<br clear="none" class="ydpe431f14dyiv0019821171">
> > selects and no matter what identical time points one selects there<br clear="none" class="ydpe431f14dyiv0019821171">
> > are always large differences in the phase difference between<br clear="none" class="ydpe431f14dyiv0019821171">
> > channels in all frequency bands. I also computed the average phase<br clear="none" class="ydpe431f14dyiv0019821171">
> > difference in the artifact free parts of the record and the averages<br clear="none" class="ydpe431f14dyiv0019821171">
> > were statistically significantly different at P < 0.0001 and the same for the FFT.<br clear="none" class="ydpe431f14dyiv0019821171">
> ><br clear="none" class="ydpe431f14dyiv0019821171">
> > Proof of phase difference adulteration is in the power points. I am<br clear="none" class="ydpe431f14dyiv0019821171">
> > again copying the hyperlink here:<br clear="none" class="ydpe431f14dyiv0019821171">
> ><br clear="none" class="ydpe431f14dyiv0019821171">
> ><br clear="none" class="ydpe431f14dyiv0019821171">
> ><br clear="none" class="ydpe431f14dyiv0019821171">
> <a shape="rect" class="ydpe431f14dyiv0019821171" href="http://www.appliedneuroscience.com/Phase_Diff-Original_&_Delorme-Post-" rel="nofollow" target="_blank">http://www.appliedneuroscience .com/Phase_Diff-Original_&_ Delorme-Post-</a><br clear="none" class="ydpe431f14dyiv0019821171">
> ICA-4_time_points.zip<br clear="none" class="ydpe431f14dyiv0019821171">
> ><br clear="none" class="ydpe431f14dyiv0019821171">
> ><br clear="none" class="ydpe431f14dyiv0019821171">
> > This cannot be explained by a low quality ICA reconstruction because<br clear="none" class="ydpe431f14dyiv0019821171">
> > the ICA reconstruction was conducted by Arnu using EEGLab software.<br clear="none" class="ydpe431f14dyiv0019821171">
> ><br clear="none" class="ydpe431f14dyiv0019821171">
> > Robert<br clear="none" class="ydpe431f14dyiv0019821171">
> > On Thursday, June 22, 2017, 2:00:19 PM EDT, Stefan Debener<br clear="none" class="ydpe431f14dyiv0019821171">
> > <<a shape="rect" class="ydpe431f14dyiv0019821171" href="mailto:stefan.debener@uni-oldenburg.de" rel="nofollow" target="_blank">stefan.debener@uni-oldenburg. de</a><br clear="none" class="ydpe431f14dyiv0019821171">
> <mailto:<a shape="rect" class="ydpe431f14dyiv0019821171" href="mailto:stefan.debener@uni-oldenburg.de" rel="nofollow" target="_blank">stefan.debener@uni-old enburg.de</a>>> wrote:<br clear="none" class="ydpe431f14dyiv0019821171">
> ><br clear="none" class="ydpe431f14dyiv0019821171">
> ><br clear="none" class="ydpe431f14dyiv0019821171">
> > Dear Robert,<br clear="none" class="ydpe431f14dyiv0019821171">
> ><br clear="none" class="ydpe431f14dyiv0019821171">
> > I looked up some own data and find absolutely no evidence in favour<br clear="none" class="ydpe431f14dyiv0019821171">
> > of your ICA phase adulteration claim, see the attached pdf report. I<br clear="none" class="ydpe431f14dyiv0019821171">
> > guess you simply used a poor ICA implementation, and/or a poor<br clear="none" class="ydpe431f14dyiv0019821171">
> > component selection. The attached example is in full accordance with<br clear="none" class="ydpe431f14dyiv0019821171">
> > Arnos reply, with the difference that I zoom into a clearly visibile<br clear="none" class="ydpe431f14dyiv0019821171">
> > alpha oscillation, to have a reference brain signal. The example<br clear="none" class="ydpe431f14dyiv0019821171">
> > shows no evidence that occipital alpha phase is biased by ICA eye<br clear="none" class="ydpe431f14dyiv0019821171">
> > blink correction. This is a very typical example and based on a<br clear="none" class="ydpe431f14dyiv0019821171">
> > quick and dirty ICA decomposition, nothing fancy, to keep this demo<br clear="none" class="ydpe431f14dyiv0019821171">
> > simple. Better preprocessing and component selection would easily<br clear="none" class="ydpe431f14dyiv0019821171">
> > further improve the signal quality.<br clear="none" class="ydpe431f14dyiv0019821171">
> ><br clear="none" class="ydpe431f14dyiv0019821171">
> > Best,<br clear="none" class="ydpe431f14dyiv0019821171">
> ><br clear="none" class="ydpe431f14dyiv0019821171">
> > Stefan<br clear="none" class="ydpe431f14dyiv0019821171">
> ><br clear="none" class="ydpe431f14dyiv0019821171">
> ><br clear="none" class="ydpe431f14dyiv0019821171">
> ><br clear="none" class="ydpe431f14dyiv0019821171">
> > Am 20.06.17 um 19:53 schrieb Robert Thatcher:<br clear="none" class="ydpe431f14dyiv0019821171">
> > ><br clear="none" class="ydpe431f14dyiv0019821171">
> > > Dear Arno,<br clear="none" class="ydpe431f14dyiv0019821171">
> > ><br clear="none" class="ydpe431f14dyiv0019821171">
> > > 1)*On Phase Differences in the Original vs the Delorme ICA<br clear="none" class="ydpe431f14dyiv0019821171">
> > > Reconstruction: *We can agree or disagree about whether or not<br clear="none" class="ydpe431f14dyiv0019821171">
> > > some small eye movement artifact was in the hand selection that I<br clear="none" class="ydpe431f14dyiv0019821171">
> > > did. But that misses the main point here. That is the ICA<br clear="none" class="ydpe431f14dyiv0019821171">
> > > reconstruction alters each and every data point in the entire<br clear="none" class="ydpe431f14dyiv0019821171">
> > > record including all artifact free portions no matter what one<br clear="none" class="ydpe431f14dyiv0019821171">
> > > selects. For example, the record is 6 minutes and 51 seconds = 411<br clear="none" class="ydpe431f14dyiv0019821171">
> > > seconds. The Mitsar sample rate was 250 samples per second =<br clear="none" class="ydpe431f14dyiv0019821171">
> > > 102,750 data samples. Phase difference for each frequency band for<br clear="none" class="ydpe431f14dyiv0019821171">
> > > each and every one of the<br clear="none" class="ydpe431f14dyiv0019821171">
> > > 102,750 data samples has been altered by your own ICA<br clear="none" class="ydpe431f14dyiv0019821171">
> > > reconstruction in the EDF file that you emailed to me. Unless you<br clear="none" class="ydpe431f14dyiv0019821171">
> > > were to sit next to me or if we do a Team Viewer it is not<br clear="none" class="ydpe431f14dyiv0019821171">
> > > possible for me to demonstrate this for all of the data points and<br clear="none" class="ydpe431f14dyiv0019821171">
> > > then create a power point for all of these data samples. However,<br clear="none" class="ydpe431f14dyiv0019821171">
> > > I can show some exemplars, for example, I have created two figures<br clear="none" class="ydpe431f14dyiv0019821171">
> > > at 4 different time points (1 sec;<br clear="none" class="ydpe431f14dyiv0019821171">
> > > 2:27 sec; 42 sec & 5:49 sec) that you can download. You can<br clear="none" class="ydpe431f14dyiv0019821171">
> > > extract each screen capture and expand them so that you can see<br clear="none" class="ydpe431f14dyiv0019821171">
> > > that the exact same time points were selected and the Hilbert<br clear="none" class="ydpe431f14dyiv0019821171">
> > > transform JTFA for the<br clear="none" class="ydpe431f14dyiv0019821171">
> > > 4 time points resulted in different phase differences in all<br clear="none" class="ydpe431f14dyiv0019821171">
> > > channel combinations with respect to O1 for all frequencies. The<br clear="none" class="ydpe431f14dyiv0019821171">
> > > same is true no matter which channel is selected to compute the<br clear="none" class="ydpe431f14dyiv0019821171">
> > > phase differences in degrees. The same is true also if one<br clear="none" class="ydpe431f14dyiv0019821171">
> > > computes averages of the instantaneous phase differences or if one<br clear="none" class="ydpe431f14dyiv0019821171">
> > > uses the FFT. Here is the download URL:<br clear="none" class="ydpe431f14dyiv0019821171">
> > ><br clear="none" class="ydpe431f14dyiv0019821171">
> > ><br clear="none" class="ydpe431f14dyiv0019821171">
> ><br clear="none" class="ydpe431f14dyiv0019821171">
> <a shape="rect" class="ydpe431f14dyiv0019821171" href="http://www.appliedneuroscience.com/Phase_Diff-Original_&_Delorme-Post-" rel="nofollow" target="_blank">http://www.appliedneuroscience .com/Phase_Diff-Original_&_ Delorme-Post-</a><br clear="none" class="ydpe431f14dyiv0019821171">
> ICA-4_time_points.zip<br clear="none" class="ydpe431f14dyiv0019821171">
> > ><br clear="none" class="ydpe431f14dyiv0019821171">
> > ><br clear="none" class="ydpe431f14dyiv0019821171">
> > > 2)*On the WinEEG ICA Reconstruction: *I agree that having access<br clear="none" class="ydpe431f14dyiv0019821171">
> > > to ICA components themselves and the topography is critical in<br clear="none" class="ydpe431f14dyiv0019821171">
> > > understanding exactly what the WinEEG software did. Unfortunately,<br clear="none" class="ydpe431f14dyiv0019821171">
> > > I personally do not have access to the WinEEG software.<br clear="none" class="ydpe431f14dyiv0019821171">
> > > Clinician/Scientists in Australia use the WinEEG software and they<br clear="none" class="ydpe431f14dyiv0019821171">
> > > were the ones that expressed concern about phase difference<br clear="none" class="ydpe431f14dyiv0019821171">
> > > distortion at a workshop in Adelaide and gave me the original and<br clear="none" class="ydpe431f14dyiv0019821171">
> > > the WinEEG ICA eye movement corrected files in EDF format. They<br clear="none" class="ydpe431f14dyiv0019821171">
> > > explained that they removed only one ICA component for eye<br clear="none" class="ydpe431f14dyiv0019821171">
> > > movement before they reconstructed a new time series. At first, I<br clear="none" class="ydpe431f14dyiv0019821171">
> > > was impressed because the eye movements were absent in the<br clear="none" class="ydpe431f14dyiv0019821171">
> > > reconstructed time series. I then was able to use JTFA (Hilbert<br clear="none" class="ydpe431f14dyiv0019821171">
> > > transform) to compare the two edf files and discovered that all of<br clear="none" class="ydpe431f14dyiv0019821171">
> > > the phase differences for all channels for all frequencies had<br clear="none" class="ydpe431f14dyiv0019821171">
> > > been altered by the ICA reconstruction including artifact free<br clear="none" class="ydpe431f14dyiv0019821171">
> > > periods. I could demonstrate this by individual time comparisons<br clear="none" class="ydpe431f14dyiv0019821171">
> > > or averages of instantaneous phase differences or by the FFT. A<br clear="none" class="ydpe431f14dyiv0019821171">
> > > user of WinEEG explained that they do not throw away the original<br clear="none" class="ydpe431f14dyiv0019821171">
> > > raw digital data, however I was told that they believe that the<br clear="none" class="ydpe431f14dyiv0019821171">
> > > ICA reconstructed times series is artifact free and therefore they<br clear="none" class="ydpe431f14dyiv0019821171">
> > > compute means and standard deviations for their normative database<br clear="none" class="ydpe431f14dyiv0019821171">
> > > using the ICA reconstructed data and not the hand edited or<br clear="none" class="ydpe431f14dyiv0019821171">
> > > artifact deleted original data samples like other commercial<br clear="none" class="ydpe431f14dyiv0019821171">
> > > companies do. Your ICA reconstructed time series is actually less<br clear="none" class="ydpe431f14dyiv0019821171">
> > > different than the original phase difference in comparison to the<br clear="none" class="ydpe431f14dyiv0019821171">
> > > WinEEG ICA. Nonetheless, both your ICA reconstruction and the<br clear="none" class="ydpe431f14dyiv0019821171">
> > > WinEEG reconstructions are significantly different than the<br clear="none" class="ydpe431f14dyiv0019821171">
> > > original<br clear="none" class="ydpe431f14dyiv0019821171">
> recording.<br clear="none" class="ydpe431f14dyiv0019821171">
> > ><br clear="none" class="ydpe431f14dyiv0019821171">
> > > Best regards,<br clear="none" class="ydpe431f14dyiv0019821171">
> > ><br clear="none" class="ydpe431f14dyiv0019821171">
> > > Robert<br clear="none" class="ydpe431f14dyiv0019821171">
> > ><br clear="none" class="ydpe431f14dyiv0019821171">
> > > Cp<br clear="none" class="ydpe431f14dyiv0019821171">
> > ><br clear="none" class="ydpe431f14dyiv0019821171">
> > ><br clear="none" class="ydpe431f14dyiv0019821171">
> > > On Tuesday, June 20, 2017, 1:12:41 AM EDT, Arnaud Delorme<br clear="none" class="ydpe431f14dyiv0019821171">
> > > <<a shape="rect" class="ydpe431f14dyiv0019821171" href="mailto:arno@ucsd.edu" rel="nofollow" target="_blank">arno@ucsd.edu</a> <mailto:<a shape="rect" class="ydpe431f14dyiv0019821171" href="mailto:arno@ucsd.edu" rel="nofollow" target="_blank">arno@ucsd.edu</a>> <mailto:<a shape="rect" class="ydpe431f14dyiv0019821171" href="mailto:arno@ucsd.edu" rel="nofollow" target="_blank">arno@ucsd.edu</a><br clear="none" class="ydpe431f14dyiv0019821171">
> <mailto:<a shape="rect" class="ydpe431f14dyiv0019821171" href="mailto:arno@ucsd.edu" rel="nofollow" target="_blank">arno@ucsd.edu</a>>>> wrote:<br clear="none" class="ydpe431f14dyiv0019821171">
> > ><br clear="none" class="ydpe431f14dyiv0019821171">
> > ><br clear="none" class="ydpe431f14dyiv0019821171">
> > > Dear Robert,<br clear="none" class="ydpe431f14dyiv0019821171">
> > ><br clear="none" class="ydpe431f14dyiv0019821171">
> > > 1) *On my ICA decomposition analysis on your data.* You have<br clear="none" class="ydpe431f14dyiv0019821171">
> > > selected a subset of the file where there is 1 minute and 41<br clear="none" class="ydpe431f14dyiv0019821171">
> > > second data of eye free data. I was only able to select 40 seconds<br clear="none" class="ydpe431f14dyiv0019821171">
> > > in the same file, and I also showed that even in this short file,<br clear="none" class="ydpe431f14dyiv0019821171">
> > > there was some residual eye movements. Jason and Stefan agreed<br clear="none" class="ydpe431f14dyiv0019821171">
> > > with me. This is the reason why ICA components power spectrum over<br clear="none" class="ydpe431f14dyiv0019821171">
> > > frontal channels (and frontal channels only) was affected below 10<br clear="none" class="ydpe431f14dyiv0019821171">
> > > Hz frequency band in my data analysis. So on my ICA decomposition,<br clear="none" class="ydpe431f14dyiv0019821171">
> > > our disagreement comes from the interpretation. You feel that the<br clear="none" class="ydpe431f14dyiv0019821171">
> > > power we remove at low frequency in frontal channel is not eye<br clear="none" class="ydpe431f14dyiv0019821171">
> > > movement. In an attempt to convince you, I have picked up a clean<br clear="none" class="ydpe431f14dyiv0019821171">
> > > region from your EDF dataset, and did some dipole localization at<br clear="none" class="ydpe431f14dyiv0019821171">
> > > this latency. We see that in the clean data, the best dipolar fit<br clear="none" class="ydpe431f14dyiv0019821171">
> > > (with 2 symmetrical dipoles) ends up near the eye balls with a<br clear="none" class="ydpe431f14dyiv0019821171">
> > > residual variance of 6.9%. Hopefully this convinces you that your<br clear="none" class="ydpe431f14dyiv0019821171">
> > > data is not free of eye movement artifacts. If you are willing to<br clear="none" class="ydpe431f14dyiv0019821171">
> > > take a step further you might contemplate the idea that ICA can remove this residual spurious activity.<br clear="none" class="ydpe431f14dyiv0019821171">
> > ><br clear="none" class="ydpe431f14dyiv0019821171">
> > > 2) *On the WinEEG ICA decomposition analysis.* It is critical for<br clear="none" class="ydpe431f14dyiv0019821171">
> > > us to see the scalp topography (and if possible continuous<br clear="none" class="ydpe431f14dyiv0019821171">
> > > activity) of the components the people at the Australia workshop<br clear="none" class="ydpe431f14dyiv0019821171">
> > > selected. Without this, it is not possible for us to comment on<br clear="none" class="ydpe431f14dyiv0019821171">
> > > the cleaned data. I agree with you that there was some phase<br clear="none" class="ydpe431f14dyiv0019821171">
> > > distortion in alpha (visible directly in the raw data in the first<br clear="none" class="ydpe431f14dyiv0019821171">
> > > email you sent) and that this should not be the case. However,<br clear="none" class="ydpe431f14dyiv0019821171">
> > > without seing the ICA decomposition, it is not possible for us to<br clear="none" class="ydpe431f14dyiv0019821171">
> > > conclude as to wether people selected the wrong ICA components or<br clear="none" class="ydpe431f14dyiv0019821171">
> > > if the ICA decomposition implemented in this software is buggy<br clear="none" class="ydpe431f14dyiv0019821171">
> > > (ICA is not a simple algorithm and it is sensitive to numerical<br clear="none" class="ydpe431f14dyiv0019821171">
> > > imprecision and a lot of other parameters - a suboptimal<br clear="none" class="ydpe431f14dyiv0019821171">
> > > implementation could easily explain the WinEEG results). Also, you<br clear="none" class="ydpe431f14dyiv0019821171">
> > > seem to imply that the WinEEG people were running ICA on their<br clear="none" class="ydpe431f14dyiv0019821171">
> > > data then throwing away the raw data (which is why their ICA<br clear="none" class="ydpe431f14dyiv0019821171">
> > > biased neurofeedback database is useless for practical purposes).<br clear="none" class="ydpe431f14dyiv0019821171">
> > > Is that correct? One should never throw away the raw data. If they<br clear="none" class="ydpe431f14dyiv0019821171">
> > > did throw away the raw data, it is an indication that the WinEEG<br clear="none" class="ydpe431f14dyiv0019821171">
> > > are not rigorous in their approach and therefore might not have<br clear="none" class="ydpe431f14dyiv0019821171">
> > > implemented ICA in an optimal way. If it is not the case, one may<br clear="none" class="ydpe431f14dyiv0019821171">
> > > easily reconstruct the database of measures with or without ICA<br clear="none" class="ydpe431f14dyiv0019821171">
> > > decomposition (assuming ICA is done right which does not seem to<br clear="none" class="ydpe431f14dyiv0019821171">
> > > be the case) then assess data measure distoritions (power, phase index, etc…) in a statistical fashion.<br clear="none" class="ydpe431f14dyiv0019821171">
> > ><br clear="none" class="ydpe431f14dyiv0019821171">
> > > Best wishes,<br clear="none" class="ydpe431f14dyiv0019821171">
> > ><br clear="none" class="ydpe431f14dyiv0019821171">
> > > Arno<br clear="none" class="ydpe431f14dyiv0019821171">
> > ><br clear="none" class="ydpe431f14dyiv0019821171">
> > > <a shape="rect" class="ydpe431f14dyiv0019821171" href="http://sccn.ucsd.edu/%7Earno/download/clean_edf_file_analysis2.pdf" rel="nofollow" target="_blank">http://sccn.ucsd.edu/~arno/dow nload/clean_edf_file_analysis2 .pdf</a><br clear="none" class="ydpe431f14dyiv0019821171">
> <<a shape="rect" class="ydpe431f14dyiv0019821171" href="http://sccn.ucsd.edu/%7Earno/download/clean_edf_file_analysis2.pdf%20" rel="nofollow" target="_blank">http://sccn.ucsd.edu/%7Earno/ download/clean_edf_file_analys is2.pdf%20</a><br clear="none" class="ydpe431f14dyiv0019821171">
> ><br clear="none" class="ydpe431f14dyiv0019821171">
> > <<a shape="rect" class="ydpe431f14dyiv0019821171" href="http://sccn.ucsd.edu/%7Earno/download/clean_edf_file_analysis2.pdf%" rel="nofollow" target="_blank">http://sccn.ucsd.edu/%7Earno/ download/clean_edf_file_analys is2.pdf%</a><br clear="none" class="ydpe431f14dyiv0019821171">
> > 20><br clear="none" class="ydpe431f14dyiv0019821171">
> > > <<a shape="rect" class="ydpe431f14dyiv0019821171" href="http://sccn.ucsd.edu/%7Earno/download/clean_edf_file_analysis2.pd" rel="nofollow" target="_blank">http://sccn.ucsd.edu/%7Earno/ download/clean_edf_file_analys is2.pd</a><br clear="none" class="ydpe431f14dyiv0019821171">
> > > f><br clear="none" class="ydpe431f14dyiv0019821171">
> > ><br clear="none" class="ydpe431f14dyiv0019821171">
> > >> On Jun 18, 2017, at 11:44 AM, Robert Thatcher<br clear="none" class="ydpe431f14dyiv0019821171">
> > <<a shape="rect" class="ydpe431f14dyiv0019821171" href="mailto:rwthatcher2@yahoo.com" rel="nofollow" target="_blank">rwthatcher2@yahoo.com</a> <mailto:<a shape="rect" class="ydpe431f14dyiv0019821171" href="mailto:rwthatcher2@yahoo.com" rel="nofollow" target="_blank">rwthatcher2@yahoo.com</a>><br clear="none" class="ydpe431f14dyiv0019821171">
> <mailto:<a shape="rect" class="ydpe431f14dyiv0019821171" href="mailto:rwthatcher2@yahoo.com" rel="nofollow" target="_blank">rwthatcher2@yahoo.com</a> <mailto:<a shape="rect" class="ydpe431f14dyiv0019821171" href="mailto:rwthatcher2@yahoo.com" rel="nofollow" target="_blank">rwthatcher2@yahoo.com</a>> ><br clear="none" class="ydpe431f14dyiv0019821171">
> > >> <mailto:<a shape="rect" class="ydpe431f14dyiv0019821171" href="mailto:rwthatcher2@yahoo.com" rel="nofollow" target="_blank">rwthatcher2@yahoo.com</a> <mailto:<a shape="rect" class="ydpe431f14dyiv0019821171" href="mailto:rwthatcher2@yahoo.com" rel="nofollow" target="_blank">rwthatcher2@yahoo.com</a>><br clear="none" class="ydpe431f14dyiv0019821171">
> <mailto:<a shape="rect" class="ydpe431f14dyiv0019821171" href="mailto:rwthatcher2@yahoo.com" rel="nofollow" target="_blank">rwthatcher2@yahoo.com</a> <mailto:<a shape="rect" class="ydpe431f14dyiv0019821171" href="mailto:rwthatcher2@yahoo.com" rel="nofollow" target="_blank">rwthatcher2@yahoo.com</a>> >>> wrote:<br clear="none" class="ydpe431f14dyiv0019821171">
><br clear="none" class="ydpe431f14dyiv0019821171">
> ><br clear="none" class="ydpe431f14dyiv0019821171">
> > >><br clear="none" class="ydpe431f14dyiv0019821171">
> > >> <Pre-ICA-Hand Artifact free selections.edf><br clear="none" class="ydpe431f14dyiv0019821171">
> ><br clear="none" class="ydpe431f14dyiv0019821171">
> > ><br clear="none" class="ydpe431f14dyiv0019821171">
> > ><br clear="none" class="ydpe431f14dyiv0019821171">
> > ><br clear="none" class="ydpe431f14dyiv0019821171">
> > > Dieser Nachrichteninhalt wird auf Anfrage komplett heruntergeladen.<br clear="none" class="ydpe431f14dyiv0019821171">
> ><br clear="none" class="ydpe431f14dyiv0019821171">
> ><br clear="none" class="ydpe431f14dyiv0019821171">
><br clear="none" class="ydpe431f14dyiv0019821171">
<br clear="none" class="ydpe431f14dyiv0019821171">
______________________________ _________________<br clear="none" class="ydpe431f14dyiv0019821171">
Eeglablist page: <a shape="rect" class="ydpe431f14dyiv0019821171" href="http://sccn.ucsd.edu/eeglab/eeglabmail.html" rel="nofollow" target="_blank">http://sccn.ucsd.edu/eeglab/ee glabmail.html</a><br clear="none" class="ydpe431f14dyiv0019821171">
To unsubscribe, send an empty email to <a shape="rect" class="ydpe431f14dyiv0019821171" href="mailto:eeglablist-unsubscribe@sccn.ucsd.edu" rel="nofollow" target="_blank">eeglablist-unsubscribe@sccn.uc sd.edu</a><br clear="none" class="ydpe431f14dyiv0019821171">
For digest mode, send an email with the subject "set digest mime" to <a shape="rect" class="ydpe431f14dyiv0019821171" href="mailto:eeglablist-request@sccn.ucsd.edu" rel="nofollow" target="_blank">eeglablist-request@sccn.ucsd.e du</a><br clear="none" class="ydpe431f14dyiv0019821171">
<br clear="none" class="ydpe431f14dyiv0019821171">
______________________________ _________________<br clear="none" class="ydpe431f14dyiv0019821171">
Eeglablist page: <a shape="rect" class="ydpe431f14dyiv0019821171" href="http://sccn.ucsd.edu/eeglab/eeglabmail.html" rel="nofollow" target="_blank">http://sccn.ucsd.edu/eeglab/ee glabmail.html</a><br clear="none" class="ydpe431f14dyiv0019821171">
To unsubscribe, send an empty email to <a shape="rect" class="ydpe431f14dyiv0019821171" href="mailto:eeglablist-unsubscribe@sccn.ucsd.edu" rel="nofollow" target="_blank">eeglablist-unsubscribe@sccn.uc sd.edu</a><br clear="none" class="ydpe431f14dyiv0019821171">
For digest mode, send an email with the subject "set digest mime" to <a shape="rect" class="ydpe431f14dyiv0019821171" href="mailto:eeglablist-request@sccn.ucsd.edu" rel="nofollow" target="_blank">eeglablist-request@sccn.ucsd.e du</a></div></div></blockquote></div><br clear="none" class="ydpe431f14dyiv0019821171"></div></div></div></div>
<br clear="none" class="ydpe431f14dyiv0019821171"><br clear="none" class="ydpe431f14dyiv0019821171">_______________________________________________<br clear="none" class="ydpe431f14dyiv0019821171">eeglablist mailing list <a shape="rect" class="ydpe431f14dyiv0019821171" href="mailto:eeglablist@sccn.ucsd.edu" rel="nofollow" target="_blank">eeglablist@sccn.ucsd.edu</a><br clear="none" class="ydpe431f14dyiv0019821171">Eeglablist page: <a shape="rect" class="ydpe431f14dyiv0019821171" href="http://www.sccn.ucsd.edu/eeglab/eeglabmail.html" rel="nofollow" target="_blank">http://www.sccn.ucsd.edu/eeglab/eeglabmail.html</a><br clear="none" class="ydpe431f14dyiv0019821171">To unsubscribe, send an empty email to <a shape="rect" class="ydpe431f14dyiv0019821171" href="mailto:eeglablist-unsub@sccn.ucsd.edu" rel="nofollow" target="_blank">eeglablist-unsub@sccn.ucsd.edu</a><br clear="none" class="ydpe431f14dyiv0019821171">To switch to non-digest mode, send an empty email to <a shape="rect" class="ydpe431f14dyiv0019821171" href="mailto:eeglablist-nodigest@sccn.ucsd.edu" rel="nofollow" target="_blank">eeglablist-nodigest@sccn.ucsd.edu</a></div></blockquote></div><br clear="none" class="ydpe431f14dyiv0019821171"></div>_______________________________________________<br clear="none" class="ydpe431f14dyiv0019821171">Eeglablist page: <a shape="rect" class="ydpe431f14dyiv0019821171" href="http://sccn.ucsd.edu/eeglab/eeglabmail.html" rel="nofollow" target="_blank">http://sccn.ucsd.edu/eeglab/eeglabmail.html</a><br clear="none" class="ydpe431f14dyiv0019821171">To unsubscribe, send an empty email to <a shape="rect" class="ydpe431f14dyiv0019821171" href="mailto:eeglablist-unsubscribe@sccn.ucsd.edu" rel="nofollow" target="_blank">eeglablist-unsubscribe@sccn.ucsd.edu</a><br clear="none" class="ydpe431f14dyiv0019821171">For digest mode, send an email with the subject "set digest mime" to <a shape="rect" class="ydpe431f14dyiv0019821171" href="mailto:eeglablist-request@sccn.ucsd.edu" rel="nofollow" target="_blank">eeglablist-request@sccn.ucsd.edu</a><br clear="none" class="ydpe431f14dyiv0019821171"><div class="ydpe431f14dyiv0019821171">
<div class="ydpe431f14dyiv0019821171" style="color:rgb(0, 0, 0);letter-spacing:normal;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word;"><div class="ydpe431f14dyiv0019821171" style="color:rgb(0, 0, 0);letter-spacing:normal;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word;"><div class="ydpe431f14dyiv0019821171" style="color:rgb(0, 0, 0);letter-spacing:normal;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word;"><div class="ydpe431f14dyiv0019821171" style="color:rgb(0, 0, 0);letter-spacing:normal;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word;">--------------------------------------------------------------------------------<br clear="none" class="ydpe431f14dyiv0019821171"><br clear="none" class="ydpe431f14dyiv0019821171">Joseph Dien, PhD</div><div class="ydpe431f14dyiv0019821171" style="color:rgb(0, 0, 0);letter-spacing:normal;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word;">Senior Research Scientist</div><div class="ydpe431f14dyiv0019821171" style="color:rgb(0, 0, 0);letter-spacing:normal;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word;">Maryland Neuroimaging Center</div><div class="ydpe431f14dyiv0019821171" style="color:rgb(0, 0, 0);letter-spacing:normal;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word;">University of Maryland, College Park</div><div class="ydpe431f14dyiv0019821171" style="color:rgb(0, 0, 0);letter-spacing:normal;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word;"><a shape="rect" class="ydpe431f14dyiv0019821171" href="http://joedien.com" rel="nofollow" target="_blank">http://joedien.com</a></div></div></div></div>
</div>
<br clear="none" class="ydpe431f14dyiv0019821171"></div></div></div><div class="ydpe431f14dyqt0940008496" id="ydpe431f14dyqt54627">_______________________________________________<br clear="none">Eeglablist page: <a shape="rect" href="http://sccn.ucsd.edu/eeglab/eeglabmail.html" rel="nofollow" target="_blank">http://sccn.ucsd.edu/eeglab/eeglabmail.html</a><br clear="none">To unsubscribe, send an empty email to <a shape="rect" href="mailto:eeglablist-unsubscribe@sccn.ucsd.edu" rel="nofollow" target="_blank">eeglablist-unsubscribe@sccn.ucsd.edu</a><br clear="none">For digest mode, send an email with the subject "set digest mime" to <a shape="rect" href="mailto:eeglablist-request@sccn.ucsd.edu" rel="nofollow" target="_blank">eeglablist-request@sccn.ucsd.edu</a></div></div></div></div></div></body></html>