ENH - use mad() function from stats toolbox

fieldtrip · Feb 14, 2024 · 7ba0239 · 7ba0239
1 parent 6cebb4b
commit 7ba0239
Showing 1 changed file with 1 addition and 2 deletions.
diff --git a/private/artifact_level.m b/private/artifact_level.m
@@ -82,8 +82,7 @@
   case 'maxzvalue'
     level = nanmax((dat-repmat(mval, 1, nsample))./repmat(sd, 1, nsample), [], 2);
   case 'mad'
-    mval = median(dat, 2, 'omitmissing'); % the median per channel, over time
-    level = median(abs(dat-repmat(mval, 1, nsample)), 2, 'omitmissing'); % median absolute deviation
+    level = mad(dat, 1, 2);
 
   case 'neighbexpvar'
     % this results in a Nx1 vector