Test d’hypothèse : les valeurs de b sont-elles différentes ?
Réalisez le test d’hypothèse esquissé dans l’exercice précédent. Les variables mags_pre et mags_post sont déjà chargées dans votre espace de travail, ainsi que mt = 3.
Cet exercice fait partie du cours
Études de cas en pensée statistique
Instructions
- Extrayez les magnitudes des séismes avant 2010 dont la magnitude est supérieure ou égale au seuil de complétude et écrasez
mags_preavec le résultat. Faites de même pourmags_post. - Calculez la différence observée entre les magnitudes moyennes, en soustrayant les magnitudes des séismes d’avant 2010 de celles d’après 2010.
- Générez 10 000 réplicats de permutation avec
dcst.draw_perm_reps(). Utilisezdcst.diff_of_meanscomme argument pourfunc. - Calculez et affichez la p-valeur en considérant « au moins aussi extrême que » comme signifiant que la statistique de test est plus petite que celle observée.
Exercice interactif pratique
Essayez cet exercice en complétant cet exemple de code.
# Only magnitudes above completeness threshold
mags_pre = ____[____ >= ____]
mags_post = ____[____ >= ____]
# Observed difference in mean magnitudes: diff_obs
diff_obs = ____ - ____
# Generate permutation replicates: perm_reps
perm_reps = ____
# Compute and print p-value
p_val = ____(____ < ____) / 10000
print('p =', p_val)