Permutační test: divoký typ vs. heterozygot
Pomocí permutačního testu otestuj hypotézu, že délky záchvatů u heterozygota a divokého typu mají stejné rozdělení.
Toto cvičení je součástí kurzu
Případové studie ve statistickém myšlení
Pokyny k cvičení
- Vypočítej rozdíl průměrů (délky záchvatů heterozygota minus divoký typ) skutečných dat a výsledek ulož do proměnné
diff_means_exp. Polenumpys názvembout_lengths_wtabout_lengths_hetjsou již dostupná v prostředí. - Pomocí
dcst.draw_perm_reps()vygeneruj 10 000 permutačních replikátů rozdílu průměrů. Můžeš použít funkcidcst.diff_of_means()a výsledek ulož doperm_reps. - Vypočítej p-hodnotu, přičemž za „přinejmenším stejně extrémní" považuj případ, kdy je rozdíl průměrů za nulové hypotézy větší nebo roven hodnotě pozorované experimentálně.
- Vypiš p-hodnotu na obrazovku.
Interaktivní cvičení na vyzkoušení si v praxi
Vyzkoušejte si toto cvičení dokončením tohoto ukázkového kódu.
# 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)