import numpy as np
#import scipy.stats as stats
import seaborn as sns # seaborn package를 이용
import matplotlib.pyplot as plt # matplotlib.pyplot package를 이용
from numpy import random

# set the random seed:
np.random.seed(12345)

r=10000

t1 = random.standard_t(df=5, size=r)
sns.histplot(data=t1, x=None).set(title='Histogram of Student-t w/ df=5')
plt.show()

t2 = random.standard_t(df=10, size=r)
sns.histplot(data=t2, x=None).set(title='Histogram of Student-t w/ df=10')
plt.show()

t3 = random.standard_t(df=15, size=r)
sns.histplot(data=t3, x=None).set(title='Histogram of Student-t w/ df=15')
plt.show()

t4 = random.standard_t(df=30, size=r)
sns.histplot(data=t4, x=None).set(title='Histogram of Student-t w/ df=30')
plt.show()

fig = plt.figure(figsize=(14,7))
# 두 개의 그래프를 한 페이지에 그림
fig, axs = plt.subplots(ncols=2)
fig, ays = plt.subplots(ncols=2)

sns.histplot(data=t1, x=None, ax=axs[0]).set(title='Histogram of Student-t w/ df=5 & df=10)')
sns.histplot(data=t2, x=None, ax=axs[1])
fig.subplots_adjust(wspace=0.5) # 우측그림의 좌우 간격을 조정
#plt.savefig('C:/BOOK/PyBasics/PyStat/code/t-1.png')

sns.histplot(data=t3, x=None, ax=ays[0]).set(title='Histogram of Student-t w/ df=15 & df=30)')
sns.histplot(data=t4, x=None, ax=ays[1])
fig.subplots_adjust(wspace=0.5) # 우측그림의 좌우 간격을 조정
#plt.savefig('C:/BOOK/PyBasics/PyStat/code/t-2.png')
plt.show()