To find the Poisson probability of seeing exactly $x$ occurrences of some event within some well-defined interval where the mean number of occurrences in that interval is expected to be $\lambda$, we can use the following R function:
dpois(x, lambda = λ)
The probability of seeing exactly $3$ blemishes on a randomly selected piece of sheet metal, when on average one expects 1.2 blemishes, could be calculated with the following:
> dpois(3, lambda = 1.2)  0.08674393
If one needs to find the probability that $x$ or fewer occurrences are seen when $\lambda$ are expected, one can use the related cumulative probability distribution function:
ppois(x, lambda = λ)
In a call center that gets on average 13 calls every hour, the probability that in a given $15$ minute period they will receive less than $6$ calls can be calculated with the following:
> ppois(5, lambda = 13/4)  0.8888132
R also has a function that lets you simulate random variables that follow a Poisson distribution:
rpois(n, lambda = λ)
Note: unlike the previous functions, here $n$ represents the number of realizations of your random variable you wish to produce.
Over the course of 15 weeks, every Saturday -- at the same time -- an individual stands by the side of a road and tallies the number of cars going by within a 120-minute window. Simulate a set of observations for this situation under the assumption that the mean number of cars that go by per hour is actually $52$.
> rpois(15, lambda=104)  96 113 106 104 101 96 96 114 109 91 94 95 108 113 113