fix to anlz_fibmatrix if geomean calculated on NA

tbep-tech · Oct 20, 2024 · 65fc10c · 65fc10c
1 parent a32aeab
commit 65fc10c
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 2 deletions.
diff --git a/DESCRIPTION b/DESCRIPTION
@@ -1,7 +1,7 @@
 Package: tbeptools
 Title: Data and Indicators for the Tampa Bay Estuary Program
 Version: 3.0.0.9000
-Date: 2024-10-10
+Date: 2024-10-20
 Authors@R: c(
     person(given = "Marcus",
       family = "Beck",

diff --git a/R/anlz_fibmatrix.R b/R/anlz_fibmatrix.R
@@ -281,7 +281,8 @@ anlz_fibmatrix <- function(fibdata, yrrng = NULL, stas = NULL, bay_segment = NUL
       exceed_30_prob = pbinom(sumgt - 1, tot, 0.30, lower.tail = FALSE),
       exceed_50_prob = pbinom(sumgt - 1, tot, 0.50, lower.tail = FALSE),
       exceed_75_prob = pbinom(sumgt - 1, tot, 0.75, lower.tail = FALSE)
-    )
+    ) %>%
+    dplyr::filter(!is.na(sumgt))
 
   # Put stations into binomial test groups based on significant exceedances of 400cfu criteria
   dat$MWQA <- NA