#, -*-, coding:, utf-8, -*-
"""
Created, on, Mon, Jan, , 3, 17:35:29, 2022

@author:, Sim
"""
import statsmodels.stats.api as sms
import scipy.stats as st

ctrl=[89, 78, 82, 94, 81, 82, 70, 93, 79, 86, 87, 81, 80, 
      89, 79, 78, 84]
trt=[85, 89, 75, 85, 99, 74, 82, 92, 92, 81, 77, 85, 77, 
     84, 83, 88, 85, 78, 88, 82, 94, 77, 79, 84, 73, 83]

Fval, pvalF = st.levene(ctrl, trt, center='mean')
print(Fval, pvalF)

cm = sms.CompareMeans(sms.DescrStatsW(ctrl, ddof=1), 
                      sms.DescrStatsW(trt, ddof=1))
tval, pval, df = cm.ttest_ind(usevar='pooled')
print(tval, pval, df)

lbd, ubd = cm.tconfint_diff(usevar='pooled')
print (lbd, ubd)