from scipy.stats import kurtosis # list containing numbers only l = [1.8, 2, 1.2, 1.5, 1.6, 2.1, 2.8] # switch to numpy array v = np.array(l) k = kurtosis(v) # ~ -0.25