[Mne_analysis] Bipolar EEG and Epochs rejection

pie.min at libero.it pie.min at libero.it
Fri Mar 20 05:04:57 EDT 2020
Search archives:

        External Email - Use Caution        

Thank you for you reply!


> Il 19 marzo 2020 alle 1.21 Dan McCloy <dan at mccloy.info> ha scritto:
> 
> 
>             External Email - Use Caution        
> 
> 
>     I think this might be possible, if you change one of them (say, the HEOG) to type 'misc', then you could use separate reject criteria for EOG and MISC channel types.  Only works if you don't have other MISC channels already, of course. Example: https://gist.github.com/drammock/475d6a0aef16468d2013f145d7dd0151
> 
>     -- dan
>     Daniel McCloy
>     https://dan.mccloy.info
>     Research Scientist
>     Institute for Learning and Brain Sciences
>     University of Washington
> 
> 
>     ‐‐‐‐‐‐‐ Original Message ‐‐‐‐‐‐‐
>     On Wednesday, March 18, 2020 12:11 PM, Phillip Alday <phillip.alday at mpi.nl> wrote:
> 
> 
>         > > 
> >         I don't think there's a way to do this via the Epochs constructor. There are two questions I have:
> > 
> >         (1) Did BVA use  peak-to-peak and or absolute threshold? I'm guessing absolute (see next question), but the Epochs constructor uses peak-to-peak.
> > 
> >         (2) Did you mean mV (1e-3) or µV (1e-6)? You use mV in your text, but the scaling in your example is µV.
> > 
> >         If you want absolute thresholds and not peak-to-peak, checkout philistine.abs_threshold. Assuming you want mV
> > 
> >         https://philistine.readthedocs.io/en/latest/api/philistine.mne.abs_threshold.html
> > 
> >         epochs = ....
> > 
> >         hmask = abs_threshold(epochs.copy().pick("HEOG"), 30e-3, eog=True)
> > 
> >         vmask = abs_threshold(epochs.copy().pick("VEOG"), 120e-3, eog=True)
> > 
> >         eye_mask = np.logical_or(hmask, vmask)
> > 
> >         epochs.drop(eye_mask, reason="EOG")
> > 
> >         (note that I'm the author of philistine)
> > 
> >         Best,
> > 
> >         Phillip
> > 
> >         On 18/03/2020 16:53, pie.min at libero.it mailto:pie.min at libero.it wrote:
> > 
> >             > > > 
> > >                     External Email - Use Caution        
> > > 
> > >             Goodmorning,
> > > 
> > >             I am new to python and MNE.
> > > 
> > >             I am trying to replicate a brain vision analyzer pre-processing in MNE. In this pre-processing authors rejected epochs in which horizontal EOG exceeded 30 mV and/or vertical EOG exceeded 60 mV. I'm trying to achieve this result using the command mne.Epochs. However, in the reject dictionary I don't know how to specify this "double condition".
> > > 
> > >             Ideally I would to write somenthing like that:
> > > 
> > >             reject = dict(eeg = 160e-6, ["VEOG"] = 120e-6, ["HEOG"] = 30e-6)
> > > 
> > >             epochs = mne.Epochs(raw, events=events[0], event_id=[71, 72], tmin=-0.1, tmax=0.6, proj=True, baseline=(-0.1,0), reject=reject, preload=True)
> > > 
> > >             However in reject dictionary, if I understood correctly, I can only use one key 'eog'.
> > > 
> > >             Thank you very much for your help,
> > >             Any advice will be usefull,
> > >             Bests,
> > >             Pietro
> > > 
> > > 
> > >             _______________________________________________
> > >             Mne_analysis mailing list
> > >             Mne_analysis at nmr.mgh.harvard.edu mailto:Mne_analysis at nmr.mgh.harvard.edu
> > >             https://mail.nmr.mgh.harvard.edu/mailman/listinfo/mne_analysis
> > > 
> > >         > > 
> >                 External Email - Use Caution        
> > 
> > 
> >     > 
>     _______________________________________________
>     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/20200320/8574b227/attachment.html 


More information about the Mne_analysis mailing list