# Create a data.frame
montecarlo <- data.frame(indep1 = rnorm(1000, sd = 50, mean = 100),
indep2 = rnorm(1000, sd = 100, mean = 100))
montecarlo[, 'output1'] <- montecarlo[, 'indep1'] * montecarlo[, 'indep2']
montecarlo[, 'output2'] <- (montecarlo[, 'indep1'] * (montecarlo[, 'indep2']) + 10)
# Run the multi_EVPI function
results_all <- multi_EVPI(montecarlo,"output1")
plot_evpi(results_all, decision_vars = c("output1", "output2"),
new_names = c("Decision option 1", "Decision option 2"))
Run the code above in your browser using DataLab