ceac {dampack}R Documentation

Cost-Effectiveness Acceptability Curve (CEAC)


ceac is used to compute and plot the cost-effectiveness acceptability curves (CEAC) from a probabilistic sensitivity analysis (PSA) dataset.


ceac(wtp, psa)



numeric vector with willingness-to-pay (WTP) thresholds


psa object from make_psa_obj


ceac computes the probability of each of the strategies being cost-effective at each wtp threshold. The returned object has classes ceac and data.frame, and has its own plot method (plot.ceac).


An object of class ceac that can be visualized with plot. The ceac object is a data.frame that shows the proportion of PSA samples for which each strategy at each WTP threshold is cost-effective. The final column indicates whether or not the strategy at a particular WTP is on the cost-efficient frontier.

See Also

plot.ceac, summary.ceac


# psa input provided with package
example_psa_obj <- make_psa_obj(example_psa$cost, example_psa$effectiveness,
                    example_psa$parameters, example_psa$strategies)

# define wtp threshold vector (can also use a single wtp)
wtp <- seq(1e4, 1e5, by = 1e4)
ceac_obj <- ceac(wtp, example_psa_obj)
plot(ceac_obj) # see ?plot.ceac for options

# this is most useful when there are many strategies
# warnings are printed to describe strategies that
# have been filtered out
plot(ceac_obj, min_prob = 0.5)

# standard ggplot layers can be used
plot(ceac_obj) +
    labs(title = "CEAC", y = "Pr(Cost-effective) at WTP")

# the ceac object is also a data frame

# summary() tells us the regions of cost-effectiveness for each strategy.
# Note that the range_max column is an open parenthesis, meaning that the
# interval over which that strategy is cost-effective goes up to but does not include
# the value in the range_max column.

[Package dampack version 1.0.1 Index]