<html><head></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">Max,<div><br></div><div>There is no reason why you can't simply use the same strategy for both dipoles. Each will produce its own leadfields, and then your inner product of the pinv'd leadfields and the waveforms will give you two separate waveforms.</div><div><br></div><div>Note that by using this pseudoinverse of the leadfield approach, depending on how far apart the sources are, you may have substantial correlation between the two waveforms that is an artifact of their leadfield correlation. This shouldn't be a big problem for things like auditory sources for the left and right hemisphere, or anterior and posterior sources within a hemisphere, but once you start getting sources closer together, you might consider an alternative approach. You can try using a beamformer to suppress correlated source activity, but you could also just adopt a beamformer-like set of weights for your dipole model instead. For that, you would ideally compute a covariance matrix on the epoched data, prior to averaging. Then, you can extend the logic of the pseudoinverse method by multiplication with inverse of the covariance matrix. You can see this in the fieldtrip functions that do beamformers, such as beamformer_lcmv. Around line 255, you can see:</div><div><br></div><div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font: normal normal normal 10px/normal Courier; ">filt = pinv(lf' * invCy * lf) * lf' * invCy;</div><div><br></div><div>where lf = leadfield and invCy = inverse of covariance matrix and filt = the resulting weight vector for your sensors for that particular source location. You can then do:</div><div><br></div><div>waveform = filt * data to get your timecourse. In the previous example, I used the dot function to remind myself it was an inner product, but it isn't strictly necessary in matlab.</div><div><br></div><div>Don</div><div><br></div><div><div>On May 23, 2012, at 7:46 AM, Taubert, Max wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite">
 
    <meta http-equiv="Content-Type" content="text/html; charset=UTF-8">
 
 <div><div style="margin-top: 0pt; margin-right: 0pt; margin-bottom: 0pt; margin-left: 0pt; ">
   <span>
    <span></span>
    Hi,
   </span>
  </div><p style="margin: 0pt;">
   <span> </span>
  </p><div style="margin-top: 0pt; margin-right: 0pt; margin-bottom: 0pt; margin-left: 0pt; ">
   <span>I'd like to use the leadfield from the ft_compute_leadfield function to compute the time course of two dipoles, given their location and orientation. So I searched on the mailing list and found a thread from  Apr-1 2011 (<a href="http://mailman.science.ru.nl/pipermail/fieldtrip/2011-April/003690.html">http://mailman.science.ru.nl/pipermail/fieldtrip/2011-April/003690.html</a>) with a nice explanation. By using this approach it is not very difficult to obtain the waveform of a single dipole, but in my case there are TWO dipoles.</span>
  </div><p style="margin: 0pt;">
   <span> </span>
  </p><div style="margin-top: 0pt; margin-right: 0pt; margin-bottom: 0pt; margin-left: 0pt; ">
   <span>This works when assuming a single dipole:</span>
  </div><p style="margin: 0pt;">
   <span> </span>
  </p><div style="margin-top: 0pt; margin-right: 0pt; margin-bottom: 0pt; margin-left: 0pt; ">
   <span>dip.pos = [0 0.5 0];</span>
  </div><div style="margin-top: 0pt; margin-right: 0pt; margin-bottom: 0pt; margin-left: 0pt; ">
   <span>dip.ori = [1 0 0]; </span>
  </div><div style="margin-top: 0pt; margin-right: 0pt; margin-bottom: 0pt; margin-left: 0pt; ">
   <span>
    lf = ft_compute_leadfield(dip.pos,elec,vol);
    <br>
    li = pinv(lf);
    <br>
    ori = dip.mom;
    <br>
    ori = repmat(ori,length(li),1)';
    <br>
    w = dot(li,ori);
    <br>
    waveform = dot(D.Data,repmat(w,size(D.Data,2),1)');
   </span>
  </div><p style="margin: 0pt;">
   <span> </span>
  </p><div style="margin-top: 0pt; margin-right: 0pt; margin-bottom: 0pt; margin-left: 0pt; ">
   <span>Is there any way to make this work with a two-dipole-model, so I get the waveforms of both of them? </span>
  </div><p style="margin: 0pt;">
   <span> </span>
  </p><div style="margin-top: 0pt; margin-right: 0pt; margin-bottom: 0pt; margin-left: 0pt; ">
   <span>Thanks in advance</span>
  </div><div style="margin-top: 0pt; margin-right: 0pt; margin-bottom: 0pt; margin-left: 0pt; ">
   <span>Max </span>
  </div><p style="margin: 0px; "></p>
 </div>
_______________________________________________<br>fieldtrip mailing list<br><a href="mailto:fieldtrip@science.ru.nl">fieldtrip@science.ru.nl</a><br>http://mailman.science.ru.nl/mailman/listinfo/fieldtrip</blockquote></div><br><div>
<div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">-----------------------<br>Don Rojas, Ph.D.<br>Associate Professor of Psychiatry<br>U. of Colorado Denver Anschutz Medical Campus<br>Director, UCD Magnetoencephalography Lab<br>13001 E. 17th Pl F546<br>Aurora, CO 80045<br>303-724-4994</div>
</div>
<br></div></body></html>