From 5c063c7cd76d26a453d9465491020f9dd424f1f5 Mon Sep 17 00:00:00 2001 From: Nicolas Barascud <10333715+nbara@users.noreply.github.com> Date: Fri, 8 Dec 2023 15:23:25 +0100 Subject: [PATCH 1/2] [FIX] Make outlier detection robust to nans --- meegkit/utils/denoise.py | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/meegkit/utils/denoise.py b/meegkit/utils/denoise.py index aa7b750f..4317eabc 100644 --- a/meegkit/utils/denoise.py +++ b/meegkit/utils/denoise.py @@ -187,11 +187,11 @@ def find_outlier_trials(X, thresh=None, show=True): else: n_chans, n_trials = X.shape - avg = np.mean(X, axis=-1, keepdims=True) # mean over trials + avg = np.nanmean(X, axis=-1, keepdims=True) # mean over trials d = X - avg # difference from mean - d = np.sum(d ** 2, axis=0) + d = np.nansum(d ** 2, axis=0) - d = d / (np.sum(X ** 2) / n_trials) + d = d / (np.nansum(X ** 2) / n_trials) idx = np.where(d < thresh[0])[0] if show: From b9752833e3eafa52c54a593cbc1fa6913490cb70 Mon Sep 17 00:00:00 2001 From: Nicolas Barascud <10333715+nbara@users.noreply.github.com> Date: Fri, 8 Dec 2023 15:30:10 +0100 Subject: [PATCH 2/2] bump version --- citation.cff | 4 ++-- meegkit/__init__.py | 2 +- 2 files changed, 3 insertions(+), 3 deletions(-) diff --git a/citation.cff b/citation.cff index e103a964..1e696eb3 100644 --- a/citation.cff +++ b/citation.cff @@ -5,7 +5,7 @@ authors: given-names: "Nicolas" orcid: "https://orcid.org/0000-0003-1495-561X" title: "MEEGkit" -version: 0.1.6 +version: 0.1.7 doi: 10.5281/zenodo.10210992 -date-released: 2023-11-28 +date-released: 2023-12-08 url: "https://github.com/nbara/python-meegkit" diff --git a/meegkit/__init__.py b/meegkit/__init__.py index c7a11d48..daa695ad 100644 --- a/meegkit/__init__.py +++ b/meegkit/__init__.py @@ -1,5 +1,5 @@ """M/EEG denoising utilities in python.""" -__version__ = "0.1.6" +__version__ = "0.1.7" from . import asr, cca, detrend, dss, lof, ress, sns, star, trca, tspca, utils