param <- define_parameters(
a = markov_cycle + 1 * 2
)
mat <- define_matrix(
1-1/a, 1/a,
.1, .9
)
mod <- define_model(
transition_matrix = mat,
A = define_state(cost = 10),
B = define_state(cost = 2)
)
heemod:::eval_model(
mod,
param,
init = c(10, 5),
cycles = 5,
method = "end"
)
Run the code above in your browser using DataLab