# -*- coding: utf-8 -*-
"""
Created on Wed Jun 23 14:30:23 2021

@author: Sim
"""
from scipy.stats import t
from scipy.stats import norm
import numpy as np
import matplotlib.pyplot as plt

x = np.linspace(-3,3, num=100)
pdf0 = norm.pdf(x,0,1)
pdf1 = t.pdf(x, 3)
pdf2 = t.pdf(x, 5)
pdf3 = t.pdf(x, 9)

plt.plot(x, pdf0, 'k-', label='N(0,1)')
plt.plot(x, pdf1, 'r-', label='t_3')
plt.plot(x, pdf2, 'g-.', label='t_5')
plt.plot(x, pdf3, 'b--', label='t_9')
plt.legend()
plt.show()


quantilet9 = t.ppf([0.005, 0.01, 0.025, 0.5, .95, .975, .99, .995], 9)
quantilen = norm.ppf([0.005, 0.01, 0.025, 0.5, .95, .975, .99, .995])

moments = t.stats(9, moments='mvsk')
median = t.median(9)
stdev = t.std(9)