nothing much

R/cwd.R

@@ -109,8 +109,15 @@ cwd <- function(df, varname_wbal, varname_date, thresh_terminate = 0.0, thresh_d
       }
 
       # record instance
-      this_inst <- tibble( idx_start = idx, len = iidx_drop-idx, iinst = iinst, date_start=df[[varname_date]][idx], date_end = df[[varname_date]][iidx_drop-1], deficit = max_deficit )
-      inst <- inst |> bind_rows(this_inst)
+      this_inst <- tibble( idx_start = idx,
+                           len = iidx_drop-idx,
+                           iinst = iinst,
+                           date_start = df[[varname_date]][idx],
+                           date_end = df[[varname_date]][iidx_drop-1],
+                           deficit = max_deficit
+                           )
+      inst <- inst |>
+        bind_rows(this_inst)
 
       # update
       iinst <- iinst + 1

vignettes/cwd_example.Rmd

@@ -164,12 +164,12 @@ df_return |>
   labs(x = "Return period (yr)", y = "Magnitude of annual CWD maximum (mm)")
 ```
 
-Visualise the estimated event size with a return period of $T = 80$ y on top of the distribution of cumulative water deficit events.
+Visualise the estimated event size with a return period of $T = 80$ y as the red line on top of the distribution of cumulative water deficit events.
 ```{r}
 ggplot() +
   geom_histogram(
     data = out_cwd$inst,
-    aes(x = deficit, y = ..density..),
+    aes(x = deficit, y = after_stat(density)),
     color = "black",
     position="identity",
     bins = 6