<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=""><div style="text-align: start; text-indent: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class=""><div style="text-align: start; text-indent: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class=""><div style="text-align: start; text-indent: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class=""><div style="text-align: start; text-indent: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class=""><div style="text-align: start; text-indent: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class=""><div style="text-align: start; text-indent: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class=""><div style="text-align: start; text-indent: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class=""><div style="text-align: start; text-indent: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class=""><div style="text-align: start; text-indent: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class=""><div style="text-align: start; text-indent: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class=""><div style="text-align: start; text-indent: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class=""><div style="text-align: start; text-indent: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class=""><div style="text-align: start; text-indent: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class=""><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class=""><div style="text-align: start; text-indent: 0px;"><p class="MsoNormal" style="color: rgb(0, 0, 0); letter-spacing: normal; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; margin-bottom: 6pt;">Dear Fieldtripper, </p><div style="color: rgb(0, 0, 0); letter-spacing: normal; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class="">I am currently trying to compute a source reconstruction on my EEG data (Biosemi, 64 electrodes).</div><div style="color: rgb(0, 0, 0); letter-spacing: normal; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><br class=""></div><div style="color: rgb(0, 0, 0); letter-spacing: normal; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class="">I use the following code:</div><div style="color: rgb(0, 0, 0); letter-spacing: normal; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><br class=""></div><div class=""><div style="color: rgb(178, 69, 243); letter-spacing: normal; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; margin: 0px; font-size: 10px; line-height: normal; font-family: Courier;" class=""><span style="color: #000000" class="">load </span>Subject01_sourcemodel_15684<span style="color: #000000" class="">;</span></div><div style="color: rgb(0, 0, 0); letter-spacing: normal; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; margin: 0px; font-size: 10px; line-height: normal; font-family: Courier;" class="">headmodel = ft_read_headmodel(<span style="color: #b245f3" class="">'headmodel/standard_bem.mat'</span>);</div><div style="margin: 0px; line-height: normal;" class=""><font face="Courier" size="1" class="">electr_pos=ft_read_sens(</font><span style="color: rgb(178, 69, 243); font-family: Courier; font-size: 10px; letter-spacing: normal; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class="">'Biosemi-Cap64.sfp</span><font color="#b245f3" face="Courier" size="1" class="">’</font><font face="Courier" size="1" class="">);</font></div></div><div style="color: rgb(0, 0, 0); letter-spacing: normal; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; margin: 0px; font-size: 10px; line-height: normal; font-family: Courier;" class=""><br class=""></div><div style="color: rgb(0, 0, 0); letter-spacing: normal; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><div style="margin: 0px; font-size: 10px; line-height: normal; font-family: Courier; color: rgb(178, 69, 243);" class=""><span style="color: #000000" class="">load(</span>'GA_GroupPerpetrators_Parents_Neutral_DOUL.mat'<span style="color: #000000" class="">)</span></div></div><div style="color: rgb(0, 0, 0); letter-spacing: normal; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><span style="color: #000000" class=""><br class=""></span></div><div style="color: rgb(0, 0, 0); letter-spacing: normal; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><div style="margin: 0px; font-size: 10px; line-height: normal; font-family: Courier;" class="">cfg = [];</div><div style="margin: 0px; font-size: 10px; line-height: normal; font-family: Courier;" class="">cfg.method = <span style="color: #b245f3" class="">'interactive'</span>;</div><div style="margin: 0px; font-size: 10px; line-height: normal; font-family: Courier;" class="">cfg.elec = electr_pos;</div><div style="margin: 0px; font-size: 10px; line-height: normal; font-family: Courier;" class="">cfg.headshape = headmodel.bnd(1); <span style="color: #25992d" class="">%1 = skin</span></div><div style="margin: 0px; font-size: 10px; line-height: normal; font-family: Courier;" class="">elec = ft_electroderealign(cfg, electr_pos);</div><div style="margin: 0px; font-size: 10px; line-height: normal; font-family: Courier; color: rgb(37, 153, 45);" class=""><br class=""></div><div style="margin: 0px; font-size: 10px; line-height: normal; font-family: Courier;" class="">cfg = [];</div><div style="margin: 0px; font-size: 10px; line-height: normal; font-family: Courier; color: rgb(37, 153, 45);" class=""><span style="color: #000000" class="">cfg.elec = elec; </span>% sensor information</div><div style="margin: 0px; font-size: 10px; line-height: normal; font-family: Courier;" class="">cfg.channel = elec.label;<span style="color: #25992d" class="">%; % the used channels</span></div><div style="margin: 0px; font-size: 10px; line-height: normal; font-family: Courier;" class="">cfg.sourcemodel = sourcemodel; <span style="color: #25992d" class="">% source points</span></div><div style="margin: 0px; font-size: 10px; line-height: normal; font-family: Courier;" class="">cfg.headmodel = headmodel; <span style="color: #25992d" class="">% volume conduction model</span></div><div style="margin: 0px; font-size: 10px; line-height: normal; font-family: Courier;" class="">cfg.singleshell.batchsize = 5000; <span style="color: #25992d" class="">% speeds up the computation</span></div><div style="margin: 0px; font-size: 10px; line-height: normal; font-family: Courier;" class="">leadfield = ft_prepare_leadfield(cfg);</div><div style="margin: 0px; font-size: 10px; line-height: normal; font-family: Courier; color: rgb(37, 153, 45);" class="">%%</div><div style="margin: 0px; font-size: 10px; line-height: normal; font-family: Courier;" class="">cfg = [];</div><div style="margin: 0px; font-size: 10px; line-height: normal; font-family: Courier;" class="">cfg.method = <span style="color: #b245f3" class="">'mne'</span>;</div><div style="margin: 0px; font-size: 10px; line-height: normal; font-family: Courier;" class="">cfg.elec = elec;</div><div style="margin: 0px; font-size: 10px; line-height: normal; font-family: Courier;" class="">cfg.grid = leadfield;</div><div style="margin: 0px; font-size: 10px; line-height: normal; font-family: Courier;" class="">cfg.headmodel = headmodel;</div><div style="margin: 0px; font-size: 10px; line-height: normal; font-family: Courier;" class="">cfg.mne.prewhiten = <span style="color: #b245f3" class="">'yes'</span>;</div><div style="margin: 0px; font-size: 10px; line-height: normal; font-family: Courier;" class="">cfg.mne.lambda = 3;</div><div style="margin: 0px; font-size: 10px; line-height: normal; font-family: Courier;" class="">cfg.mne.scalesourcecov = <span style="color: #b245f3" class="">'yes'</span>;</div><p style="margin: 0px; font-size: 10px; line-height: normal; font-family: Courier; min-height: 12px;" class=""> </p><div style="margin: 0px; font-size: 10px; line-height: normal; font-family: Courier;" class="">source_GA_GroupPerpetrators_Parents_Neutral_DOUL = ft_sourceanalysis(cfg, GA_GroupPerpetrators_Parents_Neutral_DOUL);</div><div style="margin: 0px; font-size: 12px; line-height: normal; font-family: Courier; min-height: 14px;" class=""><br class=""></div></div><div style="letter-spacing: normal; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class="">However, I got the following error when computing the ft_sourceanalysis:</div><div style="letter-spacing: normal; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><span style="font-size: 11px;" class=""><font color="#ff2600" class=""><br class=""></font></span></div><div class=""><div class=""><span style="font-size: 11px;" class=""><font color="#ff2600" class="">Error using svd</font></span></div><div class=""><span style="font-size: 11px;" class=""><font color="#ff2600" class="">Input to SVD must not contain NaN or Inf.</font></span></div><div class=""><span style="font-size: 11px;" class=""><font color="#ff2600" class=""><br class=""></font></span></div><div class=""><span style="font-size: 11px;" class=""><font color="#ff2600" class="">Error in beamformer_lcmv>pinv (line 415)</font></span></div><div class=""><span style="font-size: 11px;" class=""><font color="#ff2600" class=""> [U,S,V] = svd(A,0);</font></span></div><div class=""><span style="font-size: 11px;" class=""><font color="#ff2600" class=""><br class=""></font></span></div><div class=""><span style="font-size: 11px;" class=""><font color="#ff2600" class="">Error in beamformer_lcmv (line 297)</font></span></div><div class=""><span style="font-size: 11px;" class=""><font color="#ff2600" class=""> filt = pinv(lf' * invCy * lf) * lf' * invCy; % van Veen eqn.</font></span></div><div class=""><span style="font-size: 11px;" class=""><font color="#ff2600" class=""> 23, use PINV/SVD to cover rank deficient leadfield</font></span></div><div class=""><span style="font-size: 11px;" class=""><font color="#ff2600" class=""><br class=""></font></span></div><div class=""><span style="font-size: 11px;" class=""><font color="#ff2600" class="">Error in ft_sourceanalysis (line 950)</font></span></div><div class=""><span style="font-size: 11px;" class=""><font color="#ff2600" class=""> dip(i) = beamformer_lcmv(grid, sens, headmodel, squeeze_avg, squeeze_Cy,</font></span></div><div class=""><span style="font-size: 11px;" class=""><font color="#ff2600" class=""> optarg{:});</font></span></div></div><div style="color: rgb(0, 0, 0); letter-spacing: normal; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><br class=""></div><div style="color: rgb(0, 0, 0); letter-spacing: normal; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class="">I checked in all the datafiles I use, and the only NaN I found are in the leadfield.leadfield but all the previous steps appear to be correct and contain the correct electrode locations, data, covariance matrix etc.</div><div style="color: rgb(0, 0, 0); letter-spacing: normal; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><br class=""></div><div style="color: rgb(0, 0, 0); letter-spacing: normal; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class="">I am thus not sure how to solve the problem and why my leadfield.leadfield has NaN (or even if the problem comes from there).</div><div style="color: rgb(0, 0, 0); letter-spacing: normal; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><br class=""></div><div style="color: rgb(0, 0, 0); letter-spacing: normal; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class="">Thank you in advance for your help!</div><div style="color: rgb(0, 0, 0); letter-spacing: normal; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class=""><br class=""></div><div style="color: rgb(0, 0, 0); letter-spacing: normal; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px;" class="">Emilie</div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></body></html>