<html><head><meta http-equiv="Content-Type" content="text/html charset=us-ascii"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">Dear Elizabeth,<div class=""><br class=""></div><div class="">The RMS normalization of component only scale component activity and the ICA weight and inverse weight matrix. The back projection of component in channel space is identical wether you RMS normalize ICA components or not.</div><div class=""><br class=""></div><div class="">Arno</div><div class=""><br class=""><div><blockquote type="cite" class=""><div class="">On Aug 11, 2015, at 7:22 AM, Elizabeth Milne <<a href="mailto:e.milne@sheffield.ac.uk" class="">e.milne@sheffield.ac.uk</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><meta http-equiv="Content-Type" content="text/html; charset=utf-8" class=""><div dir="ltr" class="">Dear Scott,<div class="">Thanks - this is helpful. Dividing the IC activation by the RMS value of the its scalp map gives me very similar, but not identical, data to the IC activation [using the function <i class="">rms</i> in matlab].</div><div class=""><br class=""></div><div class="">In the majority of data sets and ICs I have looked at, the time-course of the RMS-projected activation is similar (albeit different magnitude) to the time-course of the forward-projected activation of the same IC. In these cases the RMS-projected activations are easy to interpret in the context of looking for characteristic VEP deflections such as P1 / N1. However, in a small number of cases, the time-course of the RMS-projected IC activation is orthogonal to the time-course of the forward-projected IC activation  i.e. the component appears inverted. In these cases, identifying P1 as a positive peak and N1 as a negative peak is no longer possible, although it is clear that these deflections are in the data. </div><div class=""><br class=""></div><div class="">In these cases is it legitimate to subtract the RMS-projected activation from 0 in order to invert it? Previously I was taking Makoto's advice and forward-projecting the IC activation to the channel that it projected to most strongly in order to ascertain the 'true' polarity of the activation. However, given your comments above it sounds like working with the RMS-projected IC activations is the preferable route - if I can be confident of the polarity of this activation. </div><div class=""><br class=""></div><div class="">Thanks for your input,</div><div class="">Elizabeth</div><div class=""><br class=""></div><div class=""><br class=""></div><div class=""><br class=""></div><div class=""><br class=""></div></div><div class="gmail_extra"><br class=""><div class="gmail_quote">On 10 August 2015 at 16:48, Scott Makeig <span dir="ltr" class=""><<a href="mailto:smakeig@ucsd.edu" target="_blank" class="">smakeig@ucsd.edu</a>></span> wrote:<br class=""><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr" class="">Elizabeth and all -<div class=""><br class=""></div><div class="">Makoto is correct that the term 'back projection' of an IC to the scalp can be viewed also (and more fundamentally)  as a <i class="">forward</i> projection from cortical source to the scalp sensors.</div><div class=""><br class=""></div><div class="">I blieve the only EEG phenomena of interest are the source data activities (e.g., IC activations) -- the source projections to the scalp channels (whose sum we record with scalp electrodes) are essentially epiphenomena (more or less arbitrary source sums and differences, in themselves of no fundamental interest).</div><div class=""><br class=""></div><div class="">Now, the most correct unit of measurement for cortical sources is source<i class=""> current density per cortical area </i>(for example, per mm^2).  However, to convert scalp measurements into that unit one needs a high-resolution electrical forward head model (sources to scalp electrodes) and an ideal inverse model (source scalp projection map to cortical source patch).  The former must include accurate values for tissue conductivities within the individual subject, not something we currently have access to (though see below).</div><div class=""><br class=""></div><div class="">Therefore, the most reliable unit of measurement we have for the activity of an identified source in scalp-recorded data is  'RMS uV' -- the root mean-square average projection from the source to the scalp. Even this has ambiguity -- To the <i class="">whole</i> scalp, even places where we did not place electrodes??   A simple alternative is to measure IC activity (unaveraged or averaged) as the RMS (mean) projection of the source to the actual scalp electrodes:</div><div class=""><br class=""></div><div class="">e.g., If IC is the index of an IC in an EEG data structure.</div><div class=""><br class=""></div><div class=""> <i class=""> >> RMS_uv(IC) = EEG.icaact(IC ) / RMS(EEG.icawinv(IC));</i>  % divide the IC activation by</div><div class="">                                                                                                   % the RMS value of its scalp map</div><div class=""><br class=""></div><div class="">Unfortunately, the function <i class="">RMS()</i> I used in the meta-code above is not a Matlab function (by that name)... I am away from a Matlab installation at present and forget the actual function call.</div><div class=""><br class=""></div><div class="">It may be that this normalization (making <i class="">RMS(EEG.icawinv) = ones()</i> ) is already performed in the default output of runica / binica / AMICA - To be safe, however, run the meta-code above explicitly or test the equality above ...</div><div class=""><br class=""></div><div class="">All units of measurement for IC activations (activities) should be labelled 'RMS uV/chan' -- I will work with Arno and Ramon to see that this is corrected in the current functions.</div><div class=""><br class=""></div><div class="">Note: This removes the problem, 'Which channel projection is the best one to use to measure the contribution of a source to the scalp (ERP or raw) data?'  </div><div class=""><br class=""></div><div class="">Scott Makeig</div><div class=""><br class=""></div><div class="">p.s. To see which ICs contribute most strongly, learn to use and interpret <i class="">envtopo()</i> -- called in the EEGLAB GUI by '<i class="">Plot > Component ERPs > With Scalp Maps</i>'.  Makoto has recently supervised building of a plug-in to perform this function at the STUDY level (here the question is, 'Which IC clusters contribute most strongly to the grand mean ERP or ERP difference?')</div><div class=""><br class=""></div><div class="">p.p.s  Zeynep Akalin Acar and I have developed a new method for estimating skull conductivity from the EEG data (given a subject MR head image allowing development of an individual electrical forward head model). Our report on this should be accepted soon - more details then or at the upcoming EEGLAB workshop in UK...</div></div><div class="gmail_extra"><div class=""><div class="h5"><br class=""><div class="gmail_quote">On Mon, Aug 3, 2015 at 6:09 PM, Makoto Miyakoshi <span dir="ltr" class=""><<a href="mailto:mmiyakoshi@ucsd.edu" target="_blank" class="">mmiyakoshi@ucsd.edu</a>></span> wrote:<br class=""><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr" class="">Dear Elizabeth,<div class=""><br class=""></div><div class="">A projection from ICs to channels is a 'forward projection'. We've been using the word 'backprojection' in the wrong way, so I corrected it (at least in my functions).</div><div class=""><br class=""></div><div class="">In ICA, the polarity of map and time-series signal can be flipped by applying -1 to the both. This is a fundamental uncertainty.</div><span class=""><div class=""><br class=""></div><div class=""><span style="font-size:12.8000001907349px" class="">> If I want to use this, or any other, IC to measure the amplitude / latency of particular ERP deflections it is necessary to know which is the correct channel to choose for back-projection. However, I am unsure as to how to make this decision objectively.</span><br class=""></div><div class=""> </div></span><div class="">I wrote pvaftopo() plugin for this purpose</div><div class=""><a href="http://sccn.ucsd.edu/wiki/Plugin_list_process" target="_blank" class="">http://sccn.ucsd.edu/wiki/Plugin_list_process</a><br class=""></div><div class="">pvaf is a new concept Scott proposed nearly 20 years ago which is</div><div class=""><br class=""></div><div class="">percent variance accounted for (pvaf) = 100-100*var(all_data - selected_ICs)/var(all_data)</div><div class=""><br class=""></div><div class="">We'll publish IEEE EMBC conference proceedings soon for this year about the study-level envtopo and pvaf, so you may want to refer to the paper once it's published (Lee et al., IEEE EMBC conf proc 2015)</div><div class=""><br class=""></div><div class="">Pick up the channel that is closest to the peak in the scalp maps and you can choose the 'most accounted' channel objectively.</div><div class=""><br class=""></div><div class="">I'm happy that there is finally a user asking for this plugin!</div><div class=""><br class=""></div><div class="">Makoto</div><div class="gmail_extra"><br class=""><div class="gmail_quote"><div class=""><div class="">On Fri, Jul 31, 2015 at 2:13 PM, Elizabeth Milne <span dir="ltr" class=""><<a href="mailto:e.milne@sheffield.ac.uk" target="_blank" class="">e.milne@sheffield.ac.uk</a>></span> wrote:<br class=""></div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class=""><div class=""><div dir="ltr" class=""><span style="font-size:12.8000001907349px" class="">Dear EEGLAB team,</span><div style="font-size:12.8000001907349px" class="">I have a question about resolving the "true" polarity of IC ERPs by back-projecting to channel data. Using erpimage, I have obtained the ERP of an IC, but the polarity is different (orthogonal) when back-projecting to different channels. In another example from a different IC, the polarity of the ERP does not change, but the amplitude of the signal is different depending on which channel is selected for back-projection.</div><div style="font-size:12.8000001907349px" class=""><br class=""></div><div style="font-size:12.8000001907349px" class="">If I want to use this, or any other, IC to measure the amplitude / latency of particular ERP deflections it is necessary to know which is the correct channel to choose for back-projection. However, I am unsure as to how to make this decision objectively. Are there particular features that I should look at to make this decision more systematic, i.e. dipole projection? % power accounted for?</div><div style="font-size:12.8000001907349px" class=""><br class=""></div><div style="font-size:12.8000001907349px" class="">If anyone can offer advice, or point out errors in my understanding I'd be very grateful. </div><div style="font-size:12.8000001907349px" class=""><br class=""></div><div style="font-size:12.8000001907349px" class="">Many thanks,</div><div style="font-size:12.8000001907349px" class="">Elizabeth</div><span class=""><font color="#888888" class=""><div class=""><br class=""></div>-- <br class=""><div class=""><div dir="ltr" class=""><div class="">Elizabeth Milne<br class="">Reader in Cognitive Neuroscience<br class="">University of Sheffield<br class="">Western Bank<br class="">Sheffield<br class="">S10 2TN<br class=""> <br class=""><a href="tel:%2B44%20%280%29%20114%20222%206558" value="+441142226558" target="_blank" class="">+44 (0) 114 222 6558</a><br class=""> <br class=""><a href="http://www.autismresearchlab.group.shef.ac.uk/" target="_blank" class="">http://www.autismresearchlab.group.shef.ac.uk/</a></div></div></div>
</font></span></div>
<br class=""></div></div>_______________________________________________<br class="">
Eeglablist page: <a href="http://sccn.ucsd.edu/eeglab/eeglabmail.html" rel="noreferrer" target="_blank" class="">http://sccn.ucsd.edu/eeglab/eeglabmail.html</a><br class="">
To unsubscribe, send an empty email to <a href="mailto:eeglablist-unsubscribe@sccn.ucsd.edu" target="_blank" class="">eeglablist-unsubscribe@sccn.ucsd.edu</a><br class="">
For digest mode, send an email with the subject "set digest mime" to <a href="mailto:eeglablist-request@sccn.ucsd.edu" target="_blank" class="">eeglablist-request@sccn.ucsd.edu</a><span class=""><font color="#888888" class=""><br class=""></font></span></blockquote></div><span class=""><font color="#888888" class=""><br class=""><br clear="all" class=""><div class=""><br class=""></div>-- <br class=""><div class=""><div dir="ltr" class="">Makoto Miyakoshi<br class="">Swartz Center for Computational Neuroscience<br class="">Institute for Neural Computation, University of California San Diego<br class=""></div></div>
</font></span></div></div>
<br class="">_______________________________________________<br class="">
Eeglablist page: <a href="http://sccn.ucsd.edu/eeglab/eeglabmail.html" rel="noreferrer" target="_blank" class="">http://sccn.ucsd.edu/eeglab/eeglabmail.html</a><br class="">
To unsubscribe, send an empty email to <a href="mailto:eeglablist-unsubscribe@sccn.ucsd.edu" target="_blank" class="">eeglablist-unsubscribe@sccn.ucsd.edu</a><br class="">
For digest mode, send an email with the subject "set digest mime" to <a href="mailto:eeglablist-request@sccn.ucsd.edu" target="_blank" class="">eeglablist-request@sccn.ucsd.edu</a><br class=""></blockquote></div><br class=""><br clear="all" class=""><div class=""><br class=""></div>-- <br class=""></div></div><div class="">Scott Makeig, Research Scientist and Director, Swartz Center for Computational Neuroscience, Institute for Neural Computation, University of California San Diego, La Jolla CA 92093-0961, <a href="http://sccn.ucsd.edu/~scott" target="_blank" class="">http://sccn.ucsd.edu/~scott</a></div>
</div>
</blockquote></div><br class=""><br clear="all" class=""><div class=""><br class=""></div>-- <br class=""><div class="gmail_signature"><div dir="ltr" class=""><div class="">Elizabeth Milne<br class="">Reader in Cognitive Neuroscience<br class="">University of Sheffield<br class="">Western Bank<br class="">Sheffield<br class="">S10 2TN<br class=""> <br class="">+44 (0) 114 222 6558<br class=""> <br class=""><a href="http://www.autismresearchlab.group.shef.ac.uk/" target="_blank" class="">http://www.autismresearchlab.group.shef.ac.uk/</a></div></div></div>
</div>
_______________________________________________<br class="">Eeglablist page: <a href="http://sccn.ucsd.edu/eeglab/eeglabmail.html" class="">http://sccn.ucsd.edu/eeglab/eeglabmail.html</a><br class="">To unsubscribe, send an empty email to <a href="mailto:eeglablist-unsubscribe@sccn.ucsd.edu" class="">eeglablist-unsubscribe@sccn.ucsd.edu</a><br class="">For digest mode, send an email with the subject "set digest mime" to <a href="mailto:eeglablist-request@sccn.ucsd.edu" class="">eeglablist-request@sccn.ucsd.edu</a></div></blockquote></div><br class=""></div></body></html>