Dear Fieldtrippers,<br><br><br>    I am having trouble using ft_dipolefitting. Everything looks fine in terms of the data and the headmodel but when i try to fit a dipole, i constantly get errors. Also the dipole fitting doesn't seem to work properly, because the outcome says 0 dipoles inside and, xx dipoles outside of the brain. then the error below pops up. Attached is also the headmodel and the ERF.<br><br><br><br>Error in ==> ft_prepare_sourcemodel at 550<br>  xmin_indx = find(grid.xgrid==xmin);<br><br><br>Error in ==> ft_dipolefitting at 348<br>  [grid, tmpcfg] = ft_prepare_sourcemodel(tmpcfg);<br><br><br>***************************************************<br>Here is the code:<br><br>% read header<br> hdr = ft_read_header('C:\fieldtrip_data\6100\e,rfhp1.0Hz');<br><br><br>% read data<br>cfg=[];<br>cfg.dataset='C:\fieldtrip_data\6100\e,rfhp1.0Hz';<br>data=ft_preprocessing(cfg);<br><br><br>% denoising<br><br>cfg.refchannel='MEGREF';<br>cfg.channel='MEG';<br>data=ft_denoise_pca(cfg,data);<br><br><br>% computing ERF<br><br><br><br>cfg=[];<br>cfg.trials='all';<br>cfg.removemean='yes';<br>cfg.vartrllength=0;<br>avg=ft_timelockanalysis(cfg,data);<br><br><br>% prepare headmodel<br><br>[shape] = ft_read_headshape('C:\fieldtrip_data\6100\hs_file')<br>cfg=[];<br>cfg.grad=hdr.grad;<br>cfg.headshape=shape.pnt;<br>cfg.singlesphere='yes';<br>vol = ft_prepare_localspheres(cfg);<br><br><br><br>% fit a dipole to the M50 and M100 components<br>cfg = [];<br>cfg.latency = [0.045 0.090]; <br>cfg.numdipoles = 1;<br>cfg.vol=vol;<br>cfg.feedback = 'textbar';<br>cfg.resolution = 2;<br>dip= ft_dipolefitting(cfg, avg);<br>******************************************************************<br><br><br>any help will be greatly appreciated. Thanks in advance<br>Mehmet