[FieldTrip] ICA - problem with ft_rejectcomponent

[Lei, A. (Anqi)]

Hello,

I'm new to Fieldtrip and I encountered a problem in ICA – after rejecting bad components with ft_rejectcomponent I noticed my data can't be plotted by databrowser anymore, but before this step it can be visualized just fine. This is how I am perform ICA – is there something wrong in the ft_rejectcomponent part? Thank you for your help!




    %% ICA

    % perform the independent component analysis (i.e., decompose the data), downsampling first

    cfg            = [];

    cfg.resamplefs = 150;

    cfg.detrend    = 'no';

    data           = ft_resampledata(cfg, data_all);



    cfg        = [];

    cfg.method = 'runica'; % this is the default and uses the implementation from EEGLAB

    comp = ft_componentanalysis(cfg, data);



    %Identify components reflecting eye artifacts

    % plot the components for visual inspection

    figure

    cfg = [];

    cfg.component = 1:20;       % specify the component(s) that should be plotted

    cfg.layout    = 'acticap-64ch-standard2_XZ.mat'; % specify the layout file that should be used for plotting

    cfg.comment   = 'no';

    ft_topoplotIC(cfg, comp)



    %Further inspection of time course of component

    cfg = [];

    cfg.layout = 'acticap-64ch-standard2_XZ.mat'; % specify the layout file that should be used for plotting

    cfg.viewmode = 'component';

    ft_databrowser(cfg, comp)



    %Removing bad artifacts and backprojecting data

    cfg = [];

    cfg.component = 1; % to be removed component(s)

    data_all = ft_rejectcomponent(cfg, comp, data_all);

Best,
Anqi


-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.science.ru.nl/pipermail/fieldtrip/attachments/20220327/518e2115/attachment.htm>