BCEA (version 2.2-3)

contour2: Specialised contour plot for objects in the class "bcea"

Description

Produces a scatterplot of the cost-effectiveness plane, with a contour-plot of the bivariate density of the differentials of cost (y-axis) and effectiveness (x-axis). Also adds the sustainability area (i.e. below the selected value of the willingness-to-pay threshold).

Usage

contour2(he, wtp=25000, xlim=NULL, ylim=NULL, comparison=NULL,  graph=c("base","ggplot2"),...)

Arguments

he

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

wtp

The selected value of the willingness-to-pay. Default is 25000.

xlim

Limits on the x-axis (default=NULL, so that R will select appropriate limits).

ylim

Limits on the y-axis (default=NULL, so that R will select appropriate limits).

comparison

The comparison being plotted. Default to NULL chooses the first comparison if graph="base". If graph="ggplot2" the default value will choose all the possible comparisons. Any subset of the possible comparisons can be selected (e.g., comparison=c(1,3)).

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

...

Arguments to be passed to ceplane.plot. See the relative manual page for more details.

Value

contour

A ggplot item containing the requested plot. Returned only if graph="ggplot2".

Plots the cost-effectiveness plane with a scatterplot of all the simulated values from the (posterior) bivariate distribution of (Delta_e,Delta_c), the differentials of effectiveness and costs; superimposes a contour of the distribution and prints the value of the ICER, together with the sustainability area.

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, ceplane.plot, contour.bcea

Examples

### create the bcea object m for the smoking cessation example
data(Smoking)
m=bcea(e,c,ref=4,interventions=treats,Kmax=500)
### produce the plot
contour2(m,wtp=200,graph="base")

### or use ggplot2 to plot multiple comparisons
contour2(m,wtp=200,ICER.size=2,graph="ggplot2")