Convolution: If X and Y independent rvs with densities f and g respectively and Z=X+Y then
Differentiating wrt z we get
This integral is called the convolution of densities f and g.
If iid Exponential then has a Gamma distribution. Density of is
for s> 0.
Proof:
Then
This telescopes to
Properties of Poisson Processes
Then N is a Poisson process on [0,T] with rate .
Indications of some proofs:
1) independent Poisson processes rates , . Let be the event of 2 or more points in N in the time interval (t,t+h], , the event of exactly one point in N in the time interval (t,t+h].
Let and be the corresponding events for .
Let denote the history of the processes up to time t; we condition on .
Note that
Since
we have checked one of the two infinitesimal conditions for a Poisson process.
Next note
so that
On the other hand let be the event of no points in N in the time interval (t,t+h] and the same for . Then
shows
Hence N is a Poisson process with rate .
2) The infinitesimal approach used for 1 can do part of this. See text for rest. Events defined as in 1): The event that there is one point in in (t,t+h] is the event, that there is exactly one point in any of the r processes together with a subset of where there are two or more points in N in (t,t+h] but exactly one is labeled i. Since
Similarly, is a subset of so
This shows each is Poisson with rate . 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
4): Fix for such that
Given N(T)=n we compute the probability of the event
Intersection of A and N(T)=n is ( ):
whose probability is
So
Divide by and let all go to 0 to get joint density of is
which is the density of order statistics from a Uniform[0,T] sample of size n.
5) Replace the event with . With A as before we want
Note that B is independent of and that we have already found the limit
We are left to compute the limit of
The denominator is
Thus
This gives the conditional density of given as in 4).
Inhomogeneous Poisson Processes
The idea of hazard rate can be used to extend the notion of Poisson Process. Suppose is a function of t. Suppose N is a counting process such that
and
Then N has independent increments and N(t+s)-N(t) has a Poisson distribution with mean
If we put
then mean of N(t+s)-N(T) is .
Jargon: is the intensity or instaneous intensity and the cumulative intensity.
Can use the model with 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.