Merge pull request #67 from ldecicco-USGS/master

Ignore na's for stats
DOI-USGS · Aug 11, 2020 · 786045d · 786045d
2 parents 8570ea3 + 4e7dc0c
commit 786045d
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,