<div dir="ltr">Dear Jumana,<div><br></div><div>When you do the bandpass at once, the filter order for 0.1Hz high-pass will be applied to your low-pass, which should be really excessive. It's better to do high-pass and low-pass separately. And I still don't like 0.1Hz high-pass if you use FIR (and I do not know how bad it is to use IIR... I've heard it can become 'unstable' but I've never seen it myself)</div><div><br></div><div>> I won't be using granger causality but I will be estimating phase during ITC. <br></div><div> </div><div>Should be ok.</div><div><br></div><div>Makoto</div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Jan 26, 2017 at 2:12 PM, Ahmad, Jumana <span dir="ltr"><<a href="mailto:jumana.ahmad@kcl.ac.uk" target="_blank">jumana.ahmad@kcl.ac.uk</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">
<div>
Dear Makoto,<br>
I actually switched to the pop eeg filt eeglab function and it now Really attenuated anything >40Hz, and my ERPs are cleaner. However I filtering between 0.1-40Hz at the same time in the GUI (I interned the high pass and low pass simultaneously). Is this OK
to do? The frequency response looks OK. <br>
<br>
The filter order was automatically set very high by the GUI, but it's continuous data and I have room without events at the beginning and end of the data so any edge effects can be disgusted. What do you think?<br>
<br>
Also, this is for my ERP analysis- I trained ICA on a 1Hz high pass filtered set.
<br>
I won't be using granger causality but I will be estimating phase during ITC. <br>
Best wishes, and thanks,<span class="gmail-"><br>
Jumana <br>
<br>
------------------------------<wbr>------------<br>
Jumana Ahmad<br>
Post-Doctoral Research Worker in Cognitive Neuroscience<br>
EU-AIMS Longitudinal European Autism Project (LEAP) & SynaG Study<br></span>
Room M1.26.Department of Forensic and Neurodevelopmental Sciences (PO 23) | Institute of Psychiatry, Psychology & Neuroscience | King’s College London | 16 De Crespigny Park | London SE5 8AF<br>
<br>
Phone: 0207 848 5359| Email: <a href="mailto:jumana.ahmad@kcl.ac.uk" target="_blank">jumana.ahmad@kcl.ac.uk</a> | Website: <a href="http://www.eu-aims.eu" target="_blank">www.eu-aims.eu</a> | Facebook: <a href="http://www.facebook.com/euaims" target="_blank">www.facebook.com/euaims</a><br>
<br>
<hr style="display:inline-block;width:98%">
<div id="gmail-m_-4031508253012889712divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" style="font-size:11pt" color="#000000"><b>From:</b> Makoto Miyakoshi <<a href="mailto:mmiyakoshi@ucsd.edu" target="_blank">mmiyakoshi@ucsd.edu</a>><br>
<b>Sent:</b> 26 January 2017 23:49:50<br>
<b>To:</b> Ahmad, Jumana<br>
<b>Cc:</b> <a href="mailto:eeglablist@sccn.ucsd.edu" target="_blank">eeglablist@sccn.ucsd.edu</a><br>
<b>Subject:</b> Re: [Eeglablist] high frequency oscillation- eeg advice</font>
<div> </div>
</div><div><div class="gmail-h5">
<div>
<div dir="ltr">Dear Jumana,
<div><br>
</div>
<div>It's a bad idea to perform ICA with 0.1Hz high-pass filtered data. The cutoff frequency is too low. See this page and the referenced paper.</div>
<div><br>
</div>
<div><a href="https://emea01.safelinks.protection.outlook.com/?url=https%3A%2F%2Fsccn.ucsd.edu%2Fwiki%2FMakoto%2527s_preprocessing_pipeline%23High-pass_filter_the_data_at_1-Hz_.28for_ICA.2C_ASR.2C_and_CleanLine.29&data=01%7C01%7Cjumana.ahmad%40kcl.ac.uk%7C35190d19320946ffd52408d446355d20%7C8370cf1416f34c16b83c724071654356%7C0&sdata=XQPtEZ22iFzlB2PVii2GjBWVdfPOZZnmnmOfGa8CbgI%3D&reserved=0" target="_blank">https://sccn.ucsd.edu/wiki/<wbr>Makoto%27s_preprocessing_<wbr>pipeline#High-pass_filter_the_<wbr>data_at_1-Hz_.28for_ICA.2C_<wbr>ASR.2C_and_CleanLine.29</a><br>
</div>
<div><br>
</div>
<div>> A 30Hz low pass does not help to get rid of the oscillation, which is really significant in the data.<br>
</div>
<div class="gmail_extra"><br>
</div>
<div class="gmail_extra">Check the channel frequency spectra and tell me if you see peaks in it. If necessary, you can cut it off using a designed low-pass filter (not like Butterworth...)</div>
<div class="gmail_extra"><br>
</div>
<div class="gmail_extra">> I use a butterworth filter, which is good for ERP analysis with low phase distortion.<br>
</div>
<div class="gmail_extra"><br>
</div>
<div class="gmail_extra">Do not make qualitative judgement just because something is NOT a classic Butterworth. Of course, if the attenuation is small, the phase 'distortion' is small. But if such small attenuation is not useful, it does not help at all! Also,
be careful with the word 'phase'. Particularly people who do not know basics of signal processing believe phase as some magical thing. If you are not performing Granger Causality Analysis or something, you don't need to be so worried about phase issue in practice.</div>
<div class="gmail_extra"><br>
</div>
<div class="gmail_extra">
<p class="MsoNormal">> I also already run ICA, but in some datasets there is a very significant high frequency oscillation.<u></u><u></u></p>
<div><br>
</div>
<div>Remember, to eliminate this is more important than being afraid of qualitative phase issue.</div>
</div>
<div class="gmail_extra"><br>
</div>
<div class="gmail_extra">
<p class="MsoNormal">> However, I can see the high frequency oscillations in my ERP, which is not ideal and now I need to try and get rid of it further.<u></u><u></u></p>
<p class="MsoNormal">Can I filter again on top of the data which already has already undergone ICA- I only use ICA to remove blinks?</p>
</div>
<div class="gmail_extra"><br>
</div>
<div class="gmail_extra">You'd better to filter the data on continuous state. If you need to filter the epoched data, the half of filter length from both ends becomes unreliable.</div>
<div class="gmail_extra"><br>
</div>
<div class="gmail_extra">
<p class="MsoNormal">> Should I do cleanline, although it would have to be after ICA now- I read this is not advisable.<u></u><u></u></p>
<p class="MsoNormal">> Should I use a notch filter?</p>
</div>
<div class="gmail_extra"><br>
</div>
<div class="gmail_extra">If you see > 20dB line noise, Cleanline may not help. In this case, I would simply apply a designed low-pass filter, either Hamming (-50dB) or Blackman (-70dB) using firfilt(). See 'Tools' -> 'Filter the data' -> 'Windowed sinc FIR
filter'.</div>
<div class="gmail_extra"><br>
</div>
<div class="gmail_extra">There are different guys saying different things about data preprocessing. It is confusing, I know! The only good solution for this is to become an engineer yourself...</div>
<div class="gmail_extra"><br>
</div>
<div class="gmail_extra">Makoto</div>
<div class="gmail_extra"><br>
</div>
<div class="gmail_extra"><br>
</div>
<div class="gmail_extra"><br>
<div class="gmail_quote">On Tue, Jan 24, 2017 at 7:00 AM, Ahmad, Jumana <span dir="ltr">
<<a href="mailto:jumana.ahmad@kcl.ac.uk" target="_blank">jumana.ahmad@kcl.ac.uk</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">
<div lang="EN-GB">
<div class="gmail-m_-4031508253012889712gmail-m_1381307076216940724WordSection1">
<p class="MsoNormal">Hi Everyone,<u></u><u></u></p>
<p class="MsoNormal">I am running a large scale ERP analysis. I filtered 1-40Hz (ICA AMICA), or 0.1-40Hz for the ERP dataset. A 30Hz low pass does not help to get rid of the oscillation, which is really significant in the data. I use a butterworth filter, which
is good for ERP analysis with low phase distortion. <u></u><u></u></p>
<p class="MsoNormal">I also already run ICA, but in some datasets there is a very significant high frequency oscillation.<u></u><u></u></p>
<p class="MsoNormal">I do not use cleanline, which is not typical in the literature I have been basing my pipeline on.<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">However, I can see the high frequency oscillations in my ERP, which is not ideal and now I need to try and get rid of it further.
<u></u><u></u></p>
<p class="MsoNormal">Can I filter again on top of the data which already has already undergone ICA- I only use ICA to remove blinks?<u></u><u></u></p>
<p class="MsoNormal">Should I do cleanline, although it would have to be after ICA now- I read this is not advisable.<u></u><u></u></p>
<p class="MsoNormal">Should I use a notch filter?<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">Any help would be appreciated.<u></u><u></u></p>
<p class="MsoNormal">Best wishes,<u></u><u></u></p>
<p class="MsoNormal">Jumana <u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal" style="background-image:initial;background-color:white;background-size:initial;background-origin:initial;background-clip:initial;background-position:initial;background-repeat:initial">
<b><span style="font-size:10pt;color:rgb(31,73,125)">------------------------------<wbr>------------</span></b><span style="font-family:helvetica,sans-serif;color:black"><u></u><u></u></span></p>
<p class="MsoNormal" style="background-image:initial;background-color:white;background-size:initial;background-origin:initial;background-clip:initial;background-position:initial;background-repeat:initial">
<b><span style="font-size:10pt;color:rgb(31,73,125);background-image:initial;background-color:rgb(255,238,148);background-size:initial;background-origin:initial;background-clip:initial;background-position:initial;background-repeat:initial">Jumana Ahmad</span></b><span style="font-family:helvetica,sans-serif;color:black"><u></u><u></u></span></p>
<p class="MsoNormal" style="background-image:initial;background-color:white;background-size:initial;background-origin:initial;background-clip:initial;background-position:initial;background-repeat:initial">
<span style="font-size:10pt;color:rgb(31,73,125)">Post-Doctoral Research Worker in Cognitive Neuroscience </span><span style="font-family:helvetica,sans-serif;color:black"><u></u><u></u></span></p>
<p class="MsoNormal" style="background-image:initial;background-color:white;background-size:initial;background-origin:initial;background-clip:initial;background-position:initial;background-repeat:initial">
<i><span style="font-size:10pt;color:rgb(31,73,125)">EU-AIMS Longitudinal European Autism Project (LEAP) & SynaG Study</span></i><span style="font-family:helvetica,sans-serif;color:black"><u></u><u></u></span></p>
<p class="MsoNormal" style="background-image:initial;background-color:white;background-size:initial;background-origin:initial;background-clip:initial;background-position:initial;background-repeat:initial">
<span style="font-size:10pt;color:rgb(31,73,125)">Room M1.09. Department of Forensic and Neurodevelopmental Sciences (PO 23) | Institute of Psychiatry, Psychology & Neuroscience | King’s College London | 16 De Crespigny Park | London SE5 8AF</span><span style="font-family:helvetica,sans-serif;color:black"><u></u><u></u></span></p>
<p class="MsoNormal" style="background-image:initial;background-color:white;background-size:initial;background-origin:initial;background-clip:initial;background-position:initial;background-repeat:initial">
<span style="font-size:10pt;color:rgb(31,73,125)"> </span><span style="font-family:helvetica,sans-serif;color:black"><u></u><u></u></span></p>
<p class="MsoNormal" style="background-image:initial;background-color:white;background-size:initial;background-origin:initial;background-clip:initial;background-position:initial;background-repeat:initial">
<b><span style="font-size:10pt;color:rgb(31,73,125)">Phone:</span></b><span style="font-size:10pt;color:rgb(31,73,125)"> </span><span lang="EN-US" style="font-size:10pt;color:navy">0207848 0260</span><span style="font-size:10pt;color:rgb(31,73,125)">| <b>Email:</b> <span style="background-image:initial;background-color:rgb(255,238,148);background-size:initial;background-origin:initial;background-clip:initial;background-position:initial;background-repeat:initial">jumana.ahmad</span></span><span style="color:black"><a href="mailto:antonia.sanjose@kcl.ac.uk" target="_blank"><span style="font-size:10pt;color:purple">@kcl.<wbr>ac.uk</span></a></span><span style="font-size:10pt;color:rgb(31,73,125)"> | <b>Website:</b> </span><span style="color:black"><a href="https://emea01.safelinks.protection.outlook.com/?url=http%3A%2F%2Fwww.eu-aims.eu%2F&data=01%7C01%7Cjumana.ahmad%40kcl.ac.uk%7C35190d19320946ffd52408d446355d20%7C8370cf1416f34c16b83c724071654356%7C0&sdata=vVSbBawKS2L6e9vVZXkIMrj%2B%2FwqhrreJRZRNr6zUN2s%3D&reserved=0" target="_blank"><span style="font-size:10pt;color:purple">www.eu-aims.e<wbr>u</span></a></span><span style="font-size:10pt;color:rgb(31,73,125)"> | <b>Facebook:</b> </span><span style="color:black"><a href="https://emea01.safelinks.protection.outlook.com/?url=http%3A%2F%2Fwww.facebook.com%2Feuaims&data=01%7C01%7Cjumana.ahmad%40kcl.ac.uk%7C35190d19320946ffd52408d446355d20%7C8370cf1416f34c16b83c724071654356%7C0&sdata=pOPzKLDwxDuoGPJKPTlPGwurARTFmxrwy2VEjUlBtP8%3D&reserved=0" target="_blank"><span style="font-size:10pt;color:purple">www.facebook.com<wbr>/euaims</span></a></span><span style="font-family:helvetica,sans-serif;color:black"><u></u><u></u></span></p>
<p class="MsoNormal"><span><u></u> <u></u></span></p>
<p class="MsoNormal"><span style="font-size:10pt;color:rgb(31,73,125)">******************************<wbr>******************************<wbr>***************************<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:10pt;color:rgb(31,73,125)">******************************<wbr>******************************<wbr>***************************<u></u><u></u></span></p>
<p class="MsoNormal"><b><span style="font-size:10pt;color:rgb(31,73,125)">We are currently looking for volunteers with mild intellectual disability to be part of our exciting and world-leading European project into brain development and social behaviour. Please,
do get in touch if you know of anyone who may be interested in taking part. </span>
</b><span style="font-size:10pt;color:rgb(31,73,125)"><u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:10pt;color:rgb(31,73,125)">******************************<wbr>******************************<wbr>***************************<u></u><u></u></span></p>
<p class="MsoNormal"><span style="font-size:10pt;color:rgb(31,73,125)">******************************<wbr>******************************<wbr>***************************<u></u><u></u></span></p>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
</div>
<br>
______________________________<wbr>_________________<br>
Eeglablist page: <a href="https://emea01.safelinks.protection.outlook.com/?url=http%3A%2F%2Fsccn.ucsd.edu%2Feeglab%2Feeglabmail.html&data=01%7C01%7Cjumana.ahmad%40kcl.ac.uk%7C35190d19320946ffd52408d446355d20%7C8370cf1416f34c16b83c724071654356%7C0&sdata=33Izy%2B79hiSorIQnQtazMnLGL2LcChe7GLi5GrxVclo%3D&reserved=0" rel="noreferrer" target="_blank">
http://sccn.ucsd.edu/eeglab/ee<wbr>glabmail.html</a><br>
To unsubscribe, send an empty email to <a href="mailto:eeglablist-unsubscribe@sccn.ucsd.edu" target="_blank">
eeglablist-unsubscribe@sccn.uc<wbr>sd.edu</a><br>
For digest mode, send an email with the subject "set digest mime" to <a href="mailto:eeglablist-request@sccn.ucsd.edu" target="_blank">
eeglablist-request@sccn.ucsd.e<wbr>du</a><br>
</blockquote>
</div>
<br>
<br clear="all">
<div><br>
</div>
-- <br>
<div class="gmail-m_-4031508253012889712gmail_signature">
<div dir="ltr">Makoto Miyakoshi<br>
Swartz Center for Computational Neuroscience<br>
Institute for Neural Computation, University of California San Diego<br>
</div>
</div>
</div>
</div>
</div>
</div></div></div>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature"><div dir="ltr">Makoto Miyakoshi<br>Swartz Center for Computational Neuroscience<br>Institute for Neural Computation, University of California San Diego<br></div></div>
</div></div>