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

Очевидно, что после начала закачки сточных вод частота землетрясений резко возросла. Тем не менее для практики вы проведёте проверку гипотезы. Вы не будете проверять, одинаково ли распределены промежутки между землетрясениями до и после 2010 года, поскольку закачка сточных вод могла повлиять на это распределение. Вместо этого предположите, что средние значения совпадают. Вычислите p-значение для гипотезы о том, что средние промежутки между землетрясениями до и после 2010 года равны, используя в качестве тестовой статистики разность среднего значения до 2010 года и среднего значения после 2010 года.

Это упражнение является частью курса

Примеры применения статистического мышления

Посмотреть курс

Инструкции к упражнению

Вычислите наблюдаемую тестовую статистику. Переменные mean_dt_pre и mean_dt_post из предыдущих упражнений уже доступны в вашем пространстве имён.
Сдвиньте данные за период после 2010 года так, чтобы их среднее совпало со средним значением до 2010 года. Сохраните результат в переменную dt_post_shift.
Получите по 10 000 бутстреп-реплик среднего значения для 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)

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