Define different Initial Counts for two strategies when running a model
kleanthisk10 opened this issue · 3 comments
Hello, is it possible to run different initial counts for two strategies?
Let's say we have the example provided from ?run_model
mod1 <- define_strategy( transition = define_transition( .5, .5, .1, .9 ), define_state( cost = 543, ly = 1 ), define_state( cost = 432, ly = 1 ) )
and
mod2 <- define_strategy( transition = define_transition( .5, .5, .1, .9 ), define_state( cost = 789, ly = 1 ), define_state( cost = 456, ly = 1 ) )
How can we change
res2 <- run_model( mod1, mod2, cycles = 10, cost = cost, effect = ly)
and set mod1 c(80,20) and for mod2 c(60,40) initial counts?
Thank you
Hi,
We've got the argument init
to specify the initial counts for each state, but nothing to change the counts for different strategies.
Can you please provide us with an example of the use cases for which it would be helpful to have such an option ?
For example you can see this paper, https://academic.oup.com/eurheartj/article/34/3/220/550289, where
there is a two-part cost-effectiveness model for the first strategy, which includes a short-term decision tree and a long-term Markov structure that was used to estimate long-term costs and health outcomes (Figure 1). The long-term Markov model is a logical extension of the short-term decision tree.
In the first strategy, the initial counts of the markov model is based on the probabilities derived from the short-term decision tree, meaning that the number of individuals in each state at the beginning is different from the other one where all the individuals start from a specific state.
I see.
It is for now not possible to set different initial counts. If you are comfortable with R, I think that you need to change the run_model_
function and eval_strategy
so that they can accept a list of initial counts. Once you're done, please submit a pull request.