Бутстрэп-тест гипотез

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

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

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

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

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

Создайте массив bout_lengths_concat, содержащий все длины серий как для особей дикого типа (bout_lengths_wt), так и для гетерозиготных особей (bout_lengths_het), с помощью np.concatenate().
Вычислите среднее значение всех длин серий из объединённого массива (bout_lengths_concat) и сохраните результат в переменную mean_bout_length.
Сдвиньте оба набора данных так, чтобы их средние стали равны mean_bout_length. Сохраните сдвинутые массивы в переменные wt_shifted и het_shifted.
Используйте dcst.draw_bs_reps(), чтобы получить 10 000 бутстрэп-реплик среднего для каждого из сдвинутых наборов данных. Сохраните реплики в bs_reps_wt и bs_reps_het соответственно.
Вычтите bs_reps_wt из bs_reps_het, чтобы получить бутстрэп-реплики разности средних. Сохраните результат в переменную bs_reps.
Вычислите p-значение, считая «не менее экстремальным» случай, когда разность средних при нулевой гипотезе больше или равна наблюдавшейся в эксперименте. Переменная diff_means_exp из предыдущего упражнения уже доступна в вашем пространстве имён.

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

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

# Concatenate arrays: bout_lengths_concat
bout_lengths_concat = ____((____, ____))

# Compute mean of all bout_lengths: mean_bout_length
mean_bout_length = ____

# Generate shifted arrays
wt_shifted = ____ - np.mean(____) + ____
het_shifted = ____ - ____ + ____

# Compute 10,000 bootstrap replicates from shifted arrays
bs_reps_wt = ____
bs_reps_het = ____

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

# Compute and print p-value: p
p = ____(____ >= ____) / len(____)
print('p-value =', p)

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