<meta http-equiv="Content-Type" content="text/html; charset=utf-8"><div dir="auto"><div dir="auto">Dear all,</div><div dir="auto"><br></div><div dir="auto">I'm having trouble rejecting components with just computed ICA comp. The data is raw type neuromag data. I have the impression it's a bug, so filed an issue to github <a href="https://github.com/fieldtrip/fieldtrip/issues/980">https://github.com/fieldtrip/fieldtrip/issues/980</a>.</div><div dir="auto"><br></div><div dir="auto">Any help will be much appreciated.</div><div dir="auto"><br></div><div dir="auto">Just after computing components, I try to reject components and get the following error.</div><div dir="auto"><br></div><div dir="auto">cfgica              = [];</div><div dir="auto">cfgica.method       = 'runica';</div><div dir="auto">cfgica.numcomponent = 50;</div><div dir="auto">comp             = ft_componentanalysis(cfgica, data);</div><div dir="auto"><br></div><div dir="auto">cfg                           = [];</div><div dir="auto">cfg.component                 = [1 10];</div><div dir="auto">data_corrected                = ft_rejectcomponent(cfg, comp, data );</div><div dir="auto"><br></div><div dir="auto">undoing the comp balancing for the gradiometer definition</div><div dir="auto">Error using ft_datatype_sens (line 279)</div><div dir="auto">unexpected number of coils (510) contributing to channel MEG0111 (1)</div><div dir="auto">Error in undobalancing (line 58)</div><div dir="auto">sens = ft_datatype_sens(sens);</div><div dir="auto">Error in ft_senstype (line 283)</div><div dir="auto">      sens = undobalancing(sens);</div><div dir="auto">Error in ft_datatype_sens (line 151)</div><div dir="auto">ismeg = ft_senstype(sens, 'meg');</div><div dir="auto">Error in ft_datatype_raw (line 112)</div><div dir="auto">      data.grad = ft_datatype_sens(data.grad);</div><div dir="auto">Error in ft_datatype_comp (line 110)</div><div dir="auto">      raw = ft_datatype_raw(raw, 'version', rawversion, 'hassampleinfo', hassampleinfo, 'hastrialinfo', hastrialinfo);</div><div dir="auto">Error in ft_checkdata (line 266)</div><div dir="auto">  data = ft_datatype_comp(data, 'hassampleinfo', hassampleinfo);</div><div dir="auto">Error in ft_rejectcomponent (line 90)</div><div dir="auto">  comp    = ft_checkdata(comp, 'datatype', 'comp'); </div><div dir="auto"><br></div><div dir="auto">Nota:</div><div dir="auto"><br></div><div dir="auto">data.grad</div><div dir="auto">ans = </div><div dir="auto">  struct with fields:</div><div dir="auto"><br></div><div dir="auto">     balance: [1×1 struct]</div><div dir="auto">     chanori: [306×3 double]</div><div dir="auto">     chanpos: [306×3 double]</div><div dir="auto">    chantype: {306×1 cell}</div><div dir="auto">    chanunit: {306×1 cell}</div><div dir="auto">     coilori: [510×3 double]</div><div dir="auto">     coilpos: [510×3 double]</div><div dir="auto">    coordsys: 'neuromag'</div><div dir="auto">       label: {306×1 cell}</div><div dir="auto">         tra: [306×510 double]</div><div dir="auto">        type: 'neuromag306'</div><div dir="auto">        unit: 'cm'</div><div dir="auto"><br></div><div dir="auto">comp.grad</div><div dir="auto">ans = </div><div dir="auto">  struct with fields:</div><div dir="auto"><br></div><div dir="auto">       balance: [1×1 struct]</div><div dir="auto">       chanori: [50×3 double]</div><div dir="auto">       chanpos: [50×3 double]</div><div dir="auto">       coilori: [510×3 double]</div><div dir="auto">       coilpos: [510×3 double]</div><div dir="auto">      coordsys: 'neuromag'</div><div dir="auto">         label: {50×1 cell}</div><div dir="auto">           tra: [50×510 double]</div><div dir="auto">          unit: 'cm'</div><div dir="auto">    chanposold: [306×3 double]</div><div dir="auto">      labelold: {306×1 cell}</div><div dir="auto">    chanoriold: [306×3 double]</div></div>