next up previous
Postscript version of this file

Convolution: If X and Y independent rvs with densities f and g respectively and Z=X+Y then

displaymath188

Differentiating wrt z we get

displaymath192

This integral is called the convolution of densities f and g.

If tex2html_wrap_inline198 iid Exponential tex2html_wrap_inline200 then tex2html_wrap_inline202 has a Gamma tex2html_wrap_inline204 distribution. Density of tex2html_wrap_inline206 is

displaymath208

for s> 0.

Proof:

align14

Then

align19

This telescopes to

displaymath208

Properties of Poisson Processes

1)
If tex2html_wrap_inline214 and tex2html_wrap_inline216 are independent Poisson processes with rates tex2html_wrap_inline218 and tex2html_wrap_inline220 , respectively the tex2html_wrap_inline222 is a Poisson processes with rates tex2html_wrap_inline224 .

2)
Suppose N is a Poisson process with rate tex2html_wrap_inline228 . Suppose each point is marked with a label, say one of tex2html_wrap_inline230 independently of all other occurences. Suppose tex2html_wrap_inline232 is the probability that a given point receives label tex2html_wrap_inline234 . Let tex2html_wrap_inline236 count the points with label i (so that tex2html_wrap_inline240 ). Then tex2html_wrap_inline242 are independent Poisson processes with rates tex2html_wrap_inline244 .

3)
Suppose tex2html_wrap_inline246 independent rvs, each uniformly distributed on [0,T]. Suppose M is a Poisson tex2html_wrap_inline252 random variable independent of the U's. Let

displaymath256

Then N is a Poisson process on [0,T] with rate tex2html_wrap_inline228 .

4)
Suppose N is a Poisson process with rate tex2html_wrap_inline228 . Let tex2html_wrap_inline268 be the times at which points arrive Given N(T)=n tex2html_wrap_inline272 have the same distribution as the order statistics of a sample of size n from the uniform distribution on [0,T].

5)
Given tex2html_wrap_inline278 , tex2html_wrap_inline272 have the same distribution as the order statistics of a sample of size n from the uniform distribution on [0,T].

Indications of some proofs:

1) tex2html_wrap_inline242 independent Poisson processes rates tex2html_wrap_inline288 , tex2html_wrap_inline290 . Let tex2html_wrap_inline292 be the event of 2 or more points in N in the time interval (t,t+h], tex2html_wrap_inline298 , the event of exactly one point in N in the time interval (t,t+h].

Let tex2html_wrap_inline304 and tex2html_wrap_inline306 be the corresponding events for tex2html_wrap_inline236 .

Let tex2html_wrap_inline310 denote the history of the processes up to time t; we condition on tex2html_wrap_inline310 .

Note that

displaymath316

Since

align57

we have checked one of the two infinitesimal conditions for a Poisson process.

Next note

displaymath318

so that

align65

On the other hand let tex2html_wrap_inline320 be the event of no points in N in the time interval (t,t+h] and tex2html_wrap_inline326 the same for tex2html_wrap_inline236 . Then

align69

shows

align73

Hence N is a Poisson process with rate tex2html_wrap_inline332 .

2) The infinitesimal approach used for 1 can do part of this. See text for rest. Events defined as in 1): The event tex2html_wrap_inline306 that there is one point in tex2html_wrap_inline236 in (t,t+h] is the event, tex2html_wrap_inline298 that there is exactly one point in any of the r processes together with a subset of tex2html_wrap_inline292 where there are two or more points in N in (t,t+h] but exactly one is labeled i. Since tex2html_wrap_inline352

align78

Similarly, tex2html_wrap_inline304 is a subset of tex2html_wrap_inline292 so

displaymath358

This shows each tex2html_wrap_inline236 is Poisson with rate tex2html_wrap_inline362 . To get independence requires more work; see the text for the algebraic method which is easier.

3) Fix s< t. Let N(s,t) be the number of points in (s,t]. Given N=n the conditional distribution of N(s,t) is Binomial(n,p) with p=(s-t)/T. So

align85

4): Fix tex2html_wrap_inline378 for tex2html_wrap_inline380 such that

displaymath382

Given N(T)=n we compute the probability of the event

displaymath386

Intersection of A and N(T)=n is ( tex2html_wrap_inline392 ):

multline111

whose probability is

displaymath394

So

align116

Divide by tex2html_wrap_inline396 and let all tex2html_wrap_inline398 go to 0 to get joint density of tex2html_wrap_inline272 is

displaymath402

which is the density of order statistics from a Uniform[0,T] sample of size n.

5) Replace the event tex2html_wrap_inline408 with tex2html_wrap_inline410 . With A as before we want

multline129

Note that B is independent of tex2html_wrap_inline416 and that we have already found the limit

displaymath418

We are left to compute the limit of

displaymath420

The denominator is

multline139

Thus

align142

This gives the conditional density of tex2html_wrap_inline272 given tex2html_wrap_inline408 as in 4).

Inhomogeneous Poisson Processes

The idea of hazard rate can be used to extend the notion of Poisson Process. Suppose tex2html_wrap_inline426 is a function of t. Suppose N is a counting process such that

displaymath432

and

displaymath434

Then N has independent increments and N(t+s)-N(t) has a Poisson distribution with mean

displaymath440

If we put

displaymath442

then mean of N(t+s)-N(T) is tex2html_wrap_inline446 .

Jargon: tex2html_wrap_inline228 is the intensity or instaneous intensity and tex2html_wrap_inline450 the cumulative intensity.

Can use the model with tex2html_wrap_inline450 any non-decreasing right continuous function, possibly without a derivative. This allows ties.

Continuous Time Markov Chains

Model for population growth: at time t population size is N(t). In next h time units: population goes up 1 or down 1 (or stays same). Larger movements unlikely.

Natural model has probability of birth large when population size large, small when population size small.


next up previous



Richard Lockhart
Monday October 23 13:58:14 PDT 2000