import numpy as np
import pandas as pd
import matplotlib.pyplot as plt

data = pd.read_excel('TissueMNIST.xlsx')  

DA = data.iloc[:,0]
SVD = data.iloc[:,1]

fig = plt.figure(figsize=(10, 6))
ax1 = fig.add_subplot(2, 1, 1)

SVD.hist(bins=30,alpha=0.5,ax=ax1)
SVD.plot(kind='kde', secondary_y=True,ax=ax1)     
plt.grid() 

ax2 = fig.add_subplot(2, 1, 2)
DA.hist(bins=30,alpha=0.5,ax=ax2)
DA.plot(kind='kde', secondary_y=True,ax=ax2)
plt.grid()      
 
plt.show()


from statsmodels.stats.weightstats import ztest as ztest
ztest(DA, SVD, value=0, alternative = 'larger')