inst/doc/Set_bcea_parameters.R

In BCEA: Bayesian Cost Effectiveness Analysis

## ----include = FALSE----------------------------------------------------------
knitr::opts_chunk$set(
collapse = TRUE,
comment = "#>",
fig.width = 6,
fig.height = 4
)

## ----setup--------------------------------------------------------------------
library(BCEA)

## ----eval=FALSE---------------------------------------------------------------
#  bcea(eff, cost,
#       ref = 1,
#       interventions = NULL,
#       .comparison = NULL,
#       Kmax = 50000,
#       wtp = NULL,
#       plot = FALSE)

## -----------------------------------------------------------------------------
data(Vaccine)

## -----------------------------------------------------------------------------
he_ref1 <- bcea(eff, cost,
                ref = 1,           
                interventions = treats,
                Kmax = 50000)  
str(he_ref1)

ceplane.plot(he_ref1)

## -----------------------------------------------------------------------------
he_ref2 <- bcea(eff, cost,
                ref = 2,           
                interventions = treats,
                Kmax = 50000)  

str(he_ref2[c("n_comparators", "ICER", "ref", "comp")])

## -----------------------------------------------------------------------------
setReferenceGroup(he_ref1) <- 2

str(he_ref1[c("n_comparators", "ICER", "ref", "comp")])

## -----------------------------------------------------------------------------
he_Kmax1 <- bcea(eff, cost,
                 ref = 1,           
                 interventions = treats,
                 Kmax = 50000)  

str(he_Kmax1[c("n_comparators", "ICER", "ref", "comp", "Kmax")])

## -----------------------------------------------------------------------------
he_Kmax2 <- bcea(eff, cost,
                 ref = 2,           
                 interventions = treats,
                 Kmax = 2000)  

str(he_Kmax2[c("n_comparators", "ICER", "ref", "comp", "Kmax")])

## -----------------------------------------------------------------------------
setKmax(he_Kmax1) <- 2000

str(he_Kmax1[c("n_comparators", "ICER", "ref", "comp", "Kmax")])

## -----------------------------------------------------------------------------
data(Smoking)

## -----------------------------------------------------------------------------
he_comp234 <- bcea(eff, cost,
                   ref = 1,           
                   interventions = treats,
                   Kmax = 50000)  

str(he_comp234[c("n_comparators", "ICER", "ref", "comp")])

ceplane.plot(he_comp234, wtp = 2000)

## -----------------------------------------------------------------------------
he_comp2 <- bcea(eff, cost,
                 ref = 1,
                 .comparison = 2,
                 interventions = treats,
                 Kmax = 2000)  

str(he_comp2[c("n_comparators", "ICER", "ref", "comp")])

ceplane.plot(he_comp2, wtp = 2000)

## -----------------------------------------------------------------------------
setComparisons(he_comp234) <- 2

str(he_comp234[c("n_comparators", "ICER", "ref", "comp")])

ceplane.plot(he_comp234, wtp = 2000)

## -----------------------------------------------------------------------------
he_comp24 <- bcea(eff, cost,
                  ref = 1,
                  .comparison = c(2,4),
                  interventions = treats,
                  Kmax = 2000)  
str(he_comp24[c("n_comparators", "ICER", "ref", "comp")])

ceplane.plot(he_comp24, wtp = 2000)

## -----------------------------------------------------------------------------
setComparisons(he_comp234) <- c(2,4)

str(he_comp234[c("n_comparators", "ICER", "ref", "comp")])

ceplane.plot(he_comp24, wtp = 2000)

## -----------------------------------------------------------------------------
ceplane.plot(he_comp234, comparison = 2, wtp = 2000)
ceplane.plot(he_comp234, comparison = c(2,4), wtp = 2000)

## ----echo=FALSE---------------------------------------------------------------
# create output docs
# rmarkdown::render(input = "vignettes/Set_bcea_parameters.Rmd", output_format = "pdf_document", output_dir = "vignettes")
# rmarkdown::render(input = "vignettes/Set_bcea_parameters.Rmd", output_format = "html_document", output_dir = "vignettes")