Vectorised function to calculate US KDPI percentile as published by OPTN/UNOS. Please note that this function uses creatinine measured in µmol/l by default, but can be changed to mg/dl if the optional units parameter is set to "US".

kdpi(age, height, weight, eth, htn, dm, cva, creat, hcv, dcd,
  scaling = 1, units = "SI")

Arguments

age

numeric vector of donor ages in years

height

numeric vector of donor heights in cm

weight

numeric vector of donor weights in kg

eth

character string vector of donor ethnicity, "black" or "non-black"

htn

numeric vector of donor hypertension history (1 = yes, 0 = no)

dm

numeric vector of donor diabetes history (1 = yes, 0 = no)

cva

numeric vector of whether donor death due to CVA (1 = yes, 0 = no)

creat

numeric vector of donor serum creatinine (µmol/l)

hcv

numeric vector of donor hepatitis C history (1 = yes, 0 = no)

dcd

numeric vector of type of donor (1 = DCD, 0 = DBD)

scaling

single numeric value for OPTN scaling factor (optional, defaults to 1)

units

single string value to indicate creatinine units ("SI" for µmol/l, "US" for mg/dl)

Value

numeric vector of US KDRI values

Details

The KDRI is calculated from KDRI normalised by a scaling factor based on the median KDRI in the previous year. For 2018, this was approximately 1.250609 as is published on the OPTN website. The scaling parameter in this function defaults to 1, so can be left out to calculate the KDPI without scaling.

This function requires the dplyr package to be installed.

Reference: Rao PS, Schaubel DE, Guidinger MK, et al. A Comprehensive Risk Quantification Score for Deceased Donor Kidneys: The Kidney Donor Risk Index. Transplantation 2009; 88:231-236.

Examples

# with creatinine in µmol/l (units = "SI" can be omitted) kdpi(age = 40, height = 170, weight = 80, eth = "non-black", htn = 0, dm = 0, cva = 0, creat = 120, hcv = 0, dcd = 0, scaling = 1.250609, units = "SI")
#> [1] 35
# with creatinine in mg/dl kdpi(age = 40, height = 170, weight = 80, eth = "non-black", htn = 0, dm = 0, cva = 0, creat = 1.4, hcv = 0, dcd = 0, scaling = 1.250609, units = "US")
#> [1] 36