Fixed a bug related to geometric mean calculation.

PathwayAndDataAnalysis · Jan 25, 2024 · 60d4e42 · 60d4e42
1 parent 1942ec4
commit 60d4e42
Showing 1 changed file with 5 additions and 1 deletion.
diff --git a/src/main/java/org/panda/causalpath/analyzer/ThresholdDetector.java b/src/main/java/org/panda/causalpath/analyzer/ThresholdDetector.java
@@ -88,6 +88,8 @@ public OneDataChangeDetector makeACopy()
 	 */
 	public double foldChangeGeometricMean(double[] vals)
 	{
+		if (vals.length == 1) return vals[0];
+
 		double mult = 1;
 		int cnt = 0;
 		for (double val : vals)
@@ -96,7 +98,9 @@ public double foldChangeGeometricMean(double[] vals)
 			cnt++;
 			mult *= val < 0 ? -1 / val : val;
 		}
-		return Math.pow(mult, 1D / cnt);
+		double result = Math.pow(mult, 1D / cnt);
+		if (result < 1) result = - 1 / result;
+		return result;
 	}
 
 	public double maxOfAbs(double[] vals)