Teste de hipótese: os valores de b são diferentes?
Realize o teste de hipótese esboçado no exercício anterior. As variáveis mags_pre e mags_post já estão carregadas no seu namespace, assim como mt = 3.
Este exercício faz parte do curso
Estudos de caso em pensamento estatístico
Instruções do exercício
- Faça o fatiamento das magnitudes dos terremotos antes de 2010 que têm magnitude acima (ou igual) ao limiar de completude e sobrescreva
mags_precom o resultado. Faça o mesmo paramags_post. - Calcule a diferença observada entre as magnitudes médias, subtraindo as magnitudes dos terremotos pré-2010 das dos pós-2010.
- Gere 10.000 réplicas de permutação usando
dcst.draw_perm_reps(). Usedcst.diff_of_meanscomo argumento parafunc. - Calcule e imprima o valor de p considerando "pelo menos tão extremo quanto" como significando que a estatística de teste é menor do que a observada.
Exercício interativo prático
Experimente este exercício completando este código de exemplo.
# 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)