[FieldTrip] bemcp vs dipoli

Vladimir Litvak litvak.vladimir at gmail.com
Wed Nov 4 14:51:14 CET 2015


Dear all,

I have also looked at this in a different context and found that bemcp
computation breaks down at distances smaller than mean triangle side length
from the boundary (for meshes used in SPM that's about 6mm). So you can
avoid numerical problems by making sure your sources are never closer than
that to the inner skull boundary. This can be done with cfg.inwardshift
(for grids) or cfg.moveinward (for meshes) arguments of
ft_prepare_sourcemodel, In SPM12 this is done automatically for the forward
models we generate.

Best,

Vladimir

On Wed, Nov 4, 2015 at 1:36 PM, RICHARDS, JOHN <RICHARDS at mailbox.sc.edu>
wrote:

> In addition to the answer given by Alexandre Gramfort.
>
> I have compared BEMCP, Dipoli, Sim-bio FEM, and spherical models.  My
> tests use an individual head model with that individuals data, and then
> head models based on other participants, age-appropriate MRI templates, and
> age-inappropriate MRI templates; with infants.  The BEMCP is the
> three-compartment model, the Dipoli is a four-compartment model, the Simbio
> is a full segmented head model with 10 different materials, the spherical
> is four spheres.
>
> 1—I sometimes have issues getting the forward model / source analysis with
> BEMCP.  It appears to have singularities on the borders of the compartments
> and gives the same answer as you report.  Or does not even compute.  Using
> the same (or similar) head model mesh with Dipole, it works.
>
> 2—I have been using the Dipoli with four compartments. I am comparing the
> four methods with a set of empirical data and simulated data.  I have found
> that the dipole solutions on a single individual(s) for the BEMCP and the
> spherical models are closer in solution, and the Dipoli and the SimBio-FEM
> are closer in the solution.
>
> 3—I also have been using head models from different individuals to solve
> the forward model for the EEG from a single individual, comparing the
> solutions across head models.  Both the SimBio-FEM and Dipoli show the
> greatest change in solution across individuals with different head sizes or
> ages (using age-appropriate and age-inappropriate heads), whereas the BEMCP
> and spherical model show the least change.  Interesting, the larger change
> for the two models occurs because its theoretical fit when the head comes
> from the individual is better than the fit when the head comes from an
> inappropriate age, whereas the age-appropriate head model are not as
> important for BEMCP and spherical models because they fit poorer.
>
> The results from 2 and 3 are being analyzed now, eventually will reach
> publication.
>
> John
>
>
> >------------------------------
> >
> >Message: 4
> >Date: Wed, 4 Nov 2015 00:00:35 +0100
> >From: Maris Skujevskis <icelandhouse at gmail.com>
> >To: fieldtrip at science.ru.nl
> >Subject: [FieldTrip] bemcp vs dipoli
> >Message-ID:
> >       <CAKEvDicMF9uo1G-_TVjJvGXE9vvnbXmGVg7gLZQXG-=
> dQimv7Q at mail.gmail.com>
> >Content-Type: text/plain; charset="utf-8"
> >
> >Hi Fieldtrip users,
> >
> >When constructing a volume conduction model for EEG using the boundary
> >element method (BEM), there are two methods available: 'dipoli' and
> >'bemcp'. Besides technicalities (i.e., 'dipoli' only available on Linux),
> >are there any differences that you know of/have experienced that make one
> >method better (more reliable, more accurate, or anything else that makes
> >you prefer one rather) than the other?
> >
> >For some of my EEG subjects, 'dipoli' succeeds where 'bemcp' fails. During
> >processing with 'bemcp' there is a "warning: Matrix is singular, close to
> >singular or badly scaled. Results may be inaccurate. RCOND = NaN", with
> the
> >end result being that vol.mat contains NaNs.
> >The input in both cases, besides the method chosen, is identical.
> >
> >Taking one step back in the processing pipeline, I am aware that a poor
> >segmentation outcome might be a/the cause of the warnings and eventual
> >errors when constructing a volume conduction model. But in general this
> >should hold equally for both methods. What I am wondering about is why one
> >method might deal more successfully than the other with the same input.
> >
> >Best,
> >Maris
> >-------------- next part --------------
>
> _______________________________________________
> fieldtrip mailing list
> fieldtrip at donders.ru.nl
> http://mailman.science.ru.nl/mailman/listinfo/fieldtrip
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.science.ru.nl/pipermail/fieldtrip/attachments/20151104/f441e72a/attachment-0002.html>


More information about the fieldtrip mailing list