# 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)
)
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
# Prints a summary of the results
summary(m, # uses the results of the mixed strategy analysis
wtp=25000) # (a "mixedAn" object)
# selects the relevant willingness to pay
# (default: 25,000)
Run the code above in your browser using DataLab