when_is_ramadan <- function(years){
seq_of_dates <- seq_days(
from = paste0(min(years), "-1-1"),
to = paste0(max(years), "-12-31")
)
idates <- sapply(seq_of_dates, ummalqura_algorithm)
holiday_dates <- seq_of_dates[stringr::str_detect(idates, "^1-9-....") == TRUE]
return(holiday_dates)
}
#' Create dummies for Ramadan (estimated by the Umm al-Qura algorithm)
#'
#' Create a dummy variable that takes the value 0 or 1 to indicate the
#' abscence or presence of the start of Ramadan, the ninth month of the Islamic Calender.
#' The day(s) before or after this holiday can also be recoded by using the pre and post arguments.
#' For example, if you also want to recode the days after Ramadan, set
#' post argument to 1.
#'
#' @param dates A vector with dates in the format YYYY-MM-DD.
#' @param pre non-negative integer. To encode the days before the holiday dates.
#' @param post non-negative integer. To encode the days after the holiday dates.
#' @return A dummy (binary vector).
#' @importFrom magrittr %>%
#' @export
#' @examples
#' dates <- seq_days("2021-4-10", "2021-4-29")
#' tibble::tibble(
#' dates,
#' start_of_ramadan = ramadan(dates)
#')
ramadan <- function(dates, pre=0, post=0){
years <- as.integer(unique(lubridate::year(dates)))
holiday_dates <- when_is_ramadan(years)
holiday_dates <- add_intervals(holiday_dates, pre=pre, post=post)
vector <- add_holiday_dummies(dates, holiday_dates)
return(vector)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.