<html><head><meta http-equiv="Content-Type" content="text/html; charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">I have a comprehensive automated artifact correction routine (MAAC) implemented in my EP Toolkit, which you may find helpful.  It follows the principle that different artifacts have different characteristics and so different algorithms will be most effective for each one.  I’m in the middle of writing it up.<div class=""><br class=""></div><div class=""><a href="https://sourceforge.net/projects/erppcatoolkit/" class="">https://sourceforge.net/projects/erppcatoolkit/</a><br class=""><div class=""><br class=""></div><div class="">Joe<br class=""><div><br class=""><blockquote type="cite" class=""><div class="">On Nov 30, 2018, at 12:08, Jason Taylor <jason.taylor@manchester.ac.uk> wrote:</div><br class="Apple-interchange-newline"><div class=""><div class="">Hi Aitor,<br class=""><br class="">No, not that I know of. I generally use a hacky combination of SPM, fieldtrip, and EEGLAB functions, but if you've already run ICA, you could accomplish what I suggested with some standard matlab functions. <br class=""><br class="">Best wishes,<br class="">Jason<br class=""><br class="">-----Original Message-----<br class="">From: Aitor Egurtzegi [mailto:aitor.martinezegurcegui@uzh.ch] <br class="">Sent: 30 November 2018 13:43<br class="">To: Jason Taylor; FieldTrip discussion list<br class="">Subject: Re: [FieldTrip] automatic IC rejection<br class=""><br class="">Dear Jason,<br class=""><br class="">Thanks a lot for your reply. Is there a Fieldtrip method already <br class="">implemented to run such temporal correlation? or would I have to do the <br class="">implementation in raw Matlab?<br class=""><br class="">Thanks in advance,<br class="">Aitor<br class=""><br class="">On 11/30/18 2:05 AM, Jason Taylor wrote:<br class=""><blockquote type="cite" class="">Hi Aitor,<br class=""><br class="">If you have an 'objective' measure of the artefact you're trying to remove (e.g., VEOG for blinks), a relatively straightforward method is to run a temporal correlation between each IC's activation time-course and the artefact channel's time-course. You can then reject any IC with a correlation higher than some threshold, or with a Z-score (r value relative to the distribution of IC r values) above some threshold. This tends to work very well for identifying blinks, and fairly well for eye-movements (*EOG), and can work for pulse artefact if you have recorded ECG. To avoid spurious correlations due to high-frequency noise, you can filter (e.g., 1 to 30 Hz) the component and artefact signals before correlating them (but obviously go back to the original unfiltered signals to continue with your analysis).<br class=""><br class="">Best wishes,<br class="">Jason<br class=""><br class=""><br class="">-----Original Message-----<br class="">From: fieldtrip [mailto:fieldtrip-bounces@science.ru.nl] On Behalf Of David Schubring<br class="">Sent: 28 November 2018 12:47<br class="">To: fieldtrip@science.ru.nl; aitor.martinezegurcegui@uzh.ch<br class="">Subject: Re: [FieldTrip] automatic IC rejection<br class=""><br class="">Dear Aitor,<br class=""><br class="">the closest thing I know of for a data-driven approach of selecting<br class="">independent components is COMPASS, quote:<br class=""><br class="">"COMPASS is a MATLAB and EEGLAB based algorithm with the purpose of<br class="">providing the user with a convenient technique for automatic Independent<br class="">Component (IC) selection with respect to the contributions of the ICs to<br class="">a certain ERP."<br class=""><br class="">Link to the toolbox:<br class=""><br class="">http://53450283.de.strato-hosting.eu/jrw/lab/e_compass.htm<br class=""><br class="">Paper:<br class=""><br class="">Wessel, J. R., & Ullsperger, M. (2011). Selection of independent<br class="">components representing event-related brain potentials: a data-driven<br class="">approach for greater objectivity. Neuroimage, 54(3), 2105-2115.<br class="">https://doi.org/10.1016/j.neuroimage.2010.10.033<br class=""><br class="">I have only theoretical experience with the toolbox as I only learned<br class="">about it in a workshop and did not yet have the time to test and<br class="">implement it in my personal FieldTrip workflow (even though it is on my<br class="">ever growing to-do list). So far it looked like a useful thing to try<br class="">out to me, especially as code can better be reproduced than "personal<br class="">judgement".<br class=""><br class="">Best,<br class="">David<br class=""><br class="">Am 28.11.2018 um 10:49 schrieb Aitor Egurtzegi:<br class=""><blockquote type="cite" class="">Dear researchers at Fieldtrip,<br class=""><br class=""><br class="">In order to make my work more reproducible, I would like to<br class="">automatically reject ICs instead of doing visual inspection and<br class="">rejection of the components. Unfortunately, I haven't found any<br class="">documentation for such thing. Is there a way to do it in Fieldtrip?<br class=""><br class="">Best,<br class="">Aitor<br class=""><br class="">_______________________________________________<br class="">fieldtrip mailing list<br class="">https://mailman.science.ru.nl/mailman/listinfo/fieldtrip<br class="">https://doi.org/10.1371/journal.pcbi.1002202<br class=""></blockquote></blockquote><br class="">_______________________________________________<br class="">fieldtrip mailing list<br class="">https://mailman.science.ru.nl/mailman/listinfo/fieldtrip<br class="">https://doi.org/10.1371/journal.pcbi.1002202<br class=""></div></div></blockquote></div><br class=""><div class="">
<div style="text-align: start; text-indent: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class=""><div style="text-align: start; text-indent: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class=""><div style="text-align: start; text-indent: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class=""><div style="text-align: start; text-indent: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class=""><div style="color: rgb(0, 0, 0); letter-spacing: normal; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-align: start; text-indent: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">--------------------------------------------------------------------------------<br class=""><br class="">Joseph Dien, PhD</div><div style="color: rgb(0, 0, 0); letter-spacing: normal; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-align: start; text-indent: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">Senior Research Scientist</div><div style="text-align: start; text-indent: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">Department of Human Development and Quantitative Methodology</div><div style="color: rgb(0, 0, 0); letter-spacing: normal; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-align: start; text-indent: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">University of Maryland, College Park</div><div style="color: rgb(0, 0, 0); letter-spacing: normal; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-align: start; text-indent: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">http://joedien.com</div></div></div></div></div>
</div>
<br class=""></div></div></body></html>