diff --git a/include/evalhyd/detail/probabilist/ranks.hpp b/include/evalhyd/detail/probabilist/ranks.hpp
index 3df9d3970bedef05cba405e0e420e1d0f7490a67..97bfa15e856701ad302af0523bcfecf9b4078652 100644
--- a/include/evalhyd/detail/probabilist/ranks.hpp
+++ b/include/evalhyd/detail/probabilist/ranks.hpp
@@ -170,7 +170,7 @@ namespace evalhyd
                             // compute the observed relative frequency
                             // $o_j = \sum_{k \in M_j} r_k$
                             xt::view(o_j, xt::all(), xt::all(), m, e, j) =
-                                    xt::nansum(
+                                    xt::sum(
                                             xt::equal(r_k_masked_sampled, j),
                                             -1
                                     );