Poisson Distribution

P(X=k)=limn(nk)(λn)k(1λn)nk=λkk!eλ P(X =k) = \lim_{n \rightarrow \infty}\binom{n}{k}(\frac{\lambda}{n})^k(1-\frac{\lambda}{n})^{n-k} \newline = \frac{\lambda ^k}{k!} e^{-\lambda}

λis just any given number\lambda \Longrightarrow \text{is just any given number}

Poisson Random Variables

We've already learned that we can break many problems down into terms of nn, xx, and pp and use the following formula for binomial random variables:

p(x)=(nx)px(1p)nxp(x)=\binom{n}{x} \cdot p^x \cdot (1-p)^{n-x}

But what do we do when cannot be calculated using that formula? Enter the Poisson random variable.

Poisson Experiment

A Poisson experiment is a statistical experiment that has the following properties:

  • The outcome of each trial is either success or failure.

  • The average number of successes (λ)(\lambda) that occurs in a specified region is known.

  • The probability that a success will occur is proportional to the size of the region.

  • The probability that a success will occur in an extremely small region is virtually zero.

A Poisson random variable is the number of successes that result from a Poisson experiment. The probability distribution of a Poisson random variable is called a Poisson distribution:

P(k,λ)=λkeλk!\huge P(k, \lambda) = \frac{\lambda^k e^{-\lambda}}{k!}

Here,

  • e=2.71828e=2.71828

  • λ\lambda is the average number of successes that occur in a specified region.

  • kk is the actual number of successes that occur in a specified region.

  • P(k,λ)P(k, \lambda) is the Poisson probability, which is the probability of getting exactly kk successes when the average number of successes is λ\lambda.

Example

Acme Realty company sells an average of 22 homes per day. What is the probability that exactly 33 homes will be sold tomorrow?

Here, λ=2\lambda=2 and k=3k=3, so P(k=3,λ=2)=λkeλk!=0.180P(k=3, \lambda =2)=\frac{\lambda^ke^{-\lambda}}{k!}=0.180

Example

Suppose the average number of lions seen by tourists on a one-day safari is 55. What is the probability that tourists will see fewer than lions on the next one-day safari? P(k3,λ=5)=r=03λreλr!=0.2650P(k\leq3,\lambda=5)=\sum_{r=0}^{3}\frac{\lambda^re^{-\lambda}}{r!}=0.2650

Special Case

Consider some Poisson random variable, XX. Let E[X]E[X] be the expectation of XX. Find the value of E[X2]E[X^2].

Let Var(X)Var(X) be the variance of XX. Recall that if a random variable has a Poisson distribution, then:

  • E[X]=λE[X] = \lambda

  • Var(X)=λVar(X) = \lambda

Now, we'll use the following property of expectation and variance for any random variable, XX:

Var(X)=E[X2](E[X])2E[X2]=Var(X)+(E[X])2Var(X)=E[X^2]-(E[X])^2 \Rightarrow E[X^2] = Var(X) +(E[X])^2

So, for any random variable XX having a Poisson distribution, the above result can be rewritten as:

E[X2]=λ+λ2\Rightarrow E[X^2]=\lambda + \lambda^2

Problem 1

A random variable, XX, follows Poisson distribution with mean of 2.52.5. Find the probability with which the random variable XX is equal to 55.

def factorial(n):
    if n == 0:
        return 1
    else:
        return n * factorial(n-1)

k = 5
l = 2.5
e = 2.71828

dist = factorial(k)
p = ((l**k)*(e**(l*-1)))/dist
print(round(p, 3))

Problem 2

The manager of a industrial plant is planning to buy a machine of either type AA or type BB. For each day’s operation:

  • The number of repairs, XX, that machine AA needs is a Poisson random variable with mean 0.880.88. The daily cost of operating AA is CA=160+40X2C_A = 160 + 40X^2.

  • The number of repairs, YY, that machine BB needs is a Poisson random variable with mean 1.551.55. The daily cost of operating BB is CB=128+40Y2C_B=128+40Y^2.

Assume that the repairs take a negligible amount of time and the machines are maintained nightly to ensure that they operate like new at the start of each day. Find and print the expected daily cost for each machine.

X = 0.88
Y = 1.55

print(round((160 + (40*(X + X**2))), 3))  # => 226.176
print(round((128 + (40*(Y + Y**2))), 3))  # => 286.1

Last updated