Comments (3)
Hi @balandongiv,
Currently, YASA does not have the capability to work with MNE Epochs. This would require quite a complete refactoring of the function, and I suspect it would also greatly increase computation time. Essentially, we would have to add another nested for loop here in order to loop across all the epochs.
Line 616 in 30d2501
I'm also not sure how we could align the epochs with the hypnogram (which is defined on continuous data). My recommendation would therefore be to keep using the continuous data and create your custom "epoch" mask directly with the hypno
and include
parameters. For example, hypno
could be a vector of integer where each integer refers to an epoch, e.g. hypno=[0, 0, 0, 0, 1, 1, 1, 1, 0, 0, 0, 2, 2, 2, 2, 2, 3, 3, 3, 3, ... 17, 17, 17]
and then you could use yasa.spindles_detect(raw, hypno, include=[1, 2, 3, ... 17])
.
Thanks,
Raphael
from yasa.
Hi @raphaelvallat,
Sorry for not being that clear about my question. The case I have in mind is for continuous non-sleep EEG signal that has been epochs.
One way around that I can think of is by implement something like demonstrated in tutorial .
epochs = mne.read_epochs('some-epo.fif', preload=True, verbose=False)
sp=[yasa.spindles_detect ( epoch.get_data () [0, :, :], epochs.info['sfreq'], ch_names=epochs.ch_names, multi_only=False ) for epoch in epochs]
from yasa.
I see! In that case, I think that your solution is probably the best option. Otherwise, adding support for MNE epochs directly in the function would require extensive refactoring, most probably for no or little speed gain compared to your solution.
from yasa.
Related Issues (20)
- sw_detection, How to improve its accuracy? HOT 1
- Support for MEG spindle detection HOT 1
- Epoch Fragmentation when predicting for short epochs (yasa in mice) HOT 3
- Error while using topoplot() function HOT 4
- I want to know unit of pandpower HOT 1
- extract fractal noise HOT 3
- AssertionError: hypno and include must have same dtype HOT 1
- Spindles and Slow waves summary error HOT 6
- Slow wave detection reproducibility question HOT 2
- Question about Power Spectral Analysis HOT 3
- Invalid REM latency HOT 1
- Add Latency to persistent sleep in sleep_statistics HOT 1
- plot_average plots all channels, not only indicated ones HOT 3
- InconsistentVersionWarning because LabelEncoder from version 0.24.2 of scikit-learn HOT 2
- Possibility of detect the sleep stages on less than 5min of data HOT 4
- Hypno plot errors HOT 1
- Switch to modern python packaging
- Important conversion problem - urgent help needed :) HOT 2
- Finalize evaluation module
- Please help HOT 1
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from yasa.