# -*- coding: utf-8 -*-
"""
Created on Mon Jan  3 13:21:30 2022

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

no=np.array([2.17, 0.17, 0.83, 3.04, 1.06, 3.93, 5.38, 3.53, 4.54, 
    3.98, 1.02, 3.46, 1.19, 5.64, 2.66, 1.22, 2.64, 2.31, 4.49])
yes=np.array([6.96, 13.00, 5.54, 6.19, 7.31, 14.34, 4.74, 2.50,  1.27,
    4.18,  3.10, 3.01, 6.81, 8.28,  2.33, 7.18, 5.67, 12.68,
    8.94,  3.16, 3.30, 6.99])

Fval, pvalF = st.levene(yes,no, center='mean')
print(Fval, pvalF)

cm = sms.CompareMeans(sms.DescrStatsW(no, ddof=1), sms.DescrStatsW(yes, ddof=1))
print (cm.tconfint_diff(usevar='unequal'))