plot_ceac: Plot cost effectivess efficiency frontier

In sheejamk/MarkovModel: Decision Analysis Markov Model

Description

Plot cost effectivess efficiency frontier

Usage

plot_ceac(list_markov, threshold_values, comparator, currency = "GBP")

Arguments

list_markov	markov_model objects
threshold_values	list of threshold values
comparator	the comparator
currency	currency

Value

plots

Examples

well <-  health_state("well", cost=0,utility=1)
disabled <- health_state("disabled", cost=100,utility=1)
dead <- health_state("dead", cost=0,utility=0)
tmat <- rbind(c(1, 2,3), c(NA, 4,5),c(NA,NA,6))
colnames(tmat) <- rownames(tmat) <- c("well","disabled" ,"dead")
tm <- transition_matrix(3, tmat, c(0.6,0.2,0.2,0.6,0.4,1), colnames(tmat) )
health_states <- combine_state(well,disabled,dead)
this.strategy <- strategy(tm, health_states, "control")
this_markov <-markov_model(this.strategy, 24, c(1000, 0,0),c(0,0,0))
well <-  health_state("well", cost=0,utility=1)
disabled <- health_state("disabled", cost=10,utility=0.5)
dead <- health_state("dead", cost=0,utility=0)
tmat <- rbind(c(1, 2,3), c(NA, 4,5),c(NA,NA,6))
colnames(tmat) <- rownames(tmat) <- c("well","disabled" ,"dead")
tm <- transition_matrix(3, tmat, c(0.4,0.4,0.2,0.6,0.4,1), colnames(tmat))
health_states <- combine_state(well,disabled,dead)
this.strategy <- strategy(tm, health_states, "intervention")
sec_markov <-markov_model(this.strategy, 24, c(1000, 0,0),c(0,0,0))
list_markov <- combine_markov(this_markov, sec_markov)
plot_ceac(list_markov,c(1000,2000,5000,7000,10000,150000,20000), comparator = "control")

sheejamk/MarkovModel documentation built on Jan. 23, 2020, 2:44 a.m.