# Function to compute eGFR 2021
egfr_2021 <- function(age, female, cr) {
  return(142 * min(cr / (0.9-0.2*female), 1)^(-0.302+0.061*female) *
           max(cr/(0.9-0.2*female), 1)^(-1.2) * 
           0.9938^age *
           1.012^female)
}