Quantificação dos valores b
Com base no gráfico que você gerou no exercício anterior, dá para usar com segurança um limiar de completude de mt = 3. Usando esse limiar, calcule os valores de b para o período entre 1980 e 2009 e para 2010 até meados de 2017. A função b_value() que você escreveu no capítulo anterior, que calcula o valor de b e o intervalo de confiança a partir de um conjunto de magnitudes e de um limiar de completude, está disponível no seu namespace, assim como os arrays numpy mags_pre e mags_post do último exercício, e mt.
Este exercício faz parte do curso
Estudos de caso em pensamento estatístico
Instruções do exercício
- Calcule o valor de b e o intervalo de confiança de 95% para os terremotos de 1980 até 2009 usando 10.000 réplicas de bootstrap.
- Calcule o valor de b e o intervalo de confiança de 95% para os terremotos de 2010 até meados de 2017 usando 10.000 réplicas de bootstrap.
- Clique em 'Enviar Resposta' para imprimir os resultados na tela.
Exercício interativo prático
Experimente este exercício completando este código de exemplo.
# 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))