Ignore na's for stats

DOI-USGS · Aug 11, 2020 · 4e7dc0c · 4e7dc0c
1 parent 91a8eee
commit 4e7dc0c
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/R/kalman_bootstrapping.R b/R/kalman_bootstrapping.R
@@ -86,8 +86,8 @@ makeMonths <- function(dailyBootOut, eList){
   monthFluxPcts <- matrix(data = NA, nrow = nMonths, ncol = 11)
   monthConcPcts <- matrix(data = NA, nrow = nMonths, ncol = 11)
   for(iMonth in 1:nMonths){
-    monthFluxPcts[iMonth,] <- quantile(monthsFlux[iMonth,], probs = probs, type = 6)
-    monthConcPcts[iMonth,] <- quantile(monthsConc[iMonth,], probs = probs, type = 6)
+    monthFluxPcts[iMonth,] <- quantile(monthsFlux[iMonth,], probs = probs, type = 6, na.rm = TRUE)
+    monthConcPcts[iMonth,] <- quantile(monthsConc[iMonth,], probs = probs, type = 6, na.rm = TRUE)
   }
   monthPcts <- list(monthFluxPcts = monthFluxPcts, 
                     monthConcPcts = monthConcPcts,