# Calculate NB's at a willingness-to-pay threshold of 80000 per unit of effects
data("df_pa")
df_pa$t_costs_d_comp2 <- df_pa$t_costs_d_comp * 1.09
df_pa$t_qaly_d_comp2 <- df_pa$t_qaly_d_comp * 1.01
calculate_nb_mult(df = df_pa,
outcomes = c("t_qaly_d_comp2", "t_qaly_d_int", "t_qaly_d_comp"),
costs = c("t_costs_d_int", "t_costs_d_comp2", "t_costs_d_comp"),
wtp = 50000
)
Run the code above in your browser using DataLab