R/age_codeage.R

In diegogarcilazo/deistools: Tools for Biostatistic datasets developed by DEIS.

#'This is a helper function that transform age and age code to days.
#'@param age num. Age. age will be coerced to dbl.
#'@param code_age code age Assume code_age levels 1 = years, 2 = months,
#'3 = days, 4 = hours, 5 = minutes.

rec_age2day <- function(age,code_age) {

  #checks vars
  stopifnot(is.numeric(age))
  age = as.double(age)
  if (any(age==0)) {warning('There are age equals 0 will be coerced to 1')}

  #coerce 0 values to 1.
  age = dplyr::if_else(age == 0, 1, age)

  #To days
  adays <- dplyr::case_when(
      code_age == 3 ~ age, #y
      code_age == 2 ~ age * 30, #m
      code_age == 1 ~ age * 365.25, #d
      code_age == 4 ~ age / 24 / 30, #hs
      code_age == 5 ~ age / 60 / 24 / 30, #min
      code_age == 9 ~ NaN,
      TRUE ~ NaN)

  return(adays)
}


#' Days to group age.
#' @param adays dbl: days of age.
#' @return ordered factor with age groups.


age_factor <- function(adays) {

breaks = c(0,6,27,330,365.25,730.5,1095.75,
					 1461, 3287.25, 5113.5, 6939.75, 8766, 10592.25, 12418.5, 14244.75,
					 16071, 17897.25, 19723.5, 21549.75, 23376, 25202.25, 27028.5,
					 28854.75, 30681, 43464.75)

labels = c('M1','M2','M3','01','02','03','04','05 - 09',
           '10 - 14','15 - 19','20 - 24','25 - 29',
					 '30 - 34','35 - 39','40 - 44','45 - 49','50 - 54',
					 '55 - 59','60 - 64','65 - 69','70 - 74',
					 '75 - 79','80 - 84','85 y +')

cut(adays, breaks, labels, ordered_result = T)}



#' recode age group from two variables age and code_age.
#' @param age dbl: numeric age.
#' @param code_age dbl: numeric code_age
#' @return ordered factor with age groups.


age_codeage <- function(age, code_age) {

  adays <- rec_age2day(age,code_age)

  return(age_factor(adays))

}