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

f1 = random.f(5, 10,size=r)
sns.histplot(data=f1, x=None).set(title='Histogram of F w/ df=(5,10)')
plt.show()

f2 = random.f(9, 10, size=r)
sns.histplot(data=f2, x=None).set(title='Histogram of F w/ df=(9,10)')
plt.show()

f3 = random.f(15,20, size=r)
sns.histplot(data=f3, x=None).set(title='Histogram of F / df=(15,20)')
plt.show()

f4 = random.f(38,40, size=r)
sns.histplot(data=f4, x=None).set(title='Histogram of F w/ df=(38,40)')
plt.show()

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

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

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