Back to the main page.
Bug 734 - ft_sourceanalysis: method ='sam', does not work
| Status | ASSIGNED |
| Reported | 2011-06-07 09:37:00 +0200 |
| Modified | 2011-11-14 22:33:20 +0100 |
| Product: | FieldTrip |
| Component: | core |
| Version: | unspecified |
| Hardware: | Macintosh |
| Operating System: | Mac OS |
| Importance: | P1 major |
| Assigned to: | Robert Oostenveld |
| URL: | |
| Tags: | |
| Depends on: | |
| Blocks: | |
| See also: |
Elena Orekhova - 2011-06-07 09:37:12 +0200
In my data the ft_sourceanalysis works with cfg.method = 'lcmv', but it does not work if cfg.method ='sam'. I used the yesterday's version of FT. The matlab output is below: cfg = grad: [1x1 struct] method: 'sam' grid: [1x1 struct] vol: [1x1 struct] reducerank: 2 keepfilter: 'yes' sourcepst = ft_sourceanalysis(cfg, tlckavgpst); Warning: multiple versions of SPM on your path will confuse FieldTrip > In fieldtrip-20110606/private/warning_once at 75 In ft_defaults at 45 In ft_sourceanalysis at 197 Warning: one version of SPM is found here: /Applications/Fieldtrip/fieldtrip-20110606/external/spm8/spm.m > In ft_defaults at 48 In ft_sourceanalysis at 197 Warning: one version of SPM is found here: /Applications/Fieldtrip/fieldtrip-20110606/external/spm2/spm.m > In ft_defaults at 48 In ft_sourceanalysis at 197 Warning: one version of SPM is found here: /Users/data/Programs/spm8/spm.m > In ft_defaults at 48 In ft_sourceanalysis at 197 the input is timelock data with 102 channels and 1500 timebins using headmodel specified in the configuration using gradiometers specified in the configuration computing surface normals creating dipole grid based on user specified dipole positions 1666 dipoles inside, 1904 dipoles outside brain scanning repetition 1 Warning: VOLTYPE is only a compatibility wrapper, which will soon be removed. Please instead call FT_VOLTYPE. > In voltype at 16 In beamformer_sam at 60 In ft_sourceanalysis at 868 ??? Error using ==> beamformer_sam at 66 unsupported voltype for determining the mean sphere origin Error in ==> ft_sourceanalysis at 868 dip(i) = beamformer_sam(grid, sens, vol, squeeze(avg(i,:,:)), squeeze(Cy(i,:,:)), optarg{:});
Robert Oostenveld - 2011-06-07 10:17:12 +0200
This is a limitation of Gareth's implementation, which requires the center of the volume conductor to be known. Probably you are using a singleshell model. Please first try with a multisphere volume conduction model. See FT_PREPARE_LOCALSPHERES If that works, I will consider extending the support to the singleshell.
Robert Oostenveld - 2011-06-07 10:17:29 +0200
I noticed a missing subfunction, please see attached.