R/RH_DELP_C_PUT.R

In surveyPrev: Mapping the Prevalence of Binary Indicators using Survey Data in Small Areas

##'RH_DELP_C_PUT
#' @param Rdata  data.frame from surveyPrev::getDHSdata
#'
#' @return A partially processed data.frame that will be used in surveyPrev::getDHSindicator. The whole function can be used as a parameter in surveyPrev::getDHSindicator
#'
#' @author Miaolei Bao, Yunhan Wu, Qianyu Dong
#' @examples
#' \dontrun{
#' dhsData <- getDHSdata(country = "Zambia", indicator = "RH_DELP_C_PUT", year = 2018)
#' }
#'
#' @export
RH_DELP_C_PUT<- function(Rdata){

  #Manually added function by Qianyu
  #RH_DELP_C_PUT is made up of the two year version for RH_DELP_C_PUB(five year version/ official id name)

  # /*****************************************************************************************************
  # Program: 			RH_DEL.R
  # Purpose: 			Code Delivery Care indicators
  # Data inputs: 	BR dataset
  # Data outputs:	coded variables
  # Author:			  Shireen Assaf
  # Date last modified: Sept 10, 2021 by Shireen Assaf
  # *****************************************************************************************************/
  #
  # /*----------------------------------------------------------------------------//
  # Variables created in this file:
  # rh_del_place		"Live births by place of delivery"
  # rh_del_pltype		"Live births by type of place"
  # rh_del_pv			  "Person providing assistance during birth"
  # rh_del_pvskill	"Skilled provider providing assistance during birth"
  # rh_del_ces			"Live births delivered by cesarean"
  # rh_del_cestime	"Timing of decision to have Cesarean"
  # rh_del_stay			"Duration of stay following recent birth"
  # /----------------------------------------------------------------------------*/
  #
  BRdata <- Rdata %>%
    mutate(wt = v005/1000000)

  # period and age of child
  # choose reference period, last 2 years (24 months) or last 5 years (60 months)
  # Using a period of the last 2 years will not match final report but would provide more recent information.
  BRdata <- BRdata %>%
    mutate(period = 24)

  # age of child. If b19 is not available in the data use v008 - b3
  if ("TRUE" %in% (!("b19" %in% names(BRdata))))
    BRdata [[paste("b19")]] <- NA
  if ("TRUE" %in% all(is.na(BRdata$b19)))
  { b19_included <- 0} else { b19_included <- 1}

  if (b19_included==1) {
    BRdata <- BRdata %>%
      mutate(age = b19)
  } else {
    BRdata <- BRdata %>%
      mutate(age = v008 - b3)
  }

  # ** NEW STEP: Explicitly filter the data to keep only recent births **
  # This directly changes the number of rows (the sample size).
  BRdata <- BRdata %>%
    filter(age < period)


  # //Place of delivery
  # Note: please check the categories of m15 especially for older surveys. The category values may differ.
  BRdata <- BRdata %>%
    mutate(rh_del_pltype =
             case_when(
               m15 >=20 & m15<30  ~ 1 ,
               m15 >=30 & m15<50  ~ 2 ,
               m15 >=10 & m15<20  ~ 3,
               m15 >=50 & m15<99  ~ 4 ,
               m15 == 99 ~ 9 ,
               age>=period ~ 99)) %>%
    replace_with_na(replace = list(rh_del_pltype = c(99))) %>%
    set_value_labels(rh_del_pltype = c("Health facility - public" = 1, "Health facility - private" = 2, "Home"=3, "Other"=4, "Missing"=9 )) %>%
    set_variable_labels(rh_del_pltype = "Live births by type of health facility-public")




  BRdata <- BRdata %>%
    mutate(RH_DELP_C_PUT =
             case_when(
               rh_del_pltype==1 ~ 1,
               rh_del_pltype %in% c(2,3,4,9)   ~ 0 )) %>%
    set_value_labels(RH_DELP_C_PUT = c("Yes" = 1, "No"=0)) %>%
    set_variable_labels(RH_DELP_C_PUT = "Live births by place of delivery: Health facility-public")

  colnames(BRdata)[colnames(BRdata) == 'RH_DELP_C_PUT'] <- 'value'
  return(BRdata)
}