<div dir="ltr">Hi Mathias,<div>I will suggest you try other method like eLORETA or MNE  first and compare the result.</div><div><br></div><div>Best,</div><div><br></div><div>AZ</div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Jun 23, 2016 at 9:34 AM, Dr.-Ing. Matthias Dümpelmann <span dir="ltr"><<a href="mailto:matthias.duempelmann@uniklinik-freiburg.de" target="_blank">matthias.duempelmann@uniklinik-freiburg.de</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Dear colleagues,<br>
<br>
we are trying to evaluate systematically the localization  of rhythmic activity at the onset of epileptic seizures in patients, who are subsequently evaluated with intracranial recordings.<br>
<br>
We are trying to use DICS but quite often receive results, which are obviously wrong. In the attached screenshots you can see the peace of EEG, which is evaluated and the localization result with a maximum in the cerebellum.<br>
<br>
This a typical localization for the cases with obvious miss localizations. A further seizure with a similar EEG pattern was localized in the right temporal lobe, which matches the EEG pattern and also the intracranial findings.<br>
<br>
Head model is generated with ASA, also the electrodes file is in the ASA format. The geometrical match of source points, head model and electrodes is checked using ft_plot_.... methods.<br>
<br>
Here is the code sniplet relevant for the localization algorithm:<br>
<br>
% Analysis frequency<br>
SourceFrequency = 4.0<br>
<br>
% read eeg, electrodes and volume conductor model<br>
%[dat] = ft_read_data(filename);<br>
hdr = ft_read_header(filename);<br>
elec =  ft_read_sens(elecfile,  'fileformat','asa_elc');<br>
headmodel = ft_read_vol(volumefilename,'fileformat','asa_vol');<br>
mri = ft_read_mri(mrifilename, 'fileformat','asa_mri');<br>
<br>
% eeg data preprocessing<br>
cfgeeg.dataset = filename;<br>
cfgeeg.continuous   = 'yes'<br>
cfg.covariance       = 'yes';<br>
cfg.covariancewindow = 'all';<br>
cfg.demean = 'yes'<br>
<br>
[data] = ft_preprocessing(cfgeeg);<br>
<br>
% Freqanalysis for beamformer<br>
cfg = [];<br>
cfg.channel      = {'*'};<br>
cfg.method       = 'mtmfft';<br>
cfg.foilim     = [1 100];<br>
cfg.tapsmofrq  = 4;<br>
cfg.taper    = 'hanning';<br>
cfg.output       = 'powandcsd';<br>
<br>
<br>
% for common filter over conditions<br>
powcsd_all      = ft_freqanalysis(cfg, data);<br>
<br>
% Create leadfield grid<br>
cfg                 = [];<br>
cfg.channel         = elec.label;<br>
cfg.elec            = elec;<br>
cfg.headmodel       = headmodel;<br>
cfg.dics.reducerank = 3; % default for MEG is 2, for EEG is 3<br>
cfg.grid.resolution = 0.5;   % use a 3-D grid with a 0.5 cm resolution<br>
cfg.grid.unit       = 'cm';<br>
cfg.grid.tight      = 'yes';<br>
[grid] = ft_prepare_leadfield(cfg);<br>
<br>
% source analysis<br>
cfg              = [];<br>
cfg.channel      = elec.label;<br>
cfg.method       = 'dics';<br>
cfg.frequency    = SourceFrequency;<br>
cfg.grid         = grid;<br>
cfg.headmodel    = headmodel;<br>
cfg.elec         = elec;<br>
cfg.senstype     = 'EEG';<br>
cfg.dics.keepfilter   = 'yes'; % We wish to use the calculated filter later on<br>
cfg.dics.projectnoise = 'yes';<br>
cfg.dics.lambda  = '0.3%';<br>
source_raw = ft_sourceanalysis(cfg, powcsd_all);<br>
<br>
<br>
I played with a couple of parameters like cfg.dics.lambda, but without success.<br>
<br>
Do you have ideas how to improve our localization results ?<br>
<br>
Thank you,<br>
<br>
Matthias<br>
<br>
<br>
<br>
UNIVERSITÄTSKLINIKUM FREIBURG<br>
Epileptologie<br>
Breisacher Str. 64 · 79106 Freiburg<br>
Telefon: +49 761 270-52410<br>
<a href="mailto:matthias.duempelmann@uniklinik-freiburg.de">matthias.duempelmann@uniklinik-freiburg.de</a><br>
<br>
<a href="http://www.uniklinik-freiburg.de" rel="noreferrer" target="_blank">www.uniklinik-freiburg.de</a><br>
<br>
<br>_______________________________________________<br>
fieldtrip mailing list<br>
<a href="mailto:fieldtrip@donders.ru.nl">fieldtrip@donders.ru.nl</a><br>
<a href="https://mailman.science.ru.nl/mailman/listinfo/fieldtrip" rel="noreferrer" target="_blank">https://mailman.science.ru.nl/mailman/listinfo/fieldtrip</a><br></blockquote></div><br></div></div>