Test par permutation : type sauvage versus hétérozygote
Vérifiez, à l'aide d'un test par permutation, l'hypothèse voulant que les durées de séquence du type sauvage et de l'hétérozygote soient distribuées de façon identique.
Cette activité fait partie du cours
Études de cas en pensée statistique
Instructions de l’exercice
- Calculez la différence des moyennes (hétérozygote moins type sauvage) des jeux de données réels et stockez le résultat dans la variable
diff_means_exp. Les tableauxnumpybout_lengths_wtetbout_lengths_hetsont déjà dans votre espace de noms. - Générez 10 000 répliques par permutation de la différence des moyennes avec
dcst.draw_perm_reps(). Vous pouvez aussi utiliser la fonctiondcst.diff_of_means(). Stockez le résultat dansperm_reps. - Calculez la valeur p, en définissant « au moins aussi extrême que » comme étant une différence des moyennes sous l'hypothèse nulle qui est supérieure ou égale à celle observée expérimentalement.
- Affichez la valeur p à l'écran.
Exercice interactif pratique
Essayez cet exercice en complétant ce code d’exemple.
# 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)