[FieldTrip] Number of epochs in the ft_freqanalysis output

[Pelt, S. van (Stan)]

Dear Elana,

The 50ms timesteps result in the 41 timebins. Per timebin, the average power spectrum (TFRhann.powspctrm) is computed over all trials that have data in that specific time bin. In your case, the power in the earlier time bins (e.g. 0.0s-0.05s) will be based on more data points (trials) than that of the later ones (e.g. 1.95-2.0s), because of your variable trial length.

If you specify cfg.keeptrials=’yes’, you can see that the individual trial’s TFRs will have different amount of data points.

Hope that helps.

Best,
Stan

--
Stan van Pelt, PhD
Donders Institute for Brain, Cognition and Behaviour
Radboud University
Montessorilaan 3, B.01.34
6525 HR Nijmegen, the Netherlands
tel: +31 24 3616288


From: fieldtrip-bounces at science.ru.nl [mailto:fieldtrip-bounces at science.ru.nl] On Behalf Of Harris, Elana
Sent: maandag 14 december 2015 13:12
To: FieldTrip discussion list
Subject: Re: [FieldTrip] Number of epochs in the ft_freqanalysis output


cfg.toi          = -0:0.05:2;                  % time window "slides" from -0.5 to 1.5 sec in steps of 0.05 sec (50 ms)



I just noticed this line has an error should be -0.5:0.5:2



Sorry, I am new to fieldtrip myself and can't be more helpful yet.



Cheers,



Elana

________________________________
From: fieldtrip-bounces at science.ru.nl<mailto:fieldtrip-bounces at science.ru.nl> [fieldtrip-bounces at science.ru.nl] on behalf of Federica Mauro [federica.ma at gmail.com]
Sent: Monday, December 14, 2015 6:30 AM
To: FieldTrip discussion list
Subject: [FieldTrip] Number of epochs in the ft_freqanalysis output
Dear all,

I'm having a problem in using the ft_freqanalysis function, so I'm writing to ask your help.
I'm trying to get power data for each 2-seconds epochs of my eeg files.
My eeg files are not of the same length (they don't have the same number of epochs), however the output of the function (namely TFRhann.powspctrm) always gives me power data for the same number of epochs, namely 41, unregarding of the effective number of epochs.

Here is the code I'm using:

cfg              = [];
cfg.output       = 'pow';
% cfg.channel
cfg.method       = 'mtmconvol';
cfg.taper        = 'hanning';
cfg.foi          = 4:1:40;                         % analysis 2 to 30 Hz in steps of 2 Hz
cfg.t_ftimwin    = ones(length(cfg.foi),1).*0.5;   % length of time window = 0.5 sec
cfg.toi          = -0:0.05:2;                  % time window "slides" from -0.5 to 1.5 sec in steps of 0.05 sec (50 ms)
TFRhann = ft_freqanalysis(cfg, eeg);

Do you have any suggestion about it?
Thank you very much for your attention.

Worm regards,
Federica


-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.science.ru.nl/pipermail/fieldtrip/attachments/20151214/2a99bf20/attachment-0002.html>