Создайте массив 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 из предыдущего упражнения уже доступна в вашем пространстве имён.

Доверительные интервалы методом бутстрепа

Create Your Free Account