From 7d9210b177f5048bfc0a68ef1f065b42d47418b7 Mon Sep 17 00:00:00 2001
From: Thibault Hallouin <thibault.hallouin@inrae.fr>
Date: Mon, 23 Jan 2023 11:39:24 +0100
Subject: [PATCH] add details on calculation of delta score

---
 include/evalhyd/detail/probabilist/ranks.hpp | 4 +++-
 1 file changed, 3 insertions(+), 1 deletion(-)

diff --git a/include/evalhyd/detail/probabilist/ranks.hpp b/include/evalhyd/detail/probabilist/ranks.hpp
index a594b8f..376a80c 100644
--- a/include/evalhyd/detail/probabilist/ranks.hpp
+++ b/include/evalhyd/detail/probabilist/ranks.hpp
@@ -303,7 +303,9 @@ namespace evalhyd
                         auto l = xt::sum(t_msk_sampled, -1);
 
                         // compute the Delta score
-                        // $\bar{o_j} = \frac{1}{n} \sum_{k \in M_j} r_k$
+                        // \Delta = \sum_{k=1}^{N+1} (r_k - \frac{M}{N+1})^2
+                        // \Delta_o = \frac{MN}{N+1}
+                        // \delta = $\frac{\Delta}{\Delta_o}
                         xt::view(DS, xt::all(), xt::all(), m, e) =
                                 xt::nansum(
                                         xt::square(
-- 
GitLab