[Mne_analysis] howto convert ICA components to SSP projection operators.

[Dan McCloy]

Luke,
I've created projectors "by hand" like this (caveat: my intent was just to
plot topomaps, not to actually apply the projector):

desc = '4Hz_artifact_ICA'
data = dict(data=my_component, nrow=1, ncol=len(my_component),
            row_names=desc, col_names=epochs.ch_names)
projector = mne.io.Projection(data=data, active=False, kind=1,
                              desc=desc, explained_var=None)
epochs.add_proj(projector)

I've never used it with ICA (only with DSS) but I think in principle it
could work the same way?
-- dan

Daniel McCloy
http://dan.mccloy.info/
Postdoctoral Research Associate
Institute for Learning and Brain Sciences
University of Washington


On Wed, Aug 16, 2017 at 9:53 AM, Alexandre Gramfort <
alexandre.gramfort at inria.fr> wrote:

> Hi Luke,
>
> we don't have such a function AFAIK and I see a use to it. My first
> reaction was to tall you to use a high n_pca_components in ICA but
> it's not fully equivalent to an SSP.
>
> proj = make_ica_proj(ica, picks)
>
> ?
>
> Alex
>
> On Wed, Aug 16, 2017 at 3:50 PM, Luke Bloy <luke.bloy at gmail.com> wrote:
> > Before I potentially reinvent a wheel, I was wondering if there was a
> way to
> > compute an ssp projection operator (proj) from the excluded basis set of
> an
> > ICA operator.
> >
> > The reason I'm trying to do this is that I have some data that has an
> > environmental artifact that is always around 4HZ but whose topology
> changes
> > slightly from subject to subject. ICA on a filtered dataset (3.5-4.5Hz)
> > makes it easy to ID the artifact components but the remaining ICA
> components
> > are not really applicable to the broadband signal so using ica.apply() on
> > the boradband signal doesn't really work for me.
> >
> > Thanks
> > Luke
> >
> > _______________________________________________
> > Mne_analysis mailing list
> > Mne_analysis at nmr.mgh.harvard.edu
> > https://mail.nmr.mgh.harvard.edu/mailman/listinfo/mne_analysis
> >
> >
> > The information in this e-mail is intended only for the person to whom
> it is
> > addressed. If you believe this e-mail was sent to you in error and the
> > e-mail
> > contains patient information, please contact the Partners Compliance
> > HelpLine at
> > http://www.partners.org/complianceline . If the e-mail was sent to you
> in
> > error
> > but does not contain patient information, please contact the sender and
> > properly
> > dispose of the e-mail.
> >
> _______________________________________________
> Mne_analysis mailing list
> Mne_analysis at nmr.mgh.harvard.edu
> https://mail.nmr.mgh.harvard.edu/mailman/listinfo/mne_analysis
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mail.nmr.mgh.harvard.edu/pipermail/mne_analysis/attachments/20170816/8c1430ac/attachment-0001.html