QDR2018A {QDiabetes}

QDiabetes-2018 Risk Calculator (A-Variant)

Description

Calculate the 10-year risk of developing type 2 diabetes, using the A-variant of the QDiabetes-2018 algorithm. This variant does not use fasting plasma glucose or glycated haemoglobin A1c.

Usage

QDR2018A(sex, age, bmi, ht, wt, ethn = "WhiteNA", smoke = "Non", tds = 0,
         fhdm = FALSE, htn = FALSE, cvd = FALSE, gdm = FALSE, pcos = FALSE,
         learn = FALSE, psy = FALSE, ster = FALSE, stat = FALSE, apsy = FALSE)

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-8 and \le14.

fhdm

Family history of diabetes in 1st degree relative.

htn

History of hypertension.

cvd

History of cardiovascular disease.

gdm

History of gestational diabetes.

pcos

History of polycystic ovary syndrome.

learn

History of one or more conditions conveying learning difficulties.

psy

History of schizophrenia or bipolar affective disorder.

ster

History of use of corticosteroids.

stat

History of use of statins.

apsy

History of use of 2nd generation antipsychotics.

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://qdiabetes.org/2018/src.php

References

Hippisley-Cox, J. and Coupland, C. (2017). Development and validation of QDiabetes-2018 risk prediction algorithm to estimate future risk of type 2 diabetes: cohort study. BMJ 359, j5019. doi: 10.1136/bmj.j5019

See Also

QDR2013 - For estimating the 1–10-year risk using QDiabetes-2013.

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
QDR2018A(sex = "Female", age = 76, ht = 1.65, wt = 70)
QDR2018A(sex = "Male", age = seq(25, 80, 5), bmi = 40, ethn = "Other", tds = 5)
QDR2018A(sex = "Female", age = 65, bmi = 35, smoke = c("Non", "Ex", "Light"), fhdm = TRUE)

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

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

## Plotting outputs
age <- seq(25, 80, 5)
risk_m <- QDR2018A(sex = "Male", age = age, bmi = 40)
risk_f <- QDR2018A(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-2018 (A-Variant)",
  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)