QDR2013 {QDiabetes}

QDiabetes-2013 Risk Calculator

Description

Calculate the risk of developing type 2 diabetes, using the QDiabetes-2013 algorithm.

Usage

QDR2013(sex, age, bmi, ht, wt, ethn = "WhiteNA", smoke = "Non", tds = 0,
        fhdm = FALSE, htn = FALSE, cvd = FALSE, ster = FALSE, surv = 10L)

Arguments

sex

Gender. Must be "Male" or "Female".

age

Age, in years. Must be \ge25 and <85.

bmi

Body-mass index, in kg/m^2. Must be \ge40/2.1^2 and \le180/1.4^2. Within this range, values >40 are set to 40, while values <20 are set to 20.

ht

Height, in m. Must be \ge1.4 and \le2.1.

wt

Weight, in kg. Must be \ge40 and \le180.

ethn

Ethnicity. Must be one of:

• "WhiteNA" (White or not stated)

• "Indian" (Indian)

• "Pakistani" (Pakistani)

• "Bangladeshi" (Bangladeshi)

• "OtherAsian" (Other Asian)

• "BlackCaribbean" (Black Caribbean)

• "BlackAfrican" (Black African)

• "Chinese" (Chinese)

• "Other" (Other ethnic group).

smoke

Smoking status. Must be one of:

• "Non" (Non-smoker)

• "Ex" (Ex-smoker)

• "Light" (Light smoker - less than 10 cigarettes per day)

• "Moderate" (Moderate smoker - 10 to 19 cigarettes per day)

• "Heavy" (Heavy smoker - 20 or more cigarettes per day).

tds

Townsend deprivation score. Must be \ge-7 and \le11.

fhdm

Family history of diabetes in 1st degree relative.

htn

History of hypertension.

cvd

History of cardiovascular disease.

ster

History of use of corticosteroids.

surv

Time point at which the Kaplan-Meier cumulative survival estimate is used to calculate risk. Must be an integer value between 1 and 10.

Value

Risk (%) of type 2 diabetes.

Note

The QDiabetes R package has been created as a research tool for scientific purposes only. The QDiabetes R package has not been granted Medicines and Healthcare products Regulatory Agency (MHRA) approval as a medical device, and hence, should not be used as part of any individualised risk assessment.

Author(s)

Benjamin G. Feakins [e-mail]

Source

https://svn.clinrisk.co.uk/opensource/qdiabetes/standard/

See Also

QDR2018A - For estimating the 10-year risk using QDiabetes-2018 (Model A).

QDR2018B - For estimating the 10-year risk using QDiabetes-2018 (Model B).

QDR2018C - For estimating the 10-year risk using QDiabetes-2018 (Model C).

Examples

## Simple usage
QDR2013(sex = "Female", age = 76, ht = 1.65, wt = 70)
QDR2013(sex = "Male", age = seq(25, 80, 5), bmi = 40, ethn = "Other", tds = 5)
QDR2013(sex = "Female", age = 65, bmi = 35, smoke = c("Non", "Ex", "Light"), fhdm = TRUE)

## Using postcodes to estimate Townsend deprivation scores
QDR2013(sex = "Male", age = 65, bmi = 40, tds = getTDS("OX3 7LF"))
QDR2013(sex = "Female", age = 60, bmi = 35, tds = getTDS(c("OX2 6NW", "OX2 6GG")))

## Data frame usage
data(dat_qdr)
with(dat_qdr, QDR2013(sex = sex, age = age, bmi = bmi))

## Plotting outputs
age <- seq(25, 80, 5)
risk_m <- QDR2013(sex = "Male", age = age, bmi = 40)
risk_f <- QDR2013(sex = "Female", age = age, bmi = 40)

oldpar <- par(no.readonly = TRUE)
par(cex = 0.8, cex.sub = 0.8)
plot.new()
plot.window(xlim = range(age), ylim = range(pretty(c(risk_m, risk_f))))
axis(1, at = age)
axis(2, at = pretty(c(risk_m, risk_f)))
title(
  main = "Diabetes 10-year risk with age:\nQDiabetes-2013",
  sub = expression("BMI set to"~40*kg/m^2*", other variables set to defaults"),
  xlab = "Age (years)", ylab = "Risk (%)"
)
lines(age, risk_m, type = "b", col = "navy", lwd = 1.5)
lines(age, risk_f, type = "b", col = "red3", lwd = 1.5)
legend("bottomright", legend = c("Male", "Female"), col = c("navy", "red3"), lty = 1, bty = "n")
par(oldpar)