scripts/define-decision-tree.R

In n8thangreen/VANTDET: Cost-Effectiveness Decision Tree

#
# VANTDET
# define-decision-tree.R
# N Green
#
# total pathway costs and health for each end state
# for each scenario


#################
# rule-out test #
#################

terminal_cost_ruleout <- function(cost) {

  c(cost$visit1 + cost$std.TB + cost$newtest,
    cost$visit1 + cost$std.TB + cost$newtest + cost$visit2,
    cost$visit1 + cost$std.nonTB + cost$newtest,
    cost$visit1 + cost$newtest)
}

terminal_health_ruleout <- function(health) {

  c(health$std.TB + health$newtest,
    health$std.TB + health$newtest + health$followup,
    health$std.nonTB + health$newtest,
    health$newtest)
}


################
# rule-in test #
################

terminal_cost_rulein <- function(cost) {

  c(cost$visit1 + cost$newtest,
    cost$visit1 + cost$newtest + cost$std.TB,
    cost$visit1 + cost$newtest + cost$twomonthTx + cost$visit2,
    cost$visit1 + cost$newtest + cost$std.nonTB)
}

terminal_health_rulein <- function(health) {

  c(health$newtest,
    health$newtest + health$std.TB,
    health$newtest + health$twomonthTx,
    health$newtest + health$std.nonTB)
}


terminal_cost_dual_rulein <- function(cost) {

  c(cost$visit1 + cost$newtest[1] + cost$newtest[2],
    cost$visit1 + cost$std.TB + cost$newtest[1] + cost$newtest[2],
    cost$visit1 + cost$std.TB + cost$newtest[1],
    cost$visit1 + cost$twomonthTx + cost$newtest[1] + cost$newtest[2] + cost$visit2,
    cost$visit1 + cost$std.nonTB + cost$newtest[1] + cost$newtest[2],
    cost$visit1 + cost$std.nonTB + cost$newtest[1]) %>%
    purrr::map_dbl(1)
}

terminal_health_dual_rulein <- function(health) {

  c(health$newtest[1] + health$newtest[2],
    health$std.TB + health$newtest[1] + health$newtest[2],
    health$std.TB + health$newtest[1],
    health$twomonthTx + health$newtest[1] + health$newtest[2],
    health$std.nonTB + health$newtest[1] + health$newtest[2],
    health$std.nonTB + health$newtest[1]) %>%
    purrr::map_dbl(1)
}