EmpezarEmpieza gratis

Prueba de hipótesis: ¿cambió la frecuencia de terremotos?

Es evidente que hubo un aumento enorme en la frecuencia de terremotos cuando comenzó la inyección de aguas residuales salinas. Aun así, harás una prueba de hipótesis para practicar. No vas a comprobar la hipótesis de que los tiempos entre terremotos tienen la misma distribución antes y después de 2010, ya que la inyección de aguas residuales puede afectar a la distribución. En su lugar, asumirás que tienen la misma media. Por tanto, calcula el valor p asociado a la hipótesis de que los tiempos entre terremotos de antes y después de 2010 tienen la misma media, usando como estadístico de prueba la media de los intervalos previos a 2010 menos la media de los posteriores a 2010.

Este ejercicio forma parte del curso

Casos prácticos de pensamiento estadístico

Ver curso

Instrucciones del ejercicio

  • Calcula el estadístico de prueba observado. Las variables mean_dt_pre y mean_dt_post de ejercicios anteriores están en tu espacio de nombres.
  • Desplaza los datos posteriores a 2010 para que tengan la misma media que los anteriores a 2010. Guarda el resultado como dt_post_shift.
  • Extrae 10.000 réplicas bootstrap de la media de dt_pre y de dt_post_shift. Guarda los resultados correspondientes en bs_reps_pre y bs_reps_post.
  • Calcula réplicas de la diferencia de medias restando bs_reps_post de bs_reps_pre.
  • Calcula e imprime el valor p. Considera «al menos tan extremo como» cuando el estadístico de prueba sea mayor o igual que el observado.

ejercicio interactivo práctico

Prueba este ejercicio completando este código de ejemplo.

# Compute the observed test statistic
mean_dt_diff = ____ - ____

# Shift the post-2010 data to have the same mean as the pre-2010 data
dt_post_shift = ____ - ____ + ____

# Compute 10,000 bootstrap replicates from arrays
bs_reps_pre = ____
bs_reps_post = ____

# Get replicates of difference of means
bs_reps = ____ - ____

# Compute and print the p-value
p_val = ____(____ >= ____) / 10000
print('p =', p_val)
Editar y ejecutar código