Quantificazione dei valori b
In base al grafico che hai generato nell'esercizio precedente, puoi usare in sicurezza una soglia di completezza mt = 3. Con questa soglia, calcola i valori di b per il periodo tra il 1980 e il 2009 e per il 2010 fino a metà 2017. La funzione b_value() che hai scritto nel capitolo precedente, che calcola il valore di b e l'intervallo di confidenza a partire da un insieme di magnitudini e da una soglia di completezza, è disponibile nel tuo namespace, così come gli array numpy mags_pre e mags_post dell'ultimo esercizio e mt.
Questo esercizio fa parte del corso
Casi di studio nel pensiero statistico
Istruzioni dell'esercizio
- Calcola il valore di b e l'intervallo di confidenza al 95% per i terremoti dal 1980 al 2009 usando 10.000 replicazioni bootstrap.
- Calcola il valore di b e l'intervallo di confidenza al 95% per i terremoti dal 2010 fino a metà 2017 usando 10.000 replicazioni bootstrap.
- Premi "Invia risposta" per stampare i risultati a schermo.
esercizio interattivo pratico
Prova questo esercizio completando questo codice di esempio.
# 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))