Aan de slagGa gratis aan de slag

Hypothesetoets: veranderde de aardbevingsfrequentie?

Er was duidelijk een enorme toename in de aardbevingsfrequentie zodra de injectie van afvalwater begon. Toch ga je voor de oefening een hypothesetoets uitvoeren. Je test niet de hypothese dat de tijden tussen aardbevingen vóór en na 2010 dezelfde verdeling hebben, omdat afvalwaterinjectie de verdeling kan beïnvloeden. In plaats daarvan ga je ervan uit dat ze hetzelfde gemiddelde hebben. Bereken dus de p-waarde die hoort bij de hypothese dat de tijden tussen aardbevingen vóór en na 2010 hetzelfde gemiddelde hebben, waarbij je het gemiddelde van de tijdsintervallen vóór 2010 min het gemiddelde van de tijdsintervallen na 2010 gebruikt als toetsstatistiek.

Deze oefening maakt deel uit van de cursus

Casestudies in statistisch denken

Cursus bekijken

Oefeninstructies

  • Bereken de geobserveerde toetsstatistiek. De variabelen mean_dt_pre en mean_dt_post uit vorige oefeningen staan in je namespace.
  • Verschuif de data na 2010 zodat ze hetzelfde gemiddelde hebben als de data vóór 2010. Sla het resultaat op als dt_post_shift.
  • Trek 10.000 bootstrap-replicaties van het gemiddelde van zowel dt_pre als dt_post_shift. Sla de respectieve resultaten op in bs_reps_pre en bs_reps_post.
  • Bereken replicaties van het verschil in gemiddelden door bs_reps_post af te trekken van bs_reps_pre.
  • Bereken en print de p-waarde. Beschouw "ten minste zo extreem als" als: de toetsstatistiek is groter dan of gelijk aan wat is geobserveerd.

Praktische interactieve oefening

Probeer deze oefening eens door deze voorbeeldcode in te vullen.

# 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)
Code bewerken en uitvoeren