Subject probability statistics Topic Poisson Distribution Poisson Distribution

Subject : probability & statistics Topic : Poisson Distribution

Poisson Distribution The Poisson distribution is used to model the number of events occurring within a given time interval. The formula for the Poisson probability density (mass) function is the shape parameter which indicates the average number of events in the given time interval.

Some events are rather rare - they don't happen that often. For instance, car accidents are the exception rather than the rule. Still, over a period of time, we can say something about the nature of rare events. An example is the improvement of traffic safety, where the government wants to know whether seat belts reduce the number of death in car accidents. Here, the Poisson distribution can be a useful tool to answer questions about benefits of seat belt use.

Other phenomena that often follow a Poisson distribution are death of infants, the number of misprints in a book, the number of customers arriving, and the number of activations of a Geiger counter. The distribution was derived by the French mathematician Siméon Poisson in 1837, and the first application was the description of the number of deaths by horse kicking in the Prussian army.

Example Arrivals at a bus-stop follow a Poisson distribution with an average of 4. 5 every quarter of an hour. Obtain a barplot of the distribution (assume a maximum of 20 arrivals in a quarter of an hour) and calculate the probability of fewer than 3 arrivals in a quarter of an hour.

The probabilities of 0 up to 2 arrivals can be calculated directly from the formula with =4. 5 So p(0) = 0. 01111

Similarly p(1)=0. 04999 and p(2)=0. 11248 So the probability of fewer than 3 arrivals is 0. 01111+ 0. 04999 + 0. 11248 =0. 17358

R Code As with the Binomial distribution, the codes dpois and ppois will do the calculations for you.

> x=dpois(0: 20, 4. 5) >x [1] 1. 110900 e-02 4. 999048 e-02 1. 124786 e-01 1. 687179 e-01 1. 898076 e-01 [6] 1. 708269 e-01 1. 281201 e-01 8. 236295 e-02 4. 632916 e-02 2. 316458 e-02 [11] 1. 042406 e-02 4. 264389 e-03 1. 599146 e-03 5. 535504 e-04 1. 779269 e-04 [16] 5. 337808 e-05 1. 501258 e-05 3. 973919 e-06 9. 934798 e-07 2. 352979 e-07 [21] 5. 294202 e-08 >

> barplot(x, names=0: 20)

Now check that ppois gives the same answer (ppois is a cumulative distribution). > ppois(2, 4. 5) [1] 0. 1735781 >

Consider a collection of graphs for different values of

=10

In the last case, the probability of 20 arrivals is no longer negligible, so values up to, say, 30 would have to be considered.

Properties of Poisson The mean and variance are both equal to . The sum of independent Poisson variables is a further Poisson variable with mean equal to the sum of the individual means. As well as cropping up in the situations already mentioned, the Poisson distribution provides an approximation for the Binomial distribution.

Approximation: If n is large and p is small, then the Binomial distribution with parameters n and p, ( B(n; p) ), is well approximated by the Poisson distribution with parameter np, i. e. by the Poisson distribution with the same mean

Example Binomial situation, n= 100, p=0. 075 Calculate the probability of fewer than 10 successes.

> pbinom(9, 100, 0. 075) [1] 0. 7832687 > This would have been very tricky with manual calculation as the factorials are very large and the probabilities very small

The Poisson approximation to the Binomial states that will be equal to np, i. e. 100 x 0. 075 so =7. 5 > ppois(9, 7. 5) [1] 0. 7764076 > So it is correct to 2 decimal places. Manually, this would have been much simpler to do than the Binomial.

Thank you