set.seed(1234)
data(d_sim_binom_cov_hgly2)
ermod_bin <- suppressWarnings(dev_ermod_bin(
data = d_sim_binom_cov_hgly2, var_resp = "AEFLAG",
var_exposure = "AUCss_1000", var_cov = c("BGLUC", "RACE"),
verbosity_level = 0,
# Below option to make the example run fast
chains = 2, iter = 1000
))
spec_coveff <- build_spec_coveff(ermod_bin)
spec_new_bgluc <- build_spec_coveff_one_variable(
"BGLUC", seq(4, 8, by = 0.1),
qi_width_cov = 0.8, show_ref_value = FALSE
)
spec_coveff_new <- replace_spec_coveff(spec_coveff, spec_new_bgluc)
plot_coveff(ermod_bin, spec_coveff = spec_coveff_new)
Run the code above in your browser using DataLab