3

Consider the following Markov chain (q=1p):

Markov Chain

I want to find the mean first passage time m(i,j)(i,j0), where m(i,j)denotes the expected number of steps to reach state j when the Markov chain starts from state i.

But I'm not sure where to start, as there are an infinite number of states, and couldn't deduce a finite number of equations.

2

All of the following arguments are based around the recursive formula that comes from conditioning on the first step:

m(i,j)=1+(1p)m(i1,j)+pm(i+1,j)


Suppose i>j:

Let's simplify notation by recognizing that shifting i and j equally has no effect:

m(i,j)m^(ij)

If p.5, then the expectation is infinite. (This will follow from our solution to the p<.5 case.)

If p<.5, then:

m^(k)=1+(1p)m^(k1)+pm^(k+1)
m^(k)m^(k1)=1+p(m^(k+1)m^(k1))
We can guess and verify that m^(k)=βk is linear.
βkβ(k1)=1+p(β(k+1)β(k1))
β=1+2pβ
β=112p
Note that β does not depend on k, so we guessed correctly. Our equation is satisfied for any k by:
m^(k)=k12p
In particular, we can see for k=1:

m^(1)=1+(1p)m^(0)+pm^(2)=1+2p12p=112p

Note that our expression for m^(k) has the correct boundary values: when p=0 it deterministically takes k steps, and the expression gets arbitrarily large as p goes to .5 from the left. In fact, this proves the expectation is infinite for the p.5 case since m^ is monotonically increasing in p.


Suppose i<j:

For any fixed j, there are finitely many states to consider and an equal number of equations.

m(0,j)=1+(1p)m(0,j)+pm(1,j)
pm(0,j)=1+pm(1,j)
m(0,j)=1/p+m(1,j)

Now let's compute the next one:

m(1,j)=1+(1p)m(0,j)+pm(2,j)
pm(1,j)=1+(1p)/p+pm(2,j)
m(1,j)=(1+(1p)/p)/p+m(2,j)

There is a pattern:

m(i,j)=Ai+m(i+1,j)
Where:
Ai=(1+(1p)Ai1)/p
A0=1/p
Thus:
Ai=k=0i(1p)kpk+1

Recursing forward until j1, we get:

m(j1,j)=Aj1+m(j,j)
m(j1,j)=Aj1
Now, we can recurse backwards:
m(i,j)=n=ij1An
m(i,j)=n=ij1k=0n(1p)kpk+1
This can be simplified further by noticing An is the sum of a finite geometric series. There are three cases to consider depending on the value of r=(1p)/p.


If (1p)/p=1, then p=.5 and An=2(n+1), thus:

m(i,j)=n=ij12(n+1)=j(j+1)i(i+1)


If (1p)/p<1, then p>.5 and An is a finite geometric series with r=(1p)/p, thus:

An=1p(11pp)1(1(1pp)n+1)=(2p1)1(1(1pp)n+1)

Their summation is then:

m(i,j)=n=ij1(2p1)1(1(1pp)n+1)=(2p1)1((ji)n=ij1(1pp)n+1)=(2p1)1((ji)p(2p1)1((1pp)i+1(1pp)j+1))


If (1p)/p>1, then p<.5 and An is a finite geometric series.

We can do exactly as in the previous case but by reversing the finite geometric series and using r=p/(1p).

(I'll write this up at some point, but it's not mathematically harder than the previous case, just a bit messier notationally.)

7
  • Thank you for the comment! About the conjecture that you have made on i<j: I reckon the equation m^(k)=1+(1p)m^(k1)+pm^(k+1) doesn't actually hold for k=0 (I think it would be m^(0)=1+pm^(1)), so perhaps the conjecture should be something like m^(k)=βk+α ?? Sorry if I am wrong... – Kay Jul 31 '20 at 8:10
  • The β would be as you have found, but I thought maybe α should be considered while thinking of the case k=0 separately. – Kay Jul 31 '20 at 8:11
  • I got α=p(1p)(12p)+11p – Kay Jul 31 '20 at 8:12
  • @Kay: In my solution, m^(0)=0 since we have successfully reached j from i. Now, you might want to define m(j,j) as the expected return time from j to j, but that wouldn't affect my solution for the ij cases since we still stop when we get to j. – Sherwin Lott Jul 31 '20 at 14:30
  • 1
    Exactly, and we can easily compute that case once we have a solution for ij. (I initially thought you were defining m so that m(j,j)=0, but that's clearly not what you had in mind? But yes, that doesn't affect my answers for the ij cases.) – Sherwin Lott Jul 31 '20 at 14:49

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.