# 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=e,c=c, # 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)
)
#
# Define the vector of values for the risk aversion parameter, r, eg:
r <- c(0.0000000001,0.005,0.020,0.035)
#
# Run the cost-effectiveness analysis accounting for risk aversion
cr <- CEriskav(m, # uses the results of the economic evalaution
# (a "bcea" object)
r=r, # defines the vector of values for the risk
# aversion parameter
comparison=1 # if more than 2 interventions, selects the
# pairwise comparison
)
Run the code above in your browser using DataLab