Solution CMC1
Solution CMC1
Solution CMC1
Question 1
(a) Let Wi be the time taken for atom i to decay. Then X1 = min{W1 , W2 , . . . , W1026 }. Since
Wi ∼ Exp(10−26 ) and the Wi ’s are independent, X1 ∼ Exp(1). To see that note that
After the first atom decays there are 1026 − 1 left. By the lack of memory property each decays
after a further time which is Exp(10−26 ). Thus
(b) Let p = P(Wi < 60) = 1 − exp{−60 × 10−26 }, so that N ∼ Bin(1026 , p). Recall that
the Poisson distribution with parameter λ = np approximates the binomial distribution B(n, p)
if n is sufficiently large and p is sufficiently small. In our cases, since 1026 is enormous and
p very small, the distribution of N is approximately Poisson with mean λ = 1026 p. But p =
1 − (1 − 60 × 10−26 + 12 (60 × 10−26 )2 + · · · ) ≈ 60 × 10−26 . Thus λ ≈ 60.
Question 2
(a) Suppose that
ith
1 if the decay is recorded
Wi =
0 if the ith decay isn’t recorded
PDt
Then Nt = i=1 Wi . Since the Wi ’s are independent and Dt is independent of the Wi ’s we have
PDt PDt
GNt (s) = E(sNt ) = E s i=1 Wi = E E s i=1 Wi Dt
Dt
!! Dt
!
Y Y
Wi Wi
= E E s Dt =E E s
i=1 i=1
Dt
!
Y
= E (GW (s))Dt = GDt (GW (s)).
= E GWi (s)
i=1
1
Since Dt ∼ Pois(µt), GDt (s) = eµt(s−1) and Gw (s) = (1 − p) + ps it follows
i.e. Nt ∼ Pois(µpt).
Alternatively, you could argue that this is nothing else than the thinning of a Poisson process!
(b) P(N1 = 1) = 3e−3 . Now
where we have used independent increments. Substituting the relevant probability mass functions
(using Nt ∼ P oi(3t)) completes the exercise.
Question 3
Applying the forward equations, it is clear that, for n ≥ 1
In addition
n n
Y
X λj
−λn pn (t) = − λi e−λi t (2)
i=0
λj − λi
j=0,j̸=i
n−1 n−1
X Y λj
λn−1 pn−1 (t) = λi e−λi t . (3)
i=0
λj − λi
j=0,j̸=i
2
The summation is equal to
n−1 n−1 n−1 n
λ2i
Y
X Y λj 1 X 2 −λi t λj
− e−λi t =− λi e
i=0
λn − λi λj − λi λn i=0 λj − λi
j=0,j̸=i j=0,j̸=i
n n
Y
1 X 2 −λi t λj
=− λi e
λn i=0 λj − λi
j=0,j̸=i
Question 4
The forward equations can be employed, as in lectures to obtain:
p′n+1 (t) + (n + 1)λpn+1 (t) = nλpn (t) = nλ(1 − e−λt )n−1 e−λt .
3
Note that
Z t
λn(1 − e−λu )n−1 e−λu eλ(n+1)u du
0
Z t
= λn(1 − e−λu )n−1 eλnu du
0
t n−1
X
n−1
Z
= λn (−1)k (e−λu )k eλnu du
0 k
k=0
Z t n−1
X n − 1
= λn (−1)k eλ(n−k)u du
0 k
k=0
n−1
X n − 1 Z t
= (−1)k n λeλ(n−k)u du
k 0
k=0
n−1
X n−1
n
(−1)k eλnt e−λkt − 1 .
=
k n−k
k=0
Thus
pn+1 (t) = eλnt (1 − e−λt )n e−λ(n+1)t = (1 − e−λt )n e−λt ,
which concludes the proof.
Alternative proof:
Alternatively, you could show that pn satisfies the forward equations. For n = 1 we have
p1 (t) = (1 − e−λt )1−1 e−λt = e−λt .
Then
p′1 (t) = −λe−λt = −λp1 (t),
hence p1 satisfies the forward equation. For n ≥ 2 we need to check whether
p′n (t) = −nλpn (t) + (n − 1)λpn−1 (t)
is satisfied.
Note that for
pn (t) = (1 − e−λt )n−1 e−λt ,
we have
p′n (t) = (n − 1)(1 − e−λt )n−2 λe−λt e−λt + (1 − e−λt )n−1 (−λ)e−λt
4
Question 5
(a) If there are i infectives there are a − i susceptibles and hence the birth rate is given by
λi = (a − i)iη if i = 1, . . . , a and 0 otherwise. (b)Let Xi be the time spent in state i then we have
that the time to complete epidemic is
T = X1 + · · · + Xa−1
Pa−1 a−1
Y
E[e−sT ] = E[e−s i=1 Xi
]= E[e−sXi ]
i=1
a−1
Y a−1
λi Y (a − i)iη
= = .
i=1
λi + s i=1
s + (a − i)iη
To compute the expectation, we calculate the log of the expectation and use the fact that
∂ −sT
E[T ] = − log E[e ] .
∂s s=0