<div dir="ltr"><div class="gmail_default" style="color:#333399">Hello Kelly, curious if you got some resolution from eeglablist about ASR's dynamics ?</div></div><br><div class="gmail_quote"><div dir="ltr">On Wed, Sep 12, 2018 at 11:48 AM Kelly Michaelis <<a href="mailto:kcmichaelis@gmail.com">kcmichaelis@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr">Hi Everyone,<div><br></div><div>I am using ASR to clean my datasets prior to ICA. I've been trying to figure out what level to set each of the ASR parameters to, and in doing so I noticed that the same input is producing massively different results each time I run it. </div><div><br></div><div>More specifically, here's what I did:</div><div><br></div><div>1. filter, resample</div><div>2. Manually remove bad channels (I noticed ASR wasnt' catching some)</div><div>3. ASR using the following parameters: </div><div><br></div><div> EEG = clean_rawdata(EEG, -1, -1, -1, -1, 24, .25);<br></div><div><br></div><div>4. Reference to average</div><div>5. Epoch</div><div><br></div><div>After I noticed some inconsistencies, I decided to run this exact series of steps several times on the same dataset (i.e. starting with the raw continuous dataset each time). On the first run, ASR inserted so many boundary events that I could't get any full epochs out. On the second run, I got 161 epochs (out of 200). On the third run, I it was back to zero full epochs. On the fourth run, I got 162 epochs.</div><div><br></div><div>Does anyone else have this issue?? I want to apply the same settings to multiple subject files, but if the results are this inconsistent that's a huge issue. I'll put the command window output below. It looks like on the runs where it inserts lots of boundary events it finishes running much faster and is maybe not doing the cleaning?</div><div><br></div><div>Please help! I want to be able to use ASR.</div><div><br></div><div>Thanks,</div><div>Kelly</div><div><br></div><div><b>From the run that produced no viable epochs: </b></div><div><br></div><div><div>>> EEG = clean_rawdata(EEG, -1, -1, -1, -1, 24, .25);</div><div><br></div><div>flatchan rej disabled.</div><div>highpass disabled.</div><div>badchan rej disabled.</div><div>noise-based rej disabled.</div><div>Finding a clean section of the data...</div><div>Determining time window rejection thresholds...done.</div><div>Keeping 21.1% (475 seconds) of the data.</div><div>eeg_insertbound(): 477 boundary (break) events added.</div><div>eeg_insertbound(): event latencies recomputed and 1682 events removed.</div><div>eeg_checkset note: upper time limit (xmax) adjusted so (xmax-xmin)*srate+1 = number of frames</div><div>eeg_checkset note: value format of event field 'stimname' made uniform</div><div>eeg_checkset note: value format of event field 'answer' made uniform</div><div>eeg_insertbound(): 477 boundary (break) events added.</div><div>eeg_insertbound(): event latencies recomputed and 1686 events removed.</div><div>BUG 1971 WARNING: IF YOU ARE USING A SCRIPT WITTEN FOR A PREVIOUS VERSION OF</div><div>EEGLAB TO CALL THIS FUNCTION, BECAUSE YOU ARE REJECTING THE ONSET OF THE DATA,</div><div>EVENTS WERE CORRUPTED. EVENT LATENCIES ARE NOW CORRECT (SEE <a href="https://sccn.ucsd.edu/wiki/EEGLAB_bug1971" target="_blank">https://sccn.ucsd.edu/wiki/EEGLAB_bug1971</a>);</div><div>Estimating calibration statistics; this may take a while...</div><div>Warning: A value of class "com.sun.management.UnixOperatingSystem" was indexed with no subscripts specified. Currently the result of</div><div>this operation is the indexed value itself, but in a future release, it will be an error. </div><div>> In asr_calibrate>hlp_memfree (line 266)</div><div>  In asr_calibrate (line 114)</div><div>  In clean_asr (line 164)</div><div>  In clean_artifacts (line 219)</div><div>  In clean_rawdata (line 83) </div><div>Determining per-component thresholds...done.</div><div>Warning: A value of class "com.sun.management.UnixOperatingSystem" was indexed with no subscripts specified. Currently the result of</div><div>this operation is the indexed value itself, but in a future release, it will be an error. </div><div>> In asr_process>hlp_memfree (line 232)</div><div>  In asr_process (line 105)</div><div>  In clean_asr (line 174)</div><div>  In clean_artifacts (line 219)</div><div>  In clean_rawdata (line 83) </div><div>Now doing final post-cleanup of the output.</div><div>Determining time window rejection thresholds...done.</div><div>Keeping 28.5% (642 seconds) of the data.</div><div>eeg_insertbound(): 587 boundary (break) events added.</div><div>eeg_insertbound(): event latencies recomputed and 1595 events removed.</div><div>eeg_checkset note: upper time limit (xmax) adjusted so (xmax-xmin)*srate+1 = number of frames</div><div>eeg_checkset note: value format of event field 'stimname' made uniform</div><div>eeg_checkset note: value format of event field 'answer' made uniform</div><div>eeg_insertbound(): 587 boundary (break) events added.</div><div>eeg_insertbound(): event latencies recomputed and 1599 events removed.</div><div>BUG 1971 WARNING: IF YOU ARE USING A SCRIPT WITTEN FOR A PREVIOUS VERSION OF</div><div>EEGLAB TO CALL THIS FUNCTION, BECAUSE YOU ARE REJECTING THE ONSET OF THE DATA,</div><div>EVENTS WERE CORRUPTED. EVENT LATENCIES ARE NOW CORRECT (SEE <a href="https://sccn.ucsd.edu/wiki/EEGLAB_bug1971" target="_blank">https://sccn.ucsd.edu/wiki/EEGLAB_bug1971</a>);</div><div>Use vis_artifacts to compare the cleaned data to the original.</div><div>Re-referencing data</div><div>Removing 1 channel(s)...</div><div>Saving dataset...</div></div><div><br></div><div><br></div><div><b>From the run that produced 162 epochs:</b></div><div><br></div><div><div>EEG = clean_rawdata(EEG, -1, -1, -1, -1, 24, .25);</div><div><br></div><div>flatchan rej disabled.</div><div>highpass disabled.</div><div>badchan rej disabled.</div><div>noise-based rej disabled.</div><div>Finding a clean section of the data...</div><div>Determining time window rejection thresholds...done.</div><div>Keeping 21.1% (475 seconds) of the data.</div><div>eeg_insertbound(): 477 boundary (break) events added.</div><div>eeg_insertbound(): event latencies recomputed and 1682 events removed.</div><div>eeg_checkset note: upper time limit (xmax) adjusted so (xmax-xmin)*srate+1 = number of frames</div><div>eeg_checkset note: value format of event field 'stimname' made uniform</div><div>eeg_checkset note: value format of event field 'answer' made uniform</div><div>eeg_insertbound(): 477 boundary (break) events added.</div><div>eeg_insertbound(): event latencies recomputed and 1686 events removed.</div><div>BUG 1971 WARNING: IF YOU ARE USING A SCRIPT WITTEN FOR A PREVIOUS VERSION OF</div><div>EEGLAB TO CALL THIS FUNCTION, BECAUSE YOU ARE REJECTING THE ONSET OF THE DATA,</div><div>EVENTS WERE CORRUPTED. EVENT LATENCIES ARE NOW CORRECT (SEE <a href="https://sccn.ucsd.edu/wiki/EEGLAB_bug1971" target="_blank">https://sccn.ucsd.edu/wiki/EEGLAB_bug1971</a>);</div><div>Estimating calibration statistics; this may take a while...</div><div>Warning: A value of class "com.sun.management.UnixOperatingSystem" was indexed with no subscripts specified. Currently the result of</div><div>this operation is the indexed value itself, but in a future release, it will be an error. </div><div>> In asr_calibrate>hlp_memfree (line 266)</div><div>  In asr_calibrate (line 114)</div><div>  In clean_asr (line 164)</div><div>  In clean_artifacts (line 219)</div><div>  In clean_rawdata (line 83) </div><div>Determining per-component thresholds...done.</div><div>Warning: A value of class "com.sun.management.UnixOperatingSystem" was indexed with no subscripts specified. Currently the result of</div><div>this operation is the indexed value itself, but in a future release, it will be an error. </div><div>> In asr_process>hlp_memfree (line 232)</div><div>  In asr_process (line 105)</div><div>  In clean_asr (line 174)</div><div>  In clean_artifacts (line 219)</div><div>  In clean_rawdata (line 83) </div><div><b>Now cleaning data in 1268 blocks.................................................................................................................................................................................................................................................................................................................................................................................................................................................................................................................................................................................................................................................................................................................................................................................................................................................................................................................................................................................................................................................................................................................................................................................................................................................................................................................... (NOTE: missing from previous run)</b></div><div>Now doing final post-cleanup of the output.</div><div>Determining time window rejection thresholds...done.</div><div>Keeping 90.5% (2035 seconds) of the data.</div><div>eeg_insertbound(): 83 boundary (break) events added.</div><div>eeg_insertbound(): event latencies recomputed and 132 events removed.</div><div>eeg_checkset note: value format of event field 'stimname' made uniform</div><div>eeg_checkset note: value format of event field 'answer' made uniform</div><div>eeg_insertbound(): 83 boundary (break) events added.</div><div>eeg_insertbound(): event latencies recomputed and 132 events removed.</div><div>BUG 1971 WARNING: IF YOU ARE USING A SCRIPT WITTEN FOR A PREVIOUS VERSION OF</div><div>EEGLAB TO CALL THIS FUNCTION, BECAUSE YOU ARE REJECTING THE ONSET OF THE DATA,</div><div>EVENTS WERE CORRUPTED. EVENT LATENCIES ARE NOW CORRECT (SEE <a href="https://sccn.ucsd.edu/wiki/EEGLAB_bug1971" target="_blank">https://sccn.ucsd.edu/wiki/EEGLAB_bug1971</a>);</div><div>Use vis_artifacts to compare the cleaned data to the original.</div><div>Saving dataset...</div></div></div></div></div></div>
_______________________________________________<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" target="_blank">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" target="_blank">eeglablist-request@sccn.ucsd.edu</a></blockquote></div>