# -*- coding: utf-8 -*-
"""
Created on Fri Nov 19 14:08:37 2021

@author: Sim
"""

import scipy.stats as st
import statsmodels.api as sm
import matplotlib.pyplot as plt
import numpy as np

fig, (ax1, ax2, ax3) = plt.subplots(1,3, figsize=(12, 5))
xx = [10, 20, 30, 40, 50]

res = st.probplot(xx, plot=ax1)
print(res)

sm.qqplot(np.array(xx), dist=st.norm, line='q', ax=ax2)

n = len(xx)
p = np.arange(1,n+1)/(n+1)
q = st.norm.ppf(p)
ax3.plot(q, xx, 'bo')