DiscussWk10_605

jbrnbrg

November 1, 2017

library(expm) # for %^% operator

11.1

Which of the following matrices are transition matrices for regular Markov chains?:

a)

A Markov chain is called a regular chain if some power of the transition matrix has only positive elements.

\(P= \begin{pmatrix} .5 & .5 \\ .5 & .5 \end{pmatrix}\)

Here, we’ve only got postie elements: so this is regular.

b)

\(P= \begin{pmatrix} .5 & .5 \\ 1 & 0 \end{pmatrix}\)

Here, we’ve got only positive entries with at least one diagonal not equal to zero so this is regular, too.

c)

\(P= \begin{pmatrix} 1/3 & 0 & 2/3\\ 0 & 1 & 0 \\ 0 & 1/5 & 4/5\\ \end{pmatrix}\)

For this one, I wanted to check the powers of the matrix to see if they converged:

P <- matrix(c(1/3,0,2/3,0,1,0,0,1/5,4/5), nrow=3, ncol=3, byrow = T)
P
##           [,1] [,2]      [,3]
## [1,] 0.3333333  0.0 0.6666667
## [2,] 0.0000000  1.0 0.0000000
## [3,] 0.0000000  0.2 0.8000000
P %^% 100
##              [,1] [,2]         [,3]
## [1,] 1.940325e-48    1 2.910051e-10
## [2,] 0.000000e+00    1 0.000000e+00
## [3,] 0.000000e+00    1 2.037036e-10

As shown, multiplying by itself repeatedly results in 0’s. Therefore this is not regular.