BCEA (version 2.3-1.1)

ceplane.plot: Cost-effectiveness plane plot

Description

Produces a scatter plot of the cost-effectiveness plane, together with the sustainability area, as a function of the selected willingness to pay threshold

Usage

ceplane.plot(he, comparison = NULL, wtp = 25000, pos = c(1, 1),
  size = NULL, graph = c("base", "ggplot2"), xlim = NULL,
  ylim = NULL, ...)

Arguments

he

A bcea object containing the results of the Bayesian modelling and the economic evaluation.

comparison

Selects the comparator, in case of more than two interventions being analysed. Default as NULL plots all the comparisons together. Any subset of the possible comparisons can be selected (e.g., comparison=c(1,3) or comparison=2).

wtp

The value of the willingness to pay parameter. Not used iff graph="base" for multiple comparisons.

pos

Parameter to set the position of the legend; for a single comparison plot, the ICER legend position. 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(1,1), that is the topright corner inside the plot area.

size

Value (in millimetres) of the size of the willingness to pay label. Used only if graph="ggplot2", otherwise is ignored with a message.

graph

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".

xlim

The range of the plot along the x-axis. If NULL (default) it is determined by the range of the simulated values for delta.e

ylim

The range of the plot along the y-axis. If NULL (default) it is determined by the range of the simulated values for delta.c

...

If graph="ggplot2" and a named theme object is supplied, it will be added to the ggplot object. Additional graphical arguments:

• label.pos=FALSE will place the willingness to pay label in a different position at the bottom of the graph - base and ggplot2 only (no label in plotly);

• point_colors: a vector of colours specifying the colour(s) associated to the cloud of points. Should be of length 1 or equal to the number of comparisons.

• point_sizes: a vector of colours specifying the size(s) of the points. Should be of length 1 or equal to the number of comparisons.

• ICER_colors: a vector of colours specifying the colour(s) of the ICER points. Should be of length 1 or equal to the number of comparisons.

• ICER_sizes: a vector of colours specifying the size(s) of the ICER points. Should be of length 1 or equal to the number of comparisons.

• area_include: logical, include or exclude the cost-effectiveness acceptability area (default is TRUE).

• area_color: a color specifying the colour of the cost-effectiveness acceptability area

Value

ceplane

If graph="ggplot2" a ggplot object, or if graph="plotly" a plotly object containing the requested plot. Nothing is returned when graph="base", the default.

The function produces a plot of the cost-effectiveness plane. Grey dots show the simulated values for the joint distribution of the effectiveness and cost differentials. The larger red dot shows the ICER and the grey area identifies the sustainability area, i.e. the part of the plan for which the simulated values are below the willingness to pay threshold. The proportion of points in the sustainability area effectively represents the CEAC for a given value of the willingness to pay. If the comparators are more than 2 and no pairwise comparison is specified, all scatterplots are graphed using different colors.

Details

In the plotly version, point_colors, ICER_colors and area_color can also be specified as rgba colours using either the toRGB{plotly::toRGB} function or a rgba colour string, e.g. 'rgba(1, 1, 1, 1)'.

References

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

Examples

# NOT RUN {
### create the bcea object m for the smoking cessation example
data(Smoking)
m <- bcea(e,c,ref=4,Kmax=500,interventions=treats)
### produce the plot
ceplane.plot(m,wtp=200,graph="base")
### select only one comparator
ceplane.plot(m,wtp=200,graph="base",comparator=3)
### or use ggplot2 instead
if(requireNamespace("ggplot2")){
ceplane.plot(m,wtp=200,pos="right",ICER_sizes=2,graph="ggplot2")
}

# }