plot.CEriskav {BCEA}R Documentation

Summary plot of the health economic analysis when risk aversion is included


Plots the EIB and the EVPI when risk aversion is included in the utility function


## S3 method for class 'CEriskav'
plot(x, pos=c(0,1), graph=c("base","ggplot2"), ...)



An object of the class CEriskav, containing the results of the economic analysis performed accounting for a risk aversion parameter (obtained as output of the function CEriskav).


Parameter to set the position of the legend. Can be given in form of a string (bottom|top)(right|left) for base graphics and bottom|top|left|right for ggplot2. It can be a two-elements vector, which specifies the relative position on the x and y axis respectively, or alternatively it can be in form of a logical variable, with FALSE indicating to use the default position and TRUE to place it on the bottom of the plot. Default value is c(0,1), that is in the topleft corner inside the plot area.


A string used to select the graphical engine to use for plotting. Should (partial-)match the two options "base" or "ggplot2". Default value is "base".


Arguments to be passed to methods, such as graphical parameters (see par).



A two-elements named list of the ggplot objects containing the requested plots. Returned only if graph="ggplot2".

The function produces two plots for the risk aversion analysis. The first one is the EIB as a function of the discrete grid approximation of the willingness parameter for each of the possible values of the risk aversion parameter, r. The second one is a similar plot for the EVPI.


Gianluca Baio, Andrea Berardi


Baio, G., Dawid, A. P. (2011). Probabilistic Sensitivity Analysis in Health Economics. Statistical Methods in Medical Research doi:10.1177/0962280211419832.

Baio G. (2012). Bayesian Methods in Health Economics. CRC/Chapman Hall, London

See Also

bcea, CEriskav


# 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
# 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)
      plot=FALSE            # inhibits graphical output
# Define the vector of values for the risk aversion parameter, r, eg:
r <- c(0.000000000001,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 

# Now produce the plots

plot(cr # uses the results of the risk aversion 
        #  analysis (a "CEriskav" object)

### Alternative options, using ggplot2

  plot="ask"          # plot option can be specified as
                      #  "" (default), "x11" or "ask"

[Package BCEA version 2.3-1.1 Index]