[FieldTrip] Error in ft_mvaranalysis

[Schoffelen, J.M. (Jan Mathijs)]

Hi Manuela,

The code in ft_mvaranalysis (at least the part that is supposed to support the time-resolved models) is a bit messy. I patched it for now, making the function a bit more robust with respect to rounding off errors in the conversion between sample indices and time windows. I hope that the newer version now works for you.

Best wishes,
Jan-Mathijs


On 4 Sep 2018, at 10:42, Manuela Costa <manuela.costa at ctb.upm.es<mailto:manuela.costa at ctb.upm.es>> wrote:

Dear Jan-Mathijs,

thank you very much for fixing the typo in the code.
I do not get anymore the same error but as you predicted I get a new one.
I hope you can help me with that as well.

Best wishes
Manuela


Matrix dimensions must agree.

Error in ft_mvaranalysis>catnan (line 607)
    datamatrix(:,begsmp:endsmp) = datacells{trials(k)}(chanindx,:).*taper(ones(nchan,1),:);

Error in ft_mvaranalysis (line 425)
dat = catnan(tmpdata.trial, chanindx, rpt{rlop}, tap(m,:), nnans, dobvar);

Error in Step2_Calculate_Granger_time_freq_updated_MC (line 124)
    mdata_emotional_Rem{subj} = ft_mvaranalysis(cfg, preprocessing_emotional_Rem_ds{subj});


I copy here the information about the structure (it's the same as before)

The data is downsampled at 250 Hz. The fields of my data structure are fsample, trial, time, label and cfg. The size of preprocessing_emotional_Rem_ds{1}.trial{1} is 3x3751 (channels, data points).

Configuration structure:

cfg = [];
    cfg.toolbox = 'bsmart';
    cfg.order   = 10;
    cfg.toi       = -0.5:0.01:3; %% the time points at which the windows are centered
    cfg.t_ftimwin = 0.25; %% length of time window in second
    mdata_emotional_Rem{subj} = ft_mvaranalysis(cfg, preprocessing_emotional_Rem_ds{subj});


2018-08-31 11:32 GMT+02:00 Schoffelen, J.M. (Jan Mathijs) <jan.schoffelen at donders.ru.nl<mailto:jan.schoffelen at donders.ru.nl>>:
Dear Manuela,

This looks like a typo in the code, which I will fix (and once fixed will become available in the daily release version of the code). To give myself some time: probably not tonight, but perhaps after the weekend :). However, the fact that this part of the code clearly doesn’t work, I am not entirely sure whether the rest will work smoothly as well, as soon as the current hurdle is fixed. But let’s see

Best wishes,
Jan-Mathijs


On 29 Aug 2018, at 13:34, Manuela Costa <manuela.costa at ctb.upm.es<mailto:manuela.costa at ctb.upm.es>> wrote:

Dear community,

I'm planning to do granger causality analysis on human iEEG data. I'm using ft_mvaranalysis to get time-dependent coefficient sets based on a sliding window, but I get an error.

The data is downsampled at 250 Hz. The fields of my data structure are fsample, trial, time, label and cfg. The size of preprocessing_emotional_Rem_ds{1}.trial{1} is 3x3751 (channels, data points). Below is the configuration structure and the error message. The ft_mvaranalysis function works if I do not use the sliding window approach. I tried to use different cfg.toi and cfg.t_ftimwin but I could not find a solution. I hope you can help me.

Best wishes,

Manuela Costa

Configuration structure:

cfg = [];
    cfg.toolbox = 'bsmart';
    cfg.order   = 10;
    cfg.toi       = -0.5:0.01:3; %% the time points at which the windows are centered
    cfg.t_ftimwin = 0.25; %% length of time window in second
    mdata_emotional_Rem{subj} = ft_mvaranalysis(cfg, preprocessing_emotional_Rem_ds{subj});




Error message:


Matrix dimensions must agree.

Error in ft_mvaranalysis (line 185)
latency(k,:) = cfg.toi + cfg.t_ftimwin.*[-0.5 0.5];

Error in Step2_Calculate_Granger_time_freq(line 108)
    mdata_emotional_Rem{subj} = ft_mvaranalysis(cfg,preprocessing_emotional_Rem_ds{subj});



_______________________________________________
fieldtrip mailing list
https://mailman.science.ru.nl/mailman/listinfo/fieldtrip
https://doi.org/10.1371/journal.pcbi.1002202


_______________________________________________
fieldtrip mailing list
https://mailman.science.ru.nl/mailman/listinfo/fieldtrip
https://doi.org/10.1371/journal.pcbi.1002202


_______________________________________________
fieldtrip mailing list
https://mailman.science.ru.nl/mailman/listinfo/fieldtrip
https://doi.org/10.1371/journal.pcbi.1002202

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.science.ru.nl/pipermail/fieldtrip/attachments/20180904/21010322/attachment-0001.html>