Hi guys, <br><br>   Thank you for your discussion. However, I did not see a direct answer to my question. My problem is, when the number of component is smaller than the number of channels, something wrong would happen. It seems that this problem appears after computing ICA. I would appreciate to to know how to solve this problem.<br>
   The ICA would converge in my situation. For Joe's question on how long, it's about one night for one dataset (122 channels, 240 trials, 3000 data points for each), if I take the same number of IC as the channels. I have 288 datasets, so it is too long to be acceptable. Probably, I can downsample the dataset to reduce the time for computing ICA, but then I still need to pick the noise components (eg. blinks) out of ~120 ones, this manual work also takes a lot of time. Thus I wish I can have less ICs.<br>
   Thanks again.<br>btw:   I am using matlab R2009a on Vista SP1, both 64-bit.<br><br><div class="gmail_quote">On Tue, Feb 2, 2010 at 12:56 AM, Arnaud Delorme <span dir="ltr"><<a href="mailto:arno@ucsd.edu">arno@ucsd.edu</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;"><div class="im">> Arno -- There was a problem with the Matlab rank() function on 64-<br>
> bit machines, I believe. Has this been solved and Is the auto rank<br>
> detection -> PCA option currently implemented in runica/binica?<br>
> Perhaps we could add a 'toy' rank() function pre-test (e.g. finding<br>
> the rank of a small full-rank matrix to detect if rank() is<br>
> working...) ? If so, run the rank test; if not, then warn the user<br>
> or build a work-around rank function that will work properly?<br>
<br>
</div>This is correct. This was a problem with Matlab 64-bit in 2007. One of<br>
the EEGLAB user, Sven Hoffman, implemented an alternate computation of<br>
the rank which we used whenever the Matlab rank is deficient. However,<br>
the new Matlab rank is fine now.<br>
<br>
Arno<br>
<div><div></div><div class="h5"><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>
</div></div></blockquote></div><br>