<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
</head>
<body style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">
Hi Stephanie,
<div class=""><br class="">
</div>
<div class="">I don’t know about the referenced website from natmeg, but my point still is that you are incorrectly using ft_freqstatistics. You cannot ask the function to do a statistical comparison between two conditions if you only input a single condition.
You say that you want to identify clusters, as it turns out in this single condition. What is the threshold heuristic that you wish to use in order to distinguish data points that don’t and do belong to a cluster? The fact that cfg.neighbours = [] does not
give you what you want at this point to me still seems to be caused by an inappropriate conceptual use of ft_freqstatistics.</div>
<div class=""><br class="">
</div>
<div class="">Best wishes,</div>
<div class="">Jan-Mathijs</div>
<div class=""><br class="">
<div><br class="">
<blockquote type="cite" class="">
<div class="">On 8 Aug 2023, at 23:31, Brener, Stephanie <<a href="mailto:Stephanie.Brener@ucsf.edu" class="">Stephanie.Brener@ucsf.edu</a>> wrote:</div>
<div class=""><br class="">
<div style="caret-color: rgb(0, 0, 0); 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; -webkit-text-stroke-width: 0px; text-decoration: none;" class="">
<div class="elementToProof" style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt;">
Hi Jan-Mathijs,</div>
<div class="elementToProof" style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt;">
<br class="">
</div>
<div class="elementToProof ContentPasted0" style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt;">
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!5gfoG1LrHSOAbgEFOKGDOajACJz2kz_CZ_HJtrPbyjvZ4Qy0zHcYrXbfucGxf_fTcWrAxHNcCEPghjGYd0ODyvYGt232RdHPmlo75A$" originalsrc="https://natmeg.se/ft_statistics/statistics.html" shash="BOjwMPWUHHofbWEa9y7QcvxArMYUswMvLAtfK+icRQaPHnqPqkjyqtUlNgmNTPSjH4iGFEKhs0SX7C8pyOzpZ2v0s5UmbMbodh9R8Ed2cAtozhjxzMXZZCpyxUBogUPmQtIOpgVLTD1wc3BNRDWvhbyegURTal/ZV5Etsv7q3HU=" id="LPlnk186767" class="">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 class="elementToProof ContentPasted0" style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt;">
<br class="">
</div>
<div class="elementToProof ContentPasted0" style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt;">
Thank you,</div>
<div class="elementToProof ContentPasted0" style="font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt;">
Stephanie</div>
<div id="appendonsend" class=""></div>
<hr tabindex="-1" style="display: inline-block; width: 897.671875px;" class="">
<div id="divRplyFwdMsg" dir="ltr" class=""><font face="Calibri, sans-serif" style="font-size: 11pt;" class=""><b class="">From:</b><span class="Apple-converted-space"> </span>fieldtrip <<a href="mailto:fieldtrip-bounces@science.ru.nl" class="">fieldtrip-bounces@science.ru.nl</a>>
on behalf of Schoffelen, J.M. (Jan Mathijs) via fieldtrip <<a href="mailto:fieldtrip@science.ru.nl" class="">fieldtrip@science.ru.nl</a>><br class="">
<b class="">Sent:</b><span class="Apple-converted-space"> </span>Tuesday, August 8, 2023 10:36 AM<br class="">
<b class="">To:</b><span class="Apple-converted-space"> </span>FieldTrip discussion list <<a href="mailto:fieldtrip@science.ru.nl" class="">fieldtrip@science.ru.nl</a>><br class="">
<b class="">Cc:</b><span class="Apple-converted-space"> </span>Schoffelen, J.M. (Jan Mathijs) <<a href="mailto:janmathijs.schoffelen@donders.ru.nl" class="">janmathijs.schoffelen@donders.ru.nl</a>><br class="">
<b class="">Subject:</b><span class="Apple-converted-space"> </span>Re: [FieldTrip] Cluster based permutation test for single channel, ft_freqstatistics</font></div>
</div>
</div>
</blockquote>
<br class="">
<blockquote type="cite" class="">
<div class="">
<div style="caret-color: rgb(0, 0, 0); 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; -webkit-text-stroke-width: 0px; text-decoration: none;" class="">
<div id="divRplyFwdMsg" dir="ltr" class="">
<div class="">Hi Stephanie,</div>
</div>
<div class="">
<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 class=""><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_x_elementToProof x_x_ContentPasted0 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_x_elementToProof x_x_ContentPasted0 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_x_elementToProof x_x_ContentPasted0 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_x_elementToProof x_x_ContentPasted0 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_ContentPasted0 x_x_x_elementToProof" style="margin: 0px; background-color: rgb(255, 255, 255);">
<br class="x_x_ContentPasted0 x_ContentPasted0">
</div>
<div class="x_ContentPasted0 x_x_x_ContentPasted2 x_x_ContentPasted0 x_x_x_elementToProof x_x_x_ContentPasted0" style="margin: 0px; background-color: rgb(255, 255, 255);">
data_ft = [];
<div class="x_x_x_ContentPasted0 x_x_ContentPasted0 x_ContentPasted0" style="margin: 0px;">
data_ft.dimord = 'rpt_chan_freq_time';</div>
<div class="x_x_x_ContentPasted0 x_x_ContentPasted0 x_ContentPasted0" style="margin: 0px;">
data_ft.label = {'Onset','Choice','Outcome'}';</div>
<div class="x_x_x_ContentPasted0 x_x_ContentPasted0 x_ContentPasted0" style="margin: 0px;">
data_ft.freq = linspace(1, 50, 50);</div>
<div class="x_x_x_ContentPasted0 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_ContentPasted0 x_x_x_ContentPasted2">
<br class="x_x_ContentPasted0 x_ContentPasted0">
</div>
<div class="x_ContentPasted0 x_x_elementToProof x_x_ContentPasted0 x_x_x_elementToProof" 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_x_elementToProof x_x_ContentPasted0 x_ContentPasted0" style="margin: 0px; background-color: rgb(255, 255, 255);">
The cfg for ft_freqstatistics looks like this:</div>
<div class="x_ContentPasted0 x_x_ContentPasted0 x_x_x_elementToProof x_x_x_ContentPasted1" style="margin: 0px; background-color: rgb(255, 255, 255);">
cfg = [];
<div class="x_ContentPasted0 x_x_elementToProof x_x_ContentPasted0 x_x_x_ContentPasted1" 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_elementToProof x_x_x_ContentPasted1" style="margin: 0px; background-color: rgb(255, 255, 255);">
<br class="x_x_ContentPasted0 x_ContentPasted0">
</div>
<div class="x_ContentPasted0 x_x_ContentPasted0 x_x_x_elementToProof x_x_x_ContentPasted1" 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_elementToProof x_x_x_ContentPasted1" style="margin: 0px; background-color: rgb(255, 255, 255);">
<br class="x_x_ContentPasted0 x_ContentPasted0">
</div>
<div class="x_ContentPasted0 x_x_ContentPasted0 x_x_x_elementToProof x_x_x_ContentPasted1 x_x_x_ContentPasted3" style="margin: 0px; background-color: rgb(255, 255, 255);">
Error using findcluster
<div class="x_x_x_ContentPasted3 x_x_ContentPasted0 x_ContentPasted0" style="margin: 0px;">
invalid dimension of spatdimneighbstructmat</div>
<div class="" style="margin: 0px;"><br class="x_x_x_ContentPasted3 x_x_ContentPasted0 x_ContentPasted0">
</div>
<div class="x_x_x_ContentPasted3 x_x_ContentPasted0 x_ContentPasted0" style="margin: 0px;">
Error in clusterstat (line 214)</div>
<div class="x_x_x_ContentPasted3 x_x_ContentPasted0 x_ContentPasted0" style="margin: 0px;">
posclusobs = findcluster(tmp, connmat, cfg.minnbchan);</div>
<div class="" style="margin: 0px;"><br class="x_x_x_ContentPasted3 x_x_ContentPasted0 x_ContentPasted0">
</div>
<div class="x_x_x_ContentPasted3 x_x_ContentPasted0 x_ContentPasted0" style="margin: 0px;">
Error in ft_statistics_montecarlo (line 364)</div>
<div class="x_x_x_ContentPasted3 x_x_ContentPasted0 x_ContentPasted0" style="margin: 0px;">
[stat, cfg] = clusterstat(cfg, statrand, statobs);</div>
<div class="" style="margin: 0px;"><br class="x_x_x_ContentPasted3 x_x_ContentPasted0 x_ContentPasted0">
</div>
<div class="x_x_x_ContentPasted3 x_x_ContentPasted0 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_ContentPasted3 x_x_x_elementToProof x_x_x_ContentPasted1" 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_ContentPasted0 x_x_ContentPasted0 x_x_x_elementToProof x_x_x_ContentPasted1 x_x_x_ContentPasted3" 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_ContentPasted1 x_x_x_ContentPasted3" 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_ContentPasted0 x_x_ContentPasted0 x_x_x_elementToProof x_x_x_ContentPasted1 x_x_x_ContentPasted3" 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_ContentPasted0 x_x_ContentPasted0 x_x_x_elementToProof x_x_x_ContentPasted1 x_x_x_ContentPasted3" 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_ContentPasted0 x_x_ContentPasted0 x_x_x_elementToProof x_x_x_ContentPasted1 x_x_x_ContentPasted3" 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$" originalsrc="https://urldefense.com/v3/__https://mailman.science.ru.nl/mailman/listinfo/fieldtrip__;!!LQC6Cpwp!qZKfpBRRtXD-MhJNgK8Umx7bPEupxVdrVdeMrETazqQDeha0_MbS_da176LrjM7A3I0fRDP8yENOe5wA5VfOGjRg_4fr9g$" shash="O9ylv98LALhZJVEfs6pdBKXewGXtK9xArmJ3883HxfnFmGdjNYlS0FPwqkyydGqbua/fyfeO/0QKpTkGZHzm0w2j8+WMruYIaBo/MDtlS3XP2qjn6/1PmAmUuYeWyP01LZi68iJERhWC0la1IkoID9+ra2gx7o0qeh+wElpRH+0=" 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;">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$" originalsrc="https://urldefense.com/v3/__https://doi.org/10.1371/journal.pcbi.1002202__;!!HJOPV4FYYWzcc1jazlU!6lvoYv-5AUVQgVqdkunLQW47GM9q2ZAfISkCDDLA_izYMc74Be5DxwaonX_MSWy1cdx5IUlflE2Mdlb_DmEgVJkDp15TAJ724jS-WQ$" shash="rEUEBeYH4uLvUFLA5oWPX5hyoOtWgZhmSGmSPU+bi/h4VFkVg0MIA+lo5kYhBOproBUHQvLI3KAWn3aYywQCek/tDhOPCQNa3pYVxy6qtHJoEJU/Ae+Velc6kjzvOR3TKxe4xsETlYngwC6ZNHsvsa+SSj3lmF6cDy61koCgJPA=" 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;">https://doi.org/10.1371/journal.pcbi.1002202</a></div>
</blockquote>
</div>
</div>
</div>
</div>
</div>
</blockquote>
</div>
<br class="">
</div>
</body>
</html>