Description Usage Arguments Value Note Author(s) Source See Also Examples
Calculate the risk of developing type 2 diabetes, using the QDiabetes-2013 algorithm.
1 2 |
sex |
Gender. Must be |
age |
Age, in years. Must be ≥25 and <85. |
bmi |
Body-mass index, in kg/m^2. Must be ≥40/2.1^2 and ≤180/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 ≥1.4 and ≤2.1. |
wt |
Weight, in kg. Must be ≥40 and ≤180. |
ethn |
Ethnicity. Must be one of:
|
smoke |
Smoking status. Must be one of:
|
tds |
Townsend deprivation score. Must be ≥-7 and ≤11. |
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. |
Risk (%) of type 2 diabetes.
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.
Benjamin G. Feakins benjamin.feakins@ndph.ox.ac.uk
https://svn.clinrisk.co.uk/opensource/qdiabetes/standard/
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).
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 | ## 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)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.