diff --git a/examples/advanced/plot_standardization.py b/examples/advanced/plot_standardization.py
index 34fbe26d..f7864749 100644
--- a/examples/advanced/plot_standardization.py
+++ b/examples/advanced/plot_standardization.py
@@ -128,7 +128,7 @@
 df_vgm = xdem.spatialstats.sample_empirical_variogram(
     values=z_dh.data.squeeze(),
     gsd=dh.res[0],
-    subsample=300,
+    subsample=1000,
     n_variograms=10,
     estimator="dowd",
     random_state=42,
diff --git a/examples/advanced/plot_variogram_estimation_modelling.py b/examples/advanced/plot_variogram_estimation_modelling.py
index a32eddfa..3a03deea 100644
--- a/examples/advanced/plot_variogram_estimation_modelling.py
+++ b/examples/advanced/plot_variogram_estimation_modelling.py
@@ -78,9 +78,7 @@
 # conveniently by :func:`xdem.spatialstats.sample_empirical_variogram`:
 # Dowd's variogram is used for robustness in conjunction with the NMAD (see :ref:`robuststats-corr`).
 
-df = xdem.spatialstats.sample_empirical_variogram(
-    values=dh.data, gsd=dh.res[0], subsample=100, n_variograms=10, estimator="dowd", random_state=42
-)
+df = xdem.spatialstats.sample_empirical_variogram(values=dh, subsample=1000, n_variograms=10, estimator="dowd", random_state=42)
 
 # %%
 # *Note: in this example, we add a* ``random_state`` *argument to yield a reproducible random sampling of pixels within