R/day_ana.R

In ERottler/meltimr: Shiny App to investigate GRDC runoff data with regard to changes in runoff timing

#' Analysis on daily level
#'
#' Apply functions on daily level. First, time series re-structured according to day of the year, then selected analytical function applied.
#'
#' @param data  Vector of input data
#' @param date Time vector
#' @param start_year Selected start year
#' @param end_year Seleced end year
#' @param do_ma Apply moving average filter before analysis (FALSE / TRUE)
#' @param window_width If do_ma is TRUE, then moving average filter with here selected window width applied
#' @param method_ana Selected analytical method('mean', 'median', 'sum', 'sens_slope')
#' @param break_day Define start day (e.g. 274 is 1.October, which is the start of  hydrological year in Switzerland).
#' @export
day_ana <- function(data, date, start_year, end_year, do_ma = F, window_width = 30, method_ana, break_day = 0){

  data_day <- ord_day(data_in = data, date = date,  start_y = start_year, end_y = end_year,
                        do_ma = do_ma, window_width = window_width, break_day = break_day)

  if(method_ana == "mean"){
    res <- apply(data_day, 2, mea_na)
  }

  if(method_ana == "median"){
    res <- apply(data_day, 2, med_na)
  }

  if(method_ana == "sum"){
    f_sum <- function(data_in){sum(data_in, na.rm = T)}
    res <- apply(data_day, 2, f_sum)
  }

  if(method_ana == "sens_slope"){

    res <- apply(data_day, 2, sens_slo)

  }

  return(res)

}