Как я использую распределение хи-квадрат с библиотекой Boost C++?
Я проверил примеры в веб-сайт Повышения, но они не то, что я ищу.
Для помещения его простой я хочу видеть, одобрено ли число на умирании, с помощью 600 списков, таким образом, средние появления каждого номера (1 - 6) должны быть 100.
И я хочу использовать распределение хи-квадрат, чтобы проверить, справедливо ли умирание.
Справка!, как я сделал бы это??
1 ответ
Предположим, что E [I] и O [I] являются массивами, удерживающими ожидаемое и наблюдаемое количество валков для каждого из 6 возможностей. В вашем случае E [I] составляет 100 для каждой корзины, а o [i] - это количество раз, когда я был прокатан в ваших 600 испытаний.
Затем вы рассчитываете статистику в квадрату Chi, суммируя (E [I] -O [I]) 2 / E [I] 6 банок. Позвольте сказать, что ваш o [i] массив вышел с 105, 95, 102, 98, 98 и 102 Считается после выполнения ваших 600 испытаний.
CHI 2 = 5 2 / 100 + 5 2 / 100 + 2 2 / 100 + 2 2 / 100 + 2 2 / 100 + 2 2 / 100 = 0660
У вас есть пять градусов свободы (количество банок минус 1). Так что вы собираетесь Имейте декларацию, такую как
boost::math::chi_squared mydist(5);
, чтобы создать объект Boost, представляющий распределение вашего хи-квадрата.
На данный момент вы будете использовать функцию Accessor CDF (на кумулятивной функции распределения) из библиотеки BOOST, чтобы посмотреть значение P-значения, соответствующее CHI-квадрату .660 с пятью степенями свободы. Отказ
p = boost::math::cdf(mydist,.660);
Вы должны получить что-то близкое к 0,015, который будет интерпретироваться как (1 - 0,015) = 98,5% вероятность наблюдения за квадратом Chi-квадрата, по меньшей мере, как экстрим, чем 0,660, если принять участие в нулевой гипотезе (что умереть справедливо) держит. Таким образом, для этого набора данных нулевая гипотеза не может быть отклонена любым разумным уровнем доверия. (Отказ от ответственности: непроверенный код! Но если я правильно понимаю повышенную документацию, вот как это должно работать.)