data-raw/score_coef.R

In vcastro/CVrisk: Compute Risk Scores for Cardiovascular Diseases

library(tibble)



ascvd_pooled_coef <- tribble(
  ~race,
  ~gender,
  ~ln_age,
  ~ln_age_squared,
  ~ln_totchol,
  ~ln_age_totchol,
  ~ln_hdl,
  ~ln_age_hdl,
  ~ln_treated_sbp,
  ~ln_age_treated_sbp,
  ~ln_untreated_sbp,
  ~ln_age_untreated_sbp,
  ~smoker,
  ~ln_age_smoker,
  ~diabetes,
  ~group_mean,
  ~baseline_survival,
  "white",
  "female",
  -29.799,
  4.884,
  13.54,
  -3.114,
  -13.578,
  3.149,
  2.019,
  0,
  1.957,
  0,
  7.574,
  -1.665,
  0.661,
  -29.18,
  0.9665,
  "aa",
  "female",
  17.114,
  0,
  0.94,
  0,
  -18.92,
  4.475,
  29.291,
  -6.432,
  27.82,
  -6.087,
  0.691,
  0,
  0.874,
  86.61,
  0.9533,
  "white",
  "male",
  12.344,
  0,
  11.853,
  -2.664,
  -7.99,
  1.769,
  1.797,
  0,
  1.764,
  0,
  7.837,
  -1.795,
  0.658,
  61.18,
  0.9144,
  "aa",
  "male",
  2.469,
  0,
  0.302,
  0,
  -0.307,
  0,
  1.916,
  0,
  1.809,
  0,
  0.549,
  0,
  0.645,
  19.54,
  0.8954
)


frs_simple_coef <-
  tribble(
    ~gender,
    ~ln_age,
    ~ln_bmi,
    ~ln_untreated_sbp,
    ~ln_treated_sbp,
    ~smoker,
    ~diabetes,
    ~group_mean,
    ~baseline_survival,
    "male",
    3.11296,
    0.79277,
    1.85508,
    1.92672,
    0.70953,
    0.53160,
    23.9388,
    0.88431,
    "female",
    2.72107,
    0.51125,
    2.81291,
    2.88267,
    0.61868,
    0.77763,
    26.0145,
    0.94833
  )


frs_coef <-
  tribble(
    ~gender,
    ~ln_age,
    ~ln_totchol,
    ~ln_hdl,
    ~ln_untreated_sbp,
    ~ln_treated_sbp,
    ~smoker,
    ~diabetes,
    ~group_mean,
    ~baseline_survival,
    "male",
    3.06117,
    1.12370,
    -0.93263,
    1.93303,
    1.99881,
    0.65451,
    0.57367,
    23.9802,
    0.88936,
    "female",
    2.32888,
    1.20904,
    -0.70833,
    2.76157,
    2.82263,
    0.52873,
    0.69154,
    26.1931,
    0.95012
  )


#### MESA
# Abstracted from:
# https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4603537/table/T2/

mesa_coef <-
  tribble(
    ~age,
    ~gender_male,
    ~race_chinese,
    ~race_aa,
    ~race_hispanic,
    ~diabetes,
    ~smoker,
    ~totchol,
    ~hdl,
    ~hld_med,
    ~sbp,
    ~bp_med,
    ~fh_heartattack,
    ~baseline_survival,
    0.0455,
    0.7496,
    -0.5055,
    -0.2111,
    -0.1900,
    0.5168,
    0.4732,
    0.0053,
    -0.0140,
    0.2473,
    0.0085,
    0.3381,
    0.4522,
    0.99963
  )

mesa_cac_coef <-
  tribble(
    ~age,
    ~gender_male,
    ~race_chinese,
    ~race_aa,
    ~race_hispanic,
    ~diabetes,
    ~smoker,
    ~totchol,
    ~hdl,
    ~hld_med,
    ~sbp,
    ~bp_med,
    ~fh_heartattack,
    ~log1p_cac,
    ~baseline_survival,
    0.0172,
    0.4079,
    -0.3475,
    0.0353,
    -0.0222,
    0.3892,
    0.3717,
    0.0043,
    -0.0114,
    0.1206,
    0.0066,
    0.2278,
    0.3239,
    0.2743,
    0.99833
  )


usethis::use_data(ascvd_pooled_coef, overwrite = TRUE)
usethis::use_data(frs_simple_coef, overwrite = TRUE)
usethis::use_data(frs_coef, overwrite = TRUE)
usethis::use_data(mesa_coef, overwrite = TRUE)
usethis::use_data(mesa_cac_coef, overwrite = TRUE)