<div dir="ltr">Hey Marijn,<div><br></div><div>Many thanks for this very basic but really important message, I think it makes me feel much better why select average reference.</div><div>Thanks,</div><div><br></div><div>Mengting</div></div><div class="gmail_extra"><br><div class="gmail_quote">2017-09-28 3:54 GMT-04:00 Marijn van Vliet <span dir="ltr"><<a href="mailto:w.m.vanvliet@gmail.com" target="_blank">w.m.vanvliet@gmail.com</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div text="#000000" bgcolor="#FFFFFF">
<p>Hi Mengting,</p>
<p>let me chime in here. Whenever you are doing EEG analysis, the
choice of reference is important. EEG measures voltage, which is
the difference in electric potential between two points. This
means, there's technically no such thing as "the voltage at
electrode Cz", there is only "the difference in electric potential
between electrodes Cz and the one I stuck on the nose" (or
whatever reference electrode you used, could be CMS).</p>
<p>After the recording, it is easy to "re-reference" the data, which
is choosing a difference electrode as a reference. For example, to
move from the nose reference to channel Pz as reference, you can
just do: Cz = Cz - Pz and you have a voltage at the Cz channel
with Pz as reference.</p>
<p>You can imagine your data often looks completely different,
depending on the reference you choose. If you choose Pz as a
reference, the signal at Pz will be 0 (Pz - Pz) and all channels
surrounding Pz only have a very small signal. If Pz contained
oscillatory activity (maybe alpha rhythms), you would see alpha
power at every electrode *except* surrounding Pz. For example you
would see "negative" alpha at Fpz (remember, we're looking at Fpz
- Pz).</p>
<p>So maybe Pz is not a very useful reference. We would like to have
some "neutral" reference that doesn't contain signals of interest.
This is why nose and mastoids are popular locations for the
reference. But there's another option: choose the average of all
electrodes as "virtual reference". With this reference, at any
time, mean(Cz, Pz, ...) = 0. This is the reference assumed by
forward models, which is why when doing anything that involves the
forward model (source localization, connectivity analysis at the
source level, etc.) your data needs to be in this reference.</p>
<p>One thing you need to be aware of when re-reference to an average
reference is that the reference is, you know, the average of all
electrodes. If one of the electrodes is noisy (maybe it got loose
during the experiment), this ruins the reference signal, and since
the reference is subtracted from all sensors, it ruins the signal
at *all* sensors.</p>
<p>So, first use a sensors which you know for sure has a good signal
as a reference (preferably a "neutral" reference like the noise,
mastoids or earlobe), then inspect the data to spot all sensors
that have problems. Either remove them from the data or have
MNE-Python reconstruct the signal at the bad sensor by
interpolating the neighbouring sensors. Then, when the data is
squeaky clean, move to the average reference.<br>
</p>
<p>regards,<br>
Marijn.</p><div><div class="h5">
<p><br>
</p>
<br>
<div class="m_-3253518844974229613moz-cite-prefix">On 09/28/2017 09:58 AM, Christopher
Bailey wrote:<br>
</div>
<blockquote type="cite">
Hi Mengting,
<div><br>
</div>
<div>
<blockquote type="cite">
<div dir="ltr">
<div>So as I understand, reference is only a
required step for source localization,</div>
</div>
</blockquote>
</div>
<div><br>
</div>
<div>The _average_ reference is required for source
localisation</div>
<div><br>
</div>
<div>
<blockquote type="cite">
<div dir="ltr">
<div>but not any other steps like measuring the
connectivity between brain regions. </div>
</div>
</blockquote>
</div>
<div><br>
</div>
<div>I’m not sure I understand what you mean by “brain
regions” here. If you’re going to measure connectivity in source
space, your need an (inverse) operator of some sort to project
your electrode-data “into the brain”.</div>
<div><br>
</div>
<div>If you want to do sensor-level connectivity
calculations on EEG data, the effect of the reference might
depend on the connectivity metric you use. For those
time/frequency-domain measures I’m aware of (but have little
first-hand experience with), it doesn’t matter which reference
your data is in.</div>
<div><br>
</div>
<div>
<blockquote type="cite">
<div dir="ltr">
<div>Sounds reference removed the DC componets but
remains the oscillation. If so, does it influence the
power in DC components? </div>
</div>
</blockquote>
<br>
</div>
<div>Yes and no. Relative differences within and between
channels remain the same in any (proper) reference. Note though
that the time-domain is not involved in (re-)referencing, so
possible ‘DC components’ still remain in the data.</div>
<div><br>
</div>
<div>
<blockquote type="cite">
<div dir="ltr">
<div> Also, does it matter the conventional source
localization methods such as dipole fitting?</div>
</div>
</blockquote>
<br>
</div>
<div>All localisation methods require the average
reference, otherwise the recorded and predicted (forward model)
data cannot be compared.</div>
<div><br>
</div>
<div>Best,</div>
<div><br>
</div>
<div>Chris</div>
<div><br>
<div>
<blockquote type="cite">
<div>On 27 Sep 2017, at 20.54, Liu Mengting <<a href="mailto:bigting84@gmail.com" target="_blank">bigting84@gmail.com</a>>
wrote:</div>
<br class="m_-3253518844974229613Apple-interchange-newline">
<div>
<div dir="ltr">Hi Chris,
<div><br>
</div>
<div>Thanks for the info, this really helps in
understanding the whole procedures. So as I
understand, reference is only a required step for
source localization, but not any other steps like
measuring the connectivity between brain regions.
Sounds reference removed the DC componets but remains
the oscillation. If so, does it influence the power in
DC components? Also, does it matter the conventional
source localization methods such as dipole fitting?</div>
<div><br>
</div>
<div>Really appreciate for help,</div>
<div><br>
</div>
<div>Mengting</div>
</div>
<div class="gmail_extra"><br>
<div class="gmail_quote">2017-09-25 4:55 GMT-04:00
Christopher Bailey <span dir="ltr">
<<a href="mailto:cjb@cfin.au.dk" target="_blank">cjb@cfin.au.dk</a>></span>:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div style="word-wrap:break-word">Hi
Mengting,
<div><br>
</div>
<div>The average reference is the only
valid scheme when performing ‘source
localisation’, no matter what you do afterwards.
To first order, the localisation procedure is</div>
<div><br>
</div>
<div>- have measured (EEG) data</div>
<div>- build forward model that computes
sensor-level readings for know sources</div>
<div>- compare measured and predicted
(forward-projected) data to each other</div>
<div>- minimise prediction error under
chosen prior/model</div>
<div><br>
</div>
<div>The predicted data are calculated
relative to a hypothetical absolute reference
potential of zero at infinity, whereas a real
dataset could be referenced to a number of
points on the head. To be able to compare the
two datasets, both are re-referenced to their
respective average: the average reference <i>does not depend on the location of
the on-line reference electrode</i> (as long
as it was functioning properly). After this
re-scaling, the values can be compared directly.</div>
<div><br>
</div>
<div>Note that like any (proper)
re-referencing procedure, taking the common
average only shifts the zero-point; relative
differences between electrode readings remain
unaltered.</div>
<div><br>
</div>
<div>/Chris</div>
<div><br>
</div>
<div>
<div>
<blockquote type="cite">
<div>
<div class="m_-3253518844974229613h5">
<div>On 24 Sep 2017, at 01.28,
Liu Mengting <<a href="mailto:bigting84@gmail.com" target="_blank">bigting84@gmail.com</a>>
wrote:</div>
<br class="m_-3253518844974229613m_7606626925763248917Apple-interchange-newline">
</div>
</div>
<div>
<div>
<div class="m_-3253518844974229613h5">
<div dir="ltr">Hello MNE
users,
<div><br>
</div>
<div>I noticed that in MNE
inverse operations, all EEG data
were forced by an average
referencing. </div>
<div>Does anyone has
insight about how does the average
reference influence the inverse
operations in MNE (I mean use dSPM
or sLoreta)? Especially would this
average referencing influence
functional connectivity measure in
source space (e.g. using phase
locking)?</div>
<div><br>
</div>
<div>Thanks,</div>
<div>Mengting</div>
</div>
</div>
</div>
______________________________<wbr>_________________<br>
Mne_analysis mailing list<br>
<a href="mailto:Mne_analysis@nmr.mgh.harvard.edu" target="_blank">Mne_analysis@nmr.mgh.harvard.e<wbr>du</a><br>
<a href="https://mail.nmr.mgh.harvard.edu/mailman/listinfo/mne_analysis" target="_blank">https://mail.nmr.mgh.harvard.e<wbr>du/mailman/listinfo/mne_analys<wbr>is</a><br>
<br>
<br>
The information in this e-mail is intended
only for the person to whom it is<br>
addressed. If you believe this e-mail was
sent to you in error and the e-mail<br>
contains patient information, please
contact the Partners Compliance HelpLine
at<br>
<a href="http://www.partners.org/complianceline" target="_blank">http://www.partners.org/compli<wbr>anceline</a> . If the
e-mail was sent to you in error<br>
but does not contain patient information,
please contact the sender and properly<br>
dispose of the e-mail.<br>
</div>
</blockquote>
</div>
<br>
</div>
</div>
<br>
______________________________<wbr>_________________<br>
Mne_analysis mailing list<br>
<a href="mailto:Mne_analysis@nmr.mgh.harvard.edu" target="_blank">Mne_analysis@nmr.mgh.harvard.e<wbr>du</a><br>
<a href="https://mail.nmr.mgh.harvard.edu/mailman/listinfo/mne_analysis" rel="noreferrer" target="_blank">https://mail.nmr.mgh.harvard.e<wbr>du/mailman/listinfo/mne_analys<wbr>is</a><br>
<br>
<br>
The information in this e-mail is intended only for
the person to whom it is<br>
addressed. If you believe this e-mail was sent to
you in error and the e-mail<br>
contains patient information, please contact the
Partners Compliance HelpLine at<br>
<a href="http://www.partners.org/complianceline" rel="noreferrer" target="_blank">http://www.partners.org/compli<wbr>anceline</a> . If the e-mail was
sent to you in error<br>
but does not contain patient information, please
contact the sender and properly<br>
dispose of the e-mail.<br>
<br>
</blockquote>
</div>
<br>
</div>
______________________________<wbr>_________________<br>
Mne_analysis mailing list<br>
<a href="mailto:Mne_analysis@nmr.mgh.harvard.edu" target="_blank">Mne_analysis@nmr.mgh.harvard.<wbr>edu</a><br>
<a class="m_-3253518844974229613moz-txt-link-freetext" href="https://mail.nmr.mgh.harvard.edu/mailman/listinfo/mne_analysis" target="_blank">https://mail.nmr.mgh.harvard.<wbr>edu/mailman/listinfo/mne_<wbr>analysis</a><br>
<br>
<br>
The information in this e-mail is intended only for the
person to whom it is<br>
addressed. If you believe this e-mail was sent to you in
error and the e-mail<br>
contains patient information, please contact the Partners
Compliance HelpLine at<br>
<a class="m_-3253518844974229613moz-txt-link-freetext" href="http://www.partners.org/complianceline" target="_blank">http://www.partners.org/<wbr>complianceline</a> . If the e-mail was
sent to you in error<br>
but does not contain patient information, please contact
the sender and properly<br>
dispose of the e-mail.<br>
</div>
</blockquote>
</div>
<br>
</div>
<br>
<fieldset class="m_-3253518844974229613mimeAttachmentHeader"></fieldset>
<br>
<pre>______________________________<wbr>_________________
Mne_analysis mailing list
<a class="m_-3253518844974229613moz-txt-link-abbreviated" href="mailto:Mne_analysis@nmr.mgh.harvard.edu" target="_blank">Mne_analysis@nmr.mgh.harvard.<wbr>edu</a>
<a class="m_-3253518844974229613moz-txt-link-freetext" href="https://mail.nmr.mgh.harvard.edu/mailman/listinfo/mne_analysis" target="_blank">https://mail.nmr.mgh.harvard.<wbr>edu/mailman/listinfo/mne_<wbr>analysis</a>
The information in this e-mail is intended only for the person to whom it is
addressed. If you believe this e-mail was sent to you in error and the e-mail
contains patient information, please contact the Partners Compliance HelpLine at
<a class="m_-3253518844974229613moz-txt-link-freetext" href="http://www.partners.org/complianceline" target="_blank">http://www.partners.org/<wbr>complianceline</a> . If the e-mail was sent to you in error
but does not contain patient information, please contact the sender and properly
dispose of the e-mail.
</pre>
</blockquote>
<br>
</div></div></div>
<br>______________________________<wbr>_________________<br>
Mne_analysis mailing list<br>
<a href="mailto:Mne_analysis@nmr.mgh.harvard.edu">Mne_analysis@nmr.mgh.harvard.<wbr>edu</a><br>
<a href="https://mail.nmr.mgh.harvard.edu/mailman/listinfo/mne_analysis" rel="noreferrer" target="_blank">https://mail.nmr.mgh.harvard.<wbr>edu/mailman/listinfo/mne_<wbr>analysis</a><br>
<br>
<br>
The information in this e-mail is intended only for the person to whom it is<br>
addressed. If you believe this e-mail was sent to you in error and the e-mail<br>
contains patient information, please contact the Partners Compliance HelpLine at<br>
<a href="http://www.partners.org/complianceline" rel="noreferrer" target="_blank">http://www.partners.org/<wbr>complianceline</a> . If the e-mail was sent to you in error<br>
but does not contain patient information, please contact the sender and properly<br>
dispose of the e-mail.<br>
<br></blockquote></div><br></div>