Quantification des valeurs de b
D’après le graphique généré à l’exercice précédent, vous pouvez retenir sans risque un seuil de complétude de mt = 3. En utilisant ce seuil, calculez les valeurs de b pour la période de 1980 à 2009, puis pour 2010 jusqu’à mi-2017. La fonction b_value() que vous avez écrite au chapitre précédent, qui calcule la valeur de b et l’intervalle de confiance à partir d’un ensemble de magnitudes et d’un seuil de complétude, est disponible dans votre espace de travail, tout comme les tableaux numpy mags_pre et mags_post de l’exercice précédent, ainsi que mt.
Cet exercice fait partie du cours
Études de cas en pensée statistique
Instructions
- Calculez la valeur de b et l’intervalle de confiance à 95 % pour les séismes de 1980 à 2009 en utilisant 10 000 répliques bootstrap.
- Calculez la valeur de b et l’intervalle de confiance à 95 % pour les séismes de 2010 jusqu’à mi-2017 en utilisant 10 000 répliques bootstrap.
- Cliquez sur "Soumettre la réponse" pour afficher les résultats à l’écran.
Exercice interactif pratique
Essayez cet exercice en complétant cet exemple de code.
# Compute b-value and confidence interval for pre-2010
b_pre, conf_int_pre = ____(____, ____, ____=[____, ____], ____=____)
# Compute b-value and confidence interval for post-2010
b_post, conf_int_post = ____
# Report the results
print("""
1980 through 2009
b-value: {0:.2f}
95% conf int: [{1:.2f}, {2:.2f}]
2010 through mid-2017
b-value: {3:.2f}
95% conf int: [{4:.2f}, {5:.2f}]
""".format(b_pre, *conf_int_pre, b_post, *conf_int_post))