View source: R/evi.plot.mixedAn.R
| evi.plot.mixedAn | R Documentation |
Compares the optimal scenario to the mixed case in terms of the EVPI.
## S3 method for class 'mixedAn'
evi.plot(he, y.limits = NULL, pos = c(0, 1), graph = c("base", "ggplot2"), ...)
he |
An object of class |
y.limits |
Range of the y-axis for the graph. The default value is
|
pos |
Parameter to set the position of the legend (only relevant for
multiple interventions, ie more than 2 interventions being compared).
Can be given in form
of a string |
graph |
A string used to select the graphical engine to use for
plotting. Should (partial-)match the two options |
... |
Arguments to be passed to methods, such as graphical parameters
(see |
evi |
A ggplot object containing the plot. Returned only if
|
The function produces a graph showing the difference between the ”optimal” version of the EVPI (when only the most cost-effective intervention is included in the market) and the mixed strategy one (when more than one intervention is considered in the market).
Gianluca Baio, Andrea Berardi
Baio2009BCEA
\insertRefBaio2011BCEA
\insertRefBaio2013BCEA
bcea(),
mixedAn()
# See Baio G., Dawid A.P. (2011) for a detailed description of the
# Bayesian model and economic problem
#
# Load the processed results of the MCMC simulation model
data(Vaccine)
# Runs the health economic evaluation using BCEA
m <- bcea(e=eff, c=cost, # defines the variables of
# effectiveness and cost
ref=2, # selects the 2nd row of (e,c)
# as containing the reference intervention
interventions=treats, # defines the labels to be associated
# with each intervention
Kmax=50000, # maximum value possible for the willingness
# to pay threshold; implies that k is chosen
# in a grid from the interval (0,Kmax)
plot=FALSE # inhibits graphical output
)
mixedAn(m) <- NULL # uses the results of the mixed strategy
# analysis (a "mixedAn" object)
# the vector of market shares can be defined
# externally. If NULL, then each of the T
# interventions will have 1/T market share
# produces the plots
evi.plot(m)
evi.plot(m, graph="base")
# Or with ggplot2
if (require(ggplot2)) {
evi.plot(m, graph="ggplot2")
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.