310 likes | 521 Views
Markov Chain Population Models in Medical Decision Making. Gordon Hazen Min Huang Northwestern University. Markov models (individual-level) in medical decision making. Intervention that reduces disease mortality rate. Conventional outcome measure— QALYs for an individual (or a cohort).
E N D
Markov Chain Population Models in Medical Decision Making Gordon Hazen Min Huang Northwestern University M. Huang Northwestern Univ.
Markov models (individual-level) in medical decision making Intervention that reduces disease mortality rate M. Huang Northwestern Univ.
Conventional outcome measure—QALYs for an individual (or a cohort) M. Huang Northwestern Univ.
From individual to population Motivation: To study a whole population • Equilibrium distribution of a population • 2. Equilibrium measure of effectiveness of • an intervention 1. no equilibrium 2. no births Individual-level models — M. Huang Northwestern Univ.
Augment model by allowing “births” Intervention that reduces disease mortality rate M. Huang Northwestern Univ.
Population model and its routing Population model Routing process M. Huang Northwestern Univ.
Population no longer dies out—reaches new equilibrium after intervention M. Huang Northwestern Univ.
Time-homogeneous individual-level Markov models Individual Markov model State space {0,1,2,..,J,-1}, where ‘-1’ representing ‘Death’ is an absorbing state Transition rates M. Huang Northwestern Univ.
Population models Population Markov model State space: Transition rates: — Open Jackson processes Serfozo Serfozo R. Introduction to Stochastic Networks. Springer 1999. M. Huang Northwestern Univ.
Routing processes Individual-level model State space: {0,1,2,..,J,-1}, where ‘-1’ is a source/sink node Transition rates M. Huang Northwestern Univ.
Properties is irreducible, then at equilibrium: If • are independent, ~Poisson( ) equilibrium population means • Conditional on total population size |n|, n is • multinomial equilibrium population proportions M. Huang Northwestern Univ.
Equilibrium population means is the unique collection of positive numbers that satisfy balance equations of routing process i.e. . Here Q is a submatrix of the rate matrix of the routing process, and also a submatrix of the rate matrix of the underlying individual model, corresponding to all nonabsorbing states, i.e., health states {0,1,…,J}. M. Huang Northwestern Univ.
What measures of quality are possible at the population level? Measures of health Individual QALYs : QALYs for an individual starting in state j Equilibrium population measures M. Huang Northwestern Univ.
Average Lifetime QALY ALQ Mean QALY of randomly selected individual from equilibrium population M. Huang Northwestern Univ.
Total Lifetime QALY TLQ Mean total QALYs of all individuals in equilibrium population M. Huang Northwestern Univ.
Average QALYs per Year AQ/yr One-year QALY of randomly selected individual from equilibrium population M. Huang Northwestern Univ.
Total QALYs per Year TQ/yr One-year QALY of all individuals in equilibrium population M. Huang Northwestern Univ.
Discount rate = 3% Discounted Total QALYs DTQ Mean total discounted QALYs for this and all subsequent generations of population. M. Huang Northwestern Univ.
Relationships between measures DTQ TQ/yr if the population is in equilibrium from t=0. TLQ ALQ TQ/yr AQ/yr TQ/yr AQ/yr M. Huang Northwestern Univ.
The simple illustrative example— differences among measures Intervention that reduces disease mortality rate M. Huang Northwestern Univ.
Evaluating interventions using these measures: M. Huang Northwestern Univ.
Insight • Problem: average measures do not account for population size increase due to better survival. • Caution in choosing population measures M. Huang Northwestern Univ.
Example: tamoxifen use to prevent breast cancerCol Col N.F., Orr R.K., Fortin J.M. Survival impact of tamoxifen use for breast cancer risk reduction: projections from a patient-specific Markov model, Med Decis Making 2002; 22: 386-393. M. Huang Northwestern Univ.
Non-homogeneous individual-level Markov models 1. Human background survival Background mortality rate (Gompertz) 2. The other factor: a homogeneous Markov process M. Huang Northwestern Univ.
Population models Mean density with respect to age a of the population in state j at time t: Theorem: M. Huang Northwestern Univ.
Notations: equilibrium mean density with respect to age a of the population in state j, equilibrium expected total population count in state j. Conclusions: M. Huang Northwestern Univ.
Measures of health Individual QALYs : QALYs for an individual starting from age a0 in state j Equilibrium population measures ALQ TLQ AQ/yr TQ/yr TLQ M. Huang Northwestern Univ.
Example: tamoxifen use to prevent breast cancerCol M. Huang Northwestern Univ.
Summary • Population Markov models for medical decision making. • Population measures of interventions • Age-dependency. M. Huang Northwestern Univ.