K = t.ppf((1 + p)/2, df['Mx']-1) в степень свободы взял среднее, надо - размер выборки ssw = sum(sum((x - group_mean)**2 for x in data[group]) for group, group_mean in data.mean().items()) посчитано sst, проверь - будет 30