Negative returns in AHMI

Author: ms609

DESCRIPTION

@@ -1,7 +1,7 @@
 Package: TreeDist
 Type: Package
 Title: Calculate and Map Distances Between Phylogenetic Trees
-Version: 2.11.1
+Version: 2.11.1.9000
 Authors@R: c(person("Martin R.", "Smith",
                     email = "martin.smith@durham.ac.uk", 
                     role = c("aut", "cre", "cph", "prg"), 

NEWS.md

@@ -1,3 +1,7 @@
+# TreeDist 2.11.1.9000 (2025-11-13)
+
+- `AHMI()` now returns negative values (previously zeroed in error).
+
 # TreeDist 2.11.1 (2025-10-13)
 
 - Improve robustness of `SpectralEigens()` tests.

R/hierarchical_mutual_information.R

@@ -174,6 +174,6 @@ AHMI <- function(tree1, tree2, Mean = max, precision = 0.01, minResample = 36) {
   num <- hmi - ehmi[[1]]
   denom <- M - ehmi[[1]]
   # Return:
-  structure(if (num < sqrt(.Machine$double.eps)) 0 else num / denom,
+  structure(if (abs(num) < sqrt(.Machine$double.eps)) 0 else num / denom,
             sem = .AHMISEM(hmi, M, ehmi[[1]], attr(ehmi, "sem")))
 }