Back to the main page.

Bug 1589 - LCMV with keepmom = 'no' crashes, when using a precomputed filter

Status CLOSED WORKSFORME
Reported 2012-07-03 14:54:00 +0200
Modified 2012-08-23 14:02:05 +0200
Product: FieldTrip
Component: core
Version: unspecified
Hardware: PC
Operating System: Windows
Importance: P3 normal
Assigned to: Jörn M. Horschig
URL:
Tags:
Depends on:
Blocks:
See also:

Jörn M. Horschig - 2012-07-03 14:54:03 +0200

??? Reference to non-existent field 'mom'. Error in ==> ft_sourceanalysis at 808 tmpmom = tmpdip.mom{tmpdip.inside(1)}; I think there's not more to say, pretty obvious one


Jörn M. Horschig - 2012-07-04 13:40:57 +0200

cfg options I used lambda = '5%'; cfg = []; cfg.method = 'lcmv'; cfg.projectnoise = 'yes'; cfg.keeptrials = 'yes'; cfg.grid = grid_lf; cfg.grad = data.hdr.grad; cfg.vol = hdm; cfg.lambda = lambda; cfg.keepfilter = 'yes'; % cfg.fixedori = 'yes'; source = ft_sourceanalysis(cfg, data); %% cfg.filter = source.avg.filter; cfg.rawtrial = 'yes'; cfg.keepmom = 'yes'; source = ft_sourceanalysis(cfg, data);


Johanna - 2012-07-04 14:04:37 +0200

This worked for me without errors: load dataFIC cfg=[]; cfg.covariance='yes'; cfg.covariancewindow=[0.5 1.5]; cfg.keeptrials='yes'; tl=ft_timelockanalysis(cfg,dataFIC); load */tutorial/beamformer/segmentedmri.mat cfg = []; cfg.method = 'singleshell'; vol = ft_prepare_headmodel(cfg, segmentedmri); cfg = []; cfg.grad = tl.grad; cfg.vol = vol; cfg.reducerank = 2; cfg.channel = {'MEG','-MLP31', '-MLO12'}; cfg.grid.resolution = 1; % use a 3-D grid with a 1 cm resolution [grid] = ft_prepare_leadfield(cfg); cfg = []; cfg.method = 'lcmv'; cfg.projectnoise = 'yes'; cfg.keeptrials = 'yes'; cfg.grid = grid; cfg.grad = dataFIC.hdr.grad; cfg.vol = vol; cfg.lambda = '5%'; cfg.keepfilter = 'yes'; cfg.fixedori = 'yes'; source = ft_sourceanalysis(cfg, tl); %% cfg.filter = source.avg.filter; cfg.rawtrial = 'yes'; cfg.keepmom = 'yes'; source = ft_sourceanalysis(cfg, tl);


Jörn M. Horschig - 2012-07-04 14:29:16 +0200

strange, works for me as well now... maybe I did something else, which lead to the error. If I encounter this again, I'll reopen


Jörn M. Horschig - 2012-08-23 14:02:05 +0200

bug closing time (http://www.youtube.com/watch?v=xGytDsqkQY8)