# -*- coding: utf-8 -*-
"""
Created on Thu Jan  6 13:25:47 2022

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

F0 = 0.6596/0.02115
SST = 1.9788+0.1692
print(F0, SST)
#--------------------------------------------
print(1-st.f.cdf(31.18676, 3,8))
#-------------------------------------------
lbd1 = 2.5 - st.t.ppf(0.975, 8)*np.sqrt(0.02115/3)
ubd1 = 2.5 + st.t.ppf(0.975, 8)*np.sqrt(0.02115/3)
print(lbd1, ubd1)
#--------------------------------------------------------
lbd2=2.5-2.7 - st.t.ppf(0.975, 8)*np.sqrt(0.02115*(1/3+1/4))
ubd2=2.5-2.7 + st.t.ppf(0.975, 8)*np.sqrt(0.02115*(1/3+1/4))
print(lbd2, ubd2)
#------------------------------------------------------
t0 = (2.5-2.7)/np.sqrt(0.02115*(1/3+1/4))
p = 2*(1-st.t.cdf(np.abs(t0),8))
print(t0, p)