<div dir="ltr"><div><div><div><div>Hello,<br></div><div><br>I'm somewhat new to ft, so forgive me if this is a seemingly trivial question.<br>I ran runica on ft preprocessed MEG data, but I also want to see the frequency power spectrum of the components. I've tried using freqanalysis on the components in the following manner:<br>
<br><br>comp.cfg.previous.tapsmofrq=4; <br>comp.cfg.previous.foi=[0:2:30]; <br>comp.cfg.previous.t_ftimwin = 5./comp.cfg.previous.foi; <br>comp.cfg.previous.toi=[-0.5:0.05:1.5];<br>comp.cfg.previous.channel='all';<br>
<br>[freqcomp] = ft_freqanalysis(comp.cfg.previous, comp)<br><br></div><br>// but the errors I receive are as follows:<br><br><br>Warning: Input arguments must be scalar. <br>> In specest/private/ft_preproc_polyremoval at 57<br>
  In ft_specest_mtmfft at 66<br>  In ft_freqanalysis at 526<br>??? For colon operator with char operands, first and last operands must be char.<br><br>Error in ==> ft_preproc_polyremoval at 58<br>for i = 0:order<br><br>
Error in ==> ft_specest_mtmfft at 66<br>  dat = ft_preproc_polyremoval(dat, polyorder, 1, ndatsample);<br><br>Error in ==> ft_freqanalysis at 526<br>        [spectrum,ntaper,foi] = ft_specest_mtmfft(dat, time, 'taper', cfg.taper, options{:}, 'feedback', fbopt);<br>
<br></div>// Is using freqanalysis on the components possible? If so, what needs to be modified in order to employ it? If not, is there another method to outputting the frequency power spectrum of ica components?<br><br></div>
Thanks,<br><br></div>Ashley<br><div><div><div><div><br><br></div></div></div></div></div>