add a meaningfull error message for the zero-likelihood case

bips-hb · Jan 17, 2024 · c2d7f27 · c2d7f27
1 parent f0ef208
commit c2d7f27
Showing 1 changed file with 4 additions and 0 deletions.
diff --git a/R/forge.R b/R/forge.R
@@ -106,6 +106,10 @@ forge <- function(
   }
   omega <- omega[wt > 0]
 
+  if (nrow(omega) == 0) {
+    stop("All leaves have zero likelihood.")
+  }
+
   # Draw random leaves with probability proportional to weight
   draws <- data.table(
     'f_idx' = omega[, sample(f_idx, size = n_synth, replace = TRUE, prob = wt)]