<html xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40"><head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
<meta name="Generator" content="Microsoft Word 15 (filtered medium)">
<style><!--
/* Font Definitions */
@font-face
{font-family:Helvetica;
panose-1:0 0 0 0 0 0 0 0 0 0;}
@font-face
{font-family:Courier;
panose-1:0 0 0 0 0 0 0 0 0 0;}
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri",sans-serif;}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:purple;
text-decoration:underline;}
p.msonormal0, li.msonormal0, div.msonormal0
{mso-style-name:msonormal;
mso-margin-top-alt:auto;
margin-right:0cm;
mso-margin-bottom-alt:auto;
margin-left:0cm;
font-size:11.0pt;
font-family:"Calibri",sans-serif;}
span.apple-converted-space
{mso-style-name:apple-converted-space;}
span.s1
{mso-style-name:s1;}
span.s2
{mso-style-name:s2;}
span.s3
{mso-style-name:s3;}
span.EmailStyle22
{mso-style-type:personal-reply;
font-family:"Calibri",sans-serif;
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;
font-size:10.0pt;}
@page WordSection1
{size:612.0pt 792.0pt;
margin:72.0pt 72.0pt 72.0pt 72.0pt;}
div.WordSection1
{page:WordSection1;}
--></style>
</head>
<body lang="EN-NZ" link="blue" vlink="purple">
<div class="WordSection1">
<p class="MsoNormal">Hi Jan-Mathijs,<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Thank you very much for your suggestion. I implemented this into my code and it gives indeed nicely the statistics for the whole group. However, before I want to do statistics on group level, I would like to do it on an individual level,
because I have seen with the frequency analysis I did earlier that there are considerable differences between my participants.<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">I based the code I created on the example I found on: <a href="http://www.fieldtriptoolbox.org/example/source_statistics/">
http://www.fieldtriptoolbox.org/example/source_statistics/</a>. In order to get a single source object with the source estimates for the single trials I used the rawtrials setting. Is there another, better option to get this?
<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Would a script making separate source objects for each trial (similar to what you suggest below for each participant times each condition) be a better solution?<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Would this in-participant analysis still be a depsamplesT test statistic?<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Regards,<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Xavier<o:p></o:p></p>
<p class="MsoNormal"><o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<div style="border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0cm 0cm 0cm">
<p class="MsoNormal"><b><span style="font-size:12.0pt;color:black">From: </span></b><span style="font-size:12.0pt;color:black">fieldtrip <fieldtrip-bounces@science.ru.nl> on behalf of "Schoffelen, J.M. (Jan Mathijs)" <jan.schoffelen@donders.ru.nl><br>
<b>Reply-To: </b>FieldTrip discussion list <fieldtrip@science.ru.nl><br>
<b>Date: </b>Friday, 22 February 2019 at 4:12 AM<br>
<b>To: </b>FieldTrip discussion list <fieldtrip@science.ru.nl><br>
<b>Subject: </b>Re: [FieldTrip] Error in source localisation statistics<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<p class="MsoNormal">Hi Xavier, <o:p></o:p></p>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">I see. Unfortunately you happened to end up in a corner of the code that is a bit shaky, and possibly not well tested (although once upon a time it may have worked). This may be due to the fact that you are trying to ‘fool’ fieldtrip into
treating a multiple subject design as if it’s a single subject design with multiple trials (where each of the trials represents a single condition estimate for a subject), correct? (as a side note: you would want to use ‘depsamplesT’ rather than ‘indepsamplesT’
as your test statistic).<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">I think that you should not do this. First, as you have noticed this leads to downstream problems (specifically handcrafting combined data could easily lead to a user-error causing the data to deviate according to what FieldTrip can handle,
but also using the cfg.rawtrial option of ft_sourceanalaysis historically and notoriously is a dysfunctional (or at least unpredictable in its behavior) option. (nominated to be removed once we find the time). Next to this, even though it seems your geometrical
data (i.e. electrode positions, headmodel, sourcemodel) are shared across subjects, the covariance structure of the data will for sure be subject-specific, and I don’t think that it is a good idea to combine data across subjects before sourceanalysis.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">This is what I would do:<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">for each subject do:<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">ft_timelockanalysis with the subject-specific data combined across conditions<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">ft_timelockanalysis with the subject-specific data for each of the conditions separately<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">ft_sourceanalysis with the subject-specific data combined across conditions AND cfg.lcmv.keepfilter = ‘yes’; (and don’t specify cfg.keeptrials or cfg.rawtrial)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">input the source.avg.filter from the previous step into the cfg<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">for each condition (in each subject) do:<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> ft_sourceanalysis<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">after this you should have a set of Nsubject times number of conditions source structures<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">then you should be able to use ft_sourcestatistics with cfg.statisti = ‘depsamplesT’:<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">ft_sourcestatistics(cfg, source{:}); % if in the previous for-loop you have stored the individual source-objects into a cell-array this syntax works (and is much cleaner than either ft_sourcestatistics(cfg, source1a, source2a ,…etc), and
cleaner than going through ft_sourcegrandaverage).<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">ft_sourcedescriptives should not be needed in this recipe.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">Good luck,<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">Jan-Mathijs<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<div>
<p class="MsoNormal">J.M.Schoffelen, MD PhD<br>
Senior Researcher, VIDI-fellow - PI, language in interaction<br>
Telephone: +31-24-3614793<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">Physical location: room 00.028<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">Donders Centre for Cognitive Neuroimaging, Nijmegen, The Netherlands<o:p></o:p></p>
</div>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
<div>
<div>
<blockquote style="margin-top:5.0pt;margin-bottom:5.0pt">
<div>
<p class="MsoNormal">On 20 Feb 2019, at 21:32, Xavier Vrijdag <<a href="mailto:x.vrijdag@auckland.ac.nz">x.vrijdag@auckland.ac.nz</a>> wrote:<o:p></o:p></p>
</div>
<p class="MsoNormal"><o:p> </o:p></p>
<div>
<div>
<p class="MsoNormal">Hi Jan-Mathijs,<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">Thank you for your help. I went in to have a look at what is going on inside those loops and found that the error is occurring in iteration 6, when the field source.trial(1).label is concatenated. It uses the indexes for the big data cells
(ori, mom & filter) which clearly doesn’t work on the small cell with the labels. This seems to me that this is a bug.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">I added a piece of code in front of the loop that was giving the error (before line 177 and make the if and elseif.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:7.5pt;font-family:Courier"> <span class="apple-converted-space"> </span><span style="color:blue">if</span><span class="apple-converted-space"> </span>isequal(fn{i},<span class="apple-converted-space"> </span><span style="color:#A020F0">'label'</span>)</span><o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:7.5pt;font-family:Courier"> <span class="apple-converted-space"> </span><span style="color:blue">for</span><span class="apple-converted-space"> </span>j=1:nrpt</span><o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:7.5pt;font-family:Courier"> source_new.label(:,j) = source.trial(j).label;</span><o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:7.5pt;font-family:Courier"> <span class="apple-converted-space"> </span><span style="color:blue">end</span></span><o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:7.5pt;font-family:Courier"> <span class="apple-converted-space"> </span><span style="color:blue">elseif</span><span class="apple-converted-space"> </span>iscell(dat)</span><o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">This worked to circumvent the problem. However, I do get the warning that the dimord of “ori” and “df” could not be determined. When using this cleaned source struct in source statistics I get further errors.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">It looks like an issue with how the source struct is built by the source analysis function with trials, and lateron needs to be restructured. SO I would like to know what is the needed structure for ft_sourcestatistics, so we can adjust
the code correctly.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">Regards,<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">Xavier<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> <o:p></o:p></p>
</div>
<div style="border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0cm 0cm 0cm">
<div>
<p class="MsoNormal"><b><span style="font-size:12.0pt">From:<span class="apple-converted-space"> </span></span></b><span style="font-size:12.0pt">fieldtrip <<a href="mailto:fieldtrip-bounces@science.ru.nl">fieldtrip-bounces@science.ru.nl</a>> on behalf of "Schoffelen,
J.M. (Jan Mathijs)" <<a href="mailto:jan.schoffelen@donders.ru.nl">jan.schoffelen@donders.ru.nl</a>><br>
<b>Reply-To:<span class="apple-converted-space"> </span></b>FieldTrip discussion list <<a href="mailto:fieldtrip@science.ru.nl">fieldtrip@science.ru.nl</a>><br>
<b>Date:<span class="apple-converted-space"> </span></b>Wednesday, 20 February 2019 at 2:59 AM<br>
<b>To:<span class="apple-converted-space"> </span></b>FieldTrip discussion list <<a href="mailto:fieldtrip@science.ru.nl">fieldtrip@science.ru.nl</a>><br>
<b>Subject:<span class="apple-converted-space"> </span></b>Re: [FieldTrip] Error in source localisation statistics</span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"> <o:p></o:p></p>
</div>
</div>
<div>
<p class="MsoNormal">Hi Xavier,<span class="apple-converted-space"> </span><o:p></o:p></p>
</div>
<div>
<div>
<p class="MsoNormal"> <o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal">My recommendation here would be to use the matlab debugger and inspect the cause of the error. The ‘Index exceeds array bounds’ is quite a clear error message, and suggests that there is a mismatch between the number of elements in dat/val,
and the numbers represented in indx(k).<o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"> <o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal">Good luck,<o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"> <o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal">Jan-Mathijs<o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"> <o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"> <o:p></o:p></p>
</div>
</div>
<div>
<div>
<blockquote style="margin-top:5.0pt;margin-bottom:5.0pt">
<div>
<div>
<p class="MsoNormal">On 18 Feb 2019, at 07:23, Xavier Vrijdag <<a href="mailto:x.vrijdag@auckland.ac.nz"><span style="color:purple">x.vrijdag@auckland.ac.nz</span></a>> wrote:<o:p></o:p></p>
</div>
</div>
<div>
<p class="MsoNormal"> <o:p></o:p></p>
</div>
<div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US">Hello,</span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US"> </span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US">I have a dataset from an experiment where 12 participants were exposed to 3 levels of nitrous oxide, while I measured their EEG. I also have a baseline recording. I took for each exposure (and baseline) a 1 minute sample
for further analysis. I want to do source localization to better understand which parts of the brain are involved. I don’t have MRI scans of my participants, so I used the standard MRI and BEM. I would like to thanks the persons who made the various manuals
and analysis examples on the website, as they helped me tremendously to do this analysis.</span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US"> </span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US">After cleaning and preprocessing, I combined the datasets into one for further analysis. I also chose to make 5 sec “trials” for further analysis. For this source analysis I loaded the data of 1 exposure and the baseline
data, to compare them.</span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US"> </span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span style="font-size:7.5pt;font-family:Courier;color:forestgreen"> %Calculate the channel covariance matric</span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span class="apple-converted-space"><span style="font-size:7.5pt;font-family:Courier"> </span></span><span style="font-size:7.5pt;font-family:Courier">cfg = [];</span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span class="apple-converted-space"><span style="font-size:7.5pt;font-family:Courier"> </span></span><span style="font-size:7.5pt;font-family:Courier">cfg.covariance =<span class="apple-converted-space"> </span><span class="s1"><span style="color:#A020F0">'yes'</span></span>;</span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span class="apple-converted-space"><span style="font-size:7.5pt;font-family:Courier"> </span></span><span style="font-size:7.5pt;font-family:Courier">cfg.covariancewindow =<span class="apple-converted-space"> </span><span class="s1"><span style="color:#A020F0">'all'</span></span>;</span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span class="apple-converted-space"><span style="font-size:7.5pt;font-family:Courier"> </span></span><span style="font-size:7.5pt;font-family:Courier">cfg.keeptrials<span class="apple-converted-space"> </span>=<span class="apple-converted-space"> </span><span class="s1"><span style="color:#A020F0">'yes'</span></span>;</span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span class="apple-converted-space"><span style="font-size:7.5pt;font-family:Courier"> </span></span><span style="font-size:7.5pt;font-family:Courier">timelock = ft_timelockanalysis(cfg, data_filt);</span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span class="apple-converted-space"><span style="font-size:7.5pt;font-family:Courier"> </span></span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span class="apple-converted-space"><span style="font-size:7.5pt;font-family:Courier;color:forestgreen"> </span></span><span style="font-size:7.5pt;font-family:Courier;color:forestgreen">%Global filter</span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span class="apple-converted-space"><span style="font-size:7.5pt;font-family:Courier"> </span></span><span style="font-size:7.5pt;font-family:Courier">cfg = [];</span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span class="apple-converted-space"><span style="font-size:7.5pt;font-family:Courier"> </span></span><span style="font-size:7.5pt;font-family:Courier">cfg.headmodel = vol;</span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span class="apple-converted-space"><span style="font-size:7.5pt;font-family:Courier"> </span></span><span style="font-size:7.5pt;font-family:Courier">cfg.elec = elec_aligned;</span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span class="apple-converted-space"><span style="font-size:7.5pt;font-family:Courier"> </span></span><span style="font-size:7.5pt;font-family:Courier">cfg.grid = grid;</span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span class="apple-converted-space"><span style="font-size:7.5pt;font-family:Courier"> </span></span><span style="font-size:7.5pt;font-family:Courier">cfg.method =<span class="apple-converted-space"> </span><span class="s1"><span style="color:#A020F0">'lcmv'</span></span>;</span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span class="apple-converted-space"><span style="font-size:7.5pt;font-family:Courier;color:forestgreen"> </span></span><span class="s2"><span style="font-size:7.5pt;font-family:Courier;color:forestgreen">cfg.lcmv.projectnoise=</span></span><span class="s1"><span style="font-size:7.5pt;font-family:Courier;color:#A020F0">'yes'</span></span><span class="s2"><span style="font-size:7.5pt;font-family:Courier;color:forestgreen">;</span></span><span class="apple-converted-space"><span style="font-size:7.5pt;font-family:Courier;color:forestgreen"> </span></span><span style="font-size:7.5pt;font-family:Courier;color:forestgreen">%needed
for neural activity index</span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span class="apple-converted-space"><span style="font-size:7.5pt;font-family:Courier;color:forestgreen"> </span></span><span class="s2"><span style="font-size:7.5pt;font-family:Courier;color:forestgreen">cfg.lcmv.fixedori =</span></span><span class="apple-converted-space"><span style="font-size:7.5pt;font-family:Courier;color:forestgreen"> </span></span><span class="s1"><span style="font-size:7.5pt;font-family:Courier;color:#A020F0">'yes'</span></span><span class="s2"><span style="font-size:7.5pt;font-family:Courier;color:forestgreen">;</span></span><span class="apple-converted-space"><span style="font-size:7.5pt;font-family:Courier;color:forestgreen"> </span></span><span style="font-size:7.5pt;font-family:Courier;color:forestgreen">%Project
onto largest variance orientation</span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span class="apple-converted-space"><span style="font-size:7.5pt;font-family:Courier"> </span></span><span style="font-size:7.5pt;font-family:Courier">cfg.lcmv.keepfilter =<span class="apple-converted-space"> </span><span class="s1"><span style="color:#A020F0">'yes'</span></span>;<span class="apple-converted-space"> </span><span class="s3"><span style="color:forestgreen">%Keep
the beamformer weights</span></span></span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span class="apple-converted-space"><span style="font-size:7.5pt;font-family:Courier"> </span></span><span style="font-size:7.5pt;font-family:Courier">cfg.lcmv.lambda =<span class="apple-converted-space"> </span><span class="s1"><span style="color:#A020F0">'5%'</span></span>;<span class="apple-converted-space"> </span><span class="s3"><span style="color:forestgreen">%Regularise
a little</span></span></span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span class="apple-converted-space"><span style="font-size:7.5pt;font-family:Courier;color:forestgreen"> </span></span><span class="s2"><span style="font-size:7.5pt;font-family:Courier;color:forestgreen">cfg.rawtrial</span></span><span class="apple-converted-space"><span style="font-size:7.5pt;font-family:Courier;color:forestgreen">
</span></span><span class="s2"><span style="font-size:7.5pt;font-family:Courier;color:forestgreen">=</span></span><span class="apple-converted-space"><span style="font-size:7.5pt;font-family:Courier;color:forestgreen"> </span></span><span class="s1"><span style="font-size:7.5pt;font-family:Courier;color:#A020F0">'yes'</span></span><span class="s2"><span style="font-size:7.5pt;font-family:Courier;color:forestgreen">;</span></span><span class="apple-converted-space"><span style="font-size:7.5pt;font-family:Courier;color:forestgreen">
</span></span><span style="font-size:7.5pt;font-family:Courier;color:forestgreen">% project each single trial through the filter.<span class="apple-converted-space"> </span></span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span class="apple-converted-space"><span style="font-size:7.5pt;font-family:Courier"> </span></span><span style="font-size:7.5pt;font-family:Courier">cfg.keeptrials<span class="apple-converted-space"> </span>=<span class="apple-converted-space"> </span><span class="s1"><span style="color:#A020F0">'yes'</span></span>;</span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span class="apple-converted-space"><span style="font-size:7.5pt;font-family:Courier"> </span></span><span style="font-size:7.5pt;font-family:Courier">source = ft_sourceanalysis(cfg, timelock);</span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span class="apple-converted-space"><span style="font-size:7.5pt;font-family:Courier"> </span></span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span class="apple-converted-space"><span style="font-size:7.5pt;font-family:Courier"> </span></span><span style="font-size:7.5pt;font-family:Courier">cfg = [];</span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span class="apple-converted-space"><span style="font-size:7.5pt;font-family:Courier"> </span></span><span style="font-size:7.5pt;font-family:Courier">cfg.dim<span class="apple-converted-space"> </span>= source.dim;</span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span class="apple-converted-space"><span style="font-size:7.5pt;font-family:Courier"> </span></span><span style="font-size:7.5pt;font-family:Courier">cfg.method<span class="apple-converted-space"> </span>=<span class="apple-converted-space"> </span><span class="s1"><span style="color:#A020F0">'montecarlo'</span></span>;</span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span class="apple-converted-space"><span style="font-size:7.5pt;font-family:Courier"> </span></span><span style="font-size:7.5pt;font-family:Courier">cfg.statistic<span class="apple-converted-space"> </span>=<span class="apple-converted-space"> </span><span class="s1"><span style="color:#A020F0">'ft_statfun_indepsamplesT'</span></span>;</span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span class="apple-converted-space"><span style="font-size:7.5pt;font-family:Courier"> </span></span><span style="font-size:7.5pt;font-family:Courier">cfg.parameter<span class="apple-converted-space"> </span>=<span class="apple-converted-space"> </span><span class="s1"><span style="color:#A020F0">'pow'</span></span>;</span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span class="apple-converted-space"><span style="font-size:7.5pt;font-family:Courier"> </span></span><span style="font-size:7.5pt;font-family:Courier">cfg.correctm<span class="apple-converted-space"> </span>=<span class="apple-converted-space"> </span><span class="s1"><span style="color:#A020F0">'cluster'</span></span>;</span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span class="apple-converted-space"><span style="font-size:7.5pt;font-family:Courier"> </span></span><span style="font-size:7.5pt;font-family:Courier">cfg.numrandomization = 1000;</span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span class="apple-converted-space"><span style="font-size:7.5pt;font-family:Courier;color:forestgreen"> </span></span><span class="s2"><span style="font-size:7.5pt;font-family:Courier;color:forestgreen">cfg.alpha</span></span><span class="apple-converted-space"><span style="font-size:7.5pt;font-family:Courier;color:forestgreen">
</span></span><span class="s2"><span style="font-size:7.5pt;font-family:Courier;color:forestgreen">= 0.05;</span></span><span class="apple-converted-space"><span style="font-size:7.5pt;font-family:Courier;color:forestgreen"> </span></span><span style="font-size:7.5pt;font-family:Courier;color:forestgreen">%
note that this only implies single-sided testing</span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span class="apple-converted-space"><span style="font-size:7.5pt;font-family:Courier"> </span></span><span style="font-size:7.5pt;font-family:Courier">cfg.tail<span class="apple-converted-space"> </span>= 0;</span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span class="apple-converted-space"><span style="font-size:7.5pt;font-family:Courier"> </span></span><span style="font-size:7.5pt;font-family:Courier">cfg.design(1,:) = [1:nTrials 1:nTrials];</span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span class="apple-converted-space"><span style="font-size:7.5pt;font-family:Courier"> </span></span><span style="font-size:7.5pt;font-family:Courier">cfg.design(2,:) = [ones(1,nTrials) ones(1,nTrials).*2];</span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span class="apple-converted-space"><span style="font-size:7.5pt;font-family:Courier;color:forestgreen"> </span></span><span class="s2"><span style="font-size:7.5pt;font-family:Courier;color:forestgreen">cfg.uvar</span></span><span class="apple-converted-space"><span style="font-size:7.5pt;font-family:Courier;color:forestgreen">
</span></span><span class="s2"><span style="font-size:7.5pt;font-family:Courier;color:forestgreen">= 1;</span></span><span class="apple-converted-space"><span style="font-size:7.5pt;font-family:Courier;color:forestgreen"> </span></span><span style="font-size:7.5pt;font-family:Courier;color:forestgreen">%
row of design matrix that contains unit variable (in this case: trials)</span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span class="apple-converted-space"><span style="font-size:7.5pt;font-family:Courier;color:forestgreen"> </span></span><span class="s2"><span style="font-size:7.5pt;font-family:Courier;color:forestgreen">cfg.ivar</span></span><span class="apple-converted-space"><span style="font-size:7.5pt;font-family:Courier;color:forestgreen">
</span></span><span class="s2"><span style="font-size:7.5pt;font-family:Courier;color:forestgreen">= 2;</span></span><span class="apple-converted-space"><span style="font-size:7.5pt;font-family:Courier;color:forestgreen"> </span></span><span style="font-size:7.5pt;font-family:Courier;color:forestgreen">%
row of design matrix that contains independent variable (the conditions)<span class="apple-converted-space"> </span></span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span class="apple-converted-space"><span style="font-size:7.5pt;font-family:Courier"> </span></span><span style="font-size:7.5pt;font-family:Courier">stat = ft_sourcestatistics(cfg, source);</span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US"> </span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US">After the source analysis the source statistics script calls the checkdata script and I get the following error:</span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US"> </span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US">Index exceeds array bounds.</span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US"> </span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US">Error in ft_datatype_source (line 183)</span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US"> val{indx(k)}(1,:,:,:) = dat{indx(k)};</span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US"> </span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US">Error in ft_checkdata (line 278)</span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US"> data = ft_datatype_source(data);</span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US"> </span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US">Error in ft_sourcestatistics (line 94)</span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US"> varargin{i} = ft_checkdata(varargin{i}, 'datatype', 'source', 'feedback', 'no');</span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US"> </span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US">Error in sourcelocalisation_N2O (line 114)</span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US"> stat = ft_sourcestatistics(cfg, source);</span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US"> </span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US">The source struct has the following content:</span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US"> time: [1×15360 double]</span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US"> dim: [24 31 25]</span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US"> inside: [18600×1 logical]</span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US"> pos: [18600×3 double]</span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US"> method: 'rawtrial'</span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US"> trial: [1×24 struct]</span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US"> df: 24</span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US"> cfg: [1×1 struct]</span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US"> </span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US">With source.trial:</span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US"> 1×24 struct array with fields:</span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US"> ori</span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US"> pow</span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US"> mom</span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US"> noise</span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US"> filter</span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US"> label</span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US"> filterdimord</span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US"> </span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US">I would like your help to understand the error and how to fix it.</span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US"> </span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US">Regards,</span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US"> </span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US">Xavier Vrijdag</span><o:p></o:p></p>
</div>
</div>
<div>
<div>
<p class="MsoNormal"><span lang="EN-US"> </span><o:p></o:p></p>
</div>
</div>
<div>
<p class="MsoNormal"><span style="font-size:10.5pt;font-family:Helvetica">_______________________________________________<br>
fieldtrip mailing list<br>
</span><a href="https://mailman.science.ru.nl/mailman/listinfo/fieldtrip"><span style="font-size:10.5pt;font-family:Helvetica;color:#954F72">https://mailman.science.ru.nl/mailman/listinfo/fieldtrip</span></a><span style="font-size:10.5pt;font-family:Helvetica"><br>
</span><a href="https://doi.org/10.1371/journal.pcbi.1002202"><span style="font-size:10.5pt;font-family:Helvetica;color:#954F72">https://doi.org/10.1371/journal.pcbi.1002202</span></a><o:p></o:p></p>
</div>
</div>
</blockquote>
</div>
<div>
<p class="MsoNormal"> <o:p></o:p></p>
</div>
</div>
<p class="MsoNormal"><span style="font-size:10.5pt;font-family:Helvetica">_______________________________________________<br>
fieldtrip mailing list<br>
<a href="https://mailman.science.ru.nl/mailman/listinfo/fieldtrip">https://mailman.science.ru.nl/mailman/listinfo/fieldtrip</a><br>
<a href="https://doi.org/10.1371/journal.pcbi.1002202">https://doi.org/10.1371/journal.pcbi.1002202</a></span><o:p></o:p></p>
</div>
</blockquote>
</div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
</div>
<p class="MsoNormal" style="margin-bottom:12.0pt"><o:p> </o:p></p>
<div>
<div>
<div>
<div>
<p class="MsoNormal"><span style="color:black"><o:p> </o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="color:black"><o:p> </o:p></span></p>
</div>
</div>
</div>
</div>
<p class="MsoNormal"><br>
<br>
<o:p></o:p></p>
</div>
</body>
</html>