<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
</head>
<body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">
Hi Nicole,
<div class=""><br class="">
</div>
<div class="">Lead field normalization is a different approach than Van Veen’s method, which is often called the Neural Activity Index (NAI) and closely related to the “unit noise gain” or “weight normalization” concept you might see in some literature.</div>
<div class=""><br class="">
</div>
<div class="">I have implemented the NAI in beamformer_lcmv.m, which you can run with:</div>
<div class="">cfg.method = ‘lcmv’;</div>
<div class="">cfg.lcmv.weightnorm = ‘nai’;</div>
<div class=""><br class="">
</div>
<div class="">However, the equivalent has not been implemented in the other beamformer variants yet (SAM, DICS). You can still get output equivalent to SAM using the LCMV method if you use cfg.keeptrials=‘yes’ and average the power of the resulting time series
(in source.avg.mom). This would give you a measure of induced power changes (rather than evoked), like the SAM procedure would. Unfortunately this procedure is not yet documented, but it’s not too tricky. (Please use a brand new version of FieldTrip if you’d
like to try this, as an old bug in the NAI orientation selection was inadvertently re-introduced in FieldTrip versions between September 2016 and last week).</div>
<div class=""><br class="">
</div>
<div class="">I personally find that the NAI gives more sensible results if you are contrasting something like post-stimulus activity to a pre-stimulus baseline. If you are instead contrasting two conditions against each other rather than a baseline, then the
different normalization approaches should give (almost) the same results anyway.</div>
<div class=""><br class="">
</div>
<div class="">Anyway, regarding lead field normalization: it does indeed do a voxel-by-voxel normalization since it cycles through all the voxels in a for loop ('for ii=1:Ndipoles' on the second line). It is purely based on the properties of the lead field,
and as you noticed, is unlike Van Veen’s method in that it does not use the noise estimate at all. BTW, I believe that the lead field "column normalization" approach has been more popular in the literature. This normalizes the x/y/z components of the lead
field independently, rather than all together. You can try this with cfg.normalize = ‘column’ and see how the results compare.</div>
<div class=""><br class="">
</div>
<div class="">Cheers,</div>
<div class="">Sarang</div>
<div class=""><br class="">
</div>
<div class=""><br class="">
</div>
<div class="">
<div>
<blockquote type="cite" class="">
<div class="">On 01 Mar 2017, at 11:43, Klink-3, N.E.C. van <<a href="mailto:N.vanKlink-2@umcutrecht.nl" class="">N.vanKlink-2@umcutrecht.nl</a>> wrote:</div>
<br class="Apple-interchange-newline">
<div class="">
<div style="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-stroke-width: 0px; direction: ltr; font-family: Tahoma; font-size: 10pt;" class="">
Dear all,
<div class=""><br class="">
</div>
<div class="">I want to do SAM beamformer source localization on single trial EEG data. I would like to normalize the leadfields to correct for depth, like mentioned in the lmcv beamformer tutorial: (<a href="http://www.fieldtriptoolbox.org/tutorial/beamformer_lcmv" class="">http://www.fieldtriptoolbox.org/tutorial/beamformer_lcmv</a>)</div>
<div class=""><font face="Courier New" class=""><br class="">
</font></div>
<div class="">
<div class=""><font face="Courier New" class="">cfg = [];</font></div>
<div class=""><font face="Courier New" class="">cfg.elec = hdr.elec; % electrode distances</font></div>
<div class=""><font face="Courier New" class="">cfg.headmodel = vol; % volume conduction headmodel</font></div>
<div class=""><font face="Courier New" class="">cfg.grid = grid; % normalized grid positions</font></div>
<div class=""><font face="Courier New" class="">cfg.channel = {'EEG'};</font></div>
<div class=""><font face="Courier New" class="">cfg.normalize = 'yes'; % to remove depth bias (Q in eq. 27 of van Veen et al, 1997). </font></div>
<div class=""><font face="Courier New" class="">lf = ft_prepare_leadfield(cfg);</font></div>
<div class=""><br class="">
</div>
<div class="">However when I look what happens with cfg.normalize='yes', the following script is used in ft_compute_leadfield, from line 570:</div>
<div class=""><br class="">
</div>
<div class="">
<div class=""><font face="Courier New" class=""> case 'yes'</font></div>
<div class=""><font face="Courier New" class=""> for ii=1:Ndipoles</font></div>
<div class=""><font face="Courier New" class=""> tmplf = lf(:, (3*ii-2):(3*ii));</font></div>
<div class=""><font face="Courier New" class=""> if normalizeparam==0.5</font></div>
<div class=""><font face="Courier New" class=""> % normalize the leadfield by the Frobenius norm of the matrix</font></div>
<div class=""><font face="Courier New" class=""> % this is the same as below in case normalizeparam is 0.5</font></div>
<div class=""><font face="Courier New" class=""> nrm = norm(tmplf, 'fro');</font></div>
<div class=""><font face="Courier New" class=""> else</font></div>
<div class=""><font face="Courier New" class=""> % normalize the leadfield by sum of squares of the elements of the leadfield matrix to the power "normalizeparam"</font></div>
<div class=""><font face="Courier New" class=""> % this is the same as the Frobenius norm if normalizeparam is 0.5</font></div>
<div class=""><font face="Courier New" class=""> nrm = sum(tmplf(:).^2)^normalizeparam;</font></div>
<div class=""><font face="Courier New" class=""> end</font></div>
<div class=""><font face="Courier New" class=""> if nrm>0</font></div>
<div class=""><font face="Courier New" class=""> tmplf = tmplf ./ nrm;</font></div>
<div class=""><font face="Courier New" class=""> end</font></div>
<div class=""><font face="Courier New" class=""> lf(:, (3*ii-2):(3*ii)) = tmplf;</font></div>
<div class=""><font face="Courier New" class=""> end</font></div>
</div>
</div>
<div class=""><br class="">
</div>
<div class="">This seems to me as independent of the dipole location, and does not use an estimate of the noise spectrum as in Eq 27 of van Veen et al 1997. </div>
<div class="">DICS beamformer has the option to estimate the noise spectrum with 'projectnoise', but SAM beamformer does not have that option. SAM does something with noise and a lambda, which is noise regularization I guess (beamformer_sam from line 102). <span style="font-size: 13.3333px;" class="">I
use Fieldtrip 20170212.</span></div>
<div class=""><br class="">
</div>
<div class="">My main question: how do I correct the leadfields for depth bias?</div>
<div class=""><br class="">
</div>
<div class="">Thanks in advance,</div>
<div class="">Nicole</div>
</div>
<hr style="font-family: Helvetica; font-size: 12px; 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-stroke-width: 0px;" class="">
<p style="font-family: Helvetica; font-size: 12px; 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-stroke-width: 0px;" class="">
<font size="-2" face="arial" class=""><i class="">De informatie opgenomen in dit bericht kan vertrouwelijk zijn en is uitsluitend bestemd voor de geadresseerde. Indien u dit bericht onterecht ontvangt, wordt u verzocht de inhoud niet te gebruiken en de afzender
direct te informeren door het bericht te retourneren. Het Universitair Medisch Centrum Utrecht is een publiekrechtelijke rechtspersoon in de zin van de W.H.W. (Wet Hoger Onderwijs en Wetenschappelijk Onderzoek) en staat geregistreerd bij de Kamer van Koophandel
voor Midden-Nederland onder nr. 30244197.</i></font></p>
<p style="font-family: Helvetica; font-size: 12px; 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-stroke-width: 0px;" class="">
<font size="-2" face="arial" color="green" class=""><i class="">Denk s.v.p aan het milieu voor u deze e-mail afdrukt.</i></font></p>
<hr style="font-family: Helvetica; font-size: 12px; 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-stroke-width: 0px;" class="">
<p style="font-family: Helvetica; font-size: 12px; 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-stroke-width: 0px;" class="">
<font size="-2" face="arial" class=""><i class="">This message may contain confidential information and is intended exclusively for the addressee. If you receive this message unintentionally, please do not use the contents but notify the sender immediately
by return e-mail. University Medical Center Utrecht is a legal person by public law and is registered at the Chamber of Commerce for Midden-Nederland under no. 30244197.</i></font></p>
<p style="font-family: Helvetica; font-size: 12px; 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-stroke-width: 0px;" class="">
<font size="-2" face="arial" color="green" class=""><i class="">Please consider the environment before printing this e-mail.</i></font></p>
<span style="font-family: Helvetica; font-size: 12px; 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-stroke-width: 0px; float: none; display: inline !important;" class="">_______________________________________________</span><br style="font-family: Helvetica; font-size: 12px; 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-stroke-width: 0px;" class="">
<span style="font-family: Helvetica; font-size: 12px; 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-stroke-width: 0px; float: none; display: inline !important;" class="">fieldtrip
mailing list</span><br style="font-family: Helvetica; font-size: 12px; 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-stroke-width: 0px;" class="">
<a href="mailto:fieldtrip@donders.ru.nl" style="font-family: Helvetica; font-size: 12px; 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-stroke-width: 0px;" class="">fieldtrip@donders.ru.nl</a><br style="font-family: Helvetica; font-size: 12px; 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-stroke-width: 0px;" class="">
<a href="https://mailman.science.ru.nl/mailman/listinfo/fieldtrip" style="font-family: Helvetica; font-size: 12px; 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-stroke-width: 0px;" class="">https://mailman.science.ru.nl/mailman/listinfo/fieldtrip</a></div>
</blockquote>
</div>
<br class="">
</div>
</body>
</html>