Визуализация тестовых статистик

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

Сначала вы рассмотрите две «генеральные совокупности», сгруппированные по ранним и поздним временным меткам, и построите распределение тестовой статистики. Затем перемешайте эти две совокупности так, чтобы данные больше не были упорядочены по времени и каждая содержала смесь ранних и поздних значений, — после чего пересчитайте распределение тестовой статистики.

Для начала работы мы заранее загрузили две группы длительностей, group_duration_short и group_duration_long, а также две функции: shuffle_and_split() и plot_test_statistic().

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

Введение в линейное моделирование на Python

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

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

С помощью np.random.choice() выполните ресэмплинг group_duration_short и group_duration_long, затем вычислите разность ресэмплов, чтобы получить test_statistic_unshuffled.
Примените shuffle_and_split() к исходным group_duration_short и group_duration_long (в указанном порядке), чтобы создать две новые смешанные совокупности.
Выполните ресэмплинг перемешанных совокупностей и вычтите resample_short из resample_long, чтобы получить новое значение test_statistic_shuffled.
Используйте plot_test_statistic() для построения обоих распределений тестовой статистики и сравните их визуально.

Интерактивное практическое упражнение

Попробуйте выполнить это упражнение, дополнив этот пример кода.

# From the unshuffled groups, compute the test statistic distribution
resample_short = np.random.choice(____, size=500, replace=____)
resample_long = np.random.choice(____, size=500, replace=____)
test_statistic_unshuffled = ____ - ____

# Shuffle two populations, cut in half, and recompute the test statistic
shuffled_half1, shuffled_half2 = shuffle_and_split(____, ____)
resample_half1 = np.random.choice(____, size=500, replace=____)
resample_half2 = np.random.choice(____, size=500, replace=____)
test_statistic_shuffled = resample_half2 - resample_half1

# Plot both the unshuffled and shuffled results and compare
fig = plot_test_statistic(____, label='Unshuffled')
fig = plot_test_statistic(____, label='Shuffled')

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