<div dir="ltr">Thanks Dan and Alex.<div><br></div><div>How are multiple rows of the projector used? In SSP (PCA) the basis vectors are orthogonal however ICA components may not be.</div><div><br></div><div>Would you suggest putting all of the excluded ica components into a single projector or to build separate projectors for each excluded ICA component?</div><div><br></div><div>--Luke</div></div><br><div class="gmail_quote"><div dir="ltr">On Wed, Aug 16, 2017 at 1:40 PM Dan McCloy <<a href="mailto:drmccloy@uw.edu">drmccloy@uw.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Luke,<br>I've created projectors "by hand" like this (caveat: my intent was just to plot topomaps, not to actually apply the projector):<br><div><br>desc = '4Hz_artifact_ICA'<br>data = dict(data=my_component, nrow=1, ncol=len(my_component),<br> row_names=desc, col_names=epochs.ch_names)<br>projector = mne.io.Projection(data=data, active=False, kind=1,<br> desc=desc, explained_var=None)<br>epochs.add_proj(projector)<br><br></div><div>I've never used it with ICA (only with DSS) but I think in principle it could work the same way? <br></div><div><div dir="ltr">-- dan<br><br>Daniel McCloy<br><a href="http://dan.mccloy.info/" target="_blank">http://dan.mccloy.info/</a><br>Postdoctoral Research Associate<br>Institute for Learning and Brain Sciences<br>University of Washington<br></div>
<br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Aug 16, 2017 at 9:53 AM, Alexandre Gramfort <span dir="ltr"><<a href="mailto:alexandre.gramfort@inria.fr" target="_blank">alexandre.gramfort@inria.fr</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi Luke,<br>
<br>
we don't have such a function AFAIK and I see a use to it. My first<br>
reaction was to tall you to use a high n_pca_components in ICA but<br>
it's not fully equivalent to an SSP.<br>
<br>
proj = make_ica_proj(ica, picks)<br>
<br>
?<br>
<br>
Alex<br>
<div><div class="m_-2600843247925618138h5"><br>
On Wed, Aug 16, 2017 at 3:50 PM, Luke Bloy <<a href="mailto:luke.bloy@gmail.com" target="_blank">luke.bloy@gmail.com</a>> wrote:<br>
> Before I potentially reinvent a wheel, I was wondering if there was a way to<br>
> compute an ssp projection operator (proj) from the excluded basis set of an<br>
> ICA operator.<br>
><br>
> The reason I'm trying to do this is that I have some data that has an<br>
> environmental artifact that is always around 4HZ but whose topology changes<br>
> slightly from subject to subject. ICA on a filtered dataset (3.5-4.5Hz)<br>
> makes it easy to ID the artifact components but the remaining ICA components<br>
> are not really applicable to the broadband signal so using ica.apply() on<br>
> the boradband signal doesn't really work for me.<br>
><br>
> Thanks<br>
> Luke<br>
><br>
</div></div>> _______________________________________________<br>
> Mne_analysis mailing list<br>
> <a href="mailto:Mne_analysis@nmr.mgh.harvard.edu" target="_blank">Mne_analysis@nmr.mgh.harvard.edu</a><br>
> <a href="https://mail.nmr.mgh.harvard.edu/mailman/listinfo/mne_analysis" rel="noreferrer" target="_blank">https://mail.nmr.mgh.harvard.edu/mailman/listinfo/mne_analysis</a><br>
><br>
><br>
> The information in this e-mail is intended only for the person to whom it is<br>
> addressed. If you believe this e-mail was sent to you in error and the<br>
> e-mail<br>
> contains patient information, please contact the Partners Compliance<br>
> HelpLine at<br>
> <a href="http://www.partners.org/complianceline" rel="noreferrer" target="_blank">http://www.partners.org/complianceline</a> . If the e-mail was sent to you in<br>
> error<br>
> but does not contain patient information, please contact the sender and<br>
> properly<br>
> dispose of the e-mail.<br>
><br>
_______________________________________________<br>
Mne_analysis mailing list<br>
<a href="mailto:Mne_analysis@nmr.mgh.harvard.edu" target="_blank">Mne_analysis@nmr.mgh.harvard.edu</a><br>
<a href="https://mail.nmr.mgh.harvard.edu/mailman/listinfo/mne_analysis" rel="noreferrer" target="_blank">https://mail.nmr.mgh.harvard.edu/mailman/listinfo/mne_analysis</a><br>
</blockquote></div><br></div>
_______________________________________________<br>
Mne_analysis mailing list<br>
<a href="mailto:Mne_analysis@nmr.mgh.harvard.edu" target="_blank">Mne_analysis@nmr.mgh.harvard.edu</a><br>
<a href="https://mail.nmr.mgh.harvard.edu/mailman/listinfo/mne_analysis" rel="noreferrer" target="_blank">https://mail.nmr.mgh.harvard.edu/mailman/listinfo/mne_analysis</a><br>
<br>
<br>
The information in this e-mail is intended only for the person to whom it is<br>
addressed. If you believe this e-mail was sent to you in error and the e-mail<br>
contains patient information, please contact the Partners Compliance HelpLine at<br>
<a href="http://www.partners.org/complianceline" rel="noreferrer" target="_blank">http://www.partners.org/complianceline</a> . If the e-mail was sent to you in error<br>
but does not contain patient information, please contact the sender and properly<br>
dispose of the e-mail.<br>
</blockquote></div>