Найти степени свободы для теста квадратов Чи в scipy?

Question

У меня есть наблюдение по максвелловскому распределению, что я отвечаю ожидаемому максвелловскому распределению. Затем я запускаю тест квадратного квадрата, чтобы узнать, насколько он подходит. Я получаю отличные результаты, однако, я также хочу узнать о степенях свободы, которые использует критерий хи-квадрат. Процитировать документацию chisquare

: p-значение вычисляется с использованием распределения по хи-квадрату с k-1-ddof степеней свободы, где k - количество наблюдаемых частот. Значение по умолчанию ddof равно 0.

Что такое k именно здесь? Общее количество точек данных (41000), которые у меня есть? Или это частота в ящике?

Answer 1

Документы соответствуют типичным именам статистических переменных. K-1 - степень свободы. K представляет количество выборок каждого размера n. Так что, по вашим словам, частота на бен.

Последний абзац http://statistics.about.com/od/Inferential-Statistics/a/What-Is-A-Degree-Of-Freedom.htm гласит:

Другой пример другого способа подсчета степеней свободы приходит с тестом F. При проведении теста F мы имеем k выборок, каждый из которых имеет размер n. Степени свободы в числителе k - 1, а в знаменателе k (n - 1).

Answer 2

k является размер f_obs, первый аргумент chisquare. Это количество бункеров.

Например, в следующем примере из строку документации,

>>> chisquare([16, 18, 16, 14, 12, 12]) 
(2.0, 0.84914503608460956) 

f_obs является [16, 18, 16, 14, 12, 12] и k является len(f_obs) или 6.