Перевірка гіпотези: чи змінилася частота землетрусів?

Очевидно, після початку закачування стічної води частота землетрусів різко зросла. Втім, для практики ви все одно виконаєте перевірку гіпотези. Ви не перевірятимете гіпотезу про те, що інтервали між землетрусами мають однаковий розподіл до і після 2010 року, адже закачування може впливати на розподіл. Натомість ви припустите, що їхні середні однакові. Отже, обчисліть p-значення для гіпотези, що середні інтервали між землетрусами до 2010 року та після 2010 року однакові, використовуючи як статистику перевірки різницю: середнє до 2010 мінус середнє після 2010.

Ця вправа є частиною курсу

Кейс-стаді зі статистичного мислення

Переглянути курс

Інструкції до вправи

Обчисліть спостережувану статистику перевірки. Змінні mean_dt_pre і mean_dt_post з попередніх вправ доступні у вашому просторі імен.
Зсуньте дані після 2010 року так, щоб їхнє середнє збігалося із середнім до 2010 року. Збережіть результат у dt_post_shift.
Згенеруйте по 10 000 bootstrap-реплікатів середнього для dt_pre і dt_post_shift. Збережіть відповідні результати в bs_reps_pre і bs_reps_post.
Обчисліть реплікати різниці середніх, віднявши bs_reps_post від bs_reps_pre.
Обчисліть і виведіть p-значення. Вважайте «не менш екстремальним», коли значення статистики перевірки більше або дорівнює спостережуваному.

Інтерактивна практична вправа

Спробуйте виконати цю вправу, доповнивши цей зразок коду.

# Compute the observed test statistic
mean_dt_diff = ____ - ____

# Shift the post-2010 data to have the same mean as the pre-2010 data
dt_post_shift = ____ - ____ + ____

# Compute 10,000 bootstrap replicates from arrays
bs_reps_pre = ____
bs_reps_post = ____

# Get replicates of difference of means
bs_reps = ____ - ____

# Compute and print the p-value
p_val = ____(____ >= ____) / 10000
print('p =', p_val)

Редагувати та запускати код