Перестановковий тест: wild type проти heterozygote

Перевірте гіпотезу про те, що розподіли довжин епізодів для heterozygote і wild type є однаковими, використавши перестановковий тест.

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

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

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

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

Обчисліть різницю середніх (heterozygote мінус wild type) для фактичних наборів даних і збережіть результат у змінній diff_means_exp. Масиви numpy bout_lengths_wt і bout_lengths_het уже є у вашому просторі імен.
Згенеруйте 10 000 перестановкових реплік різниці середніх за допомогою dcst.draw_perm_reps(). Також можна використати функцію dcst.diff_of_means() і зберегти результат у perm_reps.
Обчисліть p-значення, вважаючи „не менш екстремальним", що різниця середніх за нульової гіпотези є більшою або дорівнює спостереженій експериментально.
Виведіть p-значення на екран.

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

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

# Compute the difference of means: diff_means_exp
diff_means_exp = ____ - ____

# Draw permutation replicates: perm_reps
perm_reps = ____(____, ____, 
                               ____, size=____)

# Compute the p-value: p-val
p_val = ____(____ >= ____) / len(____)

# Print the result
print('p =', p_val)

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