Model structure

Define the structure of the oncology model.

model_structure(txseqs, n_states = c("four", "three"),
  dist = c("weibull", "gompertz", "fracpoly1", "fracpoly2"))

Arguments

txseqs	A txseq_list object.
n_states	Number of modeled health states.
dist	Parametric distribution used to model health state transitions. Options are "weibull" (Weibull), "gompertz" (gompertz), "fracpoly1"(2nd order fractional polynomial with \(p_1 = 0\) and \(p_2 = 0\)), and "fracpoly2"(2nd order fractional polynomial with \(p_1 = 0\) and \(p_2 = 1\)).

Value

A list containing the elements txseqs, n_states and dist.

Examples

txseq1 <- txseq(first = "erlotinib",
                second = c("osimertinib", "PBDC"),
               second_plus = c("PBDC + bevacizumab", "PBDC + bevacizumab"))
txseq2 <- txseq(first = "gefitinib",
               second = c("osimertinib", "PBDC"),
               second_plus = c("PBDC + bevacizumab", "PBDC + bevacizumab"))
txseqs <- txseq_list(seq1 = txseq1, seq2 = txseq2)

# Model with 4 health states
struct <- model_structure(txseqs)
names(struct)
#> [1] "txseqs"   "n_states" "dist"    
class(struct$txseqs)
#> [1] "txseq_list"
struct$n_states
#> [1] "four"
struct$dist
#> [1] "weibull"

# Model with 3 health states
struct <- model_structure(txseqs, n_states = "three", dist = "weibull")
struct$n_states
#> [1] "three"