<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=Windows-1252">
<style type="text/css" style="display:none;"> P {margin-top:0;margin-bottom:0;} </style>
</head>
<body dir="ltr">
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);" class="elementToProof">
Hi Jan-Mathijs,</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);" class="elementToProof">
<br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);" class="elementToProof ContentPasted0">
Thank you for your response. In my analysis, I am not trying to compare between the Onset event and the Choice event. I only want to find clusters in time/frequency across all trials in just the Onset event (then repeat for the Choice event). The goal is to
do the standard clustering method, like described here: <a href="https://urldefense.com/v3/__https://natmeg.se/ft_statistics/statistics.html__;!!HJOPV4FYYWzcc1jazlU!8vz-0y-BMtBAed_rHKBjXtLBx5_dTvhl5rch206OvsW0fQKUcFx778iHe96MJFEyagz0ziEo3P8heHy0Igib3HrN05LEKOnBYg$" id="LPlnk186767">https://natmeg.se/ft_statistics/statistics.html</a>, for the case where there is only one spatial channel. The statistical test
will be checking whether the population mean (across all trials) is different than 0, rather than different from another population. Various posts have mentioned that the way to find clusters by time/frequency instead of spatially is by leaving the cfg.neighbours
field blank. However, I still run into errors when doing that. Do you have any insight into finding clusters in time/frequency?</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);" class="elementToProof ContentPasted0">
<br>
</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);" class="elementToProof ContentPasted0">
Thank you,</div>
<div style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);" class="elementToProof ContentPasted0">
Stephanie</div>
<div id="appendonsend"></div>
<hr style="display:inline-block;width:98%" tabindex="-1">
<div id="divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" style="font-size:11pt" color="#000000"><b>From:</b> fieldtrip <fieldtrip-bounces@science.ru.nl> on behalf of Schoffelen, J.M. (Jan Mathijs) via fieldtrip <fieldtrip@science.ru.nl><br>
<b>Sent:</b> Tuesday, August 8, 2023 10:36 AM<br>
<b>To:</b> FieldTrip discussion list <fieldtrip@science.ru.nl><br>
<b>Cc:</b> Schoffelen, J.M. (Jan Mathijs) <janmathijs.schoffelen@donders.ru.nl><br>
<b>Subject:</b> Re: [FieldTrip] Cluster based permutation test for single channel, ft_freqstatistics</font>
<div> </div>
</div>
<div>
<div style="display:none!important; display:none; visibility:hidden; font-size:1px; color:#ffffff; line-height:1px; height:0px; max-height:0px; opacity:0; overflow:hidden">
Hi Stephanie, I think that there is something wrong with what you want to achieve (at least what I think you want to achieve), and the way you try to achieve it. If you want to make a comparison between the onset/choice/outcome segments of
</div>
<div style="display:none!important; display:none; visibility:hidden; font-size:1px; color:#ffffff; line-height:1px; height:0px; max-height:0px; opacity:0; overflow:hidden">
ZjQcmQRYFpfptBannerStart</div>
<div dir="ltr" lang="en" id="x_pfptBannerj290t1g" style="display:block!important; text-align:left!important; margin:16px 0px 16px 0px!important; padding:8px 16px 8px 16px!important; border-radius:4px!important; min-width:200px!important; background-color:#D0D8DC!important; background-color:#D0D8DC; border-top:4px solid #90a4ae!important; border-top:4px solid #90a4ae">
<div id="x_pfptBannerj290t1g" style="float:left!important; display:block!important; margin:0px 0px 1px 0px!important; max-width:600px!important">
<div id="x_pfptBannerj290t1g" style="display:block!important; visibility:visible!important; background-color:#D0D8DC!important; color:#000000!important; color:#000000; font-family:'Arial',sans-serif!important; font-family:'Arial',sans-serif; font-weight:bold!important; font-weight:bold; font-size:14px!important; line-height:18px!important; line-height:18px">
This Message Is From an External Sender </div>
<div id="x_pfptBannerj290t1g" style="display:block!important; visibility:visible!important; background-color:#D0D8DC!important; color:#000000!important; color:#000000; font-weight:normal; font-family:'Arial',sans-serif!important; font-family:'Arial',sans-serif; font-size:12px!important; line-height:18px!important; line-height:18px; margin-top:2px!important">
This message came from outside your organization. </div>
</div>
<div style="clear:both!important; display:block!important; visibility:hidden!important; line-height:0!important; font-size:0.01px!important; height:0px">
</div>
</div>
<div style="display:none!important; display:none; visibility:hidden; font-size:1px; color:#ffffff; line-height:1px; height:0px; max-height:0px; opacity:0; overflow:hidden">
ZjQcmQRYFpfptBannerEnd</div>
<style>
<!--
#x_pfptBannerj290t1g
{display:block!important;
visibility:visible!important;
opacity:1!important;
background-color:#D0D8DC!important;
max-width:none!important;
max-height:none!important}
-->
</style>Hi Stephanie,
<div class=""><br class="">
</div>
<div class="">I think that there is something wrong with what you want to achieve (at least what I think you want to achieve), and the way you try to achieve it.</div>
<div class=""><br class="">
</div>
<div class="">If you want to make a comparison between the onset/choice/outcome segments of your data, you should NOT artificially represent them as different channels. Rather, you’d need to represent each of the individual (sub-)segments of the individual
trials (if you want to use the trials as a unit of observation) as separate data matrices, or the average per (sub-)segment across trials for each patient (if you want to use the patient as a unit of observation). Then, if you want to use a(n in)dependent
T-statistic as the test statistic, you can only do a pairwise comparison, and this means that the specified design should contains 1’s and 2’s.</div>
<div class=""><br class="">
</div>
<div class="">I would expect you to create 3 data structures, one for each part of the trial onset/choice/outcome, which can be either compared in a pairwise fashion using (in)depsamplesT, or in a single shot, using an F-statistic. Note, that given the presence
of a time dimension, there’s a hard requirement that the time axes for each of the conditions has the same number of time points (and that the corresponding time axis should be faked to be numerically identical across the three conditions) in order for the
FT-machinery to work. I think it would make sense to collapse across time in the first instance, because my hunch is that 1) the condition that I mentioned in the previous sentence is not met in your data, and 2) I don’t know whether there is a specific physiological
reason to assume that it is meaningful to compare in a pairwise fashion -say- the power at x ms after task onset with the power at exact x ms after choice onset.</div>
<div class=""><br class="">
</div>
<div class="">Good luck,</div>
<div class="">Jan-Mathijs</div>
<div class=""><br class="">
</div>
<div class="">
<div><br class="">
<blockquote type="cite" class="">
<div class="">On 7 Aug 2023, at 20:46, Brener, Stephanie via fieldtrip <<a href="mailto:fieldtrip@science.ru.nl" class="">fieldtrip@science.ru.nl</a>> wrote:</div>
<br class="x_Apple-interchange-newline">
<div class="">
<div class="x_elementToProof" style="font-style:normal; font-variant-caps:normal; font-weight:normal; letter-spacing:normal; text-align:start; text-indent:0px; text-transform:none; white-space:normal; word-spacing:0px; text-decoration:none; font-family:Calibri,Arial,Helvetica,sans-serif; font-size:12pt">
<span class="x_x_elementToProof" style="font-size:12pt; margin:0px; background-color:rgb(255,255,255)"><span class="x_x_ContentPasted0 x_x_x_elementToProof x_ContentPasted0" style="margin:0px; background-color:rgb(255,255,255)">Hi fieldtrippers,</span></span></div>
<div class="x_elementToProof" style="font-style:normal; font-variant-caps:normal; font-weight:normal; letter-spacing:normal; text-align:start; text-indent:0px; text-transform:none; white-space:normal; word-spacing:0px; text-decoration:none; font-family:Calibri,Arial,Helvetica,sans-serif; font-size:12pt">
<div class="x_x_elementToProof" style="font-size:12pt; margin:0px; background-color:rgb(255,255,255)">
<div class="x_x_x_elementToProof" style="margin:0px; background-color:rgb(255,255,255)">
<br class="x_x_ContentPasted0 x_ContentPasted0">
</div>
<div class="x_x_ContentPasted0 x_x_x_elementToProof x_ContentPasted0" style="margin:0px; background-color:rgb(255,255,255)">
I am trying to perform cluster-based permutation testing on a single channel of DBS electrode data to cluster in time/frequency rather than spatially. I have already pre-processed, cleaned, and separated time-locked events for the data outside of fieldtrip,
leaving me with measures of baseline-corrected power in a data array with these dimensions: (<span class="x_x_ContentPasted0 x_ContentPasted0" style="margin:0px">50 frequencies x 500 time samples x 180 trials x 3 events).</span></div>
<div class="x_x_x_elementToProof" style="margin:0px; background-color:rgb(255,255,255)">
<br class="x_x_ContentPasted0 x_ContentPasted0">
</div>
<div class="x_x_ContentPasted0 x_x_x_elementToProof x_ContentPasted0" style="margin:0px; background-color:rgb(255,255,255)">
The dimensions can be described by:</div>
<div class="x_x_x_elementToProof" style="margin:0px; background-color:rgb(255,255,255)">
<ul data-editing-info="{"orderedStyleType":1,"unorderedStyleType":2}" class="">
<li class="" style="list-style-type:"- ""><span class="x_x_ContentPasted0 x_ContentPasted0" style="margin:0px">50 Morlet wavelets from 0 to 50Hz --> power at 50 frequencies</span></li><li class="" style="list-style-type:"- ""><span class="x_x_ContentPasted0 x_ContentPasted0" style="margin:0px">500 samples = 1 second before and after time locked event in the task (at 250Hz sampling rate)</span></li><li class="" style="list-style-type:"- ""><span class="x_x_ContentPasted0 x_ContentPasted0" style="margin:0px">180 trials of the task</span></li><li class="" style="list-style-type:"- ""><span class="x_x_ContentPasted0 x_ContentPasted0" style="margin:0px">3 events/epochs - onset of task, choice taken in task, outcome of task</span></li></ul>
</div>
<div class="x_x_x_elementToProof" style="margin:0px; background-color:rgb(255,255,255)">
<br class="x_x_ContentPasted0 x_ContentPasted0">
</div>
<div class="x_x_ContentPasted0 x_x_x_elementToProof x_ContentPasted0" style="margin:0px; background-color:rgb(255,255,255)">
If I were to plot the average power across all trials for the onset event, I can see how power changes on a spectrogram plot. My goal is to find the significant clusters that show up on the spectrogram as a first step to finding interesting frequency bands.
To do so for one subject (with a single DBS lead), I reorganized the data to have these dimensions (180 trials x 3 events x 50 freqs x 500 time) and put it into a fieldtrip struct as follows:</div>
<div class="x_x_x_elementToProof x_x_x_ContentPasted0" style="margin:0px; background-color:rgb(255,255,255)">
<br class="x_x_ContentPasted0 x_ContentPasted0">
</div>
<div class="x_x_ContentPasted0 x_x_x_elementToProof x_x_x_ContentPasted0 x_ContentPasted0 x_x_x_ContentPasted2" style="margin:0px; background-color:rgb(255,255,255)">
data_ft = [];
<div class="x_x_ContentPasted0 x_x_x_ContentPasted0 x_ContentPasted0" style="margin:0px">
data_ft.dimord = 'rpt_chan_freq_time';</div>
<div class="x_x_ContentPasted0 x_x_x_ContentPasted0 x_ContentPasted0" style="margin:0px">
data_ft.label = {'Onset','Choice','Outcome'}';</div>
<div class="x_x_ContentPasted0 x_x_x_ContentPasted0 x_ContentPasted0" style="margin:0px">
data_ft.freq = linspace(1, 50, 50);</div>
<div class="x_x_ContentPasted0 x_x_x_ContentPasted0 x_ContentPasted0" style="margin:0px">
data_ft.time = linspace(-1, 1, 500);<br class="x_x_ContentPasted0 x_ContentPasted0">
</div>
data_ft.powspctrm = data;<br class="x_ContentPasted0 x_x_x_ContentPasted2 x_x_ContentPasted0">
<br class="x_x_ContentPasted0 x_ContentPasted0">
</div>
<div class="x_x_ContentPasted0 x_x_x_elementToProof x_x_elementToProof x_ContentPasted0" style="margin:0px; background-color:rgb(255,255,255)">
where I am using the events as the "channels" so I can choose one at a time.</div>
<div class="x_x_x_elementToProof" style="margin:0px; background-color:rgb(255,255,255)">
<br class="x_x_ContentPasted0 x_ContentPasted0">
</div>
<div class="x_x_ContentPasted0 x_x_x_elementToProof x_ContentPasted0" style="margin:0px; background-color:rgb(255,255,255)">
The cfg for ft_freqstatistics looks like this:</div>
<div class="x_x_ContentPasted0 x_x_x_elementToProof x_x_x_ContentPasted1 x_ContentPasted0" style="margin:0px; background-color:rgb(255,255,255)">
cfg = [];
<div class="x_x_ContentPasted0 x_x_elementToProof x_x_x_ContentPasted1 x_ContentPasted0" style="margin:0px">
cfg.method = 'montecarlo';</div>
<div class="x_x_x_ContentPasted1 x_x_ContentPasted0 x_ContentPasted0" style="margin:0px">
cfg.statistic = 'indepsamplesT';</div>
<div class="x_x_x_ContentPasted1 x_x_ContentPasted0 x_ContentPasted0" style="margin:0px">
cfg.correctm = 'cluster';</div>
<div class="x_x_x_ContentPasted1 x_x_ContentPasted0 x_ContentPasted0" style="margin:0px">
cfg.clusteralpha = 0.05;</div>
<div class="x_x_x_ContentPasted1 x_x_ContentPasted0 x_ContentPasted0" style="margin:0px">
cfg.clustertail = 0;</div>
<div class="x_x_x_ContentPasted1 x_x_ContentPasted0 x_ContentPasted0" style="margin:0px">
cfg.clusterstatistic = 'maxsum';</div>
<div class="x_x_x_ContentPasted1 x_x_ContentPasted0 x_ContentPasted0" style="margin:0px">
cfg.tail = 0;</div>
<div class="x_x_x_ContentPasted1 x_x_ContentPasted0 x_ContentPasted0" style="margin:0px">
cfg.alpha = 0.025;</div>
<div class="x_x_x_ContentPasted1 x_x_ContentPasted0 x_ContentPasted0" style="margin:0px">
cfg.numrandomization = 1000; </div>
<div class="x_x_x_ContentPasted1 x_x_ContentPasted0 x_ContentPasted0" style="margin:0px">
cfg.design = ones(size(data_ft.powspctrm, 1), 1); % size = 180 x 1</div>
<div class="x_x_x_ContentPasted1 x_x_ContentPasted0 x_ContentPasted0" style="margin:0px">
cfg.ivar = 1; </div>
<div class="x_x_x_ContentPasted1 x_x_ContentPasted0 x_ContentPasted0" style="margin:0px">
cfg.channel = {'Onset'};</div>
<div class="x_x_x_ContentPasted1 x_x_ContentPasted0 x_ContentPasted0" style="margin:0px">
cfg.neighbours = []; </div>
<div class="" style="margin:0px"><br class="x_x_ContentPasted0 x_ContentPasted0">
</div>
[stat] = ft_freqstatistics(cfg, data_ft);<br class="x_x_ContentPasted0 x_ContentPasted0">
</div>
<div class="x_x_x_ContentPasted1 x_x_x_elementToProof" style="margin:0px; background-color:rgb(255,255,255)">
<br class="x_x_ContentPasted0 x_ContentPasted0">
</div>
<div class="x_x_ContentPasted0 x_x_x_elementToProof x_x_x_ContentPasted1 x_ContentPasted0" style="margin:0px; background-color:rgb(255,255,255)">
following other fieldtrip threads noting that the neighbours field should be left blank to denote clustering in time/frequency rather than spatially across channels. When I run ft_freqstatistics, I get the following error:</div>
<div class="x_x_x_ContentPasted1 x_x_x_elementToProof" style="margin:0px; background-color:rgb(255,255,255)">
<br class="x_x_ContentPasted0 x_ContentPasted0">
</div>
<div class="x_x_ContentPasted0 x_x_x_elementToProof x_x_x_ContentPasted3 x_x_x_ContentPasted1 x_ContentPasted0" style="margin:0px; background-color:rgb(255,255,255)">
Error using findcluster
<div class="x_x_ContentPasted0 x_x_x_ContentPasted3 x_ContentPasted0" style="margin:0px">
invalid dimension of spatdimneighbstructmat</div>
<div class="" style="margin:0px"><br class="x_x_ContentPasted0 x_x_x_ContentPasted3 x_ContentPasted0">
</div>
<div class="x_x_ContentPasted0 x_x_x_ContentPasted3 x_ContentPasted0" style="margin:0px">
Error in clusterstat (line 214)</div>
<div class="x_x_ContentPasted0 x_x_x_ContentPasted3 x_ContentPasted0" style="margin:0px">
posclusobs = findcluster(tmp, connmat, cfg.minnbchan);</div>
<div class="" style="margin:0px"><br class="x_x_ContentPasted0 x_x_x_ContentPasted3 x_ContentPasted0">
</div>
<div class="x_x_ContentPasted0 x_x_x_ContentPasted3 x_ContentPasted0" style="margin:0px">
Error in ft_statistics_montecarlo (line 364)</div>
<div class="x_x_ContentPasted0 x_x_x_ContentPasted3 x_ContentPasted0" style="margin:0px">
[stat, cfg] = clusterstat(cfg, statrand, statobs);</div>
<div class="" style="margin:0px"><br class="x_x_ContentPasted0 x_x_x_ContentPasted3 x_ContentPasted0">
</div>
<div class="x_x_ContentPasted0 x_x_x_ContentPasted3 x_ContentPasted0" style="margin:0px">
Error in ft_freqstatistics (line 194)</div>
[stat, cfg] = statmethod(cfg, dat, design);<br class="x_x_ContentPasted0 x_ContentPasted0">
</div>
<div class="x_x_x_ContentPasted1 x_x_x_elementToProof x_x_x_ContentPasted3" style="margin:0px; background-color:rgb(255,255,255)">
<br class="x_x_ContentPasted0 x_ContentPasted0">
</div>
</div>
<div class="x_x_elementToProof" style="font-size:12pt; margin:0px; background-color:rgb(255,255,255)">
<span class="x_x_ContentPasted0 x_x_x_elementToProof x_x_x_ContentPasted3 x_x_x_ContentPasted1 x_ContentPasted0" style="margin:0px; background-color:rgb(255,255,255)">The dimension of spatdimneighbstructmat is given by connmat, which is 0 in this case since
I didn't specify a neighboring scheme. Is there any other way to cluster in time/frequency and not spatially when looking at only a single channel of data?</span><br class="x_ContentPasted0">
</div>
<div class="x_x_elementToProof" style="font-size:12pt; margin:0px; background-color:rgb(255,255,255)">
<span class="x_x_ContentPasted0 x_x_x_elementToProof x_x_x_ContentPasted3 x_x_x_ContentPasted1" style="margin:0px; background-color:rgb(255,255,255)"><br class="x_ContentPasted0">
</span></div>
<div class="x_x_elementToProof" style="font-size:12pt; margin:0px; background-color:rgb(255,255,255)">
<span class="x_x_ContentPasted0 x_x_x_elementToProof x_x_x_ContentPasted3 x_x_x_ContentPasted1 x_ContentPasted0" style="margin:0px; background-color:rgb(255,255,255)">Thank you,</span></div>
<div class="x_x_elementToProof" style="font-size:12pt; margin:0px; background-color:rgb(255,255,255)">
<span class="x_x_ContentPasted0 x_x_x_elementToProof x_x_x_ContentPasted3 x_x_x_ContentPasted1 x_ContentPasted0" style="margin:0px; background-color:rgb(255,255,255)">Stephanie Brener</span></div>
<span class="x_x_elementToProof" style="font-size:12pt; margin:0px; background-color:rgb(255,255,255)"><span class="x_x_ContentPasted0 x_x_x_elementToProof x_x_x_ContentPasted3 x_x_x_ContentPasted1 x_ContentPasted0" style="margin:0px; background-color:rgb(255,255,255)">UC
San Francisco</span></span><br class="">
</div>
<span class="" style="font-family:Helvetica; font-size:14px; font-style:normal; font-variant-caps:normal; font-weight:normal; letter-spacing:normal; text-align:start; text-indent:0px; text-transform:none; white-space:normal; word-spacing:0px; text-decoration:none; float:none; display:inline!important">_______________________________________________</span><br class="" style="font-family:Helvetica; font-size:14px; font-style:normal; font-variant-caps:normal; font-weight:normal; letter-spacing:normal; text-align:start; text-indent:0px; text-transform:none; white-space:normal; word-spacing:0px; text-decoration:none">
<span class="" style="font-family:Helvetica; font-size:14px; font-style:normal; font-variant-caps:normal; font-weight:normal; letter-spacing:normal; text-align:start; text-indent:0px; text-transform:none; white-space:normal; word-spacing:0px; text-decoration:none; float:none; display:inline!important">fieldtrip
mailing list</span><br class="" style="font-family:Helvetica; font-size:14px; font-style:normal; font-variant-caps:normal; font-weight:normal; letter-spacing:normal; text-align:start; text-indent:0px; text-transform:none; white-space:normal; word-spacing:0px; text-decoration:none">
<a href="https://urldefense.com/v3/__https://mailman.science.ru.nl/mailman/listinfo/fieldtrip__;!!LQC6Cpwp!qZKfpBRRtXD-MhJNgK8Umx7bPEupxVdrVdeMrETazqQDeha0_MbS_da176LrjM7A3I0fRDP8yENOe5wA5VfOGjRg_4fr9g$" class="" style="font-family:Helvetica; font-size:14px; font-style:normal; font-variant-caps:normal; font-weight:normal; letter-spacing:normal; orphans:auto; text-align:start; text-indent:0px; text-transform:none; white-space:normal; widows:auto; word-spacing:0px">https://mailman.science.ru.nl/mailman/listinfo/fieldtrip</a><br class="" style="font-family:Helvetica; font-size:14px; font-style:normal; font-variant-caps:normal; font-weight:normal; letter-spacing:normal; text-align:start; text-indent:0px; text-transform:none; white-space:normal; word-spacing:0px; text-decoration:none">
<a href="https://urldefense.com/v3/__https://doi.org/10.1371/journal.pcbi.1002202__;!!HJOPV4FYYWzcc1jazlU!6lvoYv-5AUVQgVqdkunLQW47GM9q2ZAfISkCDDLA_izYMc74Be5DxwaonX_MSWy1cdx5IUlflE2Mdlb_DmEgVJkDp15TAJ724jS-WQ$" class="" style="font-family:Helvetica; font-size:14px; font-style:normal; font-variant-caps:normal; font-weight:normal; letter-spacing:normal; orphans:auto; text-align:start; text-indent:0px; text-transform:none; white-space:normal; widows:auto; word-spacing:0px">https://doi.org/10.1371/journal.pcbi.1002202</a></div>
</blockquote>
</div>
<br class="">
</div>
</div>
</body>
</html>