Отличная обзорная статья. Спасибо

Отличная статья

Про n-1 - не совсем. В этой формуле нужно именно n, а вот когда будем оценивать стандартное отклонение - там в знаменателе будет n-1 

Да, если в std отнимать -1 то потом уже не надо. Просто в numpyпо умолчанию не отнимают (см. ddof=0)

Реализация которая совпадет с scipy:

def t_stat(x,y):
    n = len(x)
    return (np.mean(x) - np.mean(y)) / (np.std(x-y, ddof=1)/np.sqrt(n))

S – дисперсия попарных разностей, n – число фолдов.

Вот тут ошибка. В формуле для t-статистики S это стандартное отклонение.

Еще в формуле вместо n надо брать n-1

Классная статья: минимум воды, понятный вывод

Отличная работа, надеюсь у вас хватит времени/сил завершить перевод :)