ガウス分布とポアソン分布をqq-plotする

Import

import numpy as np
from scipy import stats
import matplotlib.pyplot as plt

乱数を生成

ガウス分布、ポアソン分布に従う乱数をそれぞれ１００００個生成する。

norm = np.random.normal(0, 1, 10000)#平均０、標準偏差１
poisson = np.random.poisson(lam=10, size=10000)#λ=10

これをヒストグラムで図示する。

qq-plot

正規分布をqq-plotする。

stats.probplot(gaussian, dist='norm', plot=plt)
plt.show()

次にポアソン分布をqq-plotする。

stats.probplot(poisson, dist='poisson', sparams=10, plot=plt)#λ=10
plt.show()