# NOT RUN {
# load psa object provided with package
data("example_psa_obj")
# define wtp threshold vector (can also use a single wtp)
wtp <- seq(1e4, 1e5, by = 1e4)
evpi <- calc_evpi(example_psa_obj, wtp)
plot(evpi) # see ?plot.evpi for options
# can use plot options (# see ?plot.evpi for details)
plot(evpi, effect_units = "QALE")
# or can use ggplot layers
plot(evpi) + ggtitle("Expected Value of Perfect Information")
# }
Run the code above in your browser using DataLab