# -*- coding: utf-8 -*-
"""
Created on Mon Jan  3 12:09:08 2022

@author: Sim
"""
import scipy.stats as st, numpy as np

lbd, ubd = st.t.interval(0.90, 25-1, loc=30000, scale=5000/25**.5)
print(lbd, ubd)

t0 = (30000-35000)/(5000/25**.5)
p0 = 2*(1-st.t.cdf(np.abs(t0), df=25-1)) # 양측검정
print(t0, p0)