# family of distributions parameterized by lambda ("mean rate")poisson_values <-rpois(n =10000, lambda =5)hist(poisson_values, breaks =30)
# normalize distribution for QQpoisson_values <-(poisson_values -mean(poisson_values))/sd(poisson_values)qqnorm(poisson_values); qqline(poisson_values, col =2)