Тест хи-квадрат в Python

) Я использовал следующий код в R , чтобы определить, насколько хорошо наблюдаемые значения (20, 20, 0 и 0, например) соответствуют ожидаемым значениям / соотношениям (например, 25% для каждого из четырех случаев):

> chisq.test(c(20,20,0,0), p=c(0.25, 0.25, 0.25, 0.25))

    Chi-squared test for given probabilities

data:  c(20, 20, 0, 0)

X-squared = 40, df = 3, p-value = 1.066e-08

Как я могу воспроизвести это в Python? Я пробовал использовать функцию chisquare из scipy , но результаты, которые я получил, были очень разными; я не уверен, что это вообще правильная функция для использования. Я поискал в документации scipy , но это довольно устрашающе, поскольку он занимает более 1000 страниц; документация numpy почти на 50% больше.