ÎncepețiÎncepe gratuit

Test de ipoteză: s-a modificat frecvența cutremurelor?

Este evident că frecvența cutremurelor a crescut masiv după ce a început injectarea apelor reziduale. Cu toate acestea, vei efectua un test de ipoteză pentru practică. Nu vei testa ipoteza că timpii dintre cutremure au aceeași distribuție înainte și după 2010, deoarece injectarea apelor reziduale poate afecta distribuția. În schimb, vei presupune că au aceeași medie. Calculează, așadar, p-valoarea asociată ipotezei că timpii dintre cutremure de dinaintea și de după 2010 au aceeași medie, folosind ca statistică de test diferența dintre media timpilor de dinainte de 2010 și media celor de după 2010.

Acest exercițiu face parte din cursul

Studii de caz în gândirea statistică

Vezi cursul

Instrucțiuni pentru exercițiu

  • Calculează statistica de test observată. Variabilele mean_dt_pre și mean_dt_post din exercițiile anterioare sunt disponibile în spațiul tău de lucru.
  • Deplasează datele de după 2010 astfel încât să aibă aceeași medie cu datele de dinainte de 2010. Stochează rezultatul în dt_post_shift.
  • Extrage câte 10.000 de replici bootstrap ale mediei pentru dt_pre și, respectiv, dt_post_shift. Stochează rezultatele în bs_reps_pre și bs_reps_post.
  • Calculează replicile diferenței de medii scăzând bs_reps_post din bs_reps_pre.
  • Calculează și afișează p-valoarea. Consideră că „cel puțin la fel de extremă" înseamnă că statistica de test este mai mare sau egală cu valoarea observată.

Exercițiu interactiv practic

Încearcă acest exercițiu completând acest cod de exemplu.

# 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)
Editează și rulează codul