<div dir="ltr">Dear Xiaoming,<br><br>> (1) I should feed ICA with 'cleaned' epochs (less non-stereotypic noise) and with enough data points (> k*N^2). <br><br>It does not have to be epochs. You can feed long continuous data if you want.<br>Just in case, (EEG.nbchan^2) * k, where k>30.<br><br>> (2), I can calculate ICA on short epochs, then apply the ICA (weights and sphere matrices) to long epochs from the same dateset.<br><br>Sure.<br><br><span style="font-family:Tahoma;font-size:13.3333px">> I want ICA focus on this short window to identify the artifact-components.</span><br><br>My impression is that feeding all data point is good after all. If the artifact is bad enough, ICA will pick it up. It's generally hard to imagine the situation (just from my experience though, no guarantee) that the artifact is so subtle that it can be ONLY decomposed in your suggested approach... It is because what we see and ICA sees in data are different, and our intuitive understanding of data is often wrong.<div><br></div><div>Makoto<br><div class="gmail_extra"><span style="font-family:Tahoma;font-size:13.3333px"><br></span></div><div class="gmail_extra"><span style="font-family:Tahoma;font-size:13.3333px"><br></span></div><div class="gmail_extra"><span style="font-family:Tahoma;font-size:13.3333px"><br></span></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Jun 7, 2016 at 12:52 PM, Xiaoming Du <span dir="ltr"><<a href="mailto:XDu@mprc.umaryland.edu" target="_blank">XDu@mprc.umaryland.edu</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 style="margin:4px 4px 1px;font-style:normal;font-variant:normal;font-weight:normal;font-stretch:normal;font-size:10pt;line-height:normal;font-family:Tahoma">
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<div>
<div>Hi all,<br><br></div>I have a rookie question about using ICA to remove specific artifacts. Please correct me if I was wrong.<br><br></div> (1) I should feed ICA with 'cleaned' epochs (less non-stereotypic noise) and with enough data points (> k*N^2). <br> <br> (2), I can calculate ICA on short epochs, then apply the ICA (weights and sphere matrices) to long epochs from the same dateset.<br><br></div>(reference link: <a href="http://sccn.ucsd.edu/wiki/Chapter_09:_Decomposing_Data_Using_ICA#Studying_and_removing_ICA_components" target="_blank">http://sccn.ucsd.edu/wiki/Chapter_09:_Decomposing_Data_Using_ICA#Studying_and_removing_ICA_components</a>)<br><br></div>For my data set, I want to remove the known artifacts that occurs right after the events (within 100 ms). Also, there are only 13 channels (including 2 eye-channels) for each subject.<br></div>
<div>My goal is to remove this special artifacts using ICA on long epochs. However, because I already know the artifacts occurs shortly after events, I want ICA focus on this short window to identify the artifact-components. I want to minimize the effects of other artifacts or noise on ICA. <br></div>
<div><br></div>My current thoughts are: 1, cut continuous file into long epochs (4 seconds). <br></div> 2, remove epochs with non-stereotypic noise.<br></div> 3, 1-hz high-pass filter the remaining long epochs, then cut short epochs (-30 to 100 ms). In this way, the short epochs mainly contains the artifacts I want to remove. <br></div> 4, apply ICA on the 0.13-second short epochs. The sampling rate of our data is 1000 Hz. There are more than 50 trials. Therefore, The data point for ICA are around 1000*0.13*50 = 6500. This is larger than 5070 (30*13^2), so I should have enough data for ICA on those short epochs.<br></div> 5, identify the artifact-components from step 4.<br></div> 6, apply the ICA weights and sphere to long epochs (before 1-hz highpass filter), then reject artifact-components identified from step 5. <br></div> 7, for now, I should have long epochs that are artifact-corrected and not filtered. <br><br></div>Please let me know if those steps are reasonable. Any suggestions or comments are appreciated!<br><br></div>Thanks.
<div>
<div><img src="cid:ZVWRZPATLBYK.cleardot.gif"></div><span class=""><font color="#888888">
<div> </div>
<div>Xiaoming</div></font></span></div></div><br>_______________________________________________<br>
Eeglablist page: <a href="http://sccn.ucsd.edu/eeglab/eeglabmail.html" rel="noreferrer" target="_blank">http://sccn.ucsd.edu/eeglab/eeglabmail.html</a><br>
To unsubscribe, send an empty email to <a href="mailto:eeglablist-unsubscribe@sccn.ucsd.edu">eeglablist-unsubscribe@sccn.ucsd.edu</a><br>
For digest mode, send an email with the subject "set digest mime" to <a href="mailto:eeglablist-request@sccn.ucsd.edu">eeglablist-request@sccn.ucsd.edu</a><br></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="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></div>