[FieldTrip] Time frequency analysis using evoked potentials

[Herring, J.D. (Jim)]

Dear Alex,

Change your scaling to ‘maxabs’ (cfg.zlim = ‘maxabs’;) to get a better feel of the scaling of your data. Now you have a scaling from 0 to 1, which, given your data, gives you very saturated images.

Also, I would argue using a relative baseline on a TFR of an ERP (or ERF) is risky as in theory, your baseline should be zero (you even explicitly try to impose that by performing a baseline correction when calculating the ERP). This means that you are dividing your activation period by values close to zero. This will give you arbitrarily large values. An absolute baseline is preferred.

Best,

Jim


On 1 Aug 2017, at 11:52, Alexander Mentink <amentink at chdr.nl<mailto:amentink at chdr.nl>> wrote:

Hi Fieldtrip community,

I am trying to analyse time frequency representations of electrical evoked potentials, but I am not getting any reliable results at all.

Using hanning window (1-30 Hz) or (4-30 Hz), wavelet transform with absolute baselinetype and relative baselinetype gives me four totally different results (see link)

http://tinyurl.com/y7nerh8v

The EEG data is defined 1000 ms pre stimulus and 1200 ms post stimulus, artifacts removed and preprocessed with baseline correction between 200 ms and 0 ms, bandpass filtered between 0.5 - 30 Hz and referenced to earlobes. The timelocked analysis is attached to the dropbox link as well.

The code I use for hanning window is
%%%%%%%%%%%%%%%%%%

cfg = [];
cfg.output       = 'pow';
cfg.channel      = 'Cz';
cfg.method       = 'mtmconvol';
cfg.taper        = 'hanning';
cfg.foi          = 1:0.025:30;  % OR 4:0.025:30
cfg.t_ftimwin    = 3./cfg.foi;
cfg.toi          = -0.5:0.01:0.7;
TFRwave = ft_freqanalysis(cfg, trialdata);

% Visualization
cfg = [];
cfg.baseline     = [-0.5 -0.1];
cfg.baselinetype = 'absolute';
cfg.maskstyle = 'saturation';
cfg.zlim         = [0 1];
cfg.showlabels   = 'yes';
figure
ft_singleplotTFR(cfg, TFRwave)
%%%%%%%%%%%%%%%%%%

And for wavelet transform is
%%%%%%%%%%%%%%%%%%
cfg = [];
cfg.channel    = 'Cz';
cfg.method     = 'wavelet';
cfg.width      = 7;
cfg.output     = 'pow';
cfg.foi        = 1:0.5:30;
cfg.toi        = -0.5:0.05:1;
TFRwave = ft_freqanalysis(cfg, trialdata);

cfg              = [];
cfg.baseline     = [-0.5 -0.1];
cfg.baselinetype = 'relative'; % OR absolute
cfg.maskstyle    = 'saturation';
cfg.zlim         = [0 1];
cfg.channel      = 'Cz';
figure
ft_singleplotTFR(cfg, TFRwave);
%%%%%%%%%%%%%%%%%%

I have tried the tutorials online, but this doesn't give me any better results...

Does anyone know what I am doing wrong?

Thanks many times!

Best regards,

Alex (Master student




Alexander Mentink
Centre for Human Drug Research
Zernikedreef 8
2333 CL Leiden
The Netherlands
T +31 (0) 71 7517166
amentink at chdr.nl<mailto:amentink at chdr.nl>
http://www.chdr.nl<http://www.chdr.nl/>

Follow us:
<imaged4db2c.PNG><http://www.twitter.com/chdr_nl>       <imagecccaa4.PNG><http://www.linkedin.com/company/centre-for-human-drug-research>




_______________________________________________
fieldtrip mailing list
fieldtrip at donders.ru.nl<mailto:fieldtrip at donders.ru.nl>
https://mailman.science.ru.nl/mailman/listinfo/fieldtrip

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