0% found this document useful (0 votes)
81 views53 pages

Markov Chains: Scott Sheffield

The document discusses Markov chains, which are sequences of random variables where the probability of the next state depends only on the current state. It provides examples, including a simple weather model with rainy and sunny states. Markov chains can be represented using transition probability matrices, where entry Pij gives the probability of moving from state i to state j. The document explores how to use the transition matrix to determine the probability distribution over states at future times.

Uploaded by

Ed Z
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
0% found this document useful (0 votes)
81 views53 pages

Markov Chains: Scott Sheffield

The document discusses Markov chains, which are sequences of random variables where the probability of the next state depends only on the current state. It provides examples, including a simple weather model with rainy and sunny states. Markov chains can be represented using transition probability matrices, where entry Pij gives the probability of moving from state i to state j. The document explores how to use the transition matrix to determine the probability distribution over states at future times.

Uploaded by

Ed Z
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 53

18.

600: Lecture 32
Markov Chains

Scott Sheffield

MIT
Outline

Markov chains

Examples

Ergodicity and stationarity


Outline

Markov chains

Examples

Ergodicity and stationarity


Markov chains

I Consider a sequence of random variables X0 , X1 , X2 , . . . each


taking values in the same state space, which for now we take
to be a finite set that we label by {0, 1, . . . , M}.
Markov chains

I Consider a sequence of random variables X0 , X1 , X2 , . . . each


taking values in the same state space, which for now we take
to be a finite set that we label by {0, 1, . . . , M}.
I Interpret Xn as state of the system at time n.
Markov chains

I Consider a sequence of random variables X0 , X1 , X2 , . . . each


taking values in the same state space, which for now we take
to be a finite set that we label by {0, 1, . . . , M}.
I Interpret Xn as state of the system at time n.
I Sequence is called a Markov chain if we have a fixed
collection of numbers Pij (one for each pair
i, j {0, 1, . . . , M}) such that whenever the system is in state
i, there is probability Pij that system will next be in state j.
Markov chains

I Consider a sequence of random variables X0 , X1 , X2 , . . . each


taking values in the same state space, which for now we take
to be a finite set that we label by {0, 1, . . . , M}.
I Interpret Xn as state of the system at time n.
I Sequence is called a Markov chain if we have a fixed
collection of numbers Pij (one for each pair
i, j {0, 1, . . . , M}) such that whenever the system is in state
i, there is probability Pij that system will next be in state j.
I Precisely,
P{Xn+1 = j|Xn = i, Xn1 = in1 , . . . , X1 = i1 , X0 = i0 } = Pij .
Markov chains

I Consider a sequence of random variables X0 , X1 , X2 , . . . each


taking values in the same state space, which for now we take
to be a finite set that we label by {0, 1, . . . , M}.
I Interpret Xn as state of the system at time n.
I Sequence is called a Markov chain if we have a fixed
collection of numbers Pij (one for each pair
i, j {0, 1, . . . , M}) such that whenever the system is in state
i, there is probability Pij that system will next be in state j.
I Precisely,
P{Xn+1 = j|Xn = i, Xn1 = in1 , . . . , X1 = i1 , X0 = i0 } = Pij .
I Kind of an almost memoryless property. Probability
distribution for next state depends only on the current state
(and not on the rest of the state history).
Simple example

I For example, imagine a simple weather model with two states:


rainy and sunny.
Simple example

I For example, imagine a simple weather model with two states:


rainy and sunny.
I If its rainy one day, theres a .5 chance it will be rainy the
next day, a .5 chance it will be sunny.
Simple example

I For example, imagine a simple weather model with two states:


rainy and sunny.
I If its rainy one day, theres a .5 chance it will be rainy the
next day, a .5 chance it will be sunny.
I If its sunny one day, theres a .8 chance it will be sunny the
next day, a .2 chance it will be rainy.
Simple example

I For example, imagine a simple weather model with two states:


rainy and sunny.
I If its rainy one day, theres a .5 chance it will be rainy the
next day, a .5 chance it will be sunny.
I If its sunny one day, theres a .8 chance it will be sunny the
next day, a .2 chance it will be rainy.
I In this climate, sun tends to last longer than rain.
Simple example

I For example, imagine a simple weather model with two states:


rainy and sunny.
I If its rainy one day, theres a .5 chance it will be rainy the
next day, a .5 chance it will be sunny.
I If its sunny one day, theres a .8 chance it will be sunny the
next day, a .2 chance it will be rainy.
I In this climate, sun tends to last longer than rain.
I Given that it is rainy today, how many days to I expect to
have to wait to see a sunny day?
Simple example

I For example, imagine a simple weather model with two states:


rainy and sunny.
I If its rainy one day, theres a .5 chance it will be rainy the
next day, a .5 chance it will be sunny.
I If its sunny one day, theres a .8 chance it will be sunny the
next day, a .2 chance it will be rainy.
I In this climate, sun tends to last longer than rain.
I Given that it is rainy today, how many days to I expect to
have to wait to see a sunny day?
I Given that it is sunny today, how many days to I expect to
have to wait to see a rainy day?
Simple example

I For example, imagine a simple weather model with two states:


rainy and sunny.
I If its rainy one day, theres a .5 chance it will be rainy the
next day, a .5 chance it will be sunny.
I If its sunny one day, theres a .8 chance it will be sunny the
next day, a .2 chance it will be rainy.
I In this climate, sun tends to last longer than rain.
I Given that it is rainy today, how many days to I expect to
have to wait to see a sunny day?
I Given that it is sunny today, how many days to I expect to
have to wait to see a rainy day?
I Over the long haul, what fraction of days are sunny?
Matrix representation

I To describe a Markov chain, we need to define Pij for any


i, j {0, 1, . . . , M}.
Matrix representation

I To describe a Markov chain, we need to define Pij for any


i, j {0, 1, . . . , M}.
I It is convenient to represent the collection of transition
probabilities Pij as a matrix:

P00 P01 . . . P0M

P10 P11 . . . P1M


A=




PM0 PM1 . . . PMM
Matrix representation

I To describe a Markov chain, we need to define Pij for any


i, j {0, 1, . . . , M}.
I It is convenient to represent the collection of transition
probabilities Pij as a matrix:

P00 P01 . . . P0M

P10 P11 . . . P1M


A=




PM0 PM1 . . . PMM

I For this to make sense, we require Pij 0 for all i, j and


PM
j=0 Pij = 1 for each i. That is, the rows sum to one.
Transitions via matrices
I Suppose that pi is the probability that system is in state i at
time zero.
Transitions via matrices
I Suppose that pi is the probability that system is in state i at
time zero.
I What does the following product represent?

P00 P01 . . . P0M
P10 P11 . . . P1M


p0 p1 . . . pM




PM0 PM1 . . . PMM
Transitions via matrices
I Suppose that pi is the probability that system is in state i at
time zero.
I What does the following product represent?

P00 P01 . . . P0M
P10 P11 . . . P1M


p0 p1 . . . pM




PM0 PM1 . . . PMM

I Answer: the probability distribution at time one.


Transitions via matrices
I Suppose that pi is the probability that system is in state i at
time zero.
I What does the following product represent?

P00 P01 . . . P0M
P10 P11 . . . P1M


p0 p1 . . . pM




PM0 PM1 . . . PMM

I Answer: the probability distribution at time one.


I How about the following product?

An

p0 p1 . . . pM
Transitions via matrices
I Suppose that pi is the probability that system is in state i at
time zero.
I What does the following product represent?

P00 P01 . . . P0M
P10 P11 . . . P1M


p0 p1 . . . pM




PM0 PM1 . . . PMM

I Answer: the probability distribution at time one.


I How about the following product?

An

p0 p1 . . . pM

I Answer: the probability distribution at time n.


Powers of transition matrix

(n)
I We write Pij for the probability to go from state i to state j
over n steps.
Powers of transition matrix

(n)
I We write Pij for the probability to go from state i to state j
over n steps.
I From the matrix point of view
(n) (n) (n)
n
P00 P01 . . . P0M P00 P01 . . . P0M
(n) (n) (n)
P10 P11 . . . P1M 10 P11 . . . P1M
P



=






(n) (n) (n) PM0 PM1 . . . PMM
PM0 P M1 ... P MM
Powers of transition matrix

(n)
I We write Pij for the probability to go from state i to state j
over n steps.
I From the matrix point of view
(n) (n) (n)
n
P00 P01 . . . P0M P00 P01 . . . P0M
(n) (n) (n)
P10 P11 . . . P1M 10 P11 . . . P1M
P



=






(n) (n) (n) PM0 PM1 . . . PMM
PM0 P M1 ... P MM

I If A is the one-step transition matrix, then An is the n-step


transition matrix.
Questions

I What does it mean if all of the rows are identical?


Questions

I What does it mean if all of the rows are identical?


I Answer: state sequence Xi consists of i.i.d. random variables.
Questions

I What does it mean if all of the rows are identical?


I Answer: state sequence Xi consists of i.i.d. random variables.
I What if matrix is the identity?
Questions

I What does it mean if all of the rows are identical?


I Answer: state sequence Xi consists of i.i.d. random variables.
I What if matrix is the identity?
I Answer: states never change.
Questions

I What does it mean if all of the rows are identical?


I Answer: state sequence Xi consists of i.i.d. random variables.
I What if matrix is the identity?
I Answer: states never change.
I What if each Pij is either one or zero?
Questions

I What does it mean if all of the rows are identical?


I Answer: state sequence Xi consists of i.i.d. random variables.
I What if matrix is the identity?
I Answer: states never change.
I What if each Pij is either one or zero?
I Answer: state evolution is deterministic.
Outline

Markov chains

Examples

Ergodicity and stationarity


Outline

Markov chains

Examples

Ergodicity and stationarity


Simple example

I Consider the simple weather example: If its rainy one day,


theres a .5 chance it will be rainy the next day, a .5 chance it
will be sunny. If its sunny one day, theres a .8 chance it will
be sunny the next day, a .2 chance it will be rainy.
Simple example

I Consider the simple weather example: If its rainy one day,


theres a .5 chance it will be rainy the next day, a .5 chance it
will be sunny. If its sunny one day, theres a .8 chance it will
be sunny the next day, a .2 chance it will be rainy.
I Let rainy be state zero, sunny state one, and write the
transition matrix by
 
.5 .5
A=
.2 .8
Simple example

I Consider the simple weather example: If its rainy one day,


theres a .5 chance it will be rainy the next day, a .5 chance it
will be sunny. If its sunny one day, theres a .8 chance it will
be sunny the next day, a .2 chance it will be rainy.
I Let rainy be state zero, sunny state one, and write the
transition matrix by
 
.5 .5
A=
.2 .8

I Note that  
2 .64 .35
A =
.26 .74
Simple example

I Consider the simple weather example: If its rainy one day,


theres a .5 chance it will be rainy the next day, a .5 chance it
will be sunny. If its sunny one day, theres a .8 chance it will
be sunny the next day, a .2 chance it will be rainy.
I Let rainy be state zero, sunny state one, and write the
transition matrix by
 
.5 .5
A=
.2 .8

I Note that  
2 .64 .35
A =
.26 .74
 
.285719 .714281
I Can compute A10 =
.285713 .714287
Does relationship status have the Markov property?

In a relationship

Single Its complicated

Married Engaged
Does relationship status have the Markov property?

In a relationship

Single Its complicated

Married Engaged

I Can we assign a probability to each arrow?


Does relationship status have the Markov property?

In a relationship

Single Its complicated

Married Engaged

I Can we assign a probability to each arrow?


I Markov model implies time spent in any state (e.g., a
marriage) before leaving is a geometric random variable.
Does relationship status have the Markov property?

In a relationship

Single Its complicated

Married Engaged

I Can we assign a probability to each arrow?


I Markov model implies time spent in any state (e.g., a
marriage) before leaving is a geometric random variable.
I Not true... Can we make a better model with more states?
Outline

Markov chains

Examples

Ergodicity and stationarity


Outline

Markov chains

Examples

Ergodicity and stationarity


Ergodic Markov chains
I Say Markov chain is ergodic if some power of the transition
matrix has all non-zero entries.
Ergodic Markov chains
I Say Markov chain is ergodic if some power of the transition
matrix has all non-zero entries.
I Turns out that if chain has this property, then
(n)
j := limn Pij exists and the j are the unique
non-negative solutions of j = M
P
k=0 k Pkj that sum to one.
Ergodic Markov chains
I Say Markov chain is ergodic if some power of the transition
matrix has all non-zero entries.
I Turns out that if chain has this property, then
(n)
j := limn Pij exists and the j are the unique
non-negative solutions of j = M
P
k=0 k Pkj that sum to one.
I This means that the row vector

= 0 1 . . . M

is a left eigenvector of A with eigenvalue 1, i.e., A = .


Ergodic Markov chains
I Say Markov chain is ergodic if some power of the transition
matrix has all non-zero entries.
I Turns out that if chain has this property, then
(n)
j := limn Pij exists and the j are the unique
non-negative solutions of j = M
P
k=0 k Pkj that sum to one.
I This means that the row vector

= 0 1 . . . M

is a left eigenvector of A with eigenvalue 1, i.e., A = .


I We call the stationary distribution of the Markov chain.
Ergodic Markov chains
I Say Markov chain is ergodic if some power of the transition
matrix has all non-zero entries.
I Turns out that if chain has this property, then
(n)
j := limn Pij exists and the j are the unique
non-negative solutions of j = M
P
k=0 k Pkj that sum to one.
I This means that the row vector

= 0 1 . . . M

is a left eigenvector of A with eigenvalue 1, i.e., A = .


I We call the stationary distribution of the Markov chain.
I One can
P solve the system of linear equations
j = M k=0 k Pkj to compute the values j . Equivalent to
considering A fixed and solving A = . Or solving
(A I ) = 0. This determines
P up to a multiplicative
constant, and fact that j = 1 determines the constant.
Simple example
 
.5 .5
I If A = , then we know
.2 .8
 
 .5 .5 
A = 0 1 = 0 1 = .
.2 .8
Simple example
 
.5 .5
I If A = , then we know
.2 .8
 
 .5 .5 
A = 0 1 = 0 1 = .
.2 .8

I This means that .50 + .21 = 0 and .50 + .81 = 1 and


we also know that 0 + 1 = 1. Solving these equations gives
0 = 2/7 and 1 = 5/7, so = 2/7 5/7 .
Simple example
 
.5 .5
I If A = , then we know
.2 .8
 
 .5 .5 
A = 0 1 = 0 1 = .
.2 .8

I This means that .50 + .21 = 0 and .50 + .81 = 1 and


we also know that 0 + 1 = 1. Solving these equations gives
0 = 2/7 and 1 = 5/7, so = 2/7 5/7 .
I Indeed,
 
 .5 .5 
A = 2/7 5/7 = 2/7 5/7 = .
.2 .8
Simple example
 
.5 .5
I If A = , then we know
.2 .8
 
 .5 .5 
A = 0 1 = 0 1 = .
.2 .8

I This means that .50 + .21 = 0 and .50 + .81 = 1 and


we also know that 0 + 1 = 1. Solving these equations gives
0 = 2/7 and 1 = 5/7, so = 2/7 5/7 .
I Indeed,
 
 .5 .5 
A = 2/7 5/7 = 2/7 5/7 = .
.2 .8

I Recall 
that     
10 .285719 .714281 2/7 5/7
A = =
.285713 .714287 2/7 5/7

You might also like