tfanalysis - NaN

Virginie van Wassenhove vvw at CALTECH.EDU
Thu Jun 21 20:44:41 CEST 2007


Hello again.

I am finding NaN whatever method I use (multitaper and wavelet
analysis), yet using a different analysis tool with similar
parameters I obtain clean results.
epoch length    = 1.1 s
baseline length = 0.3 s
sampling rate   = 500Hz

I tried several suggestions listed in the previous message list for a
similar problem, but none seem to provide reasonable results (aka always NaN).
Also tried with different subjects, no particular error otherwise. I
am sending a sample script in case someone has a suggestion on this...
Thanks for your help once more!

Virginie

%% ================================================
% TF analysis - MULTITAPER
cfg             = [];
cfg.method      = 'mtmconvol';
cfg.taper       = 'hanning';
cfg.output      = 'powandcsd';
%cfg.channel    = 'all';
cfg.channel     = 'E75';

%win_length     = 0.128;
%win_n  = floor(2.*(1.1./win_length)-1);

cfg.foi = 1:5:45;
%cfg.t_ftimwin = ones(1,win_n).*win_length;
cfg.t_ftimwin   = 4./cfg.foi;
%cfg.toi        = -0.3:win_length/2:0.8;
cfg.toi         = -0.3:0.05:0.8;

% cfg.pad       = 'maxperlen';
% cfg.keeptrials = 'no';
% cfg.keeptapers ='no';
s08_ill_mult = freqanalysis(cfg,s08_ill);

%% ========================
% plot freq data
cfg             = [];
cfg.baseline    = [-0.3 0];
cfg.layout      = egilay;
cfg.zlim        = [-0.00001 0.000001];          % a comparer entre participants
% warning off                 % NaN
clf
%multiplotTFR(cfg,s08_ill_mult);
singleplotTFR(cfg,s08_ill_mult);

%% ================================================
% TF analysis - wavelet
cfg             = [];
cfg.method      = 'wltconvol';
cfg.width               = 4;
cfg.output       = 'pow';
cfg.foi                 = 1:2:30;
cfg.toi                 = -0.3:0.05:0.8;
s14_ill_wvlt    = freqanalysis(cfg, s14_ill);
%% ========================
cfg = [];
cfg.baseline    = [-0.3 0];
cfg.zlim        = [-1 1];
cfg.showlabels = 'yes';
cfg.layout      = egilay;
clf
multiplotTFR(cfg, s14_ill_wvlt)

%% ================================================






Virginie van Wassenhove, PhD

:::::::::::: contact info  :::::::::::::
Caltech - Division of Biology
1200 E. California Blvd M/C 139-74
Pasadena CA 91125 USA
:::::::::::::::::::::::::::::::::::::::::::::::::
vvw at caltech.edu
Virginie.van.Wassenhove at gmail.com
W: 626.395.8959
http://www.its.caltech.edu/~vvw

:::::::::::::::::: extras ::::::::::::::::::::
http://www.kiva.org
http://www.thehungersite.com/
http://www.agloco.com/r/BBBS1539
:::::::::::::::::::::::::::::::::::::::::::::::::

----------------------------------
The aim of this list is to facilitate the discussion between users of the FieldTrip  toolbox, to share experiences and to discuss new ideas for MEG and EEG analysis. See also http://listserv.surfnet.nl/archives/fieldtrip.html and http://www.ru.nl/fcdonders/fieldtrip.



More information about the fieldtrip mailing list