<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
</head>
<body style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">
Hi Casper,
<div class=""><br class="">
</div>
<div class="">These are all very good and relevant questions. As a matter of fact, a group of colleagues and friends (i.e. methods people, epilepsy researchers and core developers of fieldtrip/mne_python/spm/brainstorm)) have been working on this topic on and
off over the past years, meeting at Aston University, to try and understand the best way to do beamforming analysis on Elekta MEG data. We think that we now understand what’s the best way to approach this, and work is underway to disseminate this to the community.</div>
<div class=""><br class="">
</div>
<div class="">In principle, gradiometers and magnetometers can be combined in a single analysis, but this is unlikely to be optimal. Reason: the different order of magnitude of the signals (which may be aggravated when changing the metric units of the data
(e.g. going from T/m to T/cm for the gradiometers)) causes the covariance matrix (used for the beamformer spatial filter computation) to be differently weighing in the mags versus the grads. Indeed, if it’s deep sources you’re after, you may want to exploit
the depth sensitivity pattern of the mags.</div>
<div class=""><br class="">
</div>
<div class="">Some software packages address this by differently weighing the sensor types (by means of a baked-in default weighing scheme). A more principled way of addressing this would be to spatially whiten the data (and apply the same whitening to the
leadfields). In FieldTrip this prewhitening can be done with ft_denoise_prewhiten, using the covariance estimated in a baseline window as the prewhitener. If the balanced gradiometer array is subsequently used for the forward computation, the same prewhitener
is applied to the leadfield, so everything stays consistent.</div>
<div class=""><br class="">
</div>
<div class="">One thing that is important to take into account, is the rank-deficiency of maxfiltered data. This can cause problems in the mathematical inversion of the covariance matrix, when done inappropriately. In order to have more flexibility in optimally
inverting covariance matrices, we now have a dedicated ft_inv function (in fieldtrip/private), which allows for different ways of (regularized) matrix inversions. For maxfiltered filtered data, it would be recommended to use a ‘kappa’ parameter for inversion,
i.e. to truncate the svd of the coviance matrix at (or just below) its rank.</div>
<div class=""><br class="">
</div>
<div class="">Currently I am working on some documentation (that uses the Wakeman and Henson dataset) that demonstrates some of these issues. It will be presented at this years practicalMEEG workshop in Paris (<a href="http://practicalmeeg2019.org" class="">http://practicalmeeg2019.org</a>),
so stay tuned! This documentation will become available soon.</div>
<div class=""><br class="">
</div>
<div class="">Best wishes,</div>
<div class="">Jan-Mathijs</div>
<div class=""><br class="">
</div>
<div class=""><br class="">
</div>
<div class=""><br class="">
<div><br class="">
<blockquote type="cite" class="">
<div class="">On 15 Nov 2019, at 10:18, Casper Kerren <<a href="mailto:C.Kerren@pgr.bham.ac.uk" class="">C.Kerren@pgr.bham.ac.uk</a>> wrote:</div>
<br class="Apple-interchange-newline">
<div class="">
<div class="WordSection1" style="page: WordSection1; caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 14px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;">
<div style="margin: 0cm 0cm 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif;" class="">
Hi,<o:p class=""></o:p></div>
<div style="margin: 0cm 0cm 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif;" class="">
<o:p class=""> </o:p></div>
<div style="margin: 0cm 0cm 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif;" class="">
I am currently source localising an effect using LCMV on Elekta MEG data. I have both gradiometers and magnetometers and have until recently treated these as the same when preparing the leadfield. However, after having had discussions with my colleagues about
the different units and strength for grads and mags, we are hesitant as to what might be the best way of selecting channels. I have searched the fiedltrip mailing list, and have found some answers, but nothing that I can say is very convincing. Could someone
point me in the right direction regarding this? Perhaps someone who has done some checks on using both grads and mags and using only grads or mags. We are interested in a deeper source, and naively thought that magnetometers would be necessary to include.
However, is this just as I say: naïve?<o:p class=""></o:p></div>
<div style="margin: 0cm 0cm 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif;" class="">
<o:p class=""> </o:p></div>
<div style="margin: 0cm 0cm 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif;" class="">
In the ft_prepareleadfield function I can see that gradiometers is the only channels that are thought of being included, so does fieldtrip actually expect us to exclude magnetometers?<o:p class=""></o:p></div>
<div style="margin: 0cm 0cm 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif;" class="">
<o:p class=""> </o:p></div>
<div style="margin: 0cm 0cm 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif;" class="">
Kind regards,<o:p class=""></o:p></div>
<div style="margin: 0cm 0cm 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif;" class="">
<o:p class=""> </o:p></div>
<div style="margin: 0cm 0cm 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif; background-color: white;" class="">
<span style="font-size: 10pt; font-family: Helvetica, sans-serif;" class="">Casper Kerrén, PhD student<o:p class=""></o:p></span></div>
<div style="margin: 0cm 0cm 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif; background-color: white;" class="">
<span style="font-size: 10pt; font-family: Arial, sans-serif;" class="">School of Psychology <br class="">
College of Life and Environmental Sciences <br class="">
University of Birmingham <br class="">
Edgbaston <br class="">
Birmingham B15 2TT<br class="">
e-mail: <a href="mailto:cxk699@student.bham.ac.uk" style="color: purple; text-decoration: underline;" class="">cxk699@student.bham.ac.uk</a> <br class="">
<br class="">
</span><span style="" class=""><o:p class=""></o:p></span></div>
<div style="margin: 0cm 0cm 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif;" class="">
<o:p class=""> </o:p></div>
<div style="margin: 0cm 0cm 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif;" class="">
<o:p class=""> </o:p></div>
</div>
<span style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 14px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; float: none; display: inline !important;" class="">_______________________________________________</span><br style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 14px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;" class="">
<span style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 14px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; float: none; display: inline !important;" class="">fieldtrip
mailing list</span><br style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 14px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;" class="">
<a href="https://mailman.science.ru.nl/mailman/listinfo/fieldtrip" style="color: purple; text-decoration: underline; font-family: Helvetica; font-size: 14px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px;" class="">https://mailman.science.ru.nl/mailman/listinfo/fieldtrip</a><br style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 14px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;" class="">
<a href="https://doi.org/10.1371/journal.pcbi.1002202" style="color: purple; text-decoration: underline; font-family: Helvetica; font-size: 14px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px;" class="">https://doi.org/10.1371/journal.pcbi.1002202</a></div>
</blockquote>
</div>
<br class="">
</div>
</body>
</html>