# -*- coding: utf-8 -*-
"""
Created on Thu Nov 11 15:37:44 2021

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

x = np.array([0.005, 0.600, 1.325, 1.505, 1.520, 2.000, 2.255, 2.575, 
     3.000, 3.300, 3.500])
y = np.array([0.005, 0.005, 0.255, 0.505, 0.525, 1.000, 1.250, 1.505, 
     1.750, 2.750])

fval, pval = st.levene(x,y, center='mean')
print(fval, pval)

## 직접계산
mx = np.mean(x); my = np.mean(y)
z1 = np.abs(x-mx); z2 = np.abs(y-my)
tval, pval = st.ttest_ind(z1, z2, equal_var=True, 
                          alternative='two-sided')
fval2 = tval ** 2
print(fval2, pval)

## center = median

mx = np.median(x); my = np.median(y)
z1 = np.abs(x-mx); z2 = np.abs(y-my)
tval, pval = st.ttest_ind(z1, z2, equal_var=True, 
                          alternative='two-sided')
fval2 = tval ** 2
print(fval2, pval)