diff --git a/include/evalhyd/detail/probabilist/brier.hpp b/include/evalhyd/detail/probabilist/brier.hpp
index 5e7e73db9f65b635c57b906d1280f94e7baa9d9b..e39e523acf1b720e3e4e05557bb980c5b9cd5160 100644
--- a/include/evalhyd/detail/probabilist/brier.hpp
+++ b/include/evalhyd/detail/probabilist/brier.hpp
@@ -412,7 +412,7 @@ namespace evalhyd
                         );
 
                         // compute number of forecasts in each forecast bin $N_i$
-                        auto N_i = xt::eval(xt::nansum(msk_bins, -1));
+                        auto N_i = xt::eval(xt::sum(msk_bins, -1));
 
                         xt::view(REL_DIAG, xt::all(), xt::all(), m, e, xt::all(),
                                  xt::all(), 2) = N_i;