# -*- coding: utf-8 -*-
"""
Created on Tue Jan  4 09:20:22 2022

@author: Sim
"""

import pandas as pd
import statsmodels.stats.api as sms
import scipy.stats as st

df = pd.read_csv(r'D:\HTEX\Pythonbk\codesdata\bmi.csv')
gdf = df.groupby('gender')

mht = gdf.get_group('M')['height']
fht = gdf.get_group('F')['height']

Fval, pvalF = st.levene(mht, fht, center='mean')
print(Fval, pvalF)

cm = sms.CompareMeans(sms.DescrStatsW(mht, ddof=1), 
                      sms.DescrStatsW(fht, ddof=1))
tval, pval, degf = cm.ttest_ind(usevar='pooled')
print(tval, pval, degf)

#----------------------------------------------------------
mwt = gdf.get_group('M')['weight']
fwt = gdf.get_group('F')['weight']

Fval, pvalF = st.levene(mwt, fwt, center='mean')
print(Fval, pvalF)

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