Chapter 6

Markov chains

[190/100%]
Show
LaTeX
Problem 6.1.1

Show that any sequence of independent random variables taking values in the countable set SS is a Markov chain. Under what condition is this chain homogeneous?

?
Problem 6.1.2

A die is rolled repeatedly. Which of the following are Markov chains? For those that are, supply the transition matrix.

?
(a)

The largest number XnX_{n} shown up to the nnth roll.

(b)

The number NnN_{n} of sixes in nn rolls.

(c)

At time rr, the time CrC_{r} since the most recent six.

(d)

At time rr, the time BrB_{r} until the next six.

Problem 6.1.3

Let {Sn:n≥0}\left\{ S_{n}: n \geq 0\right\} be a simple random walk with S0=0S_{0} = 0, and show that Xn=∣Sn∣X_{n} = \left|S_{n}\right| defines a Markov chain; find the transition probabilities of this chain. Let Mn=max⁡{Sk:0≤k≤n}M_{n} = \max \left\{ S_{k}: 0 \leq k \leq n\right\}, and show that Yn=Mn−SnY_{n} = M_{n}-S_{n} defines a Markov chain. What happens if S0≠0S_{0} \neq 0?

?
Problem 6.1.4

Let XX be a Markov chain and let {nr:r≥0}\left\{ n_{r}: r \geq 0\right\} be an unbounded increasing sequence of positive integers. Show that Yr=XnrY_{r} = X_{n_{r}} constitutes a (possibly non-homogeneous) Markov chain. Find the transition matrix of YY when nr=2rn_{r} = 2 r and XX is:

?
(a)

simple random walk, and

(b)

a branching process.

Problem 6.1.5

Let XX be a Markov chain on SS, and let I:Sn→{0,1}I: S^{n} \rightarrow \left\{ 0,1\right\}. Show that the distribution of Xn,Xn+1,…X_{n}, X_{n+1}, \ldots, conditional on {I(X1,…,Xn)=1}∩{Xn=i}\left\{ I\left(X_{1}, \ldots , X_{n}\right) = 1\right\} \cap \left\{ X_{n} = i\right\}, is identical to the distribution of Xn,Xn+1,…X_{n}, X_{n+1}, \ldots conditional on {Xn=i}\left\{ X_{n} = i\right\}.

?
Problem 6.1.6

Let XX be a Markov chain on SS, and let TT be a random variable taking values in {0,1,2,…}\left\{ 0,1,2, \ldots \right\} with the property that the indicator function I{T=n}I_{\left\{ T = n\right\} }, of the event that T=nT = n, is a function of the variables X1,X2,…,XnX_{1}, X_{2}, \ldots , X_{n}. Such a random variable TT is called a stopping time, and the above definition requires that it is decidable whether or not T=nT = n with a knowledge only of the past and present, X0,X1,…,XnX_{0}, X_{1}, \ldots , X_{n}, and with no further information about the future.

Show that

P(XT+m=j∣Xk=xk for 0≤k<T,XT=i)=P(XT+m=j∣XT=i) \mathbb {P}\left(X_{T+m} = j \mid X_{k} = x_{k} \text{ for } 0 \leq k < T, X_{T} = i\right) = \mathbb {P}\left(X_{T+m} = j \mid X_{T} = i\right)

for m≥0,i,j∈Sm \geq 0, i, j \in S, and all sequences (xk)\left(x_{k}\right) of states.

?
Problem 6.1.7

Let XX be a Markov chain with state space SS, and suppose that h:S→Th: S \rightarrow T is one-one. Show that Yn=h(Xn)Y_{n} = h\left(X_{n}\right) defines a Markov chain on TT. Must this be so if hh is not one-one?

?
Problem 6.1.8

Let XX and YY be Markov chains on the set Z\mathbb {Z} of integers.

?
(a)

Is the sequence Zn=Xn+YnZ_{n} = X_{n}+Y_{n} necessarily a Markov chain?

(b)

Is ZZ a Markov chain if XX and YY are independent chains? Give a proof or a counterexample.

(c)

Show that ZZ is a Markov chain if XX and YY are independent of one another and have independent increments.

Problem 6.1.9

Let XX be a Markov chain. Which of the following are Markov chains?

?
(a)

Xm+rX_{m+r} for r≥0r \geq 0.

(b)

X2mX_{2 m} for m≥0m \geq 0.

(c)

The sequence of pairs (Xn,Xn+1)\left(X_{n}, X_{n+1}\right) for n≥0n \geq 0.

Problem 6.1.10

Let XX be a Markov chain. Show that, for 1<r<n1 < r < n,

P(Xr=k∣Xi=xi for i=1,2,…,r−1,r+1…,n)=P(Xr=k∣Xr−1=xr−1,Xr+1=xr+1) \begin{aligned} & \mathbb {P}\left(X_{r} = k \mid X_{i} = x_{i} \text{ for } i = 1,2, \ldots , r-1, r+1 \ldots , n\right) \\ & = \mathbb {P}\left(X_{r} = k \mid X_{r-1} = x_{r-1}, X_{r+1} = x_{r+1}\right) \end{aligned}
?
Problem 6.1.11

Let {Xn:n≥1}\left\{ X_{n}: n \geq 1\right\} be independent identically distributed integer-valued random variables. Let Sn=∑r=1nXrS_{n} = \sum_{r = 1}^{n} X_{r}, with S0=0,Yn=Xn+Xn−1S_{0} = 0, Y_{n} = X_{n}+X_{n-1} with X0=0X_{0} = 0, and Zn=∑r=0nSrZ_{n} = \sum_{r = 0}^{n} S_{r}. Which of the following constitute Markov chains:

?
(a)

SnS_{n},

(b)

YnY_{n},

(c)

ZnZ_{n},

(d)

the sequence of pairs (Sn,Zn)\left(S_{n}, Z_{n}\right)?

Problem 6.1.12

A stochastic matrix P\mathbf{P} is called doubly stochastic if ∑ipij=1\sum_{i} p_{i j} = 1 for all jj. It is called sub-stochastic if ∑ipij≤1\sum_{i} p_{i j} \leq 1 for all jj. Show that, if P\mathbf{P} is stochastic (respectively, doubly stochastic, sub-stochastic), then Pn\mathbf{P}^{n} is stochastic (respectively, doubly stochastic, sub-stochastic) for all nn.

?
Problem 6.1.13

Let XX be a Markov chain on the finite state space SS with transition matrix P\mathbf{P}, and let C={Cj:j∈J}\mathcal{C} = \left\{ C_{j}: j \in J\right\} be a partition of SS. Let Yn=jY_{n} = j if Xn∈CjX_{n} \in C_{j}. The chain XX is called C\mathcal{C}-lumpable if YY is a Markov chain.

Show that XX is C\mathcal{C}-lumpable if and only if, for a,b∈J,P(Xn+1∈Cb∣Xn=i)a, b \in J, \mathbb {P}\left(X_{n+1} \in C_{b} \mid X_{n} = i\right) is constant for i∈Cai \in C_{a}.

?
Problem 6.2.1

Let lij(n)=P(Xn=j,Xk≠i for 1≤k<n∣X0=i)l_{i j}(n) = \mathbb {P}\left(X_{n} = j, X_{k} \neq i \text{ for } 1 \leq k < n \mid X_{0} = i\right), the probability that the chain passes from ii to jj in nn steps without revisiting ii. Writing

Lij(s)=∑n=1∞snlij(n) L_{i j}(s) = \sum _{n = 1}^{\infty } s^{n} l_{i j}(n)

show that Pij(s)=Pii(s)Lij(s)P_{i j}(s) = P_{i i}(s) L_{i j}(s) if i≠ji \neq j. Deduce that the first passage times and last exit times have the same distribution for any Markov chain for which Pii(s)=Pjj(s)P_{i i}(s) = P_{j j}(s) for all ii and jj. Give an example of such a chain.

?
Problem 6.2.2

Let XX be a Markov chain containing an absorbing state ss with which all other states ii communicate, in the sense that pis(n)>0p_{i s}(n) > 0 for some n=n(i)n = n(i). Show that all states other than ss are transient.

?
Problem 6.2.3

Show that a state ii is recurrent if and only if the mean number of visits of the chain to ii, having started at ii, is infinite. That is to say, ii is recurrent if and only if ∑npii(n)=∞\sum_{n} p_{i i}(n) = \infty.

?
Problem 6.2.4

Let Vj=∣{n≥1:Xn=j}∣V_{j} = \left|\left\{ n \geq 1: X_{n} = j\right\} \right| be the number of visits of the Markov chain XX to jj, and define ηij=Pi(Vj=∞)\eta_{i j} = \mathbb {P}_{i}\left(V_{j} = \infty \right). Show that:

?
(a)

ηii={1 if i is recurrent, 0 if i is transient, \eta_{i i} = \begin{cases} 1 & \text{ if } i \text{ is recurrent, } \\ 0 & \text{ if } i \text{ is transient, }\end{cases}

(b)

ηij={Pi(Tj<∞) if j is recurrent, 0 if j is transient, \eta_{i j} = \begin{cases} \mathbb {P}_{i}\left(T_{j} < \infty \right) & \text{ if } j \text{ is recurrent, } \\ 0 & \text{ if } j \text{ is transient, }\end{cases} where Tj=min⁡{n≥1:Xn=j}T_{j} = \min \left\{ n \geq 1: X_{n} = j\right\}.

Problem 6.2.5

The distinct pair i,ji, j of states of a Markov chain is called symmetric if

Pi(Tj<Ti)=Pj(Ti<Tj) \mathbb {P}_{i}\left(T_{j} < T_{i}\right) = \mathbb {P}_{j}\left(T_{i} < T_{j}\right)

where Ti=min⁡{n≥1:Xn=i}T_{i} = \min \left\{ n \geq 1: X_{n} = i\right\}. Show that, if X0=iX_{0} = i and i,ji, j is symmetric, the expected number of visits to jj before the chain revisits ii is 1. [Cf. the quotation following Theorem (3.10.18).]

?
Problem 6.2.6

Let XX be a Markov chain, and let TT be a geometric random variable with P(T>n)=sn\mathbb {P}\left(T > n\right) = s^{n} for n≥0n \geq 0, which is independent of XX. By considering the expected number of visits by XX to a given state before time TT, prove Theorem (6.2.3).

?
Problem 6.2.7

Let XX be an ergodic Markov chain started from aa, and assume XX is irreducible (in that, for states i,ji, j, there exists m≥0m \geq 0 such that pij(m)>0p_{i j}(m) > 0). Let a,b,ca, b, c be distinct states, and let T(a,b,¬c)T(a, b, \neg c) be the time until XX first visits bb, with no intermediate visit to cc. That is, if XX visits bb before cc, then T(a,b,¬c)T(a, b, \neg c) is that time, and if XX visits cc before bb we set T(a,b,¬c)=∞T(a, b, \neg c) = \infty. Let

G(a,b,¬c;s)=∑n=1∞snPa(T(a,b,¬c)=n) G(a, b, \neg c ; s) = \sum _{n = 1}^{\infty } s^{n} \mathbb {P}_{a}\left(T(a, b, \neg c) = n\right)

Show that

G(a,b,¬c;s)=Fab−FacFcb1−FbcFcb G(a, b, \neg c ; s) = \frac{F_{a b}-F_{a c} F_{c b}}{1-F_{b c} F_{c b}}

where, for example, Fab(s)F_{a b}(s) is the probability generating function of the first passage time TabT_{a b} from aa to bb irrespective of intermediate visits to cc. Show further that

Pa(T(a,b,¬c)<∞)=μac+μcb−μabμbc+μcb \mathbb {P}_{a}\left(T(a, b, \neg c) < \infty \right) = \frac{\mu _{a c}+\mu _{c b}-\mu _{a b}}{\mu _{b c}+\mu _{c b}}

where, for example, μab=Ea[Tab]\mu_{a b} = \mathbb {E}_{a}\left[T_{a b}\right].

?
Problem 6.3.1

Let XX be a Markov chain on {0,1,2,…}\left\{ 0,1,2, \ldots \right\} with transition matrix given by p0j=ajp_{0 j} = a_{j} for j≥0j \geq 0, pii=rp_{i i} = r and pi,i−1=1−rp_{i, i-1} = 1-r for i≥1i \geq 1. Classify the states of the chain, and find their mean recurrence times.

?
Problem 6.3.2

Determine whether or not the random walk on the integers having transition probabilities pi,i+2=p,pi,i−1=1−pp_{i, i+2} = p, p_{i, i-1} = 1-p, for all ii, is recurrent.

?
Problem 6.3.3

Classify the states of the Markov chains with transition matrices

?
(a)
[1−2p2p0p1−2pp02p1−2p] \left[\begin{smallmatrix} 1-2 p & 2 p & 0 \\ p & 1-2 p & p \\ 0 & 2 p & 1-2 p \end{smallmatrix}\right]
(b)
[0p01−p1−p0p001−p0pp01−p0] \left[\begin{smallmatrix} 0 & p & 0 & 1-p \\ 1-p & 0 & p & 0 \\ 0 & 1-p & 0 & p \\ p & 0 & 1-p & 0 \end{smallmatrix}\right]

In each case, calculate pij(n)p_{i j}(n) and the mean recurrence times of the states.

Problem 6.3.4

A particle performs a random walk on the vertices of a cube. At each step it remains where it is with probability 14\frac{1}{4}, or moves to one of its neighbouring vertices each having probability 14\frac{1}{4}. Let vv and ww be two diametrically opposite vertices. If the walk starts at vv, find:

?
(a)

the mean number of steps until its first return to vv,

(b)

the mean number of steps until its first visit to ww,

(c)

the mean number of visits to ww before its first return to vv.

Problem 6.3.5

With the notation of Exercise (6.2.4), show that

?
(a)

if i→ji \rightarrow j and ii is recurrent, then ηij=ηji=1\eta_{i j} = \eta_{j i} = 1,

(b)

ηij=1\eta_{i j} = 1 if and only if Pi(Tj<∞)=Pj(Tj<∞)=1\mathbb {P}_{i}\left(T_{j} < \infty \right) = \mathbb {P}_{j}\left(T_{j} < \infty \right) = 1.

Problem 6.3.6

Let TA=min⁡{n≥0:Xn∈A}T_{A} = \min \left\{ n \geq 0: X_{n} \in A\right\}, where XX is a Markov chain and AA is a subset of the state space SS, and let ηj=Pj(TA<∞)\eta_{j} = \mathbb {P}_{j}\left(T_{A} < \infty \right). Show that

ηj={1 if j∈A∑k∈Spjkηk if j∉A \eta _{j} = \begin{cases} 1 & \text{ if } j \in A \\ \sum _{k \in S} p_{j k} \eta _{k} & \text{ if } j \notin A\end{cases}

Show further that if x=(xj:j∈S)\mathbf{x} = \left(x_{j}: j \in S\right) is any non-negative solution of these equations then xj≥ηjx_{j} \geq \eta_{j} for all jj.

?
Problem 6.3.7

In the notation of Exercise (6.3.6), let ρj=Ej[TA]\rho_{j} = \mathbb {E}_{j}\left[T_{A}\right]. Show that

ρj={0 if j∈A1+∑k∈Spjkρk if j∉A \rho _{j} = \begin{cases} 0 & \text{ if } j \in A \\ 1+\sum _{k \in S} p_{j k} \rho _{k} & \text{ if } j \notin A\end{cases}

and that if x=(xj:j∈S)\mathbf{x} = \left(x_{j}: j \in S\right) is any non-negative solution of these equations then xj≥ρjx_{j} \geq \rho_{j} for all jj.

?
Problem 6.3.8

Let XX be an irreducible Markov chain and let AA be a subset of the state space. Let SrS_{r} and TrT_{r} be the successive times at which the chain enters AA and visits AA respectively. Are the sequences {XSr:r≥1},{XTr:r≥1}\left\{ X_{S_{r}}: r \geq 1\right\} ,\left\{ X_{T_{r}}: r \geq 1\right\} Markov chains? What can be said about the times at which the chain exits AA?

?
Problem 6.3.9
?
(a)

Show that for each pair i,ji, j of states of an irreducible aperiodic chain, there exists N=N(i,j)N = N(i, j) such that pij(r)>0p_{i j}(r) > 0 for all r≥Nr \geq N.

(b)

Show that there exists a function ff such that, if P\mathbf{P} is the transition matrix of an irreducible aperiodic Markov chain with nn states, then pij(r)>0p_{i j}(r) > 0 for all states i,ji, j, and all r≥f(n)r \geq f(n).

(c)

Show further that f(4)≥6f(4) \geq 6 and f(n)≥(n−1)(n−2)f(n) \geq (n-1)(n-2). [Hint: The postage stamp lemma asserts that, for a,ba, b coprime, the smallest nn such that all integers strictly exceeding nn have the form αa+βb\alpha a+\beta b for some integers α,β≥0\alpha , \beta \geq 0 is (a−1)(b−1)(a-1)(b-1).]

Problem 6.3.10

An urn initially contains nn green balls and n+2n+2 red balls. A ball is picked at random: if it is green then a red ball is also removed and both are discarded; if it is red then it is replaced together with an extra red and an extra green ball. This is repeated until there are no green balls in the urn. Show that the probability the process terminates is 1/(n+1)1 /(n+1).

Now reverse the rules: if the ball is green, it is replaced together with an extra green and an extra red ball; if it is red it is discarded along with a green ball. Show that the expected number of iterations until no green balls remain is ∑j=1n(2j+1)=n(n+2)\sum_{j = 1}^{n}(2 j+1) = n(n+2). [Thus, a minor perturbation of a simple symmetric random walk can be positive recurrent, whereas the original is null recurrent.]

?
Problem 6.4.1

The proof copy of a book is read by an infinite sequence of editors checking for mistakes. Each mistake is detected with probability pp at each reading; between readings the printer corrects the detected mistakes but introduces a random number of new errors (errors may be introduced even if no mistakes were detected). Assuming as much independence as usual, and that the numbers of new errors after different readings are identically distributed, find an expression for the probability generating function of the stationary distribution of the number XnX_{n} of errors after the nnth editor-printer cycle, whenever this exists. Find it explicitly when the printer introduces a Poisson-distributed number of errors at each stage.

?
Problem 6.4.2

Do the appropriate parts of Exercises (6.3.1)-(6.3.4) again, making use of the new techniques at your disposal. In particular, for Exercise (6.3.3):

?
(a)

redo part (a);

(b)

redo part (b).

Problem 6.4.3

Let XnX_{n} be the amount of water in a reservoir at noon on day nn. During the 24 hour period beginning at this time, a quantity YnY_{n} of water flows into the reservoir, and just before noon on each day exactly one unit of water is removed (if this amount can be found). The maximum capacity of the reservoir is KK, and excessive inflows are spilled and lost. Assume that the YnY_{n} are independent and identically distributed random variables and that, by rounding off to some laughably small unit of volume, all numbers in this exercise are non-negative integers. Show that (Xn)\left(X_{n}\right) is a Markov chain, and find its transition matrix and an expression for its stationary distribution in terms of the probability generating function GG of the YnY_{n}.

Find the stationary distribution when YY has probability generating function G(s)=p(1−qs)−1G(s) = p(1-q s)^{-1}.

?
Problem 6.4.4

Show by example that chains which are not irreducible may have many different stationary distributions.

?
Problem 6.4.5

Let (xi(n):i,n≥1)\left(x_{i}(n): i, n \geq 1\right) be a bounded collection of real numbers.

?
(a)

Show that there exists an increasing sequence n1,n2,…n_{1}, n_{2}, \ldots of positive integers such that lim⁡r→∞xi(nr)\lim_{r \rightarrow \infty } x_{i}\left(n_{r}\right) exists for all ii.

(b)

Use this result to prove that, for an irreducible Markov chain, if it is not the case that pij(n)→0p_{i j}(n) \rightarrow 0 as n→∞n \rightarrow \infty for all ii and jj, then there exists a sequence (nr:r≥1)(n_{r}: r \geq 1) and a vector α(≠0)\mathbf{\alpha }( \neq \mathbf{0}) such that pij(nr)→αjp_{i j}\left(n_{r}\right) \rightarrow \alpha_{j} as r→∞r \rightarrow \infty for all ii and jj.

Problem 6.4.6

Random walk on a graph. A particle performs a random walk on the vertex set of a connected graph GG, which for simplicity we assume to have neither loops nor multiple edges. At each stage it moves to a neighbour of its current position, each such neighbour being chosen with equal probability. If GG has η(<∞)\eta ( < \infty ) edges, show that the stationary distribution is given by πv=dv/(2η)\pi_{v} = d_{v} /(2 \eta ), where dvd_{v} is the degree of vertex vv.

?
Problem 6.4.7

Show that a random walk on the infinite binary tree is transient.

?
Problem 6.4.8

At each time n=0,1,2,…n = 0,1,2, \ldots a number YnY_{n} of particles enters a chamber, where {Yn:n≥0}\left\{ Y_{n}: n \geq 0\right\} are independent and Poisson distributed with parameter λ\lambda. Lifetimes of particles are independent and geometrically distributed with parameter pp. Let XnX_{n} be the number of particles in the chamber at time nn. Show that XX is a Markov chain, and find its stationary distribution.

?
Problem 6.4.9

A random sequence of convex polygons is generated by picking two edges of the current polygon at random, joining their midpoints, and picking one of the two resulting smaller polygons at random to be the next in the sequence. Let Xn+3X_{n}+3 be the number of edges of the nnth polygon thus constructed. Find E[Xn]\mathbb {E}\left[X_{n}\right] in terms of X0X_{0}, and find the stationary distribution of the Markov chain XX.

?
Problem 6.4.10

Let ss be a state of an irreducible Markov chain on the non-negative integers. Show that the chain is recurrent if there exists a solution y\mathbf{y} to the equations yi≥∑j:j≠spijyj,i≠sy_{i} \geq \sum_{j: j \neq s} p_{i j} y_{j}, i \neq s, satisfying yi→∞y_{i} \rightarrow \infty.

?
Problem 6.4.11

A particle performs a random walk on a bow tie ABCDE drawn beneath on the left, where C is the knot. From any vertex its next step is equally likely to be to any neighbouring vertex. Initially it is at A. Find the expected value of:

?
(a)

the time of first return to A,

(b)

the number of visits to D before returning to A,

(c)

the number of visits to C before returning to A,

(d)

the time of first return to A, given no prior visit by the particle to E,

(e)

the number of visits to D before returning to A, given no prior visit by the particle to E.

Problem 6.4.12

A particle starts at A and executes a symmetric random walk on the graph drawn above on the right. Find the expected number of visits to BB before it returns to AA.

?
Problem 6.4.13

A pack contains 52 cards labelled 1,2,…,521,2, \ldots , 52, and initially they are in increasing order from top to bottom. At each stage of the shuffling process, the top card is moved to one of the 52 available places determined by the other 51 cards, this place being chosen uniformly at random, independently of all previous stages. Find the mean number of stages until card 52 is first on top.

Show that, after the moment at which the card labelled 52 is inserted at random from the top, the order of the pack is uniformly distributed over the 52! possibilities.

?
Problem 6.4.14

A pack contains 52 cards labelled 1,2,…,521,2, \ldots , 52, and initially they are in increasing order from top to bottom. At each stage, a card is picked uniformly at random from the pack and placed on top, independently of all previous stages. Find the mean number of stages until every card has been selected at least once.

Show that, after the moment at which the final card to be selected at random is placed on top, the order of the pack is uniformly distributed over the 52! possibilities.

?
Problem 6.4.15

Dick and Jim are writing exercises in sequence for inclusion in a textbook. Dick writes them and Jim checks them. Each exercise is faulty with probability pp, independently of other exercises. Jim has two modes of operation. In Mode A, he inspects every exercise as it is produced. In Mode B, he inspects each exercise with probability 1/r1 / r where r>1r > 1, independently of all other events.

Let N≥1N \geq 1. Jim operates in Mode A until he has found NN consecutive non-defective exercises, at which point he changes to Mode B. He operates in Mode B until the first defective exercise is found, and then he reverts to Mode A.

Let XX be the Markov chain which is in state ii if Jim is operating in Mode A and the last ii consecutive exercises since entering Mode A have been found to be non-defective, and is in state NN if Jim is in Mode B.

?
(a)

Write down the transition probabilities of XX and find its stationary distribution.

(b)

Show that the long run proportion of exercises that are inspected is 1/[1+(r−1)(1−pN)]1 /\left[1+(r-1)\left(1-p^{N}\right)\right].

(c)

Find an expression for the long run proportion of defective exercises which are not detected by Jim.

Problem 6.4.16

A particle performs a random walk on the non-negative integers as follows. When at position k≥0k \geq 0, its next position is uniformly distributed on the set {0,1,…,k,k+1}\left\{ 0,1, \ldots , k, k+1\right\}. Show that the sequence of positions forms an aperiodic, positive recurrent Markov chain, and find its stationary distribution.

Find the mean number μ\mu of steps required to reach position 0 for the first time from position 1.

?
Problem 6.4.17

Let {Xn:n≥0}\left\{ X_{n}: n \geq 0\right\} be an irreducible Markov chain with state space SS and transition matrix P\mathbf{P} (the chain may be either transient or recurrent). Let k∈Sk \in S and let x\mathbf{x} be a stationary measure such that xk=1x_{k} = 1. Prove that x≥ρ(k)\mathbf{x} \geq \rho (k) where ρ(k)\rho (k) is given in equation (6.4.5). If the chain is recurrent, show that x=ρ(k)\mathbf{x} = \rho (k).

?
Problem 6.5.1

A random walk on the set {0,1,2,…,b}\left\{ 0,1,2, \ldots , b\right\} has transition matrix given by p00=1−λ0,pbb=1−μb,pi,i+1=λip_{00} = 1-\lambda_{0}, p_{b b} = 1-\mu_{b}, p_{i, i+1} = \lambda_{i} and pi+1,i=μi+1p_{i+1, i} = \mu_{i+1} for 0≤i<b0 \leq i < b, where 0<λi,μi<10 < \lambda_{i}, \mu_{i} < 1 for all ii, and λi+μi=1\lambda_{i}+\mu_{i} = 1 for 1≤i<b1 \leq i < b. Show that this process is reversible in equilibrium.

?
Problem 6.5.2

Let XX be an irreducible, positive recurrent, aperiodic Markov chain on the state space SS.

?
(a)

Kolmogorov's reversibility criterion. Show that XX is reversible in equilibrium if and only if

pj1,j2pj2,j3⋯pjn−1,jnpjn,j1=pj1,jnpjn,jn−1⋯pj2,j1 p_{j_{1}, j_{2}} p_{j_{2}, j_{3}} \cdots p_{j_{n-1}, j_{n}} p_{j_{n}, j_{1}} = p_{j_{1}, j_{n}} p_{j_{n}, j_{n-1}} \cdots p_{j_{2}, j_{1}}

for all nn and all finite sequences j1,j2,…,jnj_{1}, j_{2}, \ldots , j_{n} of states.

(b)

Kelly's reversibility condition. Show that XX is reversible in equilibrium if, for all distinct triples i,j,k∈Si, j, k \in S,

pijpjkpki=pikpkjpji p_{i j} p_{j k} p_{k i} = p_{i k} p_{k j} p_{j i}

and in addition there exists c∈Sc \in S such that pic>0p_{i c} > 0 for all i≠ci \neq c.

(c)

Consider a chain with n≥3n \geq 3 states. Show that Kolmogorov's criterion, as expressed above, may require the verification of up to 12∑r=3n(nr)(r−1)\frac{1}{2} \sum_{r = 3}^{n}\binom {n}{r}(r-1) ! equations, whereas Kelly's condition, if appropriate, requires no more than (n−13)\binom {n-1}{3}.

(d)

Show that a random walk on a finite tree is reversible in equilibrium.

Problem 6.5.3

Let XX be a reversible Markov chain, and let CC be a non-empty subset of the state space SS. Define the Markov chain YY on SS by the transition matrix Q=(qij)\mathbf{Q} = \left(q_{i j}\right) where

qij={βpij if i∈C and j∉C,pij otherwise,  q_{i j} = \begin{cases} \beta p_{i j} & \text{ if } i \in C \text{ and } j \notin C, \\ p_{i j} & \text{ otherwise, }\end{cases}

for i≠ji \neq j, and where β\beta is a constant satisfying 0<β<10 < \beta < 1. The diagonal terms qiiq_{i i} are arranged so that Q\mathbf{Q} is a stochastic matrix. Show that YY is reversible in equilibrium, and find its stationary distribution. Describe the situation in the limit as β↓0\beta \downarrow 0.

?
Problem 6.5.4

Can a reversible chain be periodic?

?
Problem 6.5.5

The dog-flea model of Example (6.5.5) is a Markov chain XX on the state space {0,1,…,m}\left\{ 0,1, \ldots , m\right\} with transition probabilities

pi,i+1=1−im,pi,i−1=im, for 0≤i≤m p_{i, i+1} = 1-\frac{i}{m}, \quad p_{i, i-1} = \frac{i}{m}, \quad \text{ for } \quad 0 \leq i \leq m

Show that, if X0=iX_{0} = i,

E[Xn−m2]=(i−m2)(1−2m)n→0 as n→∞ \mathbb {E}\left[X_{n}-\frac{m}{2}\right] = \left(i-\frac{m}{2}\right)\left(1-\frac{2}{m}\right)^{n} \rightarrow 0 \quad \text{ as } n \rightarrow \infty
?
Problem 6.5.6

Which of the following (when stationary) are reversible Markov chains?

?
(a)

The chain X={Xn}X = \left\{ X_{n}\right\} having transition matrix P=(1−ααβ1−β)\mathbf{P} = \left(\begin{array}{cc}1-\alpha & \alpha \\ \beta & 1-\beta \end{array}\right) where α+β>0\alpha +\beta > 0.

(b)

The chain Y={Yn}Y = \left\{ Y_{n}\right\} having transition matrix P=(0p1−p1−p0pp1−p0)\mathbf{P} = \left(\begin{array}{ccc}0 & p & 1-p \\ 1-p & 0 & p \\ p & 1-p & 0\end{array}\right) where 0<p<10 < p < 1.

(c)

Zn=(Xn,Yn)Z_{n} = \left(X_{n}, Y_{n}\right), where XnX_{n} and YnY_{n} are independent and satisfy (a) and (b).

Problem 6.5.7

Let Xn,YnX_{n}, Y_{n} be independent simple random walks. Let ZnZ_{n} be (Xn,Yn)(X_{n}, Y_{n}) truncated to lie in the region Xn≥0,Yn≥0,Xn+Yn≤aX_{n} \geq 0, Y_{n} \geq 0, X_{n}+Y_{n} \leq a where aa is integral. Find the stationary distribution of ZnZ_{n}.

?
Problem 6.5.8

Show that an irreducible Markov chain with a finite state space and transition matrix P\mathbf{P} is reversible in equilibrium if and only if P=DS\mathbf{P} = \mathbf{D S} for some symmetric matrix S\mathbf{S} and diagonal matrix D\mathbf{D} with strictly positive diagonal entries. Show further that for reversibility in equilibrium to hold, it is necessary but not sufficient that P\mathbf{P} has real eigenvalues.

?
Problem 6.5.9

Random walk on a graph. Let GG be a finite connected graph with neither loops nor multiple edges, and let XX be a random walk on GG as in Exercise (6.4.6). Show that XX is reversible in equilibrium.

?
Problem 6.5.10

Consider a random walk on the strictly positive integers with transition probabilities

pi,i−1=12⋅i+2i+1,pi,i+1=12⋅ii+1,i≥2 p_{i, i-1} = \frac{1}{2} \cdot \frac{i+2}{i+1}, \quad p_{i, i+1} = \frac{1}{2} \cdot \frac{i}{i+1}, \quad i \geq 2

and p11=34,p12=14p_{11} = \frac{3}{4}, p_{12} = \frac{1}{4}. Show that the walk is positive recurrent and find the mean recurrence time of the state ii.

?
Problem 6.5.11

An aleatory beetle performs a random walk on five vertices comprising the principal points (labelled n,e,s,wn, e, s, w) and centre (labelled cc) of a compass, with transition probabilities

pen=pws=14,pcn=pne=pec=pce=psw=pwn=18,pnc=pes=psc=pcs=pse=pwc=116,pcw=pnw=132. \begin{aligned} p_{e n} & = p_{w s} = \frac{1}{4}, \\ p_{c n} = p_{n e} = p_{e c} & = p_{c e} = p_{s w} = p_{w n} = \frac{1}{8}, \\ p_{n c} = p_{e s} = p_{s c} & = p_{c s} = p_{s e} = p_{w c} = \frac{1}{16}, \\ p_{c w} & = p_{n w} = \frac{1}{32}. \end{aligned}

Other moves have probability 0. Show that the mean recurrence time of the centre is μc=112\mu_{c} = \frac{11}{2}.

?
Problem 6.5.12

Let XX be an irreducible (but not necessarily aperiodic) Markov chain on the countable state space SS with transition matrix P\mathbf{P}, with invariant distribution π\pi. Let a∈(0,1)a \in (0,1) and let L=aP+(1−a)I\mathbf{L} = a \mathbf{P}+(1-a) \mathbf{I} where I\mathbf{I} is the identity matrix.

?
(a)

Show that L\mathbf{L} is the transition matrix of an irreducible, aperiodic Markov chain YY with invariant distribution π\pi.

(b)

Show that, if XX is reversible in equilibrium, then so is YY.

Problem 6.6.1

The Markov-Kakutani theorem asserts that, for any convex compact subset CC of Rn\mathbb {R}^{n} and any linear continuous mapping TT of CC into C,TC, T has a fixed point (in the sense that T(x)=xT(x) = x for some x∈Cx \in C). Use this to prove that a finite stochastic matrix has a non-negative non-zero left eigenvector corresponding to the eigenvalue 1.

?
Problem 6.6.2

Let T\mathbf{T} be a m×nm \times n matrix and let v∈Rn\mathbf{v} \in \mathbb {R}^{n}. Farkas's theorem asserts that exactly one of the following holds:

(i) there exists x∈Rm\mathbf{x} \in \mathbb {R}^{m} such that x≥0\mathbf{x} \geq \mathbf{0} and xT=v\mathbf{x} \mathbf{T} = \mathbf{v},

(ii) there exists y∈Rn\mathbf{y} \in \mathbb {R}^{n} such that yv′<0\mathbf{y v}^{\prime } < 0 and Ty′≥0\mathbf{T y}^{\prime } \geq \mathbf{0}.

Use this to prove that a finite stochastic matrix has a non-negative non-zero left eigenvector corresponding to the eigenvalue 1.

?
Problem 6.6.3

Suppose you are betting on a race with mm possible outcomes. There are nn bookmakers, and a unit stake with the ii th bookmaker yields tijt_{i j} if the jj th outcome of the race occurs. A vector x=(x1,x2,…,xn)\mathbf{x} = \left(x_{1}, x_{2}, \ldots , x_{n}\right), where xr∈(−∞,∞)x_{r} \in (-\infty , \infty ) is your stake with the rr th bookmaker, is called a betting scheme. Show that exactly one of (a) and (b) holds:

(a) there exists a probability mass function p=(p1,p2,…,pm)\mathbf{p} = \left(p_{1}, p_{2}, \ldots , p_{m}\right) such that ∑j=1mtijpj=0\sum_{j = 1}^{m} t_{i j} p_{j} = 0 for all values of ii,

(b) there exists a betting scheme x\mathbf{x} for which you surely win, that is, ∑i=1nxitij>0\sum_{i = 1}^{n} x_{i} t_{i j} > 0 for all jj.

?
Problem 6.6.4

Let XX be a Markov chain with state space S={1,2,3}S = \left\{ 1,2,3\right\} and transition matrix

P=[1−pp001−ppp01−p] \mathbf{P} =\left[\begin{smallmatrix} 1-p & p & 0 \\ 0 & 1-p & p \\ p & 0 & 1-p \end{smallmatrix}\right]

where 0<p<10 < p < 1. Prove that

Pn=[a1na2na3na3na1na2na2na3na1n] \mathbf{P}^{n} =\left[\begin{smallmatrix} a_{1 n} & a_{2 n} & a_{3 n} \\ a_{3 n} & a_{1 n} & a_{2 n} \\ a_{2 n} & a_{3 n} & a_{1 n} \end{smallmatrix}\right]

where a1n+ωa2n+ω2a3n=(1−p+pω)n,ωa_{1 n}+\omega a_{2 n}+\omega^{2} a_{3 n} = (1-p+p \omega )^{n}, \omega being a complex cube root of 1.

?
Problem 6.6.5

Let P\mathbf{P} be the transition matrix of a Markov chain with finite state space. Let I\mathbf{I} be the identity matrix, U\mathbf{U} the ∣S∣×∣S∣\left|S\right| \times \left|S\right| matrix with all entries unity, and 1\mathbf{1} the row ∣S∣\left|S\right|-vector with all entries unity. Let π\mathbf{\pi } be a non-negative vector with ∑iπi=1\sum_{i} \pi_{i} = 1. Show that πP=π\mathbf{\pi } \mathbf{P} = \mathbf{\pi } if and only if π(I−P+U)=1\mathbf{\pi }(\mathbf{I}-\mathbf{P}+\mathbf{U}) = \mathbf{1}. Deduce that if P\mathbf{P} is irreducible then π=1(I−P+U)−1\mathbf{\pi } = \mathbf{1}(\mathbf{I}-\mathbf{P}+\mathbf{U})^{-1}.

?
Problem 6.6.6

A chess piece performs a random walk on a chessboard; at each step it is equally likely to make any one of the available moves. What is the mean recurrence time of a corner square if the piece is a:

?
(a)

king?

(b)

queen?

(c)

bishop?

(d)

knight?

(e)

rook?

Problem 6.6.7

A rook and a bishop perform independent symmetric random walks with synchronous steps on a 4×44 \times 4 chessboard (16 squares). If they start together at a corner, show that the expected number of steps until they meet again at the same corner is 448/3448 / 3.

?
Problem 6.6.8

Find the nn-step transition probabilities pij(n)p_{i j}(n) for the chain XX having transition matrix

P=[0121213145122314112] \mathbf{P} =\left[\begin{smallmatrix} 0 & \frac{1}{2} & \frac{1}{2} \\ \frac{1}{3} & \frac{1}{4} & \frac{5}{12} \\ \frac{2}{3} & \frac{1}{4} & \frac{1}{12} \end{smallmatrix}\right]
?
Problem 6.6.9

The pages on the worldwide web form a directed graph WW with nn vertices (representing pages) joined by directed edges (representing links). The existence of a link from ii to jj is denoted by i→ji \rightarrow j, and the graph is specified by its adjacency matrix L=(lij)L = \left(l_{i j}\right) where lij=1l_{i j} = 1 if i→ji \rightarrow j and lij=0l_{i j} = 0 otherwise. The out-degree did_{i} (respectively, in-degree cic_{i}) of vertex ii is the number of links pointing away from ii (respectively, towards ii). Vertex ii is said to dangle if di=0d_{i} = 0.

The behaviour of a swiftly bored web surfer is modelled by a random walk on WW. Let b∈(0,1)b \in (0,1). From any dangling vertex, the random walk moves to a randomly chosen vertex of WW, each vertex having probability 1/n1 / n. When at a non-dangling vertex ii, with probability b<1b < 1 the walk moves to a random linked vertex (each having probability 1/di1 / d_{i}), while with probability 1−b1-b it moves to a random vertex of WW (each having probability 1/n1 / n).

?
(a)

Show that the transition matrix P\mathbf{P} may be written in the form P=bQ+v′e\mathbf{P} = b \mathbf{Q}+\mathbf{v}^{\prime } \mathbf{e}, where Q=(qij)\mathbf{Q} = \left(q_{i j}\right) with

qij={1/di if i→j1/n if i dangles 0 otherwise  q_{i j} = \begin{cases} 1 / d_{i} & \text{ if } i \rightarrow j \\ 1 / n & \text{ if } i \text{ dangles } \\ 0 & \text{ otherwise }\end{cases}

and v=(vi)\mathbf{v} = \left(v_{i}\right) is a row vector with vi=(1−b)/nv_{i} = (1-b) / n, and e\mathbf{e} is a row vector with entries 1.

(b)

Deduce that the stationary distribution π\mathbf{\pi } is given by π={(1−b)/n}e(I−bQ)−1\mathbf{\pi } = \left\{ (1-b) / n\right\} \mathbf{e}(\mathbf{I}-b \mathbf{Q})^{-1} where I\mathbf{I} is the identity matrix.

(c)

Explain why the elements of π\mathbf{\pi }, when rearranged in decreasing order, supply a description of the relative popularities of web pages (called 'PageRank' by Google, that being their trademark for the patented algorithm).

Problem 6.6.10

Let P\mathbf{P} be the transition matrix of an irreducible Markov chain on a finite state space, and let π\mathbf{\pi } be a left eigenvector of P\mathbf{P} corresponding to the eigenvalue 1. Show from the equation π=πP\mathbf{\pi } = \mathbf{\pi } \mathbf{P} directly that the entries of π\pi are either all positive or all negative, and hence prove Theorem (6.6.1d): there exists a unique distribution π\mathbf{\pi } satisfying πP=π\mathbf{\pi } \mathbf{P} = \mathbf{\pi }, and, furthermore, all components of π\mathbf{\pi } are strictly positive.

?
Problem 6.7.1

Let ZnZ_{n} be the size of the nnth generation of a branching process with Z0=1Z_{0} = 1 and P(Z1=k)=2−k\mathbb {P}\left(Z_{1} = k\right) = 2^{-k} for k≥0k \geq 0. Show directly that, as n→∞,P(Zn≤2yn∣Zn>0)→1−e−2y,y>0n \rightarrow \infty , \mathbb {P}\left(Z_{n} \leq 2 y n \mid Z_{n} > 0\right) \rightarrow 1-e^{-2 y}, y > 0, in agreement with Theorem (6.7.8).

?
Problem 6.7.2

Let ZZ be a supercritical branching process with Z0=1Z_{0} = 1 and family-size generating function GG. Assume that the probability η\eta of extinction satisfies 0<η<10 < \eta < 1. Find a way of describing the process ZZ, conditioned on its ultimate extinction.

?
Problem 6.7.3

Let ZnZ_{n} be the size of the nnth generation of a branching process with Z0=1Z_{0} = 1 and P(Z1=k)=qpk\mathbb {P}\left(Z_{1} = k\right) = q p^{k} for k≥0k \geq 0, where p+q=1p+q = 1 and p>12p > \frac{1}{2}. Use your answer to Exercise (6.7.2) to show that, if we condition on the ultimate extinction of ZZ, then the process grows in the manner of a branching process with generation sizes Z~n\widetilde{Z}_{n} satisfying Z~0=1\widetilde{Z}_{0} = 1 and P(Z~1=k)=pqk\mathbb {P}\left(\widetilde{Z}_{1} = k\right) = p q^{k} for k≥0k \geq 0.

?
Problem 6.7.4
?
(a)

Show that E[X∣X>0]≤E[X2]/E[X]\mathbb {E}\left[X \mid X > 0\right] \leq \mathbb {E}\left[X^{2}\right] / \mathbb {E}\left[X\right] for any random variable XX taking non-negative values.

(b)

Let ZnZ_{n} be the size of the nnth generation of a branching process with Z0=1Z_{0} = 1 and P(Z1=k)=qpk\mathbb {P}\left(Z_{1} = k\right) = q p^{k} for k≥0k \geq 0, where p>12p > \frac{1}{2}. Use part (a) to show that E[Zn/μn∣Zn>0]≤2p/(p−q)\mathbb {E}\left[Z_{n} / \mu^{n} \mid Z_{n} > 0\right] \leq 2 p /(p-q), where μ=p/q\mu = p / q.

(c)

Show that, in the notation of part (b), E[Zn/μn∣Zn>0]→p/(p−q)\mathbb {E}\left[Z_{n} / \mu^{n} \mid Z_{n} > 0\right] \rightarrow p /(p-q) as n→∞n \rightarrow \infty.

Problem 6.8.1

Flies and wasps land on your dinner plate in the manner of independent Poisson processes with respective intensities λ\lambda and μ\mu. Show that the arrivals of flying objects form a Poisson process with intensity λ+μ\lambda +\mu.

?
Problem 6.8.2

Insects land in the soup in the manner of a Poisson process with intensity λ\lambda, and each such insect is green with probability pp, independently of the colours of all other insects. Show that the arrivals of green insects form a Poisson process with intensity λp\lambda p.

?
Problem 6.8.3

Let TnT_{n} be the time of the nnth arrival in a Poisson process NN with intensity λ\lambda, and define the excess lifetime process E(t)=TN(t)+1−tE(t) = T_{N(t)+1}-t, being the time one must wait subsequent to tt before the next arrival. Show by conditioning on T1T_{1} that

P(E(t)>x)=e−λ(t+x)+∫0tP(E(t−u)>x)λe−λudu \mathbb {P}\left(E(t) > x\right) = e^{-\lambda (t+x)}+\int _{0}^{t} \mathbb {P}\left(E(t-u) > x\right) \lambda e^{-\lambda u} d u

Solve this integral equation in order to find the distribution function of E(t)E(t). Explain your conclusion.

?
Problem 6.8.4

Let BB be a simple birth process of paragraph (6.8.15 b)(6.8.15 \mathrm{~ b}) with B(0)=IB(0) = I; the birth rates are λn=nλ\lambda_{n} = n \lambda. Write down the forward system of equations for the process and deduce that

P(B(t)=k)=(k−1I−1)e−Iλt(1−e−λt)k−I,k≥I \mathbb {P}\left(B(t) = k\right) = \binom {k-1}{I-1} e^{-I \lambda t}\left(1-e^{-\lambda t}\right)^{k-I}, \quad k \geq I

Show also that E[B(t)]=Ieλt\mathbb {E}\left[B(t)\right] = I e^{\lambda t} and Var⁡[B(t)]=Ie2λt(1−e−λt)\operatorname {Var}\left[B(t)\right] = I e^{2 \lambda t}\left(1-e^{-\lambda t}\right).

?
Problem 6.8.5

Let BB be a process of simple birth with immigration (6.8.11c) with parameters λ\lambda and ν\nu, and with B(0)=0B(0) = 0; the birth rates are λn=nλ+ν\lambda_{n} = n \lambda +\nu. Write down the sequence of differential-difference equations for pn(t)=P(B(t)=n)p_{n}(t) = \mathbb {P}\left(B(t) = n\right). Without solving these equations, use them to show that m(t)=E[B(t)]m(t) = \mathbb {E}\left[B(t)\right] satisfies m′(t)=λm(t)+νm^{\prime }(t) = \lambda m(t)+\nu, and solve for m(t)m(t).

?
Problem 6.8.6

Let NN be a birth process with intensities λ0,λ1,…\lambda_{0}, \lambda_{1}, \ldots, and let N(0)=0N(0) = 0. Show that pn(t)=P(N(t)=n)p_{n}(t) = \mathbb {P}\left(N(t) = n\right) is given by

pn(t)=1λn∑i=0nλie−λit∏j=0j≠inλjλj−λi p_{n}(t) = \frac{1}{\lambda _{n}} \sum _{i = 0}^{n} \lambda _{i} e^{-\lambda _{i} t} \prod _{\substack {j = 0 \\ j \neq i}}^{n} \frac{\lambda _{j}}{\lambda _{j}-\lambda _{i}}

provided that λi≠λj\lambda_{i} \neq \lambda_{j} whenever i≠ji \neq j.

?
Problem 6.8.7

Suppose that the general birth process of the previous exercise is such that ∑nλn−1<∞\sum_{n} \lambda_{n}^{-1} < \infty. Show that λnpn(t)→f(t)\lambda_{n} p_{n}(t) \rightarrow f(t) as n→∞n \rightarrow \infty where ff is the density function of the random variable T=sup⁡{t:N(t)<∞}T = \sup \left\{ t: N(t) < \infty \right\}. Deduce that E[N(t)∣N(t)<∞]\mathbb {E}\left[N(t) \mid N(t) < \infty \right] is finite or infinite depending on the convergence or divergence of ∑nnλn−1\sum_{n} n \lambda_{n}^{-1}.

Find the Laplace transform of ff in closed form for the case when λn=(n+12)2\lambda_{n} = \left(n+\frac{1}{2}\right)^{2}, and deduce an expression for ff.

?
Problem 6.8.8

A traffic light is green at time 00, and subsequently alternates between green and red at the instants of a Poisson process with intensity λ\lambda. Starting from time x>0x > 0, let W(x)W(x) be the waiting time until the light is green for the first time. Find the distribution of W(x)W(x).

?
Problem 6.8.9

Conditional property of simple birth. Let X={X(t):t≥0}X = \left\{ X(t): t \geq 0\right\} be a simple birth process with rate λ\lambda, and let X(0)=1X(0) = 1. Let b≥1b \geq 1. Show that, conditional on the event {X(t)=b+1}\left\{ X(t) = b+1\right\}, the times of the bb births have the same distribution as the order statistics of a random sample of size bb from the density function

f(x)=λe−λ(t−x)1−e−λt,0≤x≤t f(x) = \frac{\lambda e^{-\lambda (t-x)}}{1-e^{-\lambda t}}, \quad 0 \leq x \leq t
?
Problem 6.8.10

Boulders fall down a chute (or couloir) at the instants of a Poisson process with intensity λ\lambda, and mountaineers ascend the chute at the instants of a Poisson process with intensity μ\mu (the two processes are independent of one another). If a fall and an ascent occur during any interval of length cc or less, it is said that a coincidence has occurred. Show that the time TT until the first coincidence has mean

E[T]=1λ+μ{1+λ2+μ2λμ+2e−(λ+μ)c−e−2(λ+μ)c}/{1−e−2(λ+μ)c} \mathbb {E}\left[T\right] = \frac{1}{\lambda +\mu }\left\{ 1+\frac{\lambda ^{2}+\mu ^{2}}{\lambda \mu }+2 e^{-(\lambda +\mu ) c}-e^{-2(\lambda +\mu ) c}\right\} /\left\{ 1-e^{-2(\lambda +\mu ) c}\right\}

Show that cE[T]→(2λμ)−1c \mathbb {E}\left[T\right] \rightarrow (2 \lambda \mu )^{-1} as c↓0c \downarrow 0. Can you prove the last result directly?

?
Problem 6.8.11

Let S1,S2S_{1}, S_{2} be the times of the first two arrivals in a Poisson process with intensity λ\lambda, in order of arrival. Show that

P(s<S1≤t<S2)=λ(t−s)e−λt,0<s<t<∞ \mathbb {P}\left(s < S_{1} \leq t < S_{2}\right) = \lambda (t-s) e^{-\lambda t}, \quad 0 < s < t < \infty

and deduce the joint density function of S1S_{1} and S2S_{2}.

?
Problem 6.8.12

A freelance salesman is paid RR units for each sale, and commissions arrive at the instants of a Poisson process of rate λ\lambda. Living costs consume his resources at unit rate. If his initial wealth is S>0S > 0, show that the probability that he ever becomes bankrupt is aS/Ra^{S / R}, where aa is the smallest x>0x > 0 such that x=eλR(x−1)x = e^{\lambda R(x-1)}.

Prove that a∈(0,1)a \in (0,1) if λR>1\lambda R > 1, while a=1a = 1 if λR<1\lambda R < 1.

?
Problem 6.9.1

Let λμ>0\lambda \mu > 0 and let XX be a Markov chain on {1,2}\left\{ 1,2\right\} with generator

G=[−μμλ−λ] \mathbf{G} =\left[\begin{smallmatrix} -\mu & \mu \\ \lambda & -\lambda \end{smallmatrix}\right]
?
(a)

Write down the forward equations and solve them for the transition probabilities pij(t),i,j=p_{i j}(t), i, j = 1,2.

(b)

Calculate Gn\mathbf{G}^{n} and hence find ∑n=0∞(tn/n!)Gn\sum_{n = 0}^{\infty }\left(t^{n} / n!\right) \mathbf{G}^{n}. Compare your answer with that to part (a).

(c)

Solve the equation πG=0\pi \mathbf{G} = \mathbf{0} in order to find the stationary distribution. Verify that pij(t)→πjp_{i j}(t) \rightarrow \pi_{j} as t→∞t \rightarrow \infty.

Problem 6.9.2

As a continuation of the previous exercise, find:

?
(a)

P(X(t)=2∣X(0)=1,X(3t)=1)\mathbb {P}\left(X(t) = 2 \mid X(0) = 1, X(3 t) = 1\right),

(b)

P(X(t)=2∣X(0)=1,X(3t)=1,X(4t)=1)\mathbb {P}\left(X(t) = 2 \mid X(0) = 1, X(3 t) = 1, X(4 t) = 1\right).

Problem 6.9.3

Jobs arrive in a computer queue in the manner of a Poisson process with intensity λ\lambda. The central processor handles them one by one in the order of their arrival, and each has an exponentially distributed runtime with parameter μ\mu, the runtimes of different jobs being independent of each other and of the arrival process. Let X(t)X(t) be the number of jobs in the system (either running or waiting) at time tt, where X(0)=0X(0) = 0. Explain why XX is a Markov chain, and write down its generator. Show that a stationary distribution exists if and only if λ<μ\lambda < \mu, and find it in this case.

?
Problem 6.9.4

Let X={X(t):t≥0}X = \left\{ X(t): t \geq 0\right\} be a Markov chain having stationary distribution π\mathbf{\pi }. We may sample XX at the times of a Poisson process: let NN be a Poisson process with intensity λ\lambda, independent of XX, and define Yn=X(Tn)Y_{n} = X\left(T_{n}\right). Show that Y={Yn:n≥0}Y = \left\{ Y_{n}: n \geq 0\right\} is a discrete-time Markov chain with the same stationary distribution as XX. (This exemplifies the 'Pasta' property: Poisson arrivals see time averages.) [The full assumption of the independence of NN and XX is not necessary for the conclusion. It suffices that {N(s):s≥t}\left\{ N(s): s \geq t\right\} be independent of {X(s):s≤t}\left\{ X(s): s \leq t\right\}, a property known as 'lack of anticipation'. It is not even necessary that XX be Markov; the Pasta property holds for many suitable ergodic processes.]

?
Problem 6.9.5

Let XX be a continuous-time Markov chain with generator G\mathbf{G} satisfying gi=−gii>0g_{i} = -g_{i i} > 0 for all ii. Let HA=inf⁡{t≥0:X(t)∈A}H_{A} = \inf \left\{ t \geq 0: X(t) \in A\right\} be the hitting time of the set AA of states, and let ηj=Pj(HA<∞)\eta_{j} = \mathbb {P}_{j}\left(H_{A} < \infty \right) be the chance of ever reaching AA from jj. By using properties of the jump chain, which you may assume to be well behaved, show that ∑kgjkηk=0\sum_{k} g_{j k} \eta_{k} = 0 for j∉Aj \notin A.

?
Problem 6.9.6

In continuation of the preceding exercise, let μj=Ej[HA]\mu_{j} = \mathbb {E}_{j}\left[H_{A}\right]. Show that the vector μ\mathbf{\mu } is the minimal non-negative solution of the equations

μj=0 if j∈A,1+∑k∈Sgjkμk=0 if j∉A \mu _{j} = 0 \quad \text{ if } j \in A, \quad 1+\sum _{k \in S} g_{j k} \mu _{k} = 0 \quad \text{ if } j \notin A
?
Problem 6.9.7

Let XX be a continuous-time Markov chain with transition probabilities pij(t)p_{i j}(t) and define Fi=inf⁡{t>T1:X(t)=i}F_{i} = \inf \left\{ t > T_{1}: X(t) = i\right\} where T1T_{1} is the time of the first jump of XX. Show that, if gii≠0g_{i i} \neq 0, then Pi(Fi<∞)=1\mathbb {P}_{i}\left(F_{i} < \infty \right) = 1 if and only if ii is recurrent.

?
Problem 6.9.8

Let XX be the simple symmetric random walk on the integers in continuous time, so that

pi,i+1(h)=pi,i−1(h)=12λh+o(h) p_{i, i+1}(h) = p_{i, i-1}(h) = \frac{1}{2} \lambda h+\mathrm{o}(h)

Show that the walk is recurrent. Let TT be the time spent visiting mm during an excursion from 0. Find the distribution of TT.

?
Problem 6.9.9

Let ii be a transient state of a continuous-time Markov chain XX with X(0)=iX(0) = i. Show that the total time spent in state ii has an exponential distribution.

?
Problem 6.9.10

Let XX be an asymmetric simple random walk in continuous time on the non-negative integers with retention at 00, so that

pij(h)={λh+o(h) if j=i+1,i≥0μh+o(h) if j=i−1,i≥1 p_{i j}(h) = \begin{cases} \lambda h+\mathrm{o}(h) & \text{ if } j = i+1, i \geq 0 \\ \mu h+\mathrm{o}(h) & \text{ if } j = i-1, i \geq 1\end{cases}

Suppose that X(0)=0X(0) = 0 and λ>μ\lambda > \mu. Show that the total time VrV_{r} spent in state rr is exponentially distributed with parameter λ−μ\lambda -\mu.

Assume now that X(0)X(0) has some general distribution with probability generating function GG. Find the expected amount of time spent at 0 in terms of GG.

?
Problem 6.9.11

Let XX be the continuous-time Markov chain on the finite state space SS with generator G=(gij)\mathbf{G} = \left(g_{i j}\right). Show from first principles that the transition probabilities satisfy

pij(h)={1+giih+o(h) if i=jgijh+o(h) if i≠j p_{i j}(h) = \begin{cases} 1+g_{i i} h+\mathrm{o}(h) & \text{ if } i = j \\ g_{i j} h+\mathrm{o}(h) & \text{ if } i \neq j\end{cases}
?
Problem 6.9.12

Let XX be a Markov chain on the integers Z\mathbb {Z} with generator satisfying gi,i−1=gi,i+1=2ig_{i, i-1} = g_{i, i+1} = 2^{i} for i∈Zi \in \mathbb {Z}, and gi,j=0g_{i, j} = 0 for other pairs (i,j)(i, j) with i≠ji \neq j. Does XX explode?

?
Problem 6.9.13

A population grows subject to the threat of total annihilation. It is modelled as a Markov chain with generator G=(gij)\mathbf{G} = \left(g_{i j}\right) satisfying

gi,i+1=1i+2,i≥0gi,0=1(i+1)(i+2),i≥1 \begin{aligned} g_{i, i+1} & = \frac{1}{i+2}, & & i \geq 0 \\ g_{i, 0} & = \frac{1}{(i+1)(i+2)}, & & i \geq 1 \end{aligned}

the other off-diagonal elements of G\mathbf{G} being 0. Show that the chain is null recurrent.

?
Problem 6.9.14

Let Zn=X(nh)Z_{n} = X(n h) where XX is a Markov chain and h>0h > 0. Show that ii is recurrent for ZZ if and only if it is recurrent for XX. Show that ZZ is irreducible if and only if XX is irreducible.

?
Problem 6.10.1

Let N<∞N < \infty, and let Q\mathcal{Q} be the space of N×NN \times N matrices with real entries, with norm

∣Q∣=sup⁡x≠0∣Qx∣∣x∣,Q∈Q |Q| = \sup _{\mathbf{x} \neq \mathbf{0}} \frac{|Q \mathbf{x}|}{|\mathbf{x}|}, \quad Q \in \mathcal{Q}

where ∣y∣\left|\mathbf{y}\right| is the Euclidean norm of the vector y\mathbf{y}, and the supremum is over all non-zero column vectors.

?
(a)

Show for Q1,Q2∈QQ_{1}, Q_{2} \in \mathcal{Q} that

∣Q1+Q2∣≤∣Q1∣+∣Q2∣,∣Q1Q2∣≤∣Q1∣⋅∣Q2∣ \left|Q_{1}+Q_{2}\right| \leq \left|Q_{1}\right|+\left|Q_{2}\right|, \quad \left|Q_{1} Q_{2}\right| \leq \left|Q_{1}\right| \cdot \left|Q_{2}\right|
(b)

Show for Q∈QQ \in \mathcal{Q} that En:=∑k=0nQk/k!E_{n} := \sum_{k = 0}^{n} Q^{k} / k! converges with respect to ∣⋅∣\left|\cdot \right| to a limit which we denote E=E(Q)E = E(Q).

(c)

Show that E(Q1+Q2)=E(Q1)E(Q2)E(Q_{1}+Q_{2}) = E(Q_{1}) E(Q_{2}) if Q1Q_{1} and Q2Q_{2} are commuting elements of Q\mathcal{Q}.

Problem 6.10.2

Let YY be an irreducible discrete-time Markov chain on a countably infinite state space SS, having transition matrix Y=(yij)\mathbf{Y} = \left(y_{i j}\right) satisfying yii=0y_{i i} = 0 for all states ii, and with stationary distribution ν\mathbf{\nu }. Construct a continuous-time process XX on SS for which YY is the jump chain, such that XX has no stationary distribution.

?
Problem 6.10.3

Let X=(X(t):t≥0)X = (X(t): t \geq 0) be a Markov chain on Z\mathbb {Z} with generator G=(gij)\mathbf{G} = \left(g_{i j}\right) given by

gi,i−1=i2+1,gi,i=−2(i2+1),gi,i+1=i2+1,i∈Z g_{i, i-1} = i^{2}+1, \quad g_{i, i} = -2\left(i^{2}+1\right), \quad g_{i, i+1} = i^{2}+1, \quad i \in \mathbb {Z}

Show that XX is recurrent. Is XX positive recurrent?

?
Problem 6.10.4

Let XX be a Markov chain on Z\mathbb {Z} with generator G=(gij)G = \left(g_{i j}\right) satisfying

gi,i−1=3∣i∣,gi,i=−3∣i∣+1,gi,i+1=2⋅3∣i∣,i∈Z g_{i, i-1} = 3^{\left|i\right|}, \quad g_{i, i} = -3^{\left|i\right|+1}, \quad g_{i, i+1} = 2 \cdot 3^{\left|i\right|}, \quad i \in \mathbb {Z}

Show that XX is transient, but has an invariant distribution. Explain.

?
Problem 6.10.5

Let X=(Xt:t≥0)X = \left(X_{t}: t \geq 0\right) be a Markov chain with generator G=(gij)\mathbf{G} = \left(g_{i j}\right) on the finite state space SS, and let f:S→Rf: S \rightarrow \mathbb {R} be a function, which we identify with the vector f=(f(i):i∈S)f = (f(i): i \in S). Show that

Gf(i)=∑j∈Sgij(f(j)−f(i)) \mathbf{G} f(i) = \sum _{j \in S} g_{i j}(f(j)-f(i))

where Gf\mathbf{G} f denotes the standard matrix multiplication. Show that

Gf(i)=lim⁡t→01t[Ei[f(Xt)]−f(i)],i∈S \mathbf{G} f(i) = \lim _{t \rightarrow 0} \frac{1}{t}\left[\mathbb {E}_{i}\left[f\left(X_{t}\right)\right]-f(i)\right], \quad i \in S

and deduce that

Ei[f(Xt)]=f(i)+∫0tEi[Gf(Xs)]ds \mathbb {E}_{i}\left[f\left(X_{t}\right)\right] = f(i)+\int _{0}^{t} \mathbb {E}_{i}\left[\mathbf{G} f\left(X_{s}\right)\right] d s
?
Problem 6.11.1

Describe the jump chain for a birth-death process with rates λn\lambda_{n} and μn\mu_{n}.

?
Problem 6.11.2

Consider an immigration-death process XX, being a birth-death process with birth rates λn=λ\lambda_{n} = \lambda and death rates μn=nμ\mu_{n} = n \mu. Find the transition matrix of the jump chain YY, and show that it has as stationary distribution

πn=12(n!)(1+nρ)ρne−ρ \pi _{n} = \frac{1}{2(n!)}\left(1+\frac{n}{\rho }\right) \rho ^{n} e^{-\rho }

where ρ=λ/μ\rho = \lambda / \mu. Explain why this differs from the stationary distribution of XX.

?
Problem 6.11.3

Consider the birth-death process XX with λn=nλ\lambda_{n} = n \lambda and μn=nμ\mu_{n} = n \mu for all n≥0n \geq 0. Suppose X(0)=1X(0) = 1 and let η(t)=P1(X(t)=0)\eta (t) = \mathbb {P}_{1}\left(X(t) = 0\right). Show that η\eta satisfies the differential equation

η′(t)+(λ+μ)η(t)=μ+λη(t)2 \eta ^{\prime }(t)+(\lambda +\mu ) \eta (t) = \mu +\lambda \eta (t)^{2}

Hence find η(t)\eta (t), and calculate P1(X(t)=0∣X(u)=0)\mathbb {P}_{1}\left( X(t) = 0 \mid X(u) = 0 \right) for 0<t<u0 < t < u.

?
Problem 6.11.4

For the birth-death process of the previous exercise with λ<μ\lambda < \mu, show that the distribution of X(t)X(t), conditional on the event {X(t)>0}\left\{ X(t) > 0\right\}, converges as t→∞t \rightarrow \infty to a geometric distribution.

?
Problem 6.11.5

Let XX be a birth-death process with λn=nλ\lambda_{n} = n \lambda and μn=nμ\mu_{n} = n \mu, and suppose X(0)=1X(0) = 1. Show that the time TT at which X(t)X(t) first takes the value 0 satisfies

E[T∣T<∞]={1λlog⁡(μμ−λ) if λ<μ1μlog⁡(λλ−μ) if λ>μ \mathbb {E}\left[T \mid T < \infty \right] = \begin{cases} \frac{1}{\lambda } \log \left(\frac{\mu }{\mu -\lambda }\right) & \text{ if } \lambda < \mu \\ \frac{1}{\mu } \log \left(\frac{\lambda }{\lambda -\mu }\right) & \text{ if } \lambda > \mu \end{cases}

What happens when λ=μ\lambda = \mu?

?
Problem 6.11.6

Let XX be the birth-death process of Exercise (6.11.5) with λ≠μ\lambda \neq \mu, and let Vr(t)V_{r}(t) be the total amount of time the process has spent in state r≥0r \geq 0, up to time tt. Find the distribution of V1(∞)V_{1}(\infty ) and the generating function ∑rsrE[Vr(t)]\sum_{r} s^{r} \mathbb {E}\left[V_{r}(t)\right]. Hence show in two ways that E[V1(∞)]=[max⁡{λ,μ}]−1\mathbb {E}\left[V_{1}(\infty )\right] = [\max \left\{ \lambda , \mu \right\} ]^{-1}. Show further that E[Vr(∞)]=λr−1r−1[max⁡{λ,μ}]−r\mathbb {E}\left[V_{r}(\infty )\right] = \lambda^{r-1} r^{-1}[\max \left\{ \lambda , \mu \right\} ]^{-r}.

?
Problem 6.11.7

Repeat the calculations of Exercise (6.11.6) in the case λ=μ\lambda = \mu.

?
Problem 6.11.8

Consider a birth-death process XX with birth rates λn>0\lambda_{n} > 0 for n≥0n \geq 0, death rates μn>0\mu_{n} > 0 for n>0n > 0, and μ0=0\mu_{0} = 0. Let X(0)=n>0X(0) = n > 0, and let DnD_{n} be the time until the process first takes the value n−1n-1.

?
(a)

Show that dn=E[Dn]d_{n} = \mathbb {E}\left[D_{n}\right] satisfies

λndn+1=μndn−1,n≥1 \lambda _{n} d_{n+1} = \mu _{n} d_{n}-1, \quad n \geq 1
(b)

Show that the moment generating function Mn(θ)=E[eθDn]M_{n}(\theta ) = \mathbb {E}\left[e^{\theta D_{n}}\right] satisfies

(λn+μn−θ)Mn(θ)=μn+λnMn(θ)Mn+1(θ),n≥1 \left(\lambda _{n}+\mu _{n}-\theta \right) M_{n}(\theta ) = \mu _{n}+\lambda _{n} M_{n}(\theta ) M_{n+1}(\theta ), \quad n \geq 1
Problem 6.11.9

In a model for a biofilm population, we assume there are nn colonizable 'niches' (or 'food sources'). Let X(t)X(t) be the number of occupied niches at time tt, and assume XX is a Markov chain that evolves as follows. The lifetime of any colony is exponentially distributed with parameter μ\mu; if X(t)=iX(t) = i, the rate of establishment of a new colony in an empty niche is λi(n−i)\lambda i(n-i). The usual independence may be assumed.

For X(0)=1,2,…X(0) = 1,2, \ldots, find the mean time until the population is extinct, which is to say that no niches are occupied. Discuss the implications when nn is large.

?
Problem 6.12.1

Customers entering a shop are served in the order of their arrival by the single server. They arrive in the manner of a Poisson process with intensity λ\lambda, and their service times are independent exponentially distributed random variables with parameter μ\mu. By considering the jump chain, show that the expected duration of a busy period BB of the server is (μ−λ)−1(\mu -\lambda )^{-1} when λ<μ\lambda < \mu. (The busy period runs from the moment a customer arrives to find the server free until the earliest subsequent time when the server is again free.)

?
Problem 6.12.2

Immigrants arrive at the instants of a Poisson process of rate ν\nu, and each independently founds a simple birth process of rate λ\lambda. At the instants of an independent Poisson process of rate δ\delta, the population is annihilated. Find the probability generating function of the population X(t)X(t), given that X(0)=0X(0) = 0.

?
Problem 6.12.3

In the framework of Exercise (6.12.2), suppose that each immigrant gives rise to a simple birth-death process of rates λ\lambda and μ\mu. Show that the mean population size stays bounded if and only if δ>λ−μ\delta > \lambda -\mu.

?
Problem 6.12.4

The queue M/G/∞M / G / \infty. An ftp server receives clients at the times of a Poisson process with parameter λ\lambda, beginning at time 0. The ii th client remains connected for a length SiS_{i} of time, where the SiS_{i} are independent identically distributed random variables, independent of the process of arrivals. Assuming that the server has an infinite capacity, show that the number of clients being serviced at time tt has the Poisson distribution with parameter λ∫0t[1−G(x)]dx\lambda \int_{0}^{t}[1-G(x)] d x, where GG is the common distribution function of the SiS_{i}. Show that the mean of this distribution converges to λE[S]\lambda \mathbb {E}\left[S\right] as t→∞t \rightarrow \infty.

?
Problem 6.13.1

In a certain town at time t=0t = 0 there are no bears. Brown bears and grizzly bears arrive as independent Poisson processes BB and GG with respective intensities β\beta and γ\gamma.

?
(a)

Show that the first bear is brown with probability β/(β+γ)\beta /(\beta +\gamma ).

(b)

Find the probability that between two consecutive brown bears, there arrive exactly rr grizzly bears.

(c)

Given that B(1)=1B(1) = 1, find the expected value of the time at which the first bear arrived.

Problem 6.13.2

Let Π\Pi be the points of a non-homogeneous Poisson process on Rd\mathbb {R}^{d} with intensity function λ\lambda. Let S=∑x∈Πg(x)S = \sum_{\mathbf{x} \in \Pi } g(\mathbf{x}) where gg is a (measurable) function which we assume for convenience to be non-negative.

?
(a)

Show directly that E[S]=∫Rdg(x)λ(x)dx\mathbb {E}\left[S\right] = \int_{\mathbb {R}^{d}} g(\mathbf{x}) \lambda (\mathbf{x}) d \mathbf{x} and Var⁡[S]=∫Rdg(x)2λ(x)dx\operatorname {Var}\left[S\right] = \int_{\mathbb {R}^{d}} g(\mathbf{x})^{2} \lambda (\mathbf{x}) d \mathbf{x}, provided these integrals converge.

(b)

Show that

E[e−tS]=exp⁡{−∫Rd(1−e−tg(x))λ(x)dx},t>0 \mathbb {E}\left[e^{-t S}\right] = \exp \left\{ -\int _{\mathbb {R}^{d}}\left(1-e^{-t g(\mathbf{x})}\right) \lambda (\mathbf{x}) d \mathbf{x}\right\} , \quad t > 0

and deduce that P(S<∞)=1\mathbb {P}\left(S < \infty \right) = 1 if ∫Rdmin⁡{1,g(x)}λ(x)dx<∞\int_{\mathbb {R}^{d}} \min \left\{ 1, g(\mathbf{x})\right\} \lambda (\mathbf{x}) d \mathbf{x} < \infty.

(c)

If the integral condition of part (b) holds, show that the characteristic function ϕ\phi of SS satisfies

ϕ(t)=exp⁡{−∫Rd(1−eitg(x))λ(x)dx},t∈R \phi (t) = \exp \left\{ -\int _{\mathbb {R}^{d}}\left(1-e^{i t g(\mathbf{x})}\right) \lambda (\mathbf{x}) d \mathbf{x}\right\} , \quad t \in \mathbb {R}
Problem 6.13.3

Let Π\Pi be a Poisson process with constant intensity λ\lambda on the surface of the sphere of R3\mathbb {R}^{3} with radius 1. Let PP be the process given by the (X,Y)(X, Y) coordinates of the points projected on a plane passing through the centre of the sphere. Show that PP is a Poisson process, and find its intensity function.

?
Problem 6.13.4

Repeat Exercise (6.13.3), when Π\Pi is a homogeneous Poisson process on the ball {(x1,x2,x3):x12+x22+x32≤1}\left\{ \left(x_{1}, x_{2}, x_{3}\right) : x_{1}^{2}+x_{2}^{2}+x_{3}^{2} \leq 1\right\}.

?
Problem 6.13.5

You stick pins in a Mercator projection of the Earth in the manner of a Poisson process with constant intensity λ\lambda. What is the intensity function of the corresponding process on the globe? What would be the intensity function on the map if you formed a Poisson process of constant intensity λ\lambda of meteorite strikes on the surface of the Earth?

?
Problem 6.13.6

The rr th point TrT_{r} of a Poisson process NN of constant intensity λ\lambda on R+\mathbb {R}_{+}gives rise to an effect Xre−α(t−Tr)X_{r} e^{-\alpha \left(t-T_{r}\right)} at time t≥Trt \geq T_{r}, where the XrX_{r} are independent and identically distributed with finite variance. Find the mean and variance of the total effect S(t)=∑r=1N(t)Xre−α(t−Tr)S(t) = \sum_{r = 1}^{N(t)} X_{r} e^{-\alpha \left(t-T_{r}\right)} in terms of the first two moments of the XrX_{r}, and calculate Cov⁡[S(s),S(t)]\operatorname {Cov}\left[S(s), S(t)\right].

What is the behaviour of the correlation ρ(S(s),S(t))\rho (S(s), S(t)) as s→∞s \rightarrow \infty with t−st-s fixed?

?
Problem 6.13.7

Let NN be a non-homogeneous Poisson process on R+\mathbb {R}_{+}with intensity function λ\lambda. Find the joint density of the first two inter-event times, and deduce that they are not in general independent.

?
Problem 6.13.8

Let {Nr(t):r≥1}\left\{ N_{r}(t): r \geq 1\right\} be a collection of independent Poisson processes on R+\mathbb {R}_{+} with respective constant intensities {λr:r≥1}\left\{ \lambda_{r}: r \geq 1\right\}, such that ∑rλr=λ<∞\sum_{r} \lambda_{r} = \lambda < \infty. Set N(t)=∑rNr(t)N(t) = \sum_{r} N_{r}(t), and let II denote the index of the process supplying the first point in NN, occurring at time TT. Show that

P(I=i,T≥t)=P(I=i)P(T≥t)=λiλe−λt,i≥1 \mathbb {P}\left(I = i, T \geq t\right) = \mathbb {P}\left(I = i\right) \mathbb {P}\left(T \geq t\right) = \frac{\lambda _{i}}{\lambda } e^{-\lambda t}, \quad i \geq 1
?
Problem 6.13.9

Let Π\Pi be a Poisson process on R2\0\mathbb {R}^{2} \backslash \mathbf{0} with intensity function λ(u,v)=(u2+v2)−3/2\lambda (u, v) = \left(u^{2}+v^{2}\right)^{-3 / 2}. Each point (U,V)∈Π(U, V) \in \Pi gives rise to a line Ux+Vy=1U x+V y = 1 in the x/yx / y-plane L\mathbb {L}.

?
(a)

Let (p,θ)(p, \theta ) be the polar coordinates of the foot of the perpendicular from the origin onto the line ux+vy=1u x+v y = 1 in the x/yx / y-plane L\mathbb {L}. Express (p,θ)(p, \theta ) in terms of (u,v)(u, v).

(b)

Show that the Poisson line process is mapped by the map of part (a) to a uniform Poisson process on the strip S⊥=R×[0,π)S^{\perp } = \mathbb {R} \times [0, \pi ).

(c)

Show that the line process in L\mathbb {L} is invariant under translations and rotations of L\mathbb {L}.

Problem 6.13.10

A large continent is traversed by a doubly-infinite straight freeway on which lorries are parked at the points of a Poisson process with constant intensity 1. The masses of the lorries are independent, identically distributed random variables that are independent of the parking places. Let GG be the gravitational attraction due to the lorries on a pedestrian of unit mass standing beside the freeway. You may take the gravitational constant to be 1.

Show that GG has characteristic function of the form ϕ(t)=exp⁡(−c∣t∣1/2)\phi (t) = \exp \left(-c|t|^{1 / 2}\right) where c>0c > 0. Express cc in terms of the mean of a typical mass MM.

?
Problem 6.14.1

Let P\mathbf{P} be a stochastic matrix on the finite set Θ\Theta with stationary distribution π\mathbf{\pi }. Define the inner product ⟨x,y⟩=∑k∈Θxkykπk\langle \mathbf{x}, \mathbf{y}\rangle = \sum_{k \in \Theta } x_{k} y_{k} \pi_{k}, and let l2(π)={x∈RΘ:⟨x,x⟩<∞}l^{2}(\pi ) = \left\{ \mathbf{x} \in \mathbb {R}^{\Theta }:\langle \mathbf{x}, \mathbf{x}\rangle < \infty \right\}. Show, in the obvious notation, that P\mathbf{P} is reversible with respect to π\pi if and only if ⟨x,Py⟩=⟨Px,y⟩\langle \mathbf{x}, \mathbf{P y}\rangle = \langle \mathbf{P x}, \mathbf{y}\rangle for all x,y∈l2(π)\mathbf{x}, \mathbf{y} \in l^{2}(\pi ).

?
Problem 6.14.2

Show that a possible choice for the acceptance probabilities in Hastings's general algorithm is

bij=πjgjiπigij+πjgji b_{i j} = \frac{\pi _{j} g_{j i}}{\pi _{i} g_{i j}+\pi _{j} g_{j i}}

where G=(gij)\mathbf{G} = \left(g_{i j}\right) is the proposal matrix.

?
Problem 6.14.3

Let SS be a countable set. For each j∈Sj \in S, the sets Ajk,k∈SA_{j k}, k \in S, form a partition of the interval [0,1][0,1]. Let g:S×[0,1]→Sg: S \times [0,1] \rightarrow S be given by g(j,u)=kg(j, u) = k if u∈Ajku \in A_{j k}. The sequence {Xn:n≥0}\left\{ X_{n}: n \geq 0\right\} of random variables is generated recursively by Xn+1=g(Xn,Un+1),n≥0X_{n+1} = g\left(X_{n}, U_{n+1}\right), n \geq 0, where {Un:n≥1}\left\{ U_{n}: n \geq 1\right\} are independent random variables with the uniform distribution on [0,1][0,1]. Show that XX is a Markov chain, and find its transition matrix.

?
Problem 6.14.4

Let U=(ust)\mathbf{U} = \left(u_{s t}\right) be a finite ∣S∣×∣T∣\left|S\right| \times \left|T\right| stochastic matrix. Dobrushin's ergodic coefficient is defined to be

d(U)=12sup⁡i,j∈S∑t∈T∣uit−ujt∣ d(\mathbf{U}) = \frac{1}{2} \sup _{i, j \in S} \sum _{t \in T}\left|u_{i t}-u_{j t}\right|
?
(a)

Show that, if V\mathbf{V} is a finite ∣T∣×∣U∣\left|T\right| \times \left|U\right| stochastic matrix, then d(UV)≤d(U)d(V)d(\mathbf{U V}) \leq d(\mathbf{U}) d(\mathbf{V}).

(b)

Let XX and YY be discrete-time Markov chains with the same transition matrix P\mathbf{P}, and show that

∑k∣P(Xn=k)−P(Yn=k)∣≤d(P)n∑k∣P(X0=k)−P(Y0=k)∣ \sum _{k}\left|\mathbb {P}\left(X_{n} = k\right)-\mathbb {P}\left(Y_{n} = k\right)\right| \leq d(\mathbf{P})^{n} \sum _{k}\left|\mathbb {P}\left(X_{0} = k\right)-\mathbb {P}\left(Y_{0} = k\right)\right|
Problem 6.14.5

Let π\mathbf{\pi } be a positive mass function on the finite set Θ\Theta, and let P\mathbf{P} be the transition matrix of an irreducible, aperiodic Markov chain with stationary distribution π\mathbf{\pi }. Let W=(W(i):i∈Θ)W = (W(i): i \in \Theta ) be a vector of random variables such that P(W(i)=j)=pij\mathbb {P}\left(W(i) = j\right) = p_{i j} for i,j∈Θi, j \in \Theta, and use WW as an update rule in the coupling-from-the-past algorithm for sampling from π\pi.

?
(a)

If the W(i),i∈ΘW(i), i \in \Theta, are independent, show that the coalescence time is a.s. finite.

(b)

Give two examples of situations in which the coalescence time is a.s. infinite.

Problem 6.14.6

Show that the Ising distribution of Example (6.14.2) satisfies the FKG lattice condition (6.14.20).

?
Problem 6.15.1

Classify the states of the discrete-time Markov chains with state space S={1,2,3,4}S = \left\{ 1,2,3,4\right\} and the following transition matrices:

?
(a)
[13230012120014014120001] \left[\begin{smallmatrix} \frac{1}{3} & \frac{2}{3} & 0 & 0 \\ \frac{1}{2} & \frac{1}{2} & 0 & 0 \\ \frac{1}{4} & 0 & \frac{1}{4} & \frac{1}{2} \\ 0 & 0 & 0 & 1 \end{smallmatrix}\right]

Calculate f34(n)f_{34}(n), and deduce that the probability of ultimate absorption in state 44, starting from 33, equals 23\frac{2}{3}.

(b)
[01212013002310000010] \left[\begin{smallmatrix} 0 & \frac{1}{2} & \frac{1}{2} & 0 \\ \frac{1}{3} & 0 & 0 & \frac{2}{3} \\ 1 & 0 & 0 & 0 \\ 0 & 0 & 1 & 0 \end{smallmatrix}\right]

Find the mean recurrence times of the states.

Problem 6.15.2

A transition matrix is called doubly stochastic if all its column sums equal 11, that is, if ∑ipij=1\sum_{i} p_{i j} = 1 for all j∈Sj \in S.

?
(a)

Show that if a finite chain has a doubly stochastic transition matrix, then all its states are positive recurrent, and that if it is, in addition, irreducible and aperiodic then pij(n)→N−1p_{i j}(n) \rightarrow N^{-1} as n→∞n \rightarrow \infty, where NN is the number of states.

(b)

Show that, if an infinite irreducible chain has a doubly stochastic transition matrix, then its states are either all null recurrent or all transient.

Problem 6.15.3

Prove that intercommunicating states of a Markov chain have the same period.

?
Problem 6.15.4
?
(a)

Show that for each pair i,ji, j of states of an irreducible aperiodic chain, there exists N=N(i,j)N = N(i, j) such that pij(n)>0p_{i j}(n) > 0 for all n≥Nn \geq N.

(b)

Let XX and YY be independent irreducible aperiodic chains with the same state space SS and transition matrix P\mathbf{P}. Show that the bivariate chain Zn=(Xn,Yn),n≥0Z_{n} = \left(X_{n}, Y_{n}\right), n \geq 0, is irreducible and aperiodic.

(c)

Show that the bivariate chain ZZ may be reducible if XX and YY are periodic.

Problem 6.15.5

Suppose {Xn:n≥0}\left\{ X_{n}: n \geq 0\right\} is a discrete-time Markov chain with X0=iX_{0} = i. Let NN be the total number of visits made subsequently by the chain to the state jj. Show that

P(N=n)={1−fij if n=0fij(fjj)n−1(1−fjj) if n≥1 \mathbb {P}\left(N = n\right) = \begin{cases} 1-f_{i j} & \text{ if } n = 0 \\ f_{i j}\left(f_{j j}\right)^{n-1}\left(1-f_{j j}\right) & \text{ if } n \geq 1\end{cases}

and deduce that P(N=∞)=1\mathbb {P}\left(N = \infty \right) = 1 if and only if fij=fjj=1f_{i j} = f_{j j} = 1.

?
Problem 6.15.6

Let ii and jj be two states of a discrete-time Markov chain. Show that if ii communicates with jj, then there is positive probability of reaching jj from ii without revisiting ii in the meantime. Deduce that, if the chain is irreducible and recurrent, then the probability fijf_{i j} of ever reaching jj from ii equals 1 for all ii and jj.

?
Problem 6.15.7

Let {Xn:n≥0}\left\{ X_{n}: n \geq 0\right\} be a recurrent irreducible Markov chain on the state space SS with transition matrix P\mathbf{P}, and let x\mathbf{x} be a positive solution of the equation x=xP\mathbf{x} = \mathbf{x P}.

?
(a)

Show that

qij(n)=xjxipji(n),i,j∈S,n≥1 q_{i j}(n) = \frac{x_{j}}{x_{i}} p_{j i}(n), \quad i, j \in S, n \geq 1

defines the nn-step transition probabilities of a recurrent irreducible Markov chain on SS whose first-passage probabilities are given by

gij(n)=xjxilji(n),i≠j,n≥1 g_{i j}(n) = \frac{x_{j}}{x_{i}} l_{j i}(n), \quad i \neq j, n \geq 1

where lji(n)=Pj(Xn=i,T>n)l_{j i}(n) = \mathbb {P}_{j}\left(X_{n} = i, T > n\right) and T=min⁡{m>0:Xm=j}T = \min \left\{ m > 0: X_{m} = j\right\}.

(b)

Show that x\mathbf{x} is unique up to a multiplicative constant.

(c)

Let Tj=min⁡{n≥1:Xn=j}T_{j} = \min \left\{ n \geq 1: X_{n} = j\right\} and define hij=Pi(Tj≤Ti)h_{i j} = \mathbb {P}_{i}\left(T_{j} \leq T_{i}\right). Show that xihij=xjhjix_{i} h_{i j} = x_{j} h_{j i} for all i,j∈Si, j \in S.

Problem 6.15.8

The sequence u={un:n≥0}u = \left\{ u_{n}: n \geq 0\right\} is called a 'renewal sequence' if

u0=1,un=∑i=1nfiun−i for n≥1 u_{0} = 1, \quad u_{n} = \sum _{i = 1}^{n} f_{i} u_{n-i} \quad \text{ for } n \geq 1

for some collection f={fn:n≥1}f = \left\{ f_{n}: n \geq 1\right\} of non-negative numbers summing to 1.

?
(a)

Show that uu is a renewal sequence if and only if there exists a Markov chain XX on a countable state space SS such that un=P(Xn=s∣X0=s)u_{n} = \mathbb {P}\left(X_{n} = s \mid X_{0} = s\right), for some recurrent s∈Ss \in S and all n≥1n \geq 1.

(b)

Show that if uu and vv are renewal sequences then so is {unvn:n≥0}\left\{ u_{n} v_{n}: n \geq 0\right\}.

Problem 6.15.9

Consider the symmetric random walk in three dimensions on the set of points {(x,y,z):x,y,z=0,±1,±2,…}\left\{ (x, y, z): x, y, z = 0, \pm 1, \pm 2, \ldots \right\}; this process is a sequence {Xn:n≥0}\left\{ \mathbf{X}_{n}: n \geq 0\right\} of points such that P(Xn+1=Xn+ϵ)=16\mathbb {P}\left(\mathbf{X}_{n+1} = \mathbf{X}_{n}+\mathbf{\epsilon }\right) = \frac{1}{6} for ϵ=(±1,0,0),(0,±1,0),(0,0,±1)\mathbf{\epsilon } = (\pm 1,0,0),(0, \pm 1,0),(0,0, \pm 1). Suppose that X0=(0,0,0)\mathbf{X}_{0} = (0,0,0). Show that

P(X2n=(0,0,0))=(16)2n∑i+j+k=n(2n)!(i!j!k!)2=(12)2n(2nn)∑i+j+k=n(n!3ni!j!k!)2 \mathbb {P}\left(\mathbf{X}_{2 n} = (0,0,0)\right) = \left(\frac{1}{6}\right)^{2 n} \sum _{i+j+k = n} \frac{(2 n)!}{(i!j!k!)^{2}} = \left(\frac{1}{2}\right)^{2 n}\binom {2 n}{n} \sum _{i+j+k = n}\left(\frac{n!}{3^{n} i!j!k!}\right)^{2}

and deduce by Stirling's formula that the origin is a transient state.

?
Problem 6.15.10

Consider the three-dimensional version of the cancer model (6.12.12). If κ=1\kappa = 1, are the empires of Theorem (6.12.14) inevitable in this case?

?
Problem 6.15.11

Let XX be a discrete-time Markov chain with state space S={1,2}S = \left\{ 1,2\right\}, and transition matrix

P=[1−ααβ1−β] \mathbf{P} =\left[\begin{smallmatrix} 1-\alpha & \alpha \\ \beta & 1-\beta \end{smallmatrix}\right]

Classify the states of the chain. Suppose that αβ>0\alpha \beta > 0 and αβ≠1\alpha \beta \neq 1. Find the nn-step transition probabilities and show directly that they converge to the unique stationary distribution as n→∞n \rightarrow \infty. For what values of α\alpha and β\beta is the chain reversible in equilibrium?

?
Problem 6.15.12

NN black balls and NN white balls are placed in two urns so that each contains NN balls. After each unit of time one ball is selected at random from each urn, and the two balls thus selected are interchanged. Let the number of black balls in the first urn denote the state of the system. Write down the transition matrix of this Markov chain and find the unique stationary distribution. Is the chain reversible in equilibrium?

?
Problem 6.15.13

Consider a Markov chain on the set S={0,1,2,…}S = \left\{ 0,1,2, \ldots \right\} with transition probabilities pi,i+1=aip_{i, i+1} = a_{i}, pi,0=1−ai,i≥0p_{i, 0} = 1-a_{i}, i \geq 0, where (ai:i≥0)\left(a_{i}: i \geq 0\right) is a sequence of constants which satisfy 0<ai<10 < a_{i} < 1 for all ii. Let b0=1,bi=a0a1⋯ai−1b_{0} = 1, b_{i} = a_{0} a_{1} \cdots a_{i-1} for i≥1i \geq 1. Show that the chain is

?
(a)

recurrent if and only if bi→0b_{i} \rightarrow 0 as i→∞i \rightarrow \infty,

(b)

positive recurrent if and only if ∑ibi<∞\sum_{i} b_{i} < \infty, and write down the stationary distribution if the latter condition holds. Let AA and β\beta be positive constants and suppose that ai=1−Ai−βa_{i} = 1-A i^{-\beta } for all large ii. Show that the chain is

(c)

transient if β>1\beta > 1,

(d)

positive recurrent if β<1\beta < 1.

Finally, if β=1\beta = 1 show that the chain is

(e)

positive recurrent if A>1A > 1,

(f)

null recurrent if A≤1A \leq 1.

Problem 6.15.14

Let XX be a continuous-time Markov chain with countable state space SS and semigroup {Pt}\left\{ \mathbf{P}_{t}\right\}. Show that pij(t)p_{i j}(t) is a continuous function of tt. Let g(t)=−log⁡pii(t)g(t) = -\log p_{i i}(t); show that gg is a continuous function, g(0)=0g(0) = 0, and g(s+t)≤g(s)+g(t)g(s+t) \leq g(s)+g(t). We say that gg is 'subadditive', and a well known theorem gives the result that

lim⁡t↓0g(t)t=λ exists and λ=sup⁡t>0g(t)t≤∞ \lim _{t \downarrow 0} \frac{g(t)}{t} = \lambda \quad \text{ exists and } \quad \lambda = \sup _{t > 0} \frac{g(t)}{t} \leq \infty

Deduce that the limit gii=lim⁡t↓0t−1{pii(t)−1}g_{i i} = \lim_{t \downarrow 0} t^{-1}\left\{ p_{i i}(t)-1\right\} exists.

?
Problem 6.15.15

Let XX be a continuous-time Markov chain with generator G=(gij)\mathbf{G} = \left(g_{i j}\right). Show that XX is irreducible if and only if for any pair i,ji, j of distinct states there exists a sequence i,k1,k2,…,kn,ji, k_{1}, k_{2}, \ldots , k_{n}, j of distinct states such that gi,k1gk1,k2⋯gkn,j>0g_{i, k_{1}} g_{k_{1}, k_{2}} \cdots g_{k_{n}, j} > 0.

?
Problem 6.15.16
?
(a)

Let T>0T > 0 and let X={X(t):0≤t≤T}X = \left\{ X(t): 0 \leq t \leq T\right\} be an irreducible, non-explosive Markov chain with stationary distribution π\pi, and suppose that X(0)X(0) has distribution π\pi. Let Y(t)=X(T−t)Y(t) = X(T-t) for 0≤t≤T0 \leq t \leq T. We call XX reversible (in equilibrium) if XX and YY have the same joint distributions.

(i) Show that YY is a (left-continuous) Markov chain with transition probabilities p^ij(t)=(πj/πi)pji(t)\widehat{p}_{i j}(t) = \left(\pi_{j} / \pi_{i}\right) p_{j i}(t) and generator G^\widehat{\mathbf{G}} satisfying πjg^ji=πigij\pi_{j} \hat{g}_{j i} = \pi_{i} g_{i j}, where the pji(t)p_{j i}(t) and G=(gij)\mathbf{G} = \left(g_{i j}\right) are those of XX. Show that YY is irreducible and non-explosive with stationary distribution π\pi.

(ii) Show that XX is reversible in equilibrium if and only if the detailed balance equations πigij=πjgji\pi_{i} g_{i j} = \pi_{j} g_{j i} (for all ii and jj) hold.

(iii) Show that a measure ν\mathbf{\nu } satisfies νG=0\mathbf{\nu } \mathbf{G} = \mathbf{0} if it satisfies the detailed balance equations.

(b)

Let XX be irreducible and non-explosive with stationary distribution π\pi, and assume X(0)X(0) has distribution π\pi.

(i) Kolmogorov’s criterion. Show that XX is reversible if and only if, for all nn and all finite sequences k1,k2,…,knk_{1}, k_{2}, \ldots , k_{n} of states,

gk1,k2gk2,k3⋯gkn−1,kngkn,k1=gk1,kngkn,kn−1⋯gk2,k1 g_{k_{1}, k_{2}} g_{k_{2}, k_{3}} \cdots g_{k_{n-1}, k_{n}} g_{k_{n}, k_{1}} = g_{k_{1}, k_{n}} g_{k_{n}, k_{n-1}} \cdots g_{k_{2}, k_{1}}

(ii) Kelly’s criterion. Show that XX is reversible if, for all distinct triples i,j,k∈Si, j, k \in S, we have gijgjkgki=gikgkjgjig_{i j} g_{j k} g_{k i} = g_{i k} g_{k j} g_{j i}, and in addition there exists c∈Sc \in S such that gic>0g_{i c} > 0 for all i≠ci \neq c.

(c)

Show that every irreducible chain XX with exactly two states is reversible in equilibrium.

(d)

Show that every non-explosive birth-death process XX having a stationary distribution is reversible in equilibrium.

Problem 6.15.17

Show that not every discrete-time Markov chain can be imbedded in a continuous-time chain. More precisely, let

P=[α1−α1−αα] for some 0<α<1 \mathbf{P} =\left[\begin{smallmatrix} \alpha & 1-\alpha \\ 1-\alpha & \alpha \end{smallmatrix}\right] \quad \text{ for some } 0 < \alpha < 1

be a transition matrix. Show that there exists a semigroup {Pt}\left\{ \mathbf{P}_{t}\right\} of transition probabilities in continuous time such that P1=P\mathbf{P}_{1} = \mathbf{P}, if and only if 12<α<1\frac{1}{2} < \alpha < 1. In this case show that {Pt}\left\{ \mathbf{P}_{t}\right\} is unique and calculate it in terms of α\alpha.

?
Problem 6.15.18

Consider an immigration-death process X(t)X(t), being a birth-death process with rates λn=λ\lambda_{n} = \lambda, μn=nμ\mu_{n} = n \mu. Show that its generating function G(s,t)=E[sX(t)]G(s, t) = \mathbb {E}\left[s^{X(t)}\right] is given by

G(s,t)={1+(s−1)e−μt}Iexp⁡{ρ(s−1)(1−e−μt)} G(s, t) = \left\{ 1+(s-1) e^{-\mu t}\right\} ^{I} \exp \left\{ \rho (s-1)\left(1-e^{-\mu t}\right)\right\}

where ρ=λ/μ\rho = \lambda / \mu and X(0)=IX(0) = I. Deduce the limiting distribution of X(t)X(t) as t→∞t \rightarrow \infty.

?
Problem 6.15.19

Let NN be a non-homogeneous Poisson process on R+=[0,∞)\mathbb {R}_{+} = [0, \infty ) with intensity function λ\lambda.

?
(a)

Write down the forward and backward equations for NN, and solve them.

(b)

Let N(0)=0N(0) = 0, and find the density function of the time TT until the first arrival in the process. If λ(t)=c/(1+t)\lambda (t) = c /(1+t), show that E[T]<∞\mathbb {E}\left[T\right] < \infty if and only if c>1c > 1.

Problem 6.15.20

Successive offers for my house are independent identically distributed random variables X1X_{1}, X2,…X_{2}, \ldots, having density function ff and distribution function FF. Let Y1=X1Y_{1} = X_{1}, let Y2Y_{2} be the first offer exceeding Y1Y_{1}, and generally let Yn+1Y_{n+1} be the first offer exceeding YnY_{n}. Show that Y1,Y2,…Y_{1}, Y_{2}, \ldots are the times of arrivals in a non-homogeneous Poisson process with intensity function λ(t)=f(t)/(1−F(t))\lambda (t) = f(t) /(1-F(t)). The YiY_{i} are called 'record values'.

Now let Z1Z_{1} be the first offer received which is the second largest to date, and let Z2Z_{2} be the second such offer, and so on. Show that the ZiZ_{i} are the arrival times of a non-homogeneous Poisson process with intensity function λ\lambda.

?
Problem 6.15.21

Let NN be a Poisson process with constant intensity λ\lambda, and let Y1,Y2,…Y_{1}, Y_{2}, \ldots be independent random variables with common characteristic function ϕ\phi and density function ff. The process N∗(t)=Y1+Y2+⋯+YN(t)N^{*}(t) = Y_{1}+Y_{2}+\cdots +Y_{N(t)} is called a compound Poisson process. YnY_{n} is the change in the value of N∗N^{*} at the nnth arrival of the Poisson process NN. Think of it like this. A 'random alarm clock' rings at the arrival times of a Poisson process. At the nnth ring the process N∗N^{*} accumulates an extra quantity YnY_{n}. Write down a forward equation for N∗N^{*} and hence find the characteristic function of N∗(t)N^{*}(t). Can you see directly why it has the form which you have found?

?
Problem 6.15.22

If the intensity function λ\lambda of a non-homogeneous Poisson process NN is itself a random process, then NN is called a doubly stochastic Poisson process (or Cox process).

?
(a)

Consider the case when λ(t)=Λ\lambda (t) = \Lambda for all tt, and Λ\Lambda is a random variable taking either of two values λ1\lambda_{1} or λ2\lambda_{2}, each being picked with equal probability 12\frac{1}{2}. Find the probability generating function of N(t)N(t), and deduce its mean and variance.

(b)

For a doubly stochastic Poisson process NN, show that Var⁡[N(t)]≥E[N(t)]\operatorname {Var}\left[N(t)\right] \geq \mathbb {E}\left[N(t)\right].

(c)

Let MM be an ordinary Poisson process on the time-interval [0,∞)[0, \infty ) with constant rate 1. Let M∗M^{*} be obtained from MM by deleting the kk th arrival for every odd value of kk. Is M∗M^{*} either: (i) a Poisson process, or (ii) a doubly stochastic Poisson process?

Problem 6.15.23

Show that a simple birth process XX with parameter λ\lambda is a doubly stochastic Poisson process with intensity function λ(t)=λX(t)\lambda (t) = \lambda X(t).

?
Problem 6.15.24

The Markov chain X={X(t):t≥0}X = \left\{ X(t): t \geq 0\right\} is a birth process whose intensities λk(t)\lambda_{k}(t) depend also on the time tt and are given by

P(X(t+h)=k+1∣X(t)=k)=1+μk1+μth+o(h) \mathbb {P}\left(X(t+h) = k+1 \mid X(t) = k\right) = \frac{1+\mu k}{1+\mu t} h+\mathrm{o}(h)

as h↓0h \downarrow 0. Show that the probability generating function G(s,t)=E[sX(t)]G(s, t) = \mathbb {E}\left[s^{X(t)}\right] satisfies

∂G∂t=s−11+μt{G+μs∂G∂s},0<s<1 \frac{\partial G}{\partial t} = \frac{s-1}{1+\mu t}\left\{ G+\mu s \frac{\partial G}{\partial s}\right\} , \quad 0 < s < 1

Hence find the mean and variance of X(t)X(t) when X(0)=IX(0) = I.

?
Problem 6.15.25
?
(a)

Let XX be a birth-death process with strictly positive birth rates λ0,λ1,…\lambda_{0}, \lambda_{1}, \ldots and death rates μ1,μ2,…\mu_{1}, \mu_{2}, \ldots Let ηi\eta_{i} be the probability that X(t)X(t) ever takes the value 0 starting from X(0)=iX(0) = i. Show that

λjηj+1−(λj+μj)ηj+μjηj−1=0,j≥1 \lambda _{j} \eta _{j+1}-\left(\lambda _{j}+\mu _{j}\right) \eta _{j}+\mu _{j} \eta _{j-1} = 0, \quad j \geq 1

and deduce that ηi=1\eta_{i} = 1 for all ii so long as ∑1∞ej=∞\sum_{1}^{\infty } e_{j} = \infty where ej=μ1μ2⋯μj/(λ1λ2⋯λj)e_{j} = \mu_{1} \mu_{2} \cdots \mu_{j} /\left(\lambda_{1} \lambda_{2} \cdots \lambda_{j}\right).

(b)

For the discrete-time chain on the non-negative integers with

pj,j+1=(j+1)2j2+(j+1)2 and pj,j−1=j2j2+(j+1)2 p_{j, j+1} = \frac{(j+1)^{2}}{j^{2}+(j+1)^{2}} \quad \text{ and } \quad p_{j, j-1} = \frac{j^{2}}{j^{2}+(j+1)^{2}}

find the probability that the chain ever visits 00, starting from 1.

Problem 6.15.26

Find a good necessary condition and a good sufficient condition for the birth-death process XX of Problem (6.15.25a) to be honest.

?
Problem 6.15.27

Let XX be a simple symmetric birth-death process with λn=μn=nλ\lambda_{n} = \mu_{n} = n \lambda, and let TT be the time until extinction. Show that

P(T≤x∣X(0)=I)=(λx1+λx)I \mathbb {P}\left(T \leq x \mid X(0) = I\right) = \left(\frac{\lambda x}{1+\lambda x}\right)^{I}

and deduce that extinction is certain if P(X(0)<∞)=1\mathbb {P}\left(X(0) < \infty \right) = 1. Show that P(λT/I≤x∣X(0)=I)→e−1/x\mathbb {P}\left(\lambda T / I \leq x \mid X(0) = I\right) \rightarrow e^{-1 / x} as I→∞I \rightarrow \infty.

?
Problem 6.15.28

Let XX be an immigration-death-disaster process, that is, a birth-death process with parameters λi=λ,μi=iμ\lambda_{i} = \lambda , \mu_{i} = i \mu, and with the additional possibility of 'disasters' which reduce the population to 0. Disasters occur at the times of a Poisson process with intensity δ\delta, independently of all previous births and deaths.

?
(a)

Show that XX has a stationary distribution, and find an expression for the generating function of this distribution.

(b)

Show that, in equilibrium, the mean of X(t)X(t) is λ/(δ+μ)\lambda /(\delta +\mu ).

Problem 6.15.29

With any sufficiently nice (Lebesgue measurable, say) subset BB of the real line R\mathbb {R} is associated a random variable X(B)X(B) such that

(a) X(B)X(B) takes values in {0,1,2,…}\left\{ 0,1,2, \ldots \right\},

(b) if B1,B2,…,BnB_{1}, B_{2}, \ldots , B_{n} are disjoint then X(B1),X(B2),…,X(Bn)X\left(B_{1}\right), X\left(B_{2}\right), \ldots , X\left(B_{n}\right) are independent, and furthermore X(B1∪B2)=X(B1)+X(B2)X\left(B_{1} \cup B_{2}\right) = X\left(B_{1}\right)+X\left(B_{2}\right)

(c) the distribution of X(B)X(B) depends only on BB through its Lebesgue measure (’length’) ∣B∣\left|B\right|, and

P(X(B)≥1)P(X(B)=1)→1 as ∣B∣→0 \frac{\mathbb {P}\left(X(B) \geq 1\right)}{\mathbb {P}\left(X(B) = 1\right)} \rightarrow 1 \quad \text{ as }|B| \rightarrow 0

Show that XX is a Poisson process.

?
Problem 6.15.30

Let NN be a Poisson process in R2\mathbb {R}^{2} with constant intensity λ\lambda, and let R(1)<R(2)<R_{(1)} < R_{(2)} <... be the ordered distances from the origin of the points of the process.

?
(a)

Show that R(1)2,R(2)2,…R_{(1)}^{2}, R_{(2)}^{2}, \ldots are the points of a Poisson process on R+=[0,∞)\mathbb {R}_{+} = [0, \infty ) with intensity λπ\lambda \pi.

(b)

Show that R(k)R_{(k)} has density function

f(r)=2πλr(λπr2)k−1e−λπr2(k−1)!,r>0 f(r) = \frac{2 \pi \lambda r\left(\lambda \pi r^{2}\right)^{k-1} e^{-\lambda \pi r^{2}}}{(k-1)!}, \quad r > 0
Problem 6.15.31

Let XX be a nn-dimensional Poisson process with constant intensity λ\lambda. Show that the volume of the largest (nn-dimensional) sphere centred at the origin which contains no point of XX is exponentially distributed. Deduce the density function of the distance RR from the origin to the nearest point of XX. Show that E[R]=Γ(1/n)/{n(λc)1/n}\mathbb {E}\left[R\right] = \Gamma (1 / n) /\left\{ n(\lambda c)^{1 / n}\right\} where cc is the volume of the unit ball of Rn\mathbb {R}^{n} and Γ\Gamma is the gamma function.

?
Problem 6.15.32

A village of N+1N+1 people suffers an epidemic. Let X(t)X(t) be the number of ill people at time tt, and suppose that X(0)=1X(0) = 1 and XX is a birth process with rates λi=λi(N+1−i)\lambda_{i} = \lambda i(N+1-i). Let TT be the length of time required until every member of the population has succumbed to the illness. Show that

E[T]=1λ∑k=1N1k(N+1−k) \mathbb {E}\left[T\right] = \frac{1}{\lambda } \sum _{k = 1}^{N} \frac{1}{k(N+1-k)}

and deduce that

E[T]=2(log⁡N+γ)λ(N+1)+O(N−2) \mathbb {E}\left[T\right] = \frac{2(\log N+\gamma )}{\lambda (N+1)}+\mathrm{O}\left(N^{-2}\right)

where γ\gamma is Euler's constant. It is striking that E[T]\mathbb {E}\left[T\right] decreases with NN, for large NN.

?
Problem 6.15.33

A particle has velocity V(t)V(t) at time tt, where V(t)V(t) is assumed to take values in {n+12:n≥0}\left\{ n+\frac{1}{2}: n \geq 0\right\}. Transitions during (t,t+h)(t, t+h) are possible as follows:

P(V(t+h)=w∣V(t)=v)={(v+12)h+o(h) if w=v+11−2vh+o(h) if w=v(v−12)h+o(h) if w=v−1 \mathbb {P}\left(V(t+h) = w \mid V(t) = v\right) = \begin{cases} \left(v+\frac{1}{2}\right) h+\mathrm{o}(h) & \text{ if } w = v+1 \\ 1-2 v h+\mathrm{o}(h) & \text{ if } w = v \\ \left(v-\frac{1}{2}\right) h+\mathrm{o}(h) & \text{ if } w = v-1\end{cases}

Initially V(0)=12V(0) = \frac{1}{2}. Let

G(s,t)=∑n=0∞snP(V(t)=n+12) G(s, t) = \sum _{n = 0}^{\infty } s^{n} \mathbb {P}\left(V(t) = n+\frac{1}{2}\right)
?
(a)

Show that

∂G∂t=(1−s)2∂G∂s−(1−s)G \frac{\partial G}{\partial t} = (1-s)^{2} \frac{\partial G}{\partial s}-(1-s) G

and deduce that G(s,t)={1+(1−s)t}−1G(s, t) = \left\{ 1+(1-s) t\right\}^{-1}.

(b)

Show that the expected length mn(T)m_{n}(T) of time for which V=n+12V = n+\frac{1}{2} during the time interval [0,T][0, T] is given by

mn(T)=∫0TP(V(t)=n+12)dt m_{n}(T) = \int _{0}^{T} \mathbb {P}\left(V(t) = n+\frac{1}{2}\right) d t

and that, for fixed k,mk(T)−log⁡T→−∑i=1ki−1k, m_{k}(T)-\log T \rightarrow -\sum_{i = 1}^{k} i^{-1} as T→∞T \rightarrow \infty.

(c)

What is the expected velocity of the particle at time tt?

Problem 6.15.34

A sequence X0,X1,…X_{0}, X_{1}, \ldots of random integers is generated as follows. First, X0=0,X1=1X_{0} = 0, X_{1} = 1. For n≥1n \geq 1, conditional on X0,X1,…,XnX_{0}, X_{1}, \ldots , X_{n}, the next value Xn+1X_{n+1} is equally likely to be either Xn+Xn−1X_{n}+X_{n-1} or ∣Xn−Xn−1∣\left|X_{n}-X_{n-1}\right|.

?
(a)

Is XX a Markov chain?

(b)

Use the Markov chain Yn=(Xn−1,Xn)Y_{n} = \left(X_{n-1}, X_{n}\right) to find the probability that XX hits the value 3 before it revisits 0.

(c)

Show that the probability that YY ever reaches the state (1,1)(1,1), having started at (1,2)(1,2), is 12(3−5)\frac{1}{2}(3-\sqrt{5}).

Problem 6.15.35

Take a regular hexagon and join opposite corners by straight lines meeting at the point C. A particle performs a symmetric random walk on these 7 vertices, starting at A(≠C)\mathrm{A}( \neq \mathrm{C}). Find:

?
(a)

the probability of return to A without hitting C,

(b)

the expected time to return to A,

(c)

the expected number of visits to C before returning to A,

(d)

the expected time to return to A, given that there is no prior visit to C.

Problem 6.15.36

Markov chains are defined by the following procedures at any time nn :

?
(a)

Bernoulli model. Two adjacent containers A and B each contain mm particles; mm are of type I and mm are of type II. A particle is selected at random in each container. If they are of opposite types they are exchanged with probability α\alpha if the type I is in A, or with probability β\beta if the type I is in B. Let XnX_{n} be the number of type I particles in A at time nn.

(b)

Ehrenfest dog-flea model. Two adjacent containers contain mm particles in all. A particle is selected at random. If it is in A it is moved to B with probability α\alpha, if it is in B it is moved to A with probability β\beta. Let YnY_{n} be the number of particles in A at time nn. In each case find the transition matrix and stationary distribution of the chain.

Problem 6.15.37

Let XX be an irreducible continuous-time Markov chain on the state space SS with transition probabilities pjk(t)p_{j k}(t) and unique stationary distribution π\mathbf{\pi }, and write P(X(t)=j)=aj(t)\mathbb {P}\left(X(t) = j\right) = a_{j}(t). If c(x)c(x) is a concave function, show that the function d(t)=∑j∈Sπjc(aj(t)/πj)d(t) = \sum_{j \in S} \pi_{j} c\left(a_{j}(t) / \pi_{j}\right) increases to c(1)c(1) as t→∞t \rightarrow \infty.

The relative entropy (or Kullback-Leibler divergence) of two strictly positive probability mass functions f,gf, g on a subset SS of integers is defined as

D(f;g)=∑i∈Sf(i)log⁡(f(i)/g(i)) D(f ; g) = \sum _{i \in S} f(i) \log (f(i) / g(i))

Prove that, if XX has a finite state space and stationary distribution π\mathbf{\pi }, the relative entropy D(a(t);π)D(a(t) ; \mathbf{\pi }) decreases monotonely to 0 as t→∞t \rightarrow \infty.

?
Problem 6.15.38

With the notation of the preceding problem, let uk(t)=P(X(t)=k∣X(0)=0)u_{k}(t) = \mathbb {P}\left(X(t) = k \mid X(0) = 0\right), and suppose the chain is reversible in equilibrium (see Problem (6.15.16)). Show that u0(2t)=∑j(π0/πj)uj(t)2u_{0}(2 t) = \sum_{j}\left(\pi_{0} / \pi_{j}\right) u_{j}(t)^{2}, and deduce that u0(t)u_{0}(t) decreases to π0\pi_{0} as t→∞t \rightarrow \infty.

?
Problem 6.15.39

Let Π\Pi be the set of points in a Poisson process on Rd\mathbb {R}^{d} with constant intensity λ\lambda. Each point is displaced, where the displacements are independent and identically distributed. Show that the resulting point process is a Poisson process with intensity λ\lambda.

?
Problem 6.15.40

Suppose for convenience in Problem (6.15.39) that the displacements have a continuous distribution function and finite mean, and that d=1d = 1. Suppose also that you are at the origin originally, and you move to aa in the perturbed process. Let LR\mathrm{L}_{\mathrm{R}} be the number of points formerly on your left that are now on your right, and RL\mathrm{R}_{\mathrm{L}} the number of points formerly on your right that are now on your left. Show that E[LR]=E[RL]\mathbb {E}\left[\mathrm{L}_{\mathrm{R}}\right] = \mathbb {E}\left[\mathrm{R}_{\mathrm{L}}\right] if and only if a=μa = \mu where μ\mu is the mean displacement of a particle.

Deduce that if cars enter the start of a long road at the instants of a Poisson process, having independent identically distributed velocities, then, if you travel at the average speed, in the long run the rate at which you are overtaken by other cars equals the rate at which you overtake other cars.

?
Problem 6.15.41

Ants enter a kitchen at the instants of a Poisson process NN of rate λ\lambda; they each visit the pantry and then the sink, and leave. The rr th ant spends time XrX_{r} in the pantry and YrY_{r} in the sink (and Xr+YrX_{r}+Y_{r} in the kitchen altogether), where the vectors Vr=(Xr,Yr)V_{r} = \left(X_{r}, Y_{r}\right) and VsV_{s} are independent for r≠sr \neq s. At time t=0t = 0 the kitchen is free of ants. Find the joint distribution of the numbers A(t)A(t) of ants in the pantry and B(t)B(t) of ants in the sink at time tt.

Show that, as t→∞t \rightarrow \infty, the number of ants in the kitchen converges in distribution, provided E[Xr+Yr]<∞\mathbb {E}\left[X_{r}+Y_{r}\right] < \infty.

Now suppose the ants arrive in pairs at the times of the Poisson process, but then separate to behave independently as above. Find the joint distribution of the numbers of ants in the two locations.

?
Problem 6.15.42

Let {Xr:r≥1}\left\{ X_{r}: r \geq 1\right\} be independent exponential random variables with parameter λ\lambda, and set Sn=∑r=1nXrS_{n} = \sum_{r = 1}^{n} X_{r}. Show that:

?
(a)

Yk=Sk/Sn,1≤k≤n−1Y_{k} = S_{k} / S_{n}, 1 \leq k \leq n-1, have the same distribution as the order statistics of independent variables {Uk:1≤k≤n−1}\left\{ U_{k}: 1 \leq k \leq n-1\right\} which are uniformly distributed on (0,1)(0,1),

(b)

Zk=Xk/Sn,1≤k≤nZ_{k} = X_{k} / S_{n}, 1 \leq k \leq n, have the same joint distribution as the coordinates of a point (U1,…,Un)\left(U_{1}, \ldots , U_{n}\right) chosen uniformly at random on the simplex ∑r=1nur=1,ur≥0\sum_{r = 1}^{n} u_{r} = 1, u_{r} \geq 0 for all rr.

Problem 6.15.43

Let XX be a discrete-time Markov chain with a finite number of states and transition matrix P=(pij)\mathbf{P} = \left(p_{i j}\right) where pij>0p_{i j} > 0 for all i,ji, j. Show that there exists λ∈(0,1)\lambda \in (0,1) such that ∣pij(n)−πj∣<λn\left|p_{i j}(n)-\pi_{j}\right| < \lambda^{n}, where π\mathbf{\pi } is the stationary distribution.

?
Problem 6.15.44

Under the conditions of Problem (6.15.43), let Vi(n)=∑r=0n−1I{Xr=i}V_{i}(n) = \sum_{r = 0}^{n-1} I_{\left\{ X_{r} = i\right\} } be the number of visits of the chain to ii before time nn. Show that

E[∣1nVi(n)−πi∣2]→0 as n→∞ \mathbb {E}\left[\left|\frac{1}{n} V_{i}(n)-\pi _{i}\right|^{2}\right] \rightarrow 0 \quad \text{ as } n \rightarrow \infty

Show further that, if ff is any bounded function on the state space, then

E[∣1n∑r=0n−1f(Xr)−∑i∈Sf(i)πi∣2]→0 \mathbb {E}\left[\left|\frac{1}{n} \sum _{r = 0}^{n-1} f\left(X_{r}\right)-\sum _{i \in S} f(i) \pi _{i}\right|^{2}\right] \rightarrow 0
?
Problem 6.15.45

Let AA and B=(B0,B1,…,Bn)\mathbf{B} = \left(B_{0}, B_{1}, \ldots , B_{n}\right) be a discrete random variable and vector, respectively. The conditional entropy of AA with respect to B\mathbf{B} is defined as H(A∣B)=E[E[−log⁡f(A∣B)∣B]]H(A \mid \mathbf{B}) = \mathbb {E}\left[\mathbb {E}\left[-\log f(A \mid \mathbf{B}) \mid \mathbf{B}\right]\right] where f(a∣b)=P(A=a∣B=b)f(a \mid \mathbf{b}) = \mathbb {P}\left(A = a \mid \mathbf{B} = \mathbf{b}\right). Let XX be an aperiodic Markov chain on a finite state space. Show that

H(Xn+1∣X0,X1…,Xn)=H(Xn+1∣Xn) H\left(X_{n+1} \mid X_{0}, X_{1} \ldots , X_{n}\right) = H\left(X_{n+1} \mid X_{n}\right)

and that

H(Xn+1∣Xn)→−∑iπi∑jpijlog⁡pij as n→∞ H\left(X_{n+1} \mid X_{n}\right) \rightarrow -\sum _{i} \pi _{i} \sum _{j} p_{i j} \log p_{i j} \quad \text{ as } n \rightarrow \infty

if XX is aperiodic with a unique stationary distribution π\pi.

?
Problem 6.15.46

Let XX and YY be independent recurrent birth-death processes with the same parameters (and no explosions). It is not assumed that X0=Y0X_{0} = Y_{0}. Show that:

?
(a)

for any A⊆R,∣P(Xt∈A)−P(Yt∈A)∣→0A \subseteq \mathbb {R},\left|\mathbb {P}\left(X_{t} \in A\right)-\mathbb {P}\left(Y_{t} \in A\right)\right| \rightarrow 0 as t→∞t \rightarrow \infty,

(b)

if P(X0≤Y0)=1\mathbb {P}\left(X_{0} \leq Y_{0}\right) = 1, then E[g(Xt)]≤E[g(Yt)]\mathbb {E}\left[g\left(X_{t}\right)\right] \leq \mathbb {E}\left[g\left(Y_{t}\right)\right] for any increasing function gg.

Problem 6.15.47

The number of birds in a wood at time tt is a continuous-time Markov process XX. Food resources impose the constraint 0≤X(t)≤n0 \leq X(t) \leq n. Competition entails that the transition probabilities obey

pk,k+1(h)=λ(n−k)h+o(h),pk,k−1(h)=μkh+o(h) p_{k, k+1}(h) = \lambda (n-k) h+\mathrm{o}(h), \quad p_{k, k-1}(h) = \mu k h+\mathrm{o}(h)

Find E[sX(t)]\mathbb {E}\left[s^{X(t)}\right], together with the mean and variance of X(t)X(t), when X(0)=rX(0) = r. What happens as t→∞t \rightarrow \infty?

?
Problem 6.15.48

A counter performs an irreducible random walk on the vertices 0,1,20,1,2 of the triangle in the figure beneath, with transition matrix

P=[0p0q0q10p1p2q20] \mathbf{P} =\left[\begin{smallmatrix} 0 & p_{0} & q_{0} \\ q_{1} & 0 & p_{1} \\ p_{2} & q_{2} & 0 \end{smallmatrix}\right]

where pi+qi=1p_{i}+q_{i} = 1 for all ii. Show that the stationary distribution π\pi has

π0=1−q2p13−q1p0−q2p1−q0p2 \pi _{0} = \frac{1-q_{2} p_{1}}{3-q_{1} p_{0}-q_{2} p_{1}-q_{0} p_{2}}

with corresponding formulae for π1,π2\pi_{1}, \pi_{2}.

Suppose that you gain one peseta for each clockwise step of the walk, and you lose one peseta for each anticlockwise step. Show that, in equilibrium, the mean yield per step is

γ=∑i(2pi−1)πi=3(2p0p1p2−p0p1−p1p2−p2p0+p0+p1+p2−1)3−q1p0−q2p1−q0p2 \gamma = \sum _{i}\left(2 p_{i}-1\right) \pi _{i} = \frac{3\left(2 p_{0} p_{1} p_{2}-p_{0} p_{1}-p_{1} p_{2}-p_{2} p_{0}+p_{0}+p_{1}+p_{2}-1\right)}{3-q_{1} p_{0}-q_{2} p_{1}-q_{0} p_{2}}

Consider now three cases of this process: A. We have pi=12−ap_{i} = \frac{1}{2}-a for each ii, where a>0a > 0. Show that the mean yield per step satisfies γA<0\gamma_{\mathrm{A}} < 0. B. We have that p0=110−a,p1=p2=34−ap_{0} = \frac{1}{10}-a, p_{1} = p_{2} = \frac{3}{4}-a, where a>0a > 0. Show that γB<0\gamma_{\mathrm{B}} < 0 for sufficiently small aa. C. At each step the counter is equally likely to move according to the transition probabilities of case AA or case BB, the choice being made independently at every step. Show that, in this case, p0=310−a,p1=p2=58−ap_{0} = \frac{3}{10}-a, p_{1} = p_{2} = \frac{5}{8}-a. Show that γC>0\gamma_{\mathrm{C}} > 0 for sufficiently small aa. The fact that two systematically unfavourable games may be combined to make a favourable game is called Parrondo's paradox. Such bets are not available in casinos.

?
Problem 6.15.49

Cars arrive at the beginning of a long road in a Poisson stream of rate λ\lambda from time t=0t = 0 onwards. A car has a fixed velocity V>0V > 0 which is a random variable. The velocities of cars are independent and identically distributed, and independent of the arrival process. Cars can overtake each other freely. Show that the number of cars on the first xx miles of the road at time tt has the Poisson distribution with parameter λE[V−1min⁡{x,Vt}]\lambda \mathbb {E}\left[V^{-1} \min \left\{ x, V t\right\} \right].

?
Problem 6.15.50

Events occur at the times of a Poisson process with intensity λ\lambda, and you are offered a bet based on the process. Let t>0t > 0. You are required to say the word 'now' immediately after the event which you think will be the last to occur prior to time tt. You win if you succeed, otherwise you lose. If no events occur before tt you lose. If you have not selected an event before time tt you lose.

Consider the strategy in which you choose the first event to occur after a specified time ss, where 0<s<t0 < s < t.

?
(a)

Calculate an expression for the probability that you win using this strategy.

(b)

Which value of ss maximizes this probability?

(c)

If λt≥1\lambda t \geq 1, show that the probability that you win using this value of ss is e−1e^{-1}.

Problem 6.15.51

A new Oxbridge professor wishes to buy a house, and can afford to spend up to one million pounds. Declining the services of conventional estate agents, she consults her favourite internet property page on which houses are announced at the times of a Poisson process with intensity λ\lambda per day. House prices may be assumed to be independent random variables which are uniformly distributed over the interval (800,000,2,000,000)(800,000,2,000,000). She decides to view every affordable property announced during the next 30 days. The time spent viewing any given property is uniformly distributed over the range (1,2)(1,2) hours. What is the moment generating function of the total time spent viewing houses?

?
Problem 6.15.52

Let X={Xn:n≥0}X = \left\{ X_{n}: n \geq 0\right\} be an irreducible, aperiodic Markov chain on the finite state space SS, and let hij=Ei[min⁡{n≥0:Xn=j}]h_{i j} = \mathbb {E}_{i}\left[\min \left\{ n \geq 0: X_{n} = j\right\} \right] denote the mean hitting time of jj (noting that hii=0h_{i i} = 0). Let Ki=∑jhijπjK_{i} = \sum_{j} h_{i j} \pi_{j} be the mean time to hit a state ZZ chosen at random according to the stationary distribution π\pi. Show that KiK_{i} is independent of the choice of ii.

?
Problem 6.15.53

A professor travels between home and work on foot. She possesses a total of rr umbrellas, being those at home and at work. If it is raining when she leaves either home or work, she takes an umbrella with her (if there is one available). Assume that it is raining at the start of any given walk with probability pp (subject to the usual independence). Let XnX_{n} be the number of umbrellas available to her at the start of her nnth walk.

?
(a)

Explain why XX is a Markov chain, and write down its transition matrix.

(b)

Show that the chain has stationary distribution π\pi given by

πi={1−pr+1−p if i=01r+1−p if i=1,2,…,r \pi _{i} = \begin{cases} \frac{1-p}{r+1-p} & \text{ if } i = 0 \\ \frac{1}{r+1-p} & \text{ if } i = 1,2, \ldots , r\end{cases}

What proportion of walks result in her getting wet, in the long run?

(c)

Let r=1r = 1 and X1=1X_{1} = 1. Calculate the mean number of walks made before she gets wet.

Problem 6.15.54

A spider climbs a vertical spout of height hh at speed 1. At the instants of a Poisson process of constant rate λ\lambda, the spider is flushed back to the bottom of the spout. It then recommences its climb. Let TT be the time to reach the top, and NN the number of intermediate flushes. Show that

E[e−θTsN]=(λ+θ)e−(λ+θ)hλ+θ−λs(1−e−(λ+θ)h),θ,s∈R \mathbb {E}\left[e^{-\theta T} s^{N}\right] = \frac{(\lambda +\theta ) e^{-(\lambda +\theta ) h}}{\lambda +\theta -\lambda s\left(1-e^{-(\lambda +\theta ) h}\right)}, \quad \theta , s \in \mathbb {R}

By calculating E[e−θT∣N=n]\mathbb {E}\left[e^{-\theta T} \mid N = n\right] or otherwise, determine E[T∣N=n]\mathbb {E}\left[T \mid N = n\right] for n≥0n \geq 0.

?
Problem 6.15.55

Let XX be an ergodic Markov chain with transition matrix P\mathbf{P} and stationary distribution π\pi. Show for any set AA of states that

∑i∈Aj∉Aπipij=∑i∉Aj∈Aπipij \sum _{\substack {i \in A \\ j \notin A}} \pi _{i} p_{i j} = \sum _{\substack {i \notin A \\ j \in A}} \pi _{i} p_{i j}
?
Problem 6.15.56

A hiker of unit mass stands at the origin of the plain R2\mathbb {R}^{2}. Boulders with independent, identically distributed masses M1,M2,…M_{1}, M_{2}, \ldots lie scattered on the plain at the points of a Poisson process with intensity 1. Let GRG_{R} be the xx-component of the gravitational attraction on the hiker of the boulders within distance RR of the hiker. You may take the gravitational constant to be 1.

Show that, as R→∞,GRR \rightarrow \infty , G_{R} converges in distribution to a Cauchy distribution with characteristic function of the form ϕ(t)=e−c∣t∣\phi (t) = e^{-c\left|t\right|}, and express cc in terms of a typical mass MM.

?
Problem 6.15.57

Holtsmark distribution for stellar gravity. Let stars of common mass mm be positioned at the points of a Poisson process with intensity 1 in R3\mathbb {R}^{3}. Let GRG_{R} be the xx-component of the gravitational attraction due to the stars within distance RR of the origin, upon a hitchhiker of unit mass at the origin. You may take the gravitational constant to be 1.

?
(a)

Show that, as R→∞,GRR \rightarrow \infty , G_{R} converges in distribution to a symmetric distribution with characteristic function ϕ(t)=exp⁡{−c∣t∣3/2}\phi (t) = \exp \left\{ -c|t|^{3 / 2}\right\} where c>0c > 0.

(b)

What is the answer if the stars have independent, identically distributed random masses MiM_{i}?