<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
</head>
<body>
<div dir="ltr">
<div>
<div>
<div>Dear all,<br>
<br>
</div>
I am trying to look at ERPs stats with my EGI cap (I have 16 subjects - within subjects design).<br>
</div>
According to the code I use, I get in the order, an ERROR or I see something weird and I wish to have an opinion.
<br>
</div>
<div>See the script below:<br>
</div>
<div><br>
</div>
Here is my script:<br>
% ERPs grandaverage calculation for 4 conditions<br>
cfg = [];<br>
cfg.keepindividual = 'yes';<br>
ERPs_FP_grandavg = ft_timelockgrandaverage(cfg, FP_block{:});<br>
ERPs_FS_grandavg = ft_timelockgrandaverage(cfg, FS_block{:});<br>
ERPs_HP_grandavg = ft_timelockgrandaverage(cfg, HP_block{:});<br>
ERPs_HS_grandavg = ft_timelockgrandaverage(cfg, HS_block{:});<br clear="all">
<div>
<div>
<div>
<div><br>
% t-tests<br>
cfg = [];<br>
cfg.method = 'triangulation';<br>
cfg.layout = 'GSN-HydroCel-128.sfp';<br>
cfg.neighbourdist = 2;<br>
cfg.senstype = 'EEG'; <br>
neighbours_EEG = ft_prepare_neighbours(cfg, ERPs_FP_grandavg); <br>
<br>
cfg = [];<br>
cfg.channel = 'EEG';<br>
cfg.minnbchan = 2; <br>
cfg.neighbours = neighbours_EEG; <br>
cfg.latency = [1.17 1.25];<br>
cfg.avgovertime = 'no';<br>
cfg.parameter = 'avg';<br>
cfg.method = 'montecarlo';<br>
cfg.statistic = 'ft_statfun_depsamplesT';<br>
cfg.alpha = 0.05;<br>
cfg.clusteralpha = 0.05; <br>
cfg.correctm = 'cluster';<br>
cfg.correcttail = 'prob';<br>
cfg.numrandomization = 1000;<br>
cfg.tail = 0;<br>
cfg.clustertail = 0;<br>
<br>
Nsub = length(names);<br>
cfg.design(1,1:2*Nsub) = [ones(1,Nsub) 2*ones(1,Nsub)];<br>
cfg.design(2,1:2*Nsub) = [1:Nsub 1:Nsub];<br>
cfg.ivar = 1; % the 1st row in cfg.design contains the independent variable<br>
cfg.uvar = 2; % the 2nd row in cfg.design contains the subject number<br>
<br>
IF I USE THIS CODE:<br>
stat = ft_timelockstatistics(cfg,ERPs_FP_grandavg,ERPs_FS_grandavg);<br>
<br>
</div>
<div>I GET TE ERROR:<br>
Reference to non-existent field 'dat'.<br>
<br>
Error in prepare_timefreq_data>forcedimord (line 531)<br>
Nrepl = size(output.dat, repldim);<br>
<br>
Error in prepare_timefreq_data (line 87)<br>
[remember{c}, hascrsspctrm] = forcedimord(varargin{c});<br>
<br>
Error in statistics_wrapper (line 235)<br>
[cfg, data] = prepare_timefreq_data(cfg, varargin{:});<br>
<br>
Error in ft_timelockstatistics (line 113)<br>
[stat, cfg] = statistics_wrapper(cfg, varargin{:});<br>
<br>
<br>
</div>
<div>- Note that this worked in very earlier versions of FT<br>
<br>
<br>
</div>
<div><br>
IF I USE THIS CODE<br>
stat = ft_timelockstatistics(cfg,FP_block{:},FS_block{:});<br>
<br>
</div>
<div>IT RUNS, BUT I GET SOMETHING WEIRD WITH T-VALUES BELOW 1 and 4 clusters (see attched figure).<br>
<br>
<br>
</div>
<div>What I am doing wrong? Any advice would be great.<br>
<br>
<br>
</div>
<div>Many thanks,<br>
</div>
<div>Davide<br>
</div>
<div><br>
<br>
<br>
-- <br>
<div class="gmail_signature">Davide Rivolta, PhD<br>
</div>
</div>
</div>
</div>
</div>
</div>
</body>
</html>