Dear list subscribers,<div><br></div><div>I am an ERPer looking for the first time into the world of time-frequency analysis. I am trying to perform some multitaper analysis on a data set. Ideally I would like to perform similar analysis as described in a paper by Gould et al (2011) on Journal of Neurophysiology. (i.e. in their own words "Time-frequency analysis was performed on</div>
<div>all epochs using a multitaper approach (2–30 Hz, 0.5-Hz steps, 500-ms temporal smoothing window, 2-Hz frequency resolution)" )</div><div><br></div><div>I have epochs from -1100 ms to 100 time locked to the stimulus onset and I am trying to replicate the analysis they did. The problem is that I got boundary effects on all my conditions, namely I obtain NaN in the extreme 100 ms of the time window for all frequencies and for very low frequencies (e.g. 2-4Hz) I have no data in the whole time window. I thought I had time windows large enough to avoid boundary effects, however there must be something I am missing in my code. </div>
<div><br></div><div>Here's what I did following the tutorial:</div><div><br></div><div><div><div>cfg = [];</div><div>cfg.output     = 'pow';</div><div>cfg.channel    = 'EEG';</div><div>cfg.method     = 'mtmconvol'; </div>
<div>cfg.foi        = 2:2:30;</div><div>cfg.t_ftimwin  = 5./cfg.foi; </div><div>cfg.tapsmofrq  = 0.4*cfg.foi;</div><div>cfg.toi        = -1.1:0.02:0.1;</div><div>cfg.pad        = 'maxperlen'; </div><div>TFRmult = ft_freqanalysis(cfg, data); </div>
</div><div><br></div><div><br></div><div>Does anyone have a clue on this?</div><div><br></div><div><br></div><div>Thank you very much</div><div><br></div><div>Silvia</div>-- <br>__________________________________<br><br>Silvia Pagano PhD student<br>
CIMeC - Center for Mind/Brain Science<br>Corso Bettini, 31 38068 Rovereto (TN) <br>Tel. +39 0464 808709<br><br>
</div>