diff --git a/R/utils.R b/R/utils.R
index c5d4d50..a85355e 100644
--- a/R/utils.R
+++ b/R/utils.R
@@ -218,8 +218,8 @@ geometric_mean <- function(
     w_positive <- length(x_positive) / n_x
     x_negative <- x[x < 0]
     w_negative <- length(x_negative) / n_x
-    gm_positive <- exp(sum(log(x_positive)))
-    gm_negative <- -1 * exp(sum(log(abs(x_negative))))
+    gm_positive <- exp(mean(log(x_positive)))
+    gm_negative <- -1 * exp(mean(log(abs(x_negative))))
     gm <- w_positive * gm_positive + w_negative * gm_negative
   } else if (method == "optimized") {
     # The formula is: