Форум программистов, компьютерный форум, киберфорум
Статистика, теория вероятностей
Войти
Регистрация
Восстановить пароль
Блоги Сообщество Поиск Заказать работу  
 
Рейтинг 4.86/7: Рейтинг темы: голосов - 7, средняя оценка - 4.86
4 / 4 / 3
Регистрация: 28.09.2013
Сообщений: 152

Неверный результат в критерии Пирсона (критерий хи-квадрат) для проверки закона распределения

14.04.2017, 14:49. Показов 1692. Ответов 7
Метки нет (Все метки)

Студворк — интернет-сервис помощи студентам
Как я понимаю алгоритм критерия Пирсона для проверки закона распределения. На примере равномерного дискретного распределения.
Делим весь отрезок [max, min] на равные интервалы, где max - максимальное число в выборке, min - минимальное число в выборке. Подсчитываем количество попаданий в каждый интервал и делим на количество элементов в выборке. Получили эмпирические частоты. Далее считаем теоретические частоты. Для равномерного распределения - (1 / количество интервалов). Далее считаем сумму, пробегаясь по каждому интервалу по формуле - ((emp - teor) ^ 2 / teor), где emp - эмпирические частоты, teor - теоретические частоты. Далее по таблице распределения хи-квадрат ищем критическое значение и определяем уровень значимости (ошибку первого рода), тем самым мы можем судить о достоверности полученных результатов.
Теперь вопрос. Я проверяю равномерное дискретное распределение. Делю все значения на 8 интервалов, то есть степеней свободы у меня 7. Но значение хи-квадрат, полученное в результате, получается очень маленьким, что соответствует уровню значимости 0.95, а это значит, что 95% вероятности, что распределение не является равномерным, но оно точно является равномерным, так как для генерации я использую функцию, написанную и проверенную до меня другими людьми. Я использую язык R для выполнения работы.

n_size = 10000
p = 0.5
ar = irnuni_1(n_size)
size = 9
x = seq(min(ar), max(ar), length = size)
x = trunc(x)
emp = integer(size - 1)
teor = integer(size - 1)
for (i in 1:n_size) {
emp[findInterval(ar[i], x, rightmost.closed = TRUE, all.inside = TRUE)] = emp[findInterval(ar[i], x, rightmost.closed = TRUE, all.inside = TRUE)] + 1
}
for (i in 1:(size - 1)) {
emp[i] = emp[i] / n_size
}
for (i in 1:(size - 1)) {
teor[i] = 1 / (size - 1);
}
chisq_sum = double(1);
chisq_sum = 0;
for (i in 1:(size - 1)) {
chisq_sum = chisq_sum + ((emp[i] - teor[i])^2 / teor[i])
}
print(chisq_sum)

Функция irnuni_1 генерирует числа в дискретном равномерном распределении.
0
IT_Exp
Эксперт
34794 / 4073 / 2104
Регистрация: 17.06.2006
Сообщений: 32,602
Блог
14.04.2017, 14:49
Ответы с готовыми решениями:

Выбрать гипотезу о виде закона распределе-ния и проверить ее, используя критерий Пирсона.
По данным выборки выбрать гипотезу о виде закона распределе-ния и проверить ее, используя критерий Пирсона при уровне значимости α. В...

Критерий Хи-квадрат Пирсона
\chi^2 : Дана выборка X_1 \dots X_n , нужно проверить ее нормальность , т.е. из нормального распределения ли на ша выборка . ...

Объясните человеческим языком, как найти критерий Пирсона (хи-квадрат) ?
Я не могу понять суть: что и как нужно делать, чтобы сосчитать этот критерий. Я не код прошу, а помочь разобраться. У меня есть...

7
Эксперт по математике/физике
 Аватар для SSC
3390 / 1913 / 571
Регистрация: 09.04.2015
Сообщений: 5,365
14.04.2017, 15:02
Цитата Сообщение от NikitoZZZ Посмотреть сообщение
Теперь вопрос. Я проверяю равномерное дискретное распределение.
Вопрос неверный.
Нужно определить нулевую гипотезу.
Для критерия Пирсона обычно
H0: законы распределения выборок совпадают.
И если величина расчитанного критерия больше критического на заданном уровне значимости, то нулевую гипотезу отклоняют.
А утверждение
Цитата Сообщение от NikitoZZZ Посмотреть сообщение
но оно 100% является равномерным, так как для генерации я использую функцию, написанную и проверенную до меня другими людьми.
Абсолютно ничего не значит с точки зрения статистики.
0
4 / 4 / 3
Регистрация: 28.09.2013
Сообщений: 152
14.04.2017, 16:18  [ТС]
Цитата Сообщение от SSC Посмотреть сообщение
Цитата Сообщение от NikitoZZZ Посмотреть сообщение
Теперь вопрос. Я проверяю равномерное дискретное распределение.
Вопрос неверный.
Нужно определить нулевую гипотезу.
Для критерия Пирсона обычно
H0: законы распределения выборок совпадают.
И если величина расчитанного критерия больше критического на заданном уровне значимости, то нулевую гипотезу отклоняют.
А утверждение
Цитата Сообщение от NikitoZZZ Посмотреть сообщение
но оно 100% является равномерным, так как для генерации я использую функцию, написанную и проверенную до меня другими людьми.
Абсолютно ничего не значит с точки зрения статистики.
Это прекрасно, но это мне не помогает никак в данной ситуации, к сожалению.
0
 Аватар для Таланов
1965 / 1073 / 163
Регистрация: 06.12.2012
Сообщений: 4,695
14.04.2017, 16:35
Цитата Сообщение от NikitoZZZ Посмотреть сообщение
Но значение хи-квадрат, полученное в результате, получается очень маленьким, что соответствует уровню значимости 0.95, а это значит, что 95% вероятности, что распределение не является равномерным
Скорее наоборот. Гипотетическое распределение известно?
0
4 / 4 / 3
Регистрация: 28.09.2013
Сообщений: 152
14.04.2017, 16:39  [ТС]
Цитата Сообщение от Таланов Посмотреть сообщение
Скорее наоборот. Гипотетическое распределение известно?
Да, известно, то есть я могу теоретические частоты вычислить по формуле (если я правильно понял, что вы спрашиваете). Вообще мне нужно проверить дискретное логарифмическое распределение и непрерывное хи-квадрат распределение по критерию Пирсона. Просто я решил сначала разобраться на примере равномерного распределения, так как оно проще.
0
 Аватар для Таланов
1965 / 1073 / 163
Регистрация: 06.12.2012
Сообщений: 4,695
14.04.2017, 16:50
Цитата Сообщение от NikitoZZZ Посмотреть сообщение
Просто я решил сначала разобраться на примере равномерного распределения, так как оно проще.
Это правильно. Предлагаю начать с выборки объёмом 100.
0
4 / 4 / 3
Регистрация: 28.09.2013
Сообщений: 152
14.04.2017, 17:01  [ТС]
Цитата Сообщение от Таланов Посмотреть сообщение
Это правильно. Предлагаю начать с выборки объёмом 100.
Я тут нашел кое-что https://ru.wikipedia.org/wiki/... 0.B8.D1.8F
По ссылке описывается критерий Пирсона. И там есть небольшое отличие в расчете статистики критерия. Там в итоговой формуле фигурирует переменная n (количество элементов в выборке), которая умножается на сумму, которую я у себя считал. При таком раскладе критерий считается нормально, то есть значения находятся в таблице при "хорошем" уровне значимости. Я немного в замешательстве, так как в разных источниках разные формулы.
0
 Аватар для Таланов
1965 / 1073 / 163
Регистрация: 06.12.2012
Сообщений: 4,695
14.04.2017, 17:15
Цитата Сообщение от NikitoZZZ Посмотреть сообщение
Я немного в замешательстве, так как в разных источниках разные формулы.
Я тоже в замешательстве, ибо формула одна единственная.
0
Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
BasicMan
Эксперт
29316 / 5623 / 2384
Регистрация: 17.02.2009
Сообщений: 30,364
Блог
14.04.2017, 17:15
Помогаю со студенческими работами здесь

Критерии устойчивости САУ. Логарифмический критерий и критерий Найквиста
Добрый день. Возникла такая проблема: не могу точно определиться с устойчивостью/неустойчивостью замкнутой системы Передаточная функция:...

Построить графики интегральной и дифференциальной функции распределения для нормального закона распределения для ВЫБОРОК 90, 300, 1200
Помогите, пожалуйста, построить графики интегральной и дифференциальной функции распределения для нормального закона распределения для...

Применить критерий Пирсона для интервального ряда
Прошу помощи. Даже не знаю , с чего начать :( Спасибо.

Критерии Пирсона
Здравствуйте. Есть алгоритм генератора случайных чисел на заданном интервале, элементарный. Нужно всунуть сюда проверку равномерного...

Критерии Пирсона
Что такое Критерии Пирсона Для чего нужны


Искать еще темы с ответами

Или воспользуйтесь поиском по форуму:
8
Ответ Создать тему
Новые блоги и статьи
SDL3 для Web (WebAssembly): Реализация движения на Box2D v3 - трение и коллизии с повёрнутыми стенами
8Observer8 20.02.2026
Содержание блога Box2D позволяет легко создать главного героя, который не проходит сквозь стены и перемещается с заданным трением о препятствия, которые можно располагать под углом, как верхнее. . .
Конвертировать закладки radiotray-ng в m3u-плейлист
damix 19.02.2026
Это можно сделать скриптом для PowerShell. Использование . \СonvertRadiotrayToM3U. ps1 <path_to_bookmarks. json> Рядом с файлом bookmarks. json появится файл bookmarks. m3u с результатом. # Check if. . .
Семь CDC на одном интерфейсе: 5 U[S]ARTов, 1 CAN и 1 SSI
Eddy_Em 18.02.2026
Постепенно допиливаю свою "многоинтерфейсную плату". Выглядит вот так: https:/ / www. cyberforum. ru/ blog_attachment. php?attachmentid=11617&stc=1&d=1771445347 Основана на STM32F303RBT6. На борту пять. . .
Камера Toupcam IUA500KMA
Eddy_Em 12.02.2026
Т. к. у всяких "хикроботов" слишком уж мелкий пиксель, для подсмотра в ESPriF они вообще плохо годятся: уже 14 величину можно рассмотреть еле-еле лишь на экспозициях под 3 секунды (а то и больше),. . .
И ясному Солнцу
zbw 12.02.2026
И ясному Солнцу, и светлой Луне. В мире покоя нет и люди не могут жить в тишине. А жить им немного лет.
«Знание-Сила»
zbw 12.02.2026
«Знание-Сила» «Время-Деньги» «Деньги -Пуля»
SDL3 для Web (WebAssembly): Подключение Box2D v3, физика и отрисовка коллайдеров
8Observer8 12.02.2026
Содержание блога Box2D - это библиотека для 2D физики для анимаций и игр. С её помощью можно определять были ли коллизии между конкретными объектами и вызывать обработчики событий столкновения. . . .
SDL3 для Web (WebAssembly): Загрузка PNG с прозрачным фоном с помощью SDL_LoadPNG (без SDL3_image)
8Observer8 11.02.2026
Содержание блога Библиотека SDL3 содержит встроенные инструменты для базовой работы с изображениями - без использования библиотеки SDL3_image. Пошагово создадим проект для загрузки изображения. . .
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin
Copyright ©2000 - 2026, CyberForum.ru