Back to the main page.

Bug 46 - ft_connectivityanalysis coh/psi bug

Reported 2010-03-04 12:01:00 +0100
Modified 2010-03-18 10:19:17 +0100
Product: FieldTrip
Component: core
Version: unspecified
Hardware: PC
Operating System: Linux
Importance: P5 normal
Assigned to: Jan-Mathijs Schoffelen
Depends on:
See also:

Andre Bastos - 2010-03-04 12:01:19 +0100

Here is the error: ??? Subscript indices must either be real positive integers or logicals. Error in ==> checkdata>fixcsd at 882 tmpdat1 = data.fourierspctrm(indx,cmbindx(:,1),:,:); Error in ==> checkdata at 575 data = fixcsd(data, cmbrepresentation, channelcmb); Error in ==> ft_connectivityanalysis>univariate2bivariate at 992 data = checkdata(data, 'cmbrepresentation', 'sparse', 'channelcmb', cmb); Error in ==> ft_connectivityanalysis at 138 [data, powindx, hasrpt] = univariate2bivariate(data, 'fourierspctrm', 'crsspctrm', dtype, 0, cfg.channelcmb); Error in ==> run_all_measures at 43 psi1 = ft_connectivityanalysis(cfg, freq); freq is: freq = label: {2x1 cell} dimord: 'rpttap_chan_freq' freq: [1x121 double] fourierspctrm: [7500x2x121 double] cumsumcnt: [500x1 double] cumtapcnt: [500x1 double] cfg: [1x1 struct] cfg is: cfg = method: 'psi' bandwidth: 4 jackknife: 'no' The same error also occurs for cfg = method: 'coh' complex: 'abs'

Jan-Mathijs Schoffelen - 2010-03-04 13:03:01 +0100

Could reproduce the error. Uncommenting lines 124-126 fixes the problem. Don't remember having put those lines as a comment in the first place. Could it have to do with your own previous adjustment? If so, please revisit that one. For now, please update and verify that it works for me. Let me know through bugzilla so that I can close this bug.

Jan-Mathijs Schoffelen - 2010-03-04 13:04:59 +0100

You as the quality controller should verify whether the fix works for you, then you can mark the bug as verified (including me in cc list). I will then close the bug. This is the procedure we discussed in yesterday's FT-meeting

Jan-Mathijs Schoffelen - 2010-03-09 08:41:29 +0100

please follow up on bug, i.e. check whether it is OK with you and if so, mark it as verified