p-значення двома способами
Ви вже бачили, що зазвичай є два способи отримати нульовий розподіл: обчислювальний і математична апроксимація. Перевірка узгодженості за критерієм хі-квадрат не є винятком. Розподілом апроксимації знову є «розподіл хі-квадрат» зі ступенями вільності, що дорівнюють кількості категорій мінус один.
У цій вправі ви порівняєте ці два підходи, щоб обчислити p-значення, яке вимірює узгодженість розподілу перших цифр в Ірані із законом Бенфорда. Зверніть увагу, що спостережувана статистика, яку ви створили в попередній вправі, збережена у вашому робочому середовищі як chi_obs_stat.
Ця вправа є частиною курсу
Статистичні висновки для категоріальних даних у R
Інтерактивна практична вправа
Спробуйте виконати цю вправу, доповнивши цей зразок коду.
# Compute degrees of freedom
degrees_of_freedom <- ___ %>%
# Pull out first_digit vector
pull("first_digit") %>%
# Calculate n levels and subtract 1
___
# Plot both null dists
___
# Add density layer
___
# Add vertical line at obs stat
___
# Overlay chisq approx
stat_function(fun = dchisq, args = list(df = ___), color = "blue")