[Eeglablist] Error in ICLabel feature extraction

[Arnaud Delorme]

Dear Hamza,

Would you creating an issue and uploading a dataset. It might be because you have a non integer sampling rate.

https://urldefense.com/v3/__https://github.com/sccn/eeglab/issues__;!!Mih3wA!BJfnOG8fjpjavceKlwcjIsRp4ekjOoNjKC7QkkBADqi328I8-wbNeoRo_uaQqNnhR8OVQDR2abzmPh24R67GiH5u$ 

Arno

> On Jun 25, 2023, at 11:24 PM, Hamza Imtiaz via eeglablist <eeglablist at sccn.ucsd.edu> wrote:
> 
> Hi,
> 
> I am encountering an error while using the ICLabel function in EEGLab and would greatly appreciate any guidance. Automated ICA rejection is part of the preprocessing pipeline, which includes resampling > filtering > clean_artifacts > re-referencing > make the data full ranked > ICA > ICAlabel. All of this is done in the EEGLAB environment and the filetype is .edf. I've tried updating EEGLAB to the latest version (2023.0) and have reinstalled the plugin as well. Attached below are the error message and the code snippet.
> 
> Could you please provide some insights on the possible cause of this error and any recommended steps to resolve it? Thank you very much for your time and assistance.
> 
> Error message:
> Index in position 1 exceeds array bounds. Index must not exceed 101.
> Error in signal.internal.resample.uniformResample (line 135)
>    yM=ymat(indM,:);
> Error in resample (line 198)
>            signal.internal.resample.uniformResample(xIn, isDimValSet, Dim, dimIn, numericArgs{:});
> Error in eeg_autocorr_welch (line 48)
>    ac = resample(double(ac'), 100, EEG.srate)';
> Error in ICL_feature_extractor (line 77)
>            autocorr = eeg_autocorr_welch(EEG);
> Error in iclabel (line 52)
> features = ICL_feature_extractor(EEG, flag_autocorr);
> Error in pop_iclabel (line 109)
>    EEG = iclabel(EEG, icversion);
> 
> MATLAB code snippet:
>    % Make the data full-ranked is to discard (any) one channel after average referencing. If you don't do this, ICA may return invalid results.
>    EEG.data = bsxfun( @minus, EEG.data, sum( EEG.data, 1 ) / ( EEG.nbchan + 1 ) );
> 
>    % Run ICA decomposition
>    EEG = pop_runica(EEG, 'icatype', 'runica', 'extended', 1, 'interrupt', 'on');
> 
>    % Assign class labels to the ICs
>    EEG = pop_iclabel(EEG, 'default');
> 
> Kind regards,
> Hamza Imtiaz
> PhD Student at the Friedman lab (https://urldefense.com/v3/__https://www.bbbscience.com/__;!!Mih3wA!DIfUcNB8PYd4DZqKipa4kHpCrLSmL57qDS8SWisKkpvnsnr4CgojuEVI0mr8xzd6IarQ2Kln3plOMFtGI_nisRU$ )
> Dalhousie University
> _______________________________________________
> Eeglablist page: http://sccn.ucsd.edu/eeglab/eeglabmail.html
> To unsubscribe, send an empty email to eeglablist-unsubscribe at sccn.ucsd.edu
> For digest mode, send an email with the subject "set digest mime" to eeglablist-request at sccn.ucsd.edu