About the package
CVrisk is a package to calculate various cardiovascular disease risk scores.
Currently available risk scores
- Framingham 2008 10-year ASCVD risk (model with lipid labs)
- Framingham 2008 10-year ASCVD risk (model with BMI)
- ACC/AHA 2013 10-year ASCVD risk
- MESA 2015 10-year ASCVD risk (traditional risk factors)
- MESA 2015 10-year ASCVD risk (using coronary artery calcium)
- AHA PREVENT 2023 10-year ASCVD risk (base model) - New!
Installation
# Install from CRAN
install.packages("CVrisk")
# Install the development version from GitHub:
# install.packages("devtools")
devtools::install_github("vcastro/CVrisk")Usage
Calculate a single score:
library(CVrisk)
ascvd_10y_accaha(race = "aa", gender = "male", age = 55,
totchol = 213, hdl = 50, sbp = 140,
bp_med = FALSE, smoker=0, diabetes=0)
# Calculate PREVENT risk score (requires additional parameters)
ascvd_10y_prevent(gender = "female", age = 50, sbp = 160, bp_med = 1,
totchol = 200, hdl = 45, statin = 0, diabetes = 1, smoker = 0,
egfr = 90, bmi = 35)
Calculate multiple scores for a dataframe and append scores to the dataframe:
# Compute traditional risk scores
compute_CVrisk(sample_data,
age = "age", race = "race", gender = "gender", bmi = "BMI", sbp = "sbp",
hdl = "hdl", totchol = "totchol", bp_med = "bp_med", smoker = "smoker",
diabetes = "diabetes", lipid_med = "lipid_med",
fh_heartattack = "fh_heartattack", cac = "cac"
)
# To include PREVENT score, specify scores explicitly and provide egfr
compute_CVrisk(my_data,
scores = c("ascvd_10y_accaha", "ascvd_10y_prevent"),
age = "age", gender = "gender", sbp = "sbp",
hdl = "hdl", totchol = "totchol", bmi = "BMI",
bp_med = "bp_med", diabetes = "diabetes", smoker = "smoker",
statin = "lipid_med", egfr = "egfr"
)