diff --git a/scipy-stubs/ndimage/_measurements.pyi b/scipy-stubs/ndimage/_measurements.pyi
index 1f2a9302..e3854b15 100644
--- a/scipy-stubs/ndimage/_measurements.pyi
+++ b/scipy-stubs/ndimage/_measurements.pyi
@@ -88,7 +88,7 @@ def labeled_comprehension(
     labels: onp.ToComplex | onp.ToComplexND | None,
     index: onp.ToInt | onp.ToIntND | None,
     func: _ComprehensionFunc,
-    out_dtype: onp.AnyFloat64DType,
+    out_dtype: onp.AnyFloat64DType | None,
     default: onp.ToFloat,
     pass_positions: bool = False,
 ) -> onp.ArrayND[np.float64]: ...
diff --git a/scipy-stubs/stats/_stats.pyi b/scipy-stubs/stats/_stats.pyi
index e5a5e124..47142fb1 100644
--- a/scipy-stubs/stats/_stats.pyi
+++ b/scipy-stubs/stats/_stats.pyi
@@ -111,7 +111,7 @@ def gaussian_kernel_estimate(
     values: onp.Array2D[_Real],
     xi: onp.Array2D[_AsReal],
     cho_cov: onp.Array2D[_AsReal],
-    dtype: onp.AnyFloat64DType,
+    dtype: onp.AnyFloat64DType | None,
     _: _Real | float = 0.0,
 ) -> onp.Array2D[np.float64]: ...
 @overload
@@ -140,7 +140,7 @@ def gaussian_kernel_estimate_log(
     values: onp.Array2D[_Real],
     xi: onp.Array2D[_AsReal],
     cho_cov: onp.Array2D[_AsReal],
-    dtype: onp.AnyFloat64DType,
+    dtype: onp.AnyFloat64DType | None,
     _: _Real | float = 0.0,
 ) -> onp.Array2D[np.float64]: ...
 @overload