Use weighted risk

As per jxx123#63
drozzy · Sep 12, 2023 · 6b0b068 · 6b0b068
1 parent 8a9e530
commit 6b0b068
Showing 1 changed file with 4 additions and 1 deletion.
diff --git a/simglucose/analysis/risk.py b/simglucose/analysis/risk.py
@@ -13,5 +13,8 @@ def risk_index(BG, horizon):
         rh = 10 * fBG[fBG > 0]**2
         LBGI = np.nan_to_num(np.mean(rl))
         HBGI = np.nan_to_num(np.mean(rh))
-        RI = LBGI + HBGI
+        weight_rl = len(rl) / horizon
+        weight_rh = len(rh) / horizon
+
+        RI = weight_rl * LBGI + weight_rh * HBGI
     return (LBGI, HBGI, RI)