[FieldTrip] How to obtain the frequency power spectrum of ica components?

Ashley Greene agreene24 at gmail.com
Wed Aug 6 07:12:06 CEST 2014


Hello,

I'm somewhat new to ft, so forgive me if this is a seemingly trivial
question.
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:


comp.cfg.previous.tapsmofrq=4;
comp.cfg.previous.foi=[0:2:30];
comp.cfg.previous.t_ftimwin = 5./comp.cfg.previous.foi;
comp.cfg.previous.toi=[-0.5:0.05:1.5];
comp.cfg.previous.channel='all';

[freqcomp] = ft_freqanalysis(comp.cfg.previous, comp)


// but the errors I receive are as follows:


Warning: Input arguments must be scalar.
> In specest/private/ft_preproc_polyremoval at 57
  In ft_specest_mtmfft at 66
  In ft_freqanalysis at 526
??? For colon operator with char operands, first and last operands must be
char.

Error in ==> ft_preproc_polyremoval at 58
for i = 0:order

Error in ==> ft_specest_mtmfft at 66
  dat = ft_preproc_polyremoval(dat, polyorder, 1, ndatsample);

Error in ==> ft_freqanalysis at 526
        [spectrum,ntaper,foi] = ft_specest_mtmfft(dat, time, 'taper',
cfg.taper, options{:}, 'feedback', fbopt);

// 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?

Thanks,

Ashley
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.science.ru.nl/pipermail/fieldtrip/attachments/20140806/30a89463/attachment-0001.html>


More information about the fieldtrip mailing list