<p dir="ltr">Hey,</p>
<p dir="ltr">You are correct, this is a bug in dipfit that only throws an error when the number of components is fewer than the number of data channels. It wasn't that hard to fix in the code (like 3 or 4 function calls deep from pop_multifit) but I haven't submitted the bug report yet.</p>
<p dir="ltr">Michael</p>
<br><div class="gmail_quote"><div dir="ltr">On Thu, Sep 3, 2015, 13:46 <a href="mailto:lxraplikelx@qq.com">lxraplikelx@qq.com</a> <<a href="mailto:lxraplikelx@qq.com">lxraplikelx@qq.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div>
<div><span></span>Dear EEGLABLIST</div><div><br></div><div>Recently I am in trouble in using DIPFIT. I found a <a href="http://sccn.ucsd.edu/pipermail/eeglablist/2015/009530.html" target="_blank">same problem</a> in the archive of eeglablist but I couldn't find any solution. </div><div><br></div><div>The problem is "<b>atatype_raw error</b>", and the error information is "<b>inconsistent number of channels in trial 1</b>". </div><div><br></div><div>My dataset was average-referenced, and PCA to (EEG.nbchan-1) when I run ICA decomposition. I thought that maybe the discrepancy between the number of channels and components caused that (I actually know it must reduce the rank of the data when average-reference is used). Therefore I experimented not to reduce the rank of data, but directly run ICA. With the data of same numbers of channels and ICs, there is not aforementioned error any more. However, It resulted in bad dipole fitting (e.g. all dipoles are outside or converging to a same location). I do not know if my analysis pipeline is correct. My pipeline is below:</div><div><br></div><div><span>    1. import data (brain vision)</span></div><div><span><span>    2. 1-Hz high-pass filtering (considering good ICA results)</span></span></div><div><span><span><span>    3. clean line noise (CleanLine)</span></span></span></div><div><span><span><span><span>    4. downsample to 256-Hz</span></span></span></span></div>
<div><span>    5. load channel location and append online reference channel (FCz)</span></div><div><span><span style="font-family:'';font-size:10.5pt;line-height:1.5;background-color:window">    5. average reference and </span>add online reference channel back to data at the same time</span></div><div><span><span>    6. remove channels (VEO, HEOR)</span></span></div><div><span><span><span>    7. run clean_rawdata (only to identify bad channels)</span></span></span></div><div><span><span><span><span>    8. average re-reference</span></span></span></span></div><div><span><span><span><span><span>    9. epoch and baseline-correction with the mean of whole epoch data</span></span></span></span></span></div><div><span><span><span><span><span><span>    10. reject bad epoch (semi-automatic)</span></span></span></span></span></span></div><div><span><span><span><span><span><span><span>    11. run ICA (with 'pca', EEG.nbchan-1)</span></span></span></span></span></span></span></div><div><span><span><span><span><span><span><span><span>    12. DIPFIT ---> error</span></span></span></span></span></span></span></span></div><div><br></div><div><span>    PS: eeglab version is <b>13.2.2b</b>; matlab version is </span><span style="font-family:'';font-size:10.5pt;line-height:1.5;background-color:window"><b>7.13.0.564 (R2011b)</b></span></div><div><span>    PPS: I will try <b>Prep Pipeline</b> to pre-processe data in the future, so with its <b>robust average</b>, it could also run into trouble for the unbalancing number of channels and ICs.</span></div><div><span style="font-family:'';font-size:10.5pt;line-height:1.5;background-color:window">    </span></div><div><span style="font-size:10.5pt;line-height:1.5;background-color:window">If any step of my pipeline is inappropriate and causes error I came across, please tell me. I appreciate your helps!</span></div><div><br></div><div>Best regards,</div><div><br></div><div>Li Xiang</div><div>Southwest University, China</div>
</div>_______________________________________________<br>
Eeglablist page: <a href="http://sccn.ucsd.edu/eeglab/eeglabmail.html" rel="noreferrer" 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" target="_blank">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" target="_blank">eeglablist-request@sccn.ucsd.edu</a></blockquote></div>