Hi <span style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Courier; font-size: medium; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px;"><span style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Courier; font-size: medium; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px;"><span style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Courier; font-size: medium; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px;"><span style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Courier; font-size: medium; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px;"><span style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Courier; font-size: medium; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px;"><span style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Courier; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px;"><span style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Courier; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px;"><span style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Courier; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px;"><span style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Courier; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px;">Joseph</span></span></span></span></span></span></span></span></span>,<br>
<br>Most likely cause of ill-conditioned in neuromag data is maxfilter (signal space separation) which reduces the dimensionality from 306 (Channels in neuromag machine) to 64 (Spherical harmonics Multipolar expansion of order 8). Also neuromag is consist of two different type of sensors (Planar Gradiometer and magnetometer) having different sensitivities, so either ica needs to run separately on them or combined by making them equal through normalization.<br>
<br><font color="#888888">Sheraz Khan<br>Martinos Center<br>Boston, USA</font><br><br><br><br><div class="gmail_quote">On Wed, Jul 6, 2011 at 9:27 PM, Joseph Dien <span dir="ltr"><<a href="mailto:jdien07@mac.com">jdien07@mac.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;"><div style="word-wrap: break-word;"><div>I've had this same problem with EEGlab's runICA implementation when the rank of the data was less than the number of variables (i.e., it's an ill-conditioned matrix, which is to say it has multicolinearity). Including the reference channel can do this (since it's just zeroes). Using PCA to reduce the data to a subspace can help as well since the rank will now equal the number of variables (namely the number of factor scores will equal the number of retained factors). Before you go the PCA route, it's worth trying to figure out why the data is ill-conditioned. There is likely to be something wrong with one of your channels. Maybe two of them are shorted together or something along those lines. Or if you used mean mastoid reference and included both mastoids, that can do it too (since they will be exactly inversely correlated, one of them is statistically redundant, reducing the rank by one). My EP Toolkit artifact correction routine (<a href="http://sourceforge.net/projects/erppcatoolkit/" target="_blank">http://sourceforge.net/projects/erppcatoolkit/</a>) takes all these sorts of things into account when performing eyeblink correction.</div>
<div><br></div><div>Cheers!</div><div><br></div><div>Joe</div><div><div></div><div class="h5"><div><br></div><br><div><div>On Jul 6, 2011, at 8:44 PM, Sheraz Khan wrote:</div><br><blockquote type="cite">Hi Sylvana,<br><br>
Is the data SSS (maxfilter), I normally have this problem when data is SSS, try doing ICA with PCA (set at 64 components).<br><br>Sheraz Khan<br>Martinos Center<br>Boston, USA<br><br><div class="gmail_quote">
On Wed, Jul 6, 2011 at 2:09 PM, Stephen Whitmarsh <span dir="ltr"><<a href="mailto:stephen.whitmarsh@gmail.com" target="_blank">stephen.whitmarsh@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin: 0px 0px 0px 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
Dear Sylvana,<br>
<br>
>From the top op my head: When running the ICA, try only computing it<br>
for the MEG or EEG channels. (cfg.channel = 'MEG';), so it doesn't run<br>
on reference sensors in case of MEG.<br>
I have hunch that might help out.<br>
<br>
Cheers,<br>
<font color="#888888">Stephen<br>
</font><div><div></div><div><br>
<br>
On 6 July 2011 19:37, Sylvana <<a href="mailto:sylvana.schister@utah.edu" target="_blank">sylvana.schister@utah.edu</a>> wrote:<br>
> Hi all,<br>
><br>
> I am trying to use ICA to clean EOG and ECG artifacts of neuromag data.<br>
> I am running ft_topoplotIC(cfg, comp), where 'comp' is the output of<br>
> ft_componentanalysis.<br>
><br>
> I get the following error message:<br>
><br>
><br>
> ??? Error using ==> surf at 74<br>
> X, Y, Z, and C cannot be complex.<br>
><br>
> Error in ==> ft_plot_topo at 184<br>
> h = surf(Xi-deltax/2,Yi-deltay/2,zeros(size(Zi)), Zi, 'EdgeColor',<br>
> 'none', 'FaceColor', shading);<br>
><br>
> Error in ==> ft_topoplotER at 753<br>
> ft_plot_topo(chanX,chanY,datavector,'interpmethod',cfg.interpolation,...<br>
><br>
> Error in ==> ft_topoplotIC at 116<br>
> ft_topoplotER(cfg, varargin{:});<br>
><br>
><br>
> The complex variable is 'Zi', which is taken from 'comp'. I can modify<br>
> the function to take the real part of the variable (or the absolute<br>
> value), but I am not sure if this would be the correct thing to do. In<br>
> case this is a bug, I thought I should report it. Any comments?<br>
><br>
> Thanks for your help,<br>
><br>
> Sylvana Schister<br>
> Univ. of Utah<br>
> Dept. of Bioengineering<br>
><br>
> _______________________________________________<br>
> fieldtrip mailing list<br>
> <a href="mailto:fieldtrip@donders.ru.nl" target="_blank">fieldtrip@donders.ru.nl</a><br>
> <a href="http://mailman.science.ru.nl/mailman/listinfo/fieldtrip" target="_blank">http://mailman.science.ru.nl/mailman/listinfo/fieldtrip</a><br>
><br>
<br>
_______________________________________________<br>
fieldtrip mailing list<br>
<a href="mailto:fieldtrip@donders.ru.nl" target="_blank">fieldtrip@donders.ru.nl</a><br>
<a href="http://mailman.science.ru.nl/mailman/listinfo/fieldtrip" target="_blank">http://mailman.science.ru.nl/mailman/listinfo/fieldtrip</a><br>
</div></div></blockquote></div><br>
_______________________________________________<br>fieldtrip mailing list<br><a href="mailto:fieldtrip@donders.ru.nl" target="_blank">fieldtrip@donders.ru.nl</a><br><a href="http://mailman.science.ru.nl/mailman/listinfo/fieldtrip" target="_blank">http://mailman.science.ru.nl/mailman/listinfo/fieldtrip</a></blockquote>
</div><br></div></div><div>
<span style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Courier; font-size: medium; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px;"><span style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Courier; font-size: medium; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px;"><div style="word-wrap: break-word;">
<span style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Courier; font-size: medium; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px;"><div style="word-wrap: break-word;">
<span style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Courier; font-size: medium; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px;"><div style="word-wrap: break-word;">
<span style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Courier; font-size: medium; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px;"><div style="word-wrap: break-word;">
<span style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Courier; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px;"><div style="word-wrap: break-word;">
<span style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Courier; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px;"><div style="word-wrap: break-word;">
<span style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Courier; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px;"><div style="word-wrap: break-word;">
<span style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Courier; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px;"><div style="word-wrap: break-word;">
<div><div><div><div><div><br>--------------------------------------------------------------------------------</div><div><br></div><div>Joseph Dien,</div><div>Senior Research Scientist</div><div>University of Maryland </div>
<div><br></div><div>E-mail: <a href="mailto:jdien07@mac.com" target="_blank">jdien07@mac.com</a></div><div>Phone: <a href="tel:301-226-8848" value="+13012268848" target="_blank">301-226-8848</a></div><div>Fax: <a href="tel:301-226-8811" value="+13012268811" target="_blank">301-226-8811</a></div>
<div><a href="http://homepage.mac.com/jdien07/" target="_blank">http://homepage.mac.com/jdien07/</a></div></div></div><br></div><div><br></div><div><br></div><br></div></div></span></div></span></div></span></div></span></div>
</span></div></span></div></span><br></div></span><br></span><br>
</div>
<br></div><br>_______________________________________________<br>
fieldtrip mailing list<br>
<a href="mailto:fieldtrip@donders.ru.nl">fieldtrip@donders.ru.nl</a><br>
<a href="http://mailman.science.ru.nl/mailman/listinfo/fieldtrip" target="_blank">http://mailman.science.ru.nl/mailman/listinfo/fieldtrip</a><br></blockquote></div><br>