R/02-merge.R

In validmeasures/ramOP: An Implementation of Rapid Assessment Method for Older People

################################################################################
#
#' Concatenate classic and PROBIT estimates into a single data.frame
#'
#' @param x Classic estimates data frame
#' @param y Probit estimates data frame
#' @param prop2percent Logical. Should proportion type indicators be converted
#'   to percentage? Default is FALSE.
#'
#' @return Data frame of combined classic and probit estimates
#'
#' @author Ernest Guevarra
#'
#' @examples
#'   #
#'   \dontrun{
#'   test <- merge_estimates(x = classicEstimates,
#'                           y = probitEstimates)
#'   }
#'
#' @export
#'
#
################################################################################

merge_estimates <- function(x, y, prop2percent = FALSE) {

  ## Merge rows
  estimates <- rbind(x, y)

  ## Merge 'estimates' data.frame and 'language' data.frame in prepartion for
  ## reporting and maintaining the original row ordering of the 'language'
  ## data.frame ...
  temp <- subset(language, subset = language$INDICATOR %in% estimates$INDICATOR)

  temp$originalOrder <- seq_len(nrow(estimates))

  estimates <- merge(temp, estimates, by = "INDICATOR")
  estimates <- estimates[order(estimates$originalOrder), ]
  estimates <- subset(estimates, select = -originalOrder)

  if(prop2percent) {
    estimates$EST.ALL <- ifelse(estimates$TYPE == "Proportion",
                                estimates$EST.ALL * 100, estimates$EST.ALL)
    estimates$LCL.ALL <- ifelse(estimates$TYPE == "Proportion",
                                estimates$LCL.ALL * 100, estimates$LCL.ALL)
    estimates$UCL.ALL <- ifelse(estimates$TYPE == "Proportion",
                                estimates$UCL.ALL * 100, estimates$UCL.ALL)
    estimates$EST.MALES <- ifelse(estimates$TYPE == "Proportion",
                                  estimates$EST.MALES * 100, estimates$EST.MALES)
    estimates$LCL.MALES <- ifelse(estimates$TYPE == "Proportion",
                                  estimates$LCL.MALES * 100, estimates$LCL.MALES)
    estimates$UCL.MALES <- ifelse(estimates$TYPE == "Proportion",
                                  estimates$UCL.MALES * 100, estimates$UCL.MALES)
    estimates$EST.FEMALES <- ifelse(estimates$TYPE == "Proportion",
                                    estimates$EST.FEMALES * 100, estimates$EST.FEMALES)
    estimates$LCL.FEMALES <- ifelse(estimates$TYPE == "Proportion",
                                    estimates$LCL.FEMALES * 100, estimates$LCL.FEMALES)
    estimates$UCL.FEMALES <- ifelse(estimates$TYPE == "Proportion",
                                    estimates$UCL.FEMALES * 100, estimates$UCL.FEMALES)
  }

  ## Convert to tibble
  estimates <- tibble::tibble(estimates)

  ## Return
  estimates
}