352x Filetype PDF File size 1.44 MB Source: www.med.mcgill.ca
Markov Models in Medical Decision Making:
A Practical Guide
FRANK A. J. ROBERT MD
SONNENBERG, MD, BECK,
when a
are useful decision involves risk
Markov models problem that is continuous over
of events is and
when the when events more
time, timing important, important may happen
than once. such clinical with conventional decision trees is difficult
Representing settings
Markov
and unrealistic assumptions. models assume that a
may require simplifying patient
number
is in of a finite of health called
always one discrete states, Markov states. All events
are as transitions from one state to another. A Markov
model be
represented may evaluated
a
matrix as a cohort or as Monte Carlo simulation. A newer
by algebra, simulation, repre-
sentation of Markov the uses a tree of clinical
models, tree, representation
Markov-cycle
either as a cohort simulation or as a
events and be evaluated Monte Carlo
may simulation.
to events and
The of the Markov model represent repetitive the time of
ability dependence
both and utilities allows for more accurate of clinical that
probabilities representation settings
Markov
involve these issues. words: models; decision decision mak-
Key Markov-cycle tree;
Decis 1993;13:322-338)
ing. (Med Making
A decision
tree models the of a sub- of survival’ or from
prognosis patient standard life tables.’ This
paper
the
to choice of a For another method for life
sequent management
strategy. explores estimating expec-
a model the Markov
example, strategy involving surgery may the tancy, model.
events of and In Beck and Pauker described the
use of
death, 1983, Mar-
surgical surgical complications,
various outcomes of the treatment itself. For kov models for in
surgical determining medical
prognosis ap-
the must be
restricted to a Since that Markov
reasons, introduction, models
practical analysis plications.’
finite time often
referred to as the
time have been
frame, horizon with in
applied increasing frequency pub-
of the This
means aside from the lished decision
analysis. that, death, analyses.’-9 software
Microcomputer
outcomes chosen to be terminal
represented nodes has been to and eval-
by developed permit
constructing
of the tree not be
final but Markov models more For these
may outcomes, may reasons,
simply uating easily.
convenient for
represent stopping points the of a revisit of the Markov model is This
scope timely. paper
the tree
Thus, contains terminal nodes serves both as a review of the behind the
analysis. every theory Mar-
that for a kov model of and as a for
represent &dquo;subsequent
prognosis&dquo; particular prognosis practical guide
combination
of characteristics
and events. the construction of Markov
patient models microcom-
using
There are various in which a
ways decision software.
analyst puter
decision-analytic
can values to these
terminal nodes of the de- Markov models
are useful
assign particularly when a de-
cision tree. In some cases the outcome
measure
is a cision involves a risk that is over
problem ongoing time.
crude life in
expectancy; others it is a Some clinical are the risk of
quality-adjusted examples hemorrhage
life One method
expectancy.’ for life ex- while on the risk of of
estimating anticoagulant
therapy, rupture
is the an
pectancy declining exponential abdominal aortic and the risk of mor-
approximation aneurysm,
of
life which
(DEALE),2 calculates a in whether sick or There
expectancy patient- tality person,
any healthy.
rate for a combination of are two of events that
specific mortality given pa- important have
consequences
tient characteristics and comorbid
diseases. Life ex- risk. the
times at
which the will
ongoing First, events
also be obtained from models occur
are uncertain.
pectancies may Gompertz This has
important
implications
because the of an
outcome
utility often on
depends
when it occurs. For a
stroke
that im-
example, occurs
have a
mediately may different on the
Received from Division of impact patient
23, 1993, the General Internal
February than one
that occurs ten later. For economic
of Robert Wood Johnson years
Medicine, UMDNJ
Department Medicine,
both
costs and
Medical New the utilities are discounted&dquo;,&dquo;
School, Brunswick, New (FAS) and Infor- analyses,
Jersey
mation of
Medicine, such that later events have
Program, Houston, less than earlier
Technology Baylor College impact
Texas in
(JRB). Grant LM05266 from the National
Supported part
by ones. The
second is that a event
of Medicine and consequence given
Grant from the for Health
Library HS06396 Agency occur more
than once. As
may the
Care and Research. following example
Policy events that
shows, are or that
Address and to Dr. representing repetitive
correspondence reprint
requests Sonnenberg: occur with
Division of General Internal UMDNJ Robert Wood John- uncertain is difficult a
Medicine, timing using simple
son Medical 97 Paterson New NJ 08903.
School, tree
Street, Brunswick, model.
322
Downloaded from http://mdm.sagepub.com at National Institutes of Health Library on January 21, 2009
323
A Specific Example
has heart
a who a valve
Consider patient prosthetic
is Such a
and anticoagulant therapy. patient
receiving
embolic or event at time.
have an hemorrhagic any
may
Either kind of event causes and/
(short-term
morbidity
in the death. The
or and result
chronic) may patient’s
in 1 shows
tree one of
decision fragment figure way
such a The
the for patient. first
representing prognosis
chance labelled has three
node, ANTICOAG, branches,
labelled BLEED, and NO EVENT. Both BLEED and
EMBOLUS,
FATAL or NON-FATAL. If NO EVENT
EMBOLUS be either
may
the remains WELL.
occurs, patient
with
There are several this model. FIGURE 1. tree of
First, Simple fragment modeling antico-
shortcomings complications
therapy.
does not when events occur. Sec- agulant
the model specify
the structure that either or
ond, implies hemorrhage
In either
embolus occur once. fact, oc- out this for even five would
may only may carrying analysis years
more than once. at the terminal nodes
cur with hundreds of
Finally, result in a tree terminal branches.
POST POST and the
labelled EMBOLUS, BLEED, WELL, analyst Thus, a recursive model is tractable for a
only very
still is faced with the of utilities, a horizon.
problem assigning short time
task to the for each of
equivalent specifying prognosis
outcomes.
non-fatal
these The Markov Model
The first problem, specifying when events occur,
be addressed the tree structure in The Markov model a far more
may by using figure provides convenient
1 and the that either BLEED or of for
making assumption way modelling prognosis clinical with
problems
EMBOLUS occurs at the time consistent with risk. The
average ongoing model assumes that the is
patient
rate of each For if
known
the in one of a finite number of
states of
complication. example, always health
rate of is a constant 0.05
the referred to as
Markov states. All events of interest
hemorrhage per person are
the time before the occurrence
then modelled as transitions from one
state to
another.
per year, average Each
or the event
of a is 1/0.05 20 state is a and the
Thus, contribution of this
hemorrhage years. assigned utility,
of a fatal will be associated with to the overall on the
having hemorrhage utility prognosis depends length
a of 20 of survival. of
time in the In
utility years normal-quality However, spent state. our of a
example patient
the normal life be less than with a heart valve, these states are WELL,
patient’s expectancy may prosthetic
of a stroke would have
20 the occurrence and
DEAD. For the sake
Thus, DISABLED, of in this
years. simplicity
the effect of the life we assume that either a
bleed or a
paradoxical improving patient’s example, non-fatal
Other such as that embolus will result in the same
expectancy. approaches, assuming state (DISABLED) and
the stroke occurs the nor- that the is
through patient’s
halfway disability permanent.
mal life are and lessen the The time
expectancy, arbitrary may horizon of the is divided into
analysis equal
the
of increments of referred
fidelity analysis. time, to as Markov Dur-
cycles.
and the of
Both the of events each the make a
timing representation ing cycle, patient transition from
may
more than once can be ad- one
occur state
events that may to another. 3 shows a used
Figure commonly
dressed a recursive decision tree.12 In a re- of Markov called a state-
using representation
by processes,
that have
have branches
some nodes transition in
cursive which each state is
tree, ap- diagram, represented
in the tree. Each of the a circle. Arrows different
two states
previously repetition by in-
peared connecting
a convenient of time dicate
tree structure allowed transitions. Arrows from a
represents length leading state
event be considered A re- to itself indicate the remain
and that in that
any may repeatedly. may
patient
that
tree models the state in
cursive anticoagulation problem consecutive certain transitions
cycles. Only
is in 2. are allowed. For a in the WELL
depicted figure example, person state
the nodes the terminal make a transition to the DISABLED but a
Here, representing previous may state,
nodes and No EVENT transition
POST-BLEED, are re- from DISABLED to WELL is not A
POST-EMBOLUS, allowed. per-
the chance node which son in either WELL
the state or the DISABLED
by ANTICOAG, state
placed appeared may
at the
root of the tree. Each die
occurrence of and thus make a transition
to the DEAD
previously state. How-
or EMBOLUS
BLEED a distinct time so ever, a who is in the DEAD state,
represents period, person obviously,
recursive model when events
the can occur. cannot make
a transition to other state.
represent any Therefore,
this model a
and car- arrow emanates from
the
However, DEAD
despite simple state,
relatively single leading
out the recursion for two time the back to itself. It is assumed
rying only periods, that a in a
patient given
in 2 is
tree with 17 terminal state can make
branches. a state
figure &dquo;bushy,&dquo; only single transition a
during
If each level of recursion represents one then
year, cycle.
Downloaded from http://mdm.sagepub.com at National Institutes of Health Library on January 21, 2009
324
FIGURE 2. Recursive tree mod-
of
eling complications antico-
agulant therapy.
data. For
The of the is chosen to a determined the available ex-
length cycle represent by probability
time interval. For a model that if only yearly probabilities are available, there
clinically meaningful ample,
the entire life of a and is little to a
spans history patient relatively advantage using monthly cycle length.
rare events the can be one On the
length
cycle year.
other if the time frame is shorter and
hand, models
events that occur much
more the
may frequently, cycle UTILITY
INCREMENTAL
time must be for or even
shorter, example monthly
Markov the
The time Evaluation of a
weekly. cycle also must be shorter if a rate process yields average
over of the amount
time. An is the risk of number (or
changes rapidly example cycles analogously, average
each state. Seen the
of in another
infarction time)
perioperative myocardial (MI) spent way,
following pre-
to a is credit&dquo; for the time in each
vious MI that declines stable value over six months.&dquo; patient &dquo;given spent
attribute of of
The of this in state. If the interest is duration
rapidity change risk dictates a only
monthly
time. Often the choice of a time will be then one need add the
cycle cycle survival, only together average
Downloaded from http://mdm.sagepub.com at National Institutes of Health Library on January 21, 2009
325
When cost-effectiveness a
performing analyses,
incremental be for each
separate utility may specified
state, the financial cost of in that
representing being
is
state for The model evaluated
one
cycle. separately
for cost and survival. Cost-effectiveness ratios
are cal-
culated as for a standard decision tree.10,11
MARKOV
TYPES OF PROCESSES
Markov are to
processes categorized according
whether the state-transition are constant
probabilities
over In the most
time or not. general of Markov
type
the transition over
process, probabilities
may
change
time. For the transition for the
FIGURE 3. Markov-state Each circle a Markov example, probability
diagram. represents DEAD
transition from WELL to consists of two
allowed transitions. compo-
state. Arrows indicate
nents. The first is the of
component probability dying
from unrelated causes. In general, this probability
times spent in the individual states to arrive at an
as
over time the
because, older,
survival for the changes patient gets
expected process. the of from unrelated causes will
probability dying
increase The second is the
continuously. component
n of a fatal or
Expected utility = ~ ts probability suffering hemorrhage embolus
the This be
not constant
s=l i during cycle. may or may
over time.
A of Markov in
where is the time in state s. special type process which the tran-
ts spent
sition are constant over time
of probabilities is called a
the survival is consid-
Usually, however, quality
Markov chain. If it has an its
state, behavior
ered Each state is associated with a absorbing
important. quality
over
time can be determined as an exact
solution
factor the of life in that state rel- by
representing quality
matrix as discussed below.
ative to health. The that is associated simple algebra, The DEALE
perfect utility
with in a state is can be used to derive the constant mortality rates
one referred
spending cycle particular
needed to a Markov chain. the
to as the incremental Consider the Markov implement However,
utility. pro-
of software to evaluate
cess in 3. If the incremental of availability specialized Markov
depicted figure utility
and the afforded
the DISABLED state is 0.7, then the in processes greater accuracy by age-
spending cycle
specific rates have resulted in reli-
the DISABLED mortality
state contributes 0.7 greater
quality-adjusted cycles
ance on Markov with
to the expected utility. Utility accrued for the entire processes time-variant proba-
Markov is the total number of in bilities.
process cycles spent
The net of a transition from
the making one
each each incremental probability
state, multiplied by utility
state
to another a is called a
for single tran-
that state. during cycle
sition The Markov is
probability. process completely
n defined by the probability distribution among the
X the
ts Us states and the for individual
utility = ~
Expected starting probabilities
s=l i of n
allowed transitions. For a Markov model states,
there n2 When
DEAD state has an incremen- will be transition probabilities. these
Let us assume that the with to
probabilities are constant respect time, they
WELL
tal of and that the state has an in-
utility zero,* x n as shown
can be an n matrix, in
that represented by
cremental 1.0. This means for
of
utility every cycle table 1. Probabilities representing disallowed transi-
is credited
in the WELL state the with a
spent patient tions of be zero. This called the P
of to the duration of a will, course, matrix,
quantity utility equal single matrix, forms the basis for the fundamental matrix
Markov If the spends, on 2.5
cycle. patient average, described in detail
solution of Markov chains Beck
in the WELL state and 1.25 in the DISABLED by
cycles cycles and Pauker.’
DEAD the
state before the state,
entering utility assigned
would be X + X or 3.9
(2.5 1) (1.25 0.7), quality-ad-
justed This number is the life
cycles. quality-adjusted TaMe 1 . P Matrix
of the
expectancy patient.
* medical the incremental of the
For examples, utility absorbing
must
DEAD state be because the will an infinite
zero, patient spend
time in the DEAD if the
of state and incremental were
amount utility
the net for the Markov would
non-zero, utility process be infinite.
Downloaded from http://mdm.sagepub.com at National Institutes of Health Library on January 21, 2009
no reviews yet
Please Login to review.