MAX - Unfortunately, in general using two dipoles rather than one will ~always improve the fit. Even if the source is a pure single dipole, a second dipole can be used to correct for noise or errors in the forward head model. This is less always the case for the constrained spatially-symmetric dipole pairs allowed by dipfit().  However, we have not thought of an optimal way to decide between using one or (occasionally) two dipoles to fit e.g. maps of ICA brain sources.  The goal would be to decide whether the two-dipole version is fitting noise/forward model error vs actual bilateral source generation...<div>
<br></div><div>Scott Makeig<br><br><div class="gmail_quote">On Thu, Aug 9, 2012 at 1:54 AM, Maximilien Chaumon <span dir="ltr"><<a href="mailto:maximilien.chaumon@gmail.com" target="_blank">maximilien.chaumon@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hello all,<br><br>When fitting dipoles to components, we are all sooner or later puzzled by the question whether to use one or two symmetrical dipoles.<br>
<br>Would it be correct to put the problems in terms of a nested hypothesis testing?<br>

<br>We are fitting a scalp map with one or two parameters and get a residual variance after the fit.<br>Could we not use this residual variance as a measure of the SSE and compute a F statistic to decide whether to use the more complex (with two dipoles) or simpler (with one dipole) of two nested models?<br>


If yes, then how would we decide on the number of degrees of freedom? How many free parameters do we have in each case? x,y,z,and two orientations per dipole? how does the imposed symmetry affect that number? Could we really map residual variance to SSE? How many "observations" do we have in that case (see formula below)?<br>


<br>I found this formula, for F:<br>F = (SSEF-SSER)/ (kF-kR) / ((1-SSEF)/(N-kF-1))<br>where<br>SSE is sum of squared errors,<br>k is numbers of parameters,<br>N is number of observations (? what in our case?)<br>F and R indices for full and reduced model respectively (in our case two and one dipole).<br>


<br><br>Thanks a lot for any comment!<br>Best,<br>Max<br><br>dipfit<br><br><br><br><br>
<br>_______________________________________________<br>
Eeglablist page: <a href="http://sccn.ucsd.edu/eeglab/eeglabmail.html" target="_blank">http://sccn.ucsd.edu/eeglab/eeglabmail.html</a><br>
To unsubscribe, send an empty email to <a href="mailto:eeglablist-unsubscribe@sccn.ucsd.edu">eeglablist-unsubscribe@sccn.ucsd.edu</a><br>
For digest mode, send an email with the subject "set digest mime" to <a href="mailto:eeglablist-request@sccn.ucsd.edu">eeglablist-request@sccn.ucsd.edu</a><br></blockquote></div><br><br clear="all"><div><br></div>
-- <br>Scott Makeig, Research Scientist and Director, Swartz Center for Computational Neuroscience, Institute for Neural Computation; Prof. of Neurosciences (Adj.), University of California San Diego, La Jolla CA 92093-0559, <a href="http://sccn.ucsd.edu/%7Escott" target="_blank">http://sccn.ucsd.edu/~scott</a><br>

</div>