diff --git a/Cargo.toml b/Cargo.toml
index 9330045..d9f462d 100644
--- a/Cargo.toml
+++ b/Cargo.toml
@@ -44,8 +44,8 @@ itertools = "^0.13"
 lazy_static = "^1.4"
 libm = "^0.2"
 macro_const = "^0.1"
-ndarray = "^0.15"
-ndarray-stats = "^0.5"
+ndarray = "^0.16"
+ndarray-stats = "^0.6"
 num-complex = "^0.4" # the same version as used by fftw
 num-traits = "^0.2"
 paste = "1"
@@ -68,7 +68,7 @@ plotters-bitmap = "0.3.3"
 rand = "0.8"
 rand_distr = "0.4"
 rayon = "1.5"
-realfft= "3.1"
+realfft = "3.1"
 rustfft = "6.1"
 serde_json = "1.0"
 serde_test = "1.0"
diff --git a/src/features/median_absolute_deviation.rs b/src/features/median_absolute_deviation.rs
index d4b6e44..bc17000 100644
--- a/src/features/median_absolute_deviation.rs
+++ b/src/features/median_absolute_deviation.rs
@@ -62,7 +62,8 @@ where
         let sorted_deviation: SortedArray<_> =
             ts.m.sample
                 .mapv(|m| T::abs(m - m_median))
-                .into_raw_vec()
+                .as_slice_memory_order()
+                .expect("TimeSeries::m::sample::mapv(...) is supposed to be contiguous")
                 .into();
         Ok(vec![sorted_deviation.median()])
     }