Your Lecturer

profilekimberly T
week7economic.docx

Decision Analysis with Markov Models

Introduction to State-Transition Models

Limitations of Decision Trees

Flexible, analytic approach for evaluating cost effectiveness of intervention decisions as long as effects do not occur over long time course that would require:

· Discounting of outcomes and costs

· Complicated tree structures to capture multiple health-state changes over time

10 of 51

State-Transition Model (STM)

· Developed to overcome the limitations of decision trees

· Core component: Exhaustive list of mutually exclusive "states" that a participant can occupy over time

· E.g., "healthy," "sick," or "dead"

· Each associated with a health-related quality of life (HRQOL) and potentially a cost value

· Used to evaluate quality-adjusted life years (QALYs) and cost effectiveness of different decisions

STM Case Study

Published cost-effectiveness model comparing combination antiretroviral therapy to monotherapy for the treatment of HIV

· Used to explain structure and function of STMs

· Lists four different states

· State A: CD4 cell count >200 and <500

· State B: CD4 cell count <200

· State C: AIDS

· State D: Dead

STM Cycles

· Explicitly model passage of time over a series of discrete time periods

· Have lengths that are chosen depending on disease and intervention under study

· Range: Short time periods (e.g., minutes in hospital ER) to years for chronic disease

· Key consideration: Whether a participant would be likely to experience more than one event before the end of the cycle

· HIV treatment case study: One-year cycle length

36 of 51

STM Transition Probabilities

Define likelihood that an individual will move from one state to another during each cycle:

Transition from

Transition to

State A

State B

State C

State D

State A

0.721

0.202

0.067

0.010

State B

0

0.581

0.407

0.012

State C

0

0

0.750

0.250

State D

0

0

0

1

· Mutually exclusive states: Probabilities within each row must sum to 1

· Zero probability: Impossible to transition from one state to another

· Absorbing states: No possibility of leaving state (e.g., death)

39 of 51

44 of 51

Markov Property

· Markov models: Class of STMs with defining characteristic known as Markov property (alternatively, Markovian property or Markovian assumption)

· Assume that transition probabilities between states are independent of individual's history (i.e., "memory-less")

· E.g., probability of dying from AIDS not dependent on speed of transition from other two HIV states or length of time in AIDS state

· Markov chains: Models that also assume that transition probabilities remain constant over time

51 of 51

Costs and Outcomes

· Costs are generally implemented for each cycle that an individual is in a given state.

· Can also be assigned based on transitions to states (e.g., cost for transitioning to "dead" state)

· In the HIV treatment case study, cost are assigned for each year in a given state and equal for both treatment arms.

· Only cost differential: Cost of different treatments

· Only outcome measured: Life expectancy

· Each year spent in nondead state counted equally

· Can be changed to assign different health states different HRQOL in order to estimate QALYs

Markov Cohort Models

Definition of Cohort Model

· Evaluates expected outcome for average (vs. single) person

· Runs a group of identical individuals through state-transition model (STM)

· Exact size of cohort irrelevant to average per-person results

· In HIV treatment case study: 1,000 individuals run through model, all starting in State A during Cycle 0 and running for 20 cycles

13 of 56

Estimating Expected Value

1. Calculate probability of participant or proportion of group of patients being in a given state in a given cycle.

· Performed in spreadsheet or in specialized software (e.g., TreeAge)

· Called "Markov trace"

2. Estimate expected costs and outcomes from each cycle, weighted by proportion of cohort in each state.

· When life expectancy is used: Multiply those alive in each state by one.

· When quality-adjusted life years (QALYs) are used: Multiply values by health-related quality of life (HRQOL) for each state.

14 of 56

Markov Trace for HIV Case Study

Transition from

Transition to

State A

State B

State C

State D

State A

0.721

0.202

0.067

0.010

State B

0

0.581

0.407

0.012

State C

0

0

0.750

0.250

State D

0

0

0

1

Cycle

State A

State B

State C

State D

0

1,000

0

0

0

State A

1,000 × 0.721

State A

1,000 × 0.202

State A

1,000 × 0.067

State A

1,000 × 0.010

1

721

202

67

10

24 of 56

Markov Trace: Extended

Cycle

State A

State B

State C

State D

0

1,000

0

0

0

State A

1,000 × 0.721

State A

1,000 × 0.202

State A

1,000 × 0.067

State A

1,000 × 0.010

1

721

202

67

10

State A

721 × 0.721 = 520

State A

721 × 0.202 = 146

State B

202 × 0.581 = 117

State A

721 × 0.067 = 48

State B

202 × 0.407 = 83

State C

67 × 0.750 = 50

State A

721 × 0.010 = 7

State B

202 × 0.012 = 2

State C

67 × 0.250 = 17

State D

10×1 = 10

2

520

263

181

36

40 of 56

Markov Trace: Life Years

Cycle

State A

State B

State C

State D

Life Years

0

1,000

0

0

0

1

721

202

67

10

0.990

2

520

263

181

36

0.964

3

376

258

277

89

0.911

4

271

226

338

165

0.835

5

195

186

364

255

0.745

...

18

3

4

45

948

0.052

19

2

3

36

959

0.041

20

1

2

28

968

0.032

Total

7.996

Average per-person life years in Cycle 5:

(195 + 186 + 364) ÷ 1,000

745 ÷ 1,000 = 0.745

43 of 56

Evaluating Costs: HIV Case Study

· Annual cost by state (British pounds)

· State A: £2,756

· State B: £3,052

· State C: £9,007

· Additional annual cost of treatment with monotherapy: £2,278

45 of 56

Evaluating Costs in Markov Cohort

Cycle

State A

State B

State C

State D

Costs (£)

0

1,000

0

0

0

1

721

202

67

10

5,462

2

520

263

181

36

6,060

3

376

258

277

89

6,394

4

271

226

338

165

6,381

5

195

186

364

255

6,077

...

18

3

4

45

948

548

19

2

3

36

959

431

20

1

2

28

968

337

Total

63,745

Average per-person cost in Cycle 1:

(721 × (£2,756 + 2,278) + 202 × (£3,052 + 2,278) + 67 × (£9,007 + 2,278)) ÷ 1,000 = £5,462

48 of 56

Evaluating Costs: Discounting

Cycle

State A

State B

State C

State D

Costs (£)

Discounted (6%)

0

1,000

0

0

0

1

721

202

67

10

5,462

5,153

2

520

263

181

36

6,060

5,393

3

376

258

277

89

6,394

5,368

4

271

226

338

165

6,381

5,055

5

195

186

364

255

6,077

4,541

...

18

3

4

45

948

548

192

19

2

3

36

959

431

142

20

1

2

28

968

337

105

Total

63,745

44,663

Discounting in Cycle 5:

£6,077 × 1 ÷ (1 + 0.06)5 = £4,541

56 of 56

Half-Cycle Correction

· Membership in each state is counted either at beginning or end of each cycle.

· Yet transition into states occurs throughout cycle (i.e., at midpoint).

· Biased estimates result from ignoring difference between calculation and smooth transition.

· Addition of half of cycle at beginning of analysis corrects bias (i.e., half-cycle correction).

· Can be implemented using software solutions

· Debated in terms of its efficacy

· Easy to implement and favored by reviewers

Extensions to Markov Cohort Models

Despite their distinctions, both decision trees and Markov models can be used within the same model.

8 of 35

Incorporating Time Dependency Into Markov Cohort Models

Despite "memory-less" property of Markov models, time dependency can still be modeled.

· Allow transition probabilities to change depending on time cohort in model.

· View time as age-related changes in transition probabilities if cohort is all the same age.

· Increase of background mortality risk as cohort ages and faces competing risk of death.

11 of 35

Temporary States

· Added by creating a special health state only for a single cycle

· Can incorporate short-term changes in risks related to patient history without violating the "memory-less" property

13 of 35

24 of 35

Tunnel States

Special category of temporary states

· Arranged with each state leading to next (i.e., like a tunnel)

· Creates series of substates

· Allows costs, health-related quality of life (HRQOL), and other transition probabilities to vary based on time spent in given state

· E.g., for models of cancer prognosis

· Called semi-Markov processes

25 of 35

35 of 35

Limitations of Markov Cohort Models

· Become overly complicated if many different temporary states added to create memory

· Called "state explosion"

· Estimate only expected value of average participant vs. factoring in heterogeneity

· Could create separate, parallel cohorts for each type of person but would result in same state explosion

Individual Sampling Models

Microsimulation Models

· Known as first-order Monte Carlo simulation

· Retain time cycles, transition probabilities, health states from Markov cohort models

· Follow individuals, not cohorts

· Evaluate costs and effects for simulated individuals, estimate average effects

· Allow transition probabilities, other aspects to vary based on individual characteristics, history

12 of 22

Microsimulation Models (cont.)

· Require substantial data

· Insufficient data assumptions may add to complexity, detract from meaning

· Increase computational burden of running models

· Limitation declining due to increased computer power and cloud access

20 of 22

Dynamic Models: Infectious Disease, Social Processes

· Previous models assume individual outcomes do not influence others.

· This is untrue for infectious, chronic, and social-norm-component diseases.

· New infections depend on pool of infected individuals.

· Proportions change in epidemics and should be included in the model.

· This is relevant to evaluating vaccination, screening, and transmission models.

Best Practices for Model Use

Guidance on Technical Aspects of Modeling

· International Society for Pharmacoeconomics and Outcomes Research (ISPOR) and Society for Medical Decision Making (SMDM) collaborated on 2010 Good Research Practices in Modeling Task Force

· Goal: Create consensus-based guidelines for modeling studies in healthcare

· Focus: Use of state-transition models

13 of 70

State-Transition Models: Primary Prevention

· Risk reduction

· Events occurring prior to disease onset

· Cohort based on individuals free of the disease and complications

18 of 70

State-Transition Models: Screening

· Two types

· One-time: Newborns or genetic screening

· Repeated (interval): HPV or cervical cancer

· Evaluating strategies with respect to type and sequence of testing, diagnostic modes, screening interval, age at which screening begins and ends

23 of 70

State-Transition Models: Diagnosis

· Finding optimal diagnostic strategies

· Focusing on differences in test performance, sequence of different testing strategies, cutoff scores for positivity

28 of 70

State-Transition Models: Treatments

· Evaluating different treatment options

· Characterizing disease natural history, expected prognosis in the absence of treatment, and characterization of treatment effects and cost

36 of 70

Best Practice: Model Type

· If possible, cohort simulations should be chosen due to:

· Transparency

· Efficiency

· Ease of debugging

· Ability to conduct value-of-information analyses

38 of 70

Best Practice: Model Type (cont.)

· If valid representation of decision problems would lead to an unmanageable number of health states, then individual-level models are recommended.

· Do not sacrifice validity for simplicity.

44 of 70

Best Practice:

Starting Cohort

· Defined by demographic and clinical characteristics that affect transition probabilities, state values

· Outputs from a single cohort analysis allow comparison of alternative strategies for cohort.

· If optimal strategies vary by characteristics of subgroups, they should be modeled in parallel, not same cohort.

51 of 70

Best Practice:

Defining States

Specification: Biological or theoretical understanding of disease or condition

· Identify states that reflect the disease-health process in the absence of intervention.

· Capture benefit or harm of intervention, natural disease history.

· Combine decision tree and state-transition models when short-term states precede longer-term transition process.

57 of 70

Best Practice: Data Sources

· Transition probabilities and intervention effects should come from the most representative data sources.

· Natural history transition probabilities use population-based epidemiological studies.

· Control arms of intervention studies are less generalizable.

63 of 70

Best Practice: Microsimulation

Group should be large enough to generate stable estimates.

· Compare variance across model runs with same number of individuals to expected differences between strategies.

· Between-model-run variance due to random differences in sampled microsimulation populations should be smaller than smallest expected difference between strategies.

68 of 70

2010 Good Research Practices in Modeling

· 30 different guidelines

· See ISPOR site

Analytic Framework Overview

CHOICES Model

Childhood Obesity Intervention Cost Effectiveness Study

· Goal: Identify and prioritize cost-effective strategies for prevention.

· One in three adults, one in five children affected

· Higher healthcare costs and increased morbidity

· How to choose which policy or program to implement?

· Cost-effectiveness analysis (CEA) of options

· CHOICES provides menu of options for policy makers.

16 of 48

Why Use Modeling?

· Multidisciplinary approaches to public health:

· Regulatory: Tax on sugary drinks

· School policy: Increase in physical education

· Clinical approach: Gastric banding and weight loss surgeries

· Modeling offers a unified framework for evaluation and comparison.

· Structured, transparent results sensitive to uncertainty and different assumptions

· How are interventions being compared?

25 of 48

Counterfactual Comparison

· Compare intervention scenario with what would have happened without the intervention

· Example: Increasing prevalence of obesity complicates results

· Not a comparison with current rates, a comparison with what rates would have been without an intervention

· CHOICES microsimulation (Monte Carlo simulation at the individual level)

· Counterfactual comparison done person by person

· Individual's life simulated

· Intervention applied and differences noted

27 of 48

38 of 48

44 of 48

47 of 48

CHOICES Analysis Comparison

· Compare costs and effects between scenarios.

· Measure cost effectiveness.

· Determine which interventions offer value for money.

48 of 48

Next: A more detailed look at the modeling methods and basic building blocks of the CHOICES model

Cross-sectional Population

Studying Cross-sectional Populations

· Cross-sectional population: A collection of virtual individuals as they exist at one point in time

· Should incorporate all relevant domains

· Linked data systems in Netherlands and Scandinavia, fragmented in United States

· US data available but not in one place

· Statistical matching: The "borrowing" of information from different data sets and synthesizing it to create a virtual person in multiple dimensions

8 of 38

CHOICES Data Sets

National Sample Frame

State Variation

Bias Correction

US Census 2010

ACS 2010

(5-year)

BRFSS

2011

NHANES

2005-2010

CHOICES Model

State

State

State

State

Census

Tract

Census

Tract

Census

Tract

Age

Householder Age Group

Age

Age

Age

Sex

Sex

Sex

Sex

Race

Householder Race

Race

Race

Race

Ethnicity

Householder Ethnicity

Ethnicity

Ethnicity

Ethnicity

Household

Income

Income

Group

Income

Group

Income

Group

Etc.

Height/Weight (Self-report)

Height/Weight (Self-report)

Etc.

Etc.

Height/Weight (Measured)

Height/Weight (Measured)

School Attendance

Smoking Status

Etc.

Etc.

Health Insurance

Dietary Intake

Key Variable: Source of Sampled Variable; —: Matching Variable

18 of 38

Assumption: People who are similar in some respects are likely to be similar in other respects.
On average.

23 of 38

Problem: Self-reported information can be inaccurate and lead to biased estimates of BMI.

34 of 38

37 of 38

Cross-sectional Population

· Important to have the right starting place

· Make counterfactual comparison for each intervention

· Compare cost effectiveness of each intervention to the most realistic baseline

38 of 38

Next: Taking a cross-sectional population and moving it forward over time.

Longitudinal Population

Analytic Time Frame

· CHOICES model 2015–2025

· Burn-in period to 2015

· Run model for 10 years

· Long enough to measure longer-term outcomes

· Not so long that many large assumptions are made

· Attribute: Height and weight change over time

9 of 66

12 of 66

BMI Trajectories: Synthesis

· Use a statistical matching approach to pair segments with individuals with similar height and weight segments.

· Use a minimization technique to find similar segments and stitch them together over time.

· Model out how height and weight change for a virtual individual as they age.

13 of 66

21 of 66

BMI Trajectories: Calibration

· Observed data sets 10–15 years old

· Weight trajectories need adjustment to reflect current trends

· Calibration: Adjust parameters of model to realign output

· Fit regression models to mean BMI and obesity prevalence in NHANES 1999-2012

· Calibrate weight trajectories so aggregate obesity and mean BMI population trends match regression-based models

· Reasonable population-level trends

· Heterogeneity in height and weight growth maintained

22 of 66

28 of 66

Smoking: Confounding

· Look at joint distribution of smoking and obesity at the population level.

· When matching to Behavioral Risk Factor Surveillance System (BFRSS) 2011, pull in joint distribution of smoking and obesity at the individual level.

· Simulate smoking trajectories based on National Health Interview Surveys (NHIS) 1965–2009.

30 of 66

31 of 66

35 of 66

40 of 66

Obesity-Related Healthcare Costs

· The counterfactual scenario assumes obesity costs would be lower after an intervention.

· However, non-obesity-related costs might go up as individuals live long enough to have other health-related expenditures.

· Competing risks should stay in the model to avoid overestimating effect of an intervention.

44 of 66

Dietary Intake

· National Health and Nutrition Examination Survey (NHANES) data on calories consumed from different food groups

· No longitudinal (repeated) measures

· Make assumptions on how these change over time

53 of 66

Dietary Intake: Soda Example

· Use the model to calculate percentile of consumption, ranked from low to high consumption.

· Assume behavior tracks over time.

· Put a bandwidth around the percentile, and the next year sample a new percentile within it.

· Use this percentile to look at what the implied consumption is at the new age.

· Repeat process as people in the model age.

· Maintains distribution of consumption observed in the population

· Represents heterogeneity at the individual level

· Lets behaviors track over time in a realistic way

56 of 66

Open Population

· 2020 population will be different from the 2010 population

· More people, more different kinds of people

59 of 66

64 of 66

Baseline Uncertainty

· Repeat the process 50 times to capture uncertainty in the matching process.

· Create 50 populations from cross-sectional forward, independently and together.

· Explore the range of possible scenarios.

· Get a realistic comparator with which to evaluate potential interventions.

65 of 66

Baseline Uncertainty

· Repeat the process 50 times to capture uncertainty in the matching process.

· Create 50 populations from cross-sectional forward, independently and together.

· Explore the range of possible scenarios.

· Get a realistic comparator with which to evaluate potential interventions.

· Findings more robust to changes in underlying parameters

66 of 66

Next: Defining an intervention to apply to a baseline population

Defining an Intervention

Intervention: Something introduced into a counterfactual scenario that will alter an aspect of the baseline scenario

6 of 44

Intervention Components

1. Reach

2. Cost

3. Effect

11 of 44

Intervention Reach

· Who is exposed to the intervention? Who is the target?

· Example: States or counties; children in public schools; families on the Special Supplemental Nutrition Program for Women, Infants, and Children (WIC) or the Supplemental Nutrition Assistance Program (SNAP)

· Define interventions at different levels, e.g., geographic and individual characteristics from different domains.

· Create tree of recruitment probabilities to select individuals into the intervention.

19 of 44

Intervention Cost

· Societal perspective

· Account for opportunity cost of resources

· Geographic vs. individual

· Geographic: E.g., passing policy at state or county level (cost is the same regardless of number of people in intervention)

· Individual: E.g., training a teacher, sending a mailing, buying equipment

· Start-up vs. recurrent

· Doesn't change cost effectiveness but useful for budgeting and schedules

21 of 44

Intervention Cost (cont.)

· Capital is amortized over useful life.

· Costs are standardized to 2014 USD using Consumer Price Index.

26 of 44

Intervention Effect

· Modeled through individual weight trajectories

· Literature may report effect using different metrics

· BMI

· Behavior, e.g., purchasing patterns, diet, physical activity

33 of 44

Intervention Effect: Time Period

· Weight changes nonlinearly and can take a long time.

· Adults

· One year for half of the effect

· Up to three years for full effect

· Children

· 12–18 months

34 of 44

41 of 44

Intervention Parameters

· There is uncertainty about reach, cost, and effect.

· Account for this by specifying distributions about uncertain parameters.

· Example: Information from a clinical trial using confidence intervals about effect size on weight and/or BMI

· Result: A best guess with a range of uncertainty

· Model: Sample from normal distribution using different values to explore the effect on results

· Costs: Typically skew right, sample from gamma or log-normal distribution

44 of 44

Intervention Parameters (cont.)

· Specify interventions as realistically as possible.

· Incorporate relevant sources of uncertainty.

· Inform policy makers of best-guess results of intervention implications.

Running the Model

Uncertainty

· First-order uncertainty (stochastic)

· Some probabilities are known but not what happens to individuals in the model

· Eliminated when large enough sample size is simulated (1 million)

· Second-order uncertainty (parameter or systematic)

· Uncertainty about model parameters

· Model run repeatedly with different parameters

· Millions of individuals, thousands of repetitions, sampling at random

· Computationally intensive

18 of 44

Sampling

· Decide how parameters relate to each other.

· Specify if parameters are independent or linked within intervention.

· Example: Point estimate for effect is different according to gender.

· Sampling the high end of one population and low end of another cancels out the difference.

· Link parameters to avoid underestimating variance.

· Convey appropriate uncertainty; don't imply realistically precise results.

25 of 44

Point Estimate

· Distribution of outcomes needs to be distilled into a simple message.

· Take expected value (mean) as point estimate.

· Rank each outcome and take 2.5 and 97.5 percentiles to represent uncertainty intervals.

· Interpretation: 95% of the time the outcome will fall within this interval.

· It may be more complicated to report point estimates and uncertainty for ratio measures.

· Example: ICER = Cost/QALY

31 of 44

Point Estimate:

Ratio Measures Approaches

1. Expected value of ratios calculated within each iteration

· Calculate ratio in each iteration, order them, report 95% uncertainty interval based on the distribution.

2. Take a ratio of expected value across all iterations

· If results were plotted in 2D on cost-effectiveness plane, with a line drawn to center of mass, the slope of the line would be the ICER.

38 of 44

Recap

· Overview of analytic framework

· Importance of counterfactual comparisons

· Statistical matching and synthesizing data to create cross-sectional populations

· Other data sets and assumptions about population changes over time (longitudinal)

· Running models and distilling results into digestible indicators

· Conveying uncertainty

44 of 44

Conclusion

· You should:

· Understand the CHOICES model and how microsimulation modeling is used in economic evaluation

· Appreciate the challenges and assumptions

· Understand the right questions to ask of models

· See the value of a consistent framework that allows us explore scenarios in many domains

· Be able to provide guidance as to what are likely to be good value investments in health