<div dir="auto"><div dir="auto">Hi,</div><div dir="auto">I have two questions about correctly specifying the design matrix for ft_freqstatistics.</div><div dir="auto"><br></div><div dir="auto">It is a within-subject experimental design, where I want to test the (linear) effect of time. So, per subject, I have 5 consecutive sessions, with ~170 trials each. In the time-frequency domain, I would like to test if power is changing with session. For group statistics, I did the following:</div><div dir="auto"><br></div><div dir="auto">1) Compute per subject and session the mean with ft_freqdescriptives, because I don't have enough memory to work with single trial data.</div><div dir="auto"><br></div><div dir="auto">2) Testing the session effect with a dependent samples regression</div><div dir="auto"><br></div><div dir="auto">    % parameters for test</div><div dir="auto">    cfg = [];</div><div dir="auto">    cfg.parameter           = 'powspctrm';</div><div dir="auto">    cfg.latency             = 'all';  </div><div dir="auto">    cfg.method              = 'montecarlo';</div><div dir="auto">    cfg.statistic           = 'depsamplesregrT';</div><div dir="auto">    cfg.correcttail         = 'prob';  </div><div dir="auto">    cfg.tail                = 0;  </div><div dir="auto">    cfg.alpha               = 0.05;  </div><div dir="auto">    cfg.numrandomization    = 500; </div><div dir="auto">    cfg.correctm            = 'cluster';</div><div dir="auto">    cfg.clusteralpha        = 0.05;</div><div dir="auto">    cfg.clusterstatistic    = 'maxsum';</div><div dir="auto">    cfg.clusterthreshold    = 'nonparametric_common'; % method for single-sample threshold</div><div dir="auto">    cfg.clustertail         = 0; </div><div dir="auto">    cfg.minnbchan           = 3;</div><div dir="auto">    cfg.wcm_weight          = 1;</div><div dir="auto">     </div><div dir="auto">    </div><div dir="auto">    %design matrix</div><div dir="auto">    design       = zeros(2, n_sessions*n_subjects);</div><div dir="auto">    design(1, :) = repmat(1:n_subjects, 1, n_sessions);  % subject</div><div dir="auto">    design(2, :) = repelem(1:n_sessions, n_subjects);  % session</div><div dir="auto"><br></div><div dir="auto"><br></div><div dir="auto">    cfg.design              = design;</div><div dir="auto">    cfg.uvar                = 1;  % which row is the unit variable, here subjects</div><div dir="auto">    cfg.ivar                = 2;  % which row is the independent variable, here sessions</div><div dir="auto">    </div><div dir="auto">    %generate stats from two cell arrays of full fieldtrip "data" structures</div><div dir="auto">    stat = ft_freqstatistics(cfg, session1{:}, session2{:}, session3{:}, session4{:}, session5{:});</div><div dir="auto"><br></div><div dir="auto"><br></div><div dir="auto"><br></div><div dir="auto">This leads to a significant cluster. </div><div dir="auto"><br></div><div dir="auto">Now, I would like to test if the effect is significant in every single subject. So I thought I would just set n_subjects to 1 and could use the same approach. However, I get an error saying that I need as least two subjects (uvar). </div><div dir="auto"><br></div><div dir="auto">I'm wondering what is the right way to set this up. An indepsamplesregrT is running fine, but I think the data is still depended, as it is from one subject. </div><div dir="auto"><br></div><div dir="auto"><br></div><div dir="auto">Second, I would like to test the conjunction, i.e. finding a cluster that is overlapping between subjects and significant within *every* single subject. How could I do this in Fieldtrip?</div><div dir="auto"><br></div><div dir="auto">Thanks a lot for your help in advance!</div><div dir="auto"><br></div><div dir="auto">Best,</div><div dir="auto">Antonius</div><div dir="auto"><br></div></div>