Aan de slagGa gratis aan de slag

Hypothesetoets: zijn de b-waarden verschillend?

Voer de hypothesetoets uit die in de vorige oefening is uiteengezet. De variabelen mags_pre en mags_post zijn al in je namespace geladen, evenals mt = 3.

Deze oefening maakt deel uit van de cursus

Casestudies in statistisch denken

Cursus bekijken

Oefeninstructies

  • Selecteer de magnitudes van aardbevingen van vóór 2010 die groter zijn dan (of gelijk aan) de volledigheidsdrempel en overschrijf mags_pre met het resultaat. Doe hetzelfde voor mags_post.
  • Bereken het geobserveerde verschil in gemiddelde magnitudes, waarbij je de magnitudes van pre-2010-aardbevingen aftrekt van die van post-2010-aardbevingen.
  • Genereer 10.000 permutatiereplicaties met dcst.draw_perm_reps(). Gebruik dcst.diff_of_means als argument voor func.
  • Bereken en print de p-waarde, waarbij "ten minste zo extreem als" betekent dat de toetsingsgrootheid kleiner is dan wat werd geobserveerd.

Praktische interactieve oefening

Probeer deze oefening eens door deze voorbeeldcode in te vullen.

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