# -*- coding: utf-8 -*-
"""
Created on Thu Jun 24 10:42:42 2021

@author: Sim
"""
from scipy.stats import poisson
import pandas as pd
import matplotlib.pyplot as plt

x = list(range(15))
pmf = poisson.pmf(x, 4)
cdf = poisson.cdf(x, 4)
quantile = poisson.ppf([0, 0.25, 0.5, 0.75, 1], 4, 0.5)
moments = poisson.stats(4, moments='mvsk')
median = poisson.median(4)
stdev = poisson.std(4)

xx = poisson.rvs(4, size=1000)
xx2 = pd.Series(xx)
freqx = xx2.value_counts(sort=False) /1000
plt.bar(freqx.index, freqx[freqx.index])
ax2 = plt.twinx()
ax2.plot(pmf, 'r-')
plt.show()

ecdf = freqx.cumsum()
plt.bar(freqx.index, ecdf[freqx.index])
ax2 = plt.twinx()
ax2.plot(freqx.index, cdf[freqx.index], 'r-')
plt.show()