<div dir="ltr">Dear Gaëlle,<div><br></div><div>You are using default settings for downsampling, which will then try to downsample to 256Hz. Your original samplerate is 1000Hz, so that isn't a neat integer division.</div><div>I tried the following, and could reproduce your problem:</div><div><br></div><div>Fs = 1000;<br>data_eeg = [];<br>data_eeg.label = {'a','b','c'};<br>data_eeg.time{1} = 0:1/Fs:10;<br>data_eeg.trial{1} = rand(length(data_eeg.label),length(data_eeg.time{1}));<br>data_eeg.Fs = Fs;<br>data_eeg.sampleinfo = [1, length(data_eeg.time{1})];<br><br><div>% Doesn't work:<br>cfg = [];<br>[data_eeg] = ft_resampledata(cfg, data_eeg); <br></div><div><br></div><div></div>% Does work:</div><div>cfg = [];<br>cfg.resamplefs = 250;<br>[data_eeg] = ft_resampledata(cfg, data_eeg); <br></div><div><br></div><div>So that deals with your first problem. Not sure why it has a problem with the downsampling to non-integer values, but it's good practice and much faster to stick to that anyway.</div><div><br></div><div>Not sure what to do with the artefact detection though, or what you expect;  it seems you are running it with some contradictory settings/ways perhaps:</div><div>- use cfg.dataset if reading from disk directly OR give datastructure to function is you already have data in memory</div><div>- use cfg.continuous when data is not segmented</div><div>- the comment "jump" suggests you use settings for SQUID jumps, i.e. optimized for MEG, but that should not give an error :-)</div><div></div><div><br></div><div>I suspect you probably want to run artefact detection on your segmented data, so remove those references to the data on disk (datafile, headerfile), and the trl from the cfg, and give your segmented data as an argument as you did here:</div><div>%[cfg, artifact_jump] = ft_artifact_zvalue(cfg, data_eeg); %no plot error, but each trial has the same graph, also error using zeros<br></div><div><br></div><div>Let us know how it goes,</div><div><br></div><div>All the best,</div><div>Stephen</div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">Op ma 22 jun. 2020 om 17:02 schreef Gaëlle Leys <<a href="mailto:gaelle.leys@kuleuven.be">gaelle.leys@kuleuven.be</a>>:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Dear Stephen, <br>
<br>
I have included the data structures and cfg's involved in the functions.<br>
<br>
1)      A. data structure for resampling data (raw data after reading it in):<br>
data_eeg = <br>
  struct with fields:<br>
hdr: [1×1 struct]<br>
label: {129×1 cell}<br>
time: {[1×245726 double]}<br>
trial: {[129×245726 double]}<br>
fsample: 1000<br>
sampleinfo: [1 245726]<br>
elec: [1×1 struct]<br>
cfg: [1×1 struct]<br>
<br>
1)      B. cfg for resampling data (I’m using all the default options for resampling):<br>
<br>
Cfg=[];<br>
[data_eeg] = ft_resampledata(cfg, data_eeg); % defaults are all OK<br>
<br>
2)      A. data structure for artifact rejection (my data already epoched):<br>
data_eeg = <br>
struct with fields:<br>
hdr: [1×1 struct]<br>
label: {129×1 cell}<br>
 time: {1×108 cell}<br>
 trial: {1×108 cell}<br>
 fsample: 1000<br>
 sampleinfo: [108×2 double]<br>
 elec: [1×1 struct]<br>
 cfg: [1×1 struct]<br>
<br>
<br>
2)      B. cfg for artifact rejection (taken from the automatic artifact rejection example):<br>
% jump<br>
cfg = [];<br>
cfg.trl = trl;<br>
cfg.datafile = '100_1_20200305.mff';<br>
cfg.headerfile = '100_1_20200305.mff';<br>
cfg.continuous = 'yes';<br>
<br>
% channel selection, cutoff and padding<br>
cfg.artfctdef.zvalue.channel = 'EEG';<br>
cfg.artfctdef.zvalue.cutoff = 20;<br>
cfg.artfctdef.zvalue.trlpadding = 0;<br>
cfg.artfctdef.zvalue.artpadding = 0;<br>
cfg.artfctdef.zvalue.fltpadding = 0;<br>
<br>
% algorithmic parameters<br>
cfg.artfctdef.zvalue.cumulative = 'yes';<br>
cfg.artfctdef.zvalue.medianfilter = 'yes';<br>
cfg.artfctdef.zvalue.medianfiltord = 9;<br>
cfg.artfctdef.zvalue.absdiff = 'yes';<br>
<br>
% make the process interactive<br>
cfg.artfctdef.zvalue.interactive = 'yes';<br>
<br>
[cfg, artifact_jump] = ft_artifact_zvalue(cfg); %plot error and error using zeros<br>
%[cfg, artifact_jump] = ft_artifact_zvalue(cfg, data_eeg); %no plot error, but each trial has the same graph, also error using zeros<br>
<br>
<br>
Kind regards, <br>
Gaëlle Leys<br>
<br>
<br>
<br>
-----Original Message-----<br>
From: fieldtrip <<a href="mailto:fieldtrip-bounces@science.ru.nl" target="_blank">fieldtrip-bounces@science.ru.nl</a>> On Behalf Of <a href="mailto:fieldtrip-request@science.ru.nl" target="_blank">fieldtrip-request@science.ru.nl</a><br>
Sent: zaterdag 20 juni 2020 12:00<br>
To: <a href="mailto:fieldtrip@science.ru.nl" target="_blank">fieldtrip@science.ru.nl</a><br>
Subject: fieldtrip Digest, Vol 115, Issue 22<br>
<br>
Send fieldtrip mailing list submissions to<br>
        <a href="mailto:fieldtrip@science.ru.nl" target="_blank">fieldtrip@science.ru.nl</a><br>
<br>
To subscribe or unsubscribe via the World Wide Web, visit<br>
        <a href="https://mailman.science.ru.nl/mailman/listinfo/fieldtrip" rel="noreferrer" target="_blank">https://mailman.science.ru.nl/mailman/listinfo/fieldtrip</a><br>
or, via email, send a message with subject or body 'help' to<br>
        <a href="mailto:fieldtrip-request@science.ru.nl" target="_blank">fieldtrip-request@science.ru.nl</a><br>
<br>
You can reach the person managing the list at<br>
        <a href="mailto:fieldtrip-owner@science.ru.nl" target="_blank">fieldtrip-owner@science.ru.nl</a><br>
<br>
When replying, please edit your Subject line so it is more specific than "Re: Contents of fieldtrip digest..."<br>
<br>
<br>
Today's Topics:<br>
<br>
   1. Errors with downsampling and artifact rejection (Gaëlle Leys)<br>
   2. Re: Errors with downsampling and artifact rejection<br>
      (Stephen Whitmarsh)<br>
   3. MNI coordinates for the EEG easycap-M1 layout<br>
      (SILVA PEREIRA, SILVANA)<br>
<br>
<br>
----------------------------------------------------------------------<br>
<br>
Message: 1<br>
Date: Fri, 19 Jun 2020 12:48:26 +0000<br>
From: Gaëlle Leys <<a href="mailto:gaelle.leys@kuleuven.be" target="_blank">gaelle.leys@kuleuven.be</a>><br>
To: "<a href="mailto:fieldtrip@science.ru.nl" target="_blank">fieldtrip@science.ru.nl</a>" <<a href="mailto:fieldtrip@science.ru.nl" target="_blank">fieldtrip@science.ru.nl</a>><br>
Subject: [FieldTrip] Errors with downsampling and artifact rejection<br>
Message-ID:<br>
        <<a href="mailto:ccb89b4d31af42d3b4cd22485438f4ef@ICTS-S-EXMBX17.luna.kuleuven.be" target="_blank">ccb89b4d31af42d3b4cd22485438f4ef@ICTS-S-EXMBX17.luna.kuleuven.be</a>><br>
Content-Type: text/plain; charset="iso-8859-1"<br>
<br>
Dear Fieldtrip community,<br>
<br>
I have recently started with Fieldtrip for my EEG data processing and I am encountering some problems which I can't seem to resolve.<br>
I'm working with mff files from an EGI EEG system. I have no problem reading in and preprocessing the data. However, when trying to downsample and setup the automatic artifact removal, I'm getting some errors.<br>
<br>
<br>
1)      When trying to downsample with the default options, these are the errors that I get:<br>
Error using resample>getSamples (line 357) The number of elements of Tx must match the number of rows of X when X is a matrix<br>
<br>
Error in resample>nonUniformResample (line 123) [x, tx] = getSamples(varargin{1:2});<br>
<br>
Error in resample (line 116)<br>
  [varargout{1:max(1,nargout)}] = nonUniformResample(method,varargin{:});<br>
<br>
Error in ft_resampledata (line 228)<br>
        newdat = transpose(resample(transpose(olddat),fsres,fsorig));<br>
<br>
                I have tried downsampling before filtering, after filtering, on continuous data, or already segmented data. None of these things seem to resolve the error.<br>
<br>
<br>
2)      When trying to set up the automatic artifact removal process and choosing yes for the interactive plot, I get a plotting error:<br>
<br>
Error using plot<br>
<br>
Vectors must be the same length.<br>
<br>
<br>
<br>
Error in ft_artifact_zvalue>redraw_cb (line 973)<br>
<br>
    plot(opt.h1, xval, yval, 'linestyle', '-', 'color', 'b', 'displayname', 'data');<br>
<br>
<br>
<br>
Error in ft_artifact_zvalue (line 512)<br>
<br>
  redraw_cb(h);<br>
<br>
        I've noticed that the xval and yval are different numbers, which of course doesn't allow plotting. It seems like the xval changes according to the epoch length, however the yval remains the same and equals the amount of samples in my continuous data.<br>
<br>
<br>
3)      When choosing no for the interactive plot, I also get an error:<br>
Error using zeros<br>
Size inputs must be integers.<br>
<br>
Error in ft_artifact_zvalue (line 568)<br>
dum = zeros(1,max(opt.trl(:,2)));<br>
<br>
I've been looking into the functions and where it goes wrong, but I can't seem to pinpoint what exactly is causing these issues. Anyone here that can help or provide me some tips?<br>
<br>
Thanks in advance!<br>
<br>
Kind regards,<br>
<br>
Gaëlle Leys<br>
<br>
Research Associate<br>
Brain & Cognition<br>
Tiensestraat 102 box 3714<br>
3000 Leuven<br>
Room PSI 02.33<br>
tel. +32 16 37 61 98<br>
<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://mailman.science.ru.nl/pipermail/fieldtrip/attachments/20200619/e7e1bde8/attachment-0001.htm" rel="noreferrer" target="_blank">http://mailman.science.ru.nl/pipermail/fieldtrip/attachments/20200619/e7e1bde8/attachment-0001.htm</a>><br>
<br>
------------------------------<br>
<br>
Message: 2<br>
Date: Fri, 19 Jun 2020 15:08:21 +0200<br>
From: Stephen Whitmarsh <<a href="mailto:stephen.whitmarsh@gmail.com" target="_blank">stephen.whitmarsh@gmail.com</a>><br>
To: FieldTrip discussion list <<a href="mailto:fieldtrip@science.ru.nl" target="_blank">fieldtrip@science.ru.nl</a>><br>
Subject: Re: [FieldTrip] Errors with downsampling and artifact<br>
        rejection<br>
Message-ID:<br>
        <CAFrxm=<a href="mailto:w44QZ0C7xB4cpsqdy0D0ou5UcpHNM-q3bjyaqwwJxRwg@mail.gmail.com" target="_blank">w44QZ0C7xB4cpsqdy0D0ou5UcpHNM-q3bjyaqwwJxRwg@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
Dear Gaëlle,<br>
<br>
To understand the problem it would be helpful to have 1) description of the data structures going in ft_resampledata (just type in command line and copy-paste the output) and the artefact detection function you are using, and 2) the configuration (cfg) you are using for those functions.<br>
<br>
Cheers,<br>
Stephen<br>
<br>
<br>
Op vr 19 jun. 2020 om 14:56 schreef Gaëlle Leys <<a href="mailto:gaelle.leys@kuleuven.be" target="_blank">gaelle.leys@kuleuven.be</a>>:<br>
<br>
> Dear Fieldtrip community,<br>
><br>
><br>
><br>
> I have recently started with Fieldtrip for my EEG data processing and <br>
> I am encountering some problems which I can’t seem to resolve.<br>
><br>
> I’m working with mff files from an EGI EEG system. I have no problem <br>
> reading in and preprocessing the data. However, when trying to <br>
> downsample and setup the automatic artifact removal, I’m getting some errors.<br>
><br>
><br>
><br>
> 1)      When trying to downsample with the default options, these are the<br>
> errors that I get:<br>
><br>
> Error using resample>getSamples (line 357)<br>
><br>
> The number of elements of Tx must match the number of rows of X when X <br>
> is a matrix<br>
><br>
><br>
><br>
> Error in resample>nonUniformResample (line 123)<br>
><br>
> [x, tx] = getSamples(varargin{1:2});<br>
><br>
><br>
><br>
> Error in resample (line 116)<br>
><br>
>   [varargout{1:max(1,nargout)}] = <br>
> nonUniformResample(method,varargin{:});<br>
><br>
><br>
><br>
> Error in ft_resampledata (line 228)<br>
><br>
>         newdat = transpose(resample(transpose(olddat),fsres,fsorig));<br>
><br>
><br>
><br>
>                 I have tried downsampling before filtering, after <br>
> filtering, on continuous data, or already segmented data. None of <br>
> these things seem to resolve the error.<br>
><br>
><br>
><br>
> 2)      When trying to set up the automatic artifact removal process and<br>
> choosing yes for the interactive plot, I get a plotting error:<br>
><br>
> Error using plot<br>
><br>
> Vectors must be the same length.<br>
><br>
><br>
><br>
> Error in ft_artifact_zvalue>redraw_cb (line 973)<br>
><br>
>     plot(opt.h1, xval, yval, 'linestyle', '-', 'color', 'b', <br>
> 'displayname', 'data');<br>
><br>
><br>
><br>
> Error in ft_artifact_zvalue (line 512)<br>
><br>
>   redraw_cb(h);<br>
><br>
><br>
><br>
>         I’ve noticed that the xval and yval are different numbers, <br>
> which of course doesn’t allow plotting. It seems like the xval changes <br>
> according to the epoch length, however the yval remains the same and <br>
> equals the amount of samples in my continuous data.<br>
><br>
><br>
><br>
> 3)      When choosing no for the interactive plot, I also get an error:<br>
><br>
> Error using zeros<br>
><br>
> Size inputs must be integers.<br>
><br>
><br>
><br>
> Error in ft_artifact_zvalue (line 568)<br>
><br>
> dum = zeros(1,max(opt.trl(:,2)));<br>
><br>
><br>
><br>
> I’ve been looking into the functions and where it goes wrong, but I <br>
> can’t seem to pinpoint what exactly is causing these issues. Anyone <br>
> here that can help or provide me some tips?<br>
><br>
><br>
><br>
> Thanks in advance!<br>
><br>
><br>
><br>
> Kind regards,<br>
><br>
><br>
><br>
> Gaëlle Leys<br>
><br>
><br>
><br>
> Research Associate<br>
> Brain & Cognition<br>
> Tiensestraat 102 box 3714<br>
> 3000 Leuven<br>
><br>
> Room PSI 02.33<br>
> tel. +32 16 37 61 98<br>
><br>
><br>
> _______________________________________________<br>
> fieldtrip mailing list<br>
> <a href="https://mailman.science.ru.nl/mailman/listinfo/fieldtrip" rel="noreferrer" target="_blank">https://mailman.science.ru.nl/mailman/listinfo/fieldtrip</a><br>
> <a href="https://doi.org/10.1371/journal.pcbi.1002202" rel="noreferrer" target="_blank">https://doi.org/10.1371/journal.pcbi.1002202</a><br>
><br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://mailman.science.ru.nl/pipermail/fieldtrip/attachments/20200619/137e31b1/attachment-0001.htm" rel="noreferrer" target="_blank">http://mailman.science.ru.nl/pipermail/fieldtrip/attachments/20200619/137e31b1/attachment-0001.htm</a>><br>
<br>
------------------------------<br>
<br>
Message: 3<br>
Date: Fri, 19 Jun 2020 18:58:36 +0200<br>
From: "SILVA PEREIRA, SILVANA" <<a href="mailto:silvana.silva@upf.edu" target="_blank">silvana.silva@upf.edu</a>><br>
To: <a href="mailto:fieldtrip@science.ru.nl" target="_blank">fieldtrip@science.ru.nl</a><br>
Subject: [FieldTrip] MNI coordinates for the EEG easycap-M1 layout<br>
Message-ID:<br>
        <<a href="mailto:CA%2B_MWm1NufDWa-fjR71OzQ5ab88ioyUT8OJ-QFYX9ABr7gNQrA@mail.gmail.com" target="_blank">CA+_MWm1NufDWa-fjR71OzQ5ab88ioyUT8OJ-QFYX9ABr7gNQrA@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
Dear community,<br>
<br>
I've been struggling to find the 3D positions (montages) for the<br>
Acticap-64ch-standard2 on the web, which is a modified version of the easycap-M1, but I have not succeeded. In the Fieldtrip page <a href="http://www.fieldtriptoolbox.org/template/electrode/" rel="noreferrer" target="_blank">http://www.fieldtriptoolbox.org/template/electrode/</a> you have the following:<br>
<br>
You can find the template 3-D electrode sets included in FieldTrip here <<a href="https://github.com/fieldtrip/fieldtrip/tree/master/template/electrode" rel="noreferrer" target="_blank">https://github.com/fieldtrip/fieldtrip/tree/master/template/electrode</a>>.<br>
<br>
But in that link you only find a txt file with theta-phi coordinates (easycap-M1.txt). Is there any place where I could find these positions in MNI space?<br>
<br>
Thank you and best regards,<br>
Silvana<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <<a href="http://mailman.science.ru.nl/pipermail/fieldtrip/attachments/20200619/c0adf7df/attachment-0001.htm" rel="noreferrer" target="_blank">http://mailman.science.ru.nl/pipermail/fieldtrip/attachments/20200619/c0adf7df/attachment-0001.htm</a>><br>
<br>
------------------------------<br>
<br>
Subject: Digest Footer<br>
<br>
_______________________________________________<br>
fieldtrip mailing list<br>
<a href="https://mailman.science.ru.nl/mailman/listinfo/fieldtrip" rel="noreferrer" target="_blank">https://mailman.science.ru.nl/mailman/listinfo/fieldtrip</a><br>
_______________________________________________<br>
fieldtrip mailing list<br>
<a href="https://mailman.science.ru.nl/mailman/listinfo/fieldtrip" rel="noreferrer" target="_blank">https://mailman.science.ru.nl/mailman/listinfo/fieldtrip</a><br>
<a href="https://doi.org/10.1371/journal.pcbi.1002202" rel="noreferrer" target="_blank">https://doi.org/10.1371/journal.pcbi.1002202</a><br>
<br>
<br>
------------------------------<br>
<br>
End of fieldtrip Digest, Vol 115, Issue 22<br>
******************************************<br>
<br>
_______________________________________________<br>
fieldtrip mailing list<br>
<a href="https://mailman.science.ru.nl/mailman/listinfo/fieldtrip" rel="noreferrer" target="_blank">https://mailman.science.ru.nl/mailman/listinfo/fieldtrip</a><br>
<a href="https://doi.org/10.1371/journal.pcbi.1002202" rel="noreferrer" target="_blank">https://doi.org/10.1371/journal.pcbi.1002202</a><br>
</blockquote></div>