################################################################################
# This script conducts the value of information (VOI) analysis of the #
# cost-effectiveness analysis (CEA) of a hypothetical treatment for the #
# simulated cohort of the Sick-Sicker state-transition model (STM) #
# (PSA) dataset #
# #
# Authors: #
# - Fernando Alarid-Escudero, PhD, <falarid@stanford.edu> #
# - Eline Krijkamp, MS #
# - Petros Pechlivanoglou, PhD #
# - Hawre Jalal, MD, PhD #
# - Eva A. Enns, PhD #
################################################################################
# The structure of this code is according to the DARTH framework #
# https://github.com/DARTH-git/Decision-Modeling-Framework #
################################################################################
rm(list = ls()) # to clean the workspace
#### 05c.1 Load packages, functions and data ####
#### 05c.1.1 Load packages ####
library(dampack) # decision-analytic modeling visualization tool
#### 05c.1.2 Load functions ####
#### 05c.1.3 Load PSA dataset ####
data("l_psa")
#### 05c.2 Define VOI inputs ####
### Vector with willingness-to-pay (WTP) thresholds
v_wtp <- seq(0, 200000, by = 10000)
#### 05c.3 Expected value of perfect information (EVPI) ####
evpi <- calc_evpi(wtp = v_wtp, psa = l_psa)
plot(evpi, effect_units = "QALY")
ggsave("figs/05c_evpi.png", width = 8, height = 6)
#### 05c.4 Expected value of partial perfect information (EVPPI) ####
#### 05c.5 Expected value of sample information (EVSI) ####
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.