Back to the main page.
Bug 46 - ft_connectivityanalysis coh/psi bug
| Status | CLOSED FIXED |
| 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 |
| URL: | |
| Tags: | |
| Depends on: | |
| Blocks: | |
| 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