Test d'hypothèse : les valeurs 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 noms, de même que mt = 3.
Cette activité fait partie du cours
Études de cas en pensée statistique
Instructions de l’exercice
- Extrayez les magnitudes des séismes antérieurs à 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 des séismes d'après 2010.
- Générez 10 000 réplicats par permutation avec
dcst.draw_perm_reps(). Utilisezdcst.diff_of_meanscomme argument pourfunc. - Calculez et affichez la valeur p en interprétant « 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 ce code d’exemple.
# 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)