Test hipotezy: czy wartości b są różne?
Przeprowadź test hipotezy nakreślony w poprzednim ćwiczeniu. Zmienne mags_pre i mags_post są już załadowane do przestrzeni nazw, podobnie jak mt = 3.
To ćwiczenie jest częścią kursu
Studia przypadków w myśleniu statystycznym
Instrukcje do ćwiczenia
- Wyodrębnij magnitudy trzęsień ziemi sprzed 2010 roku, które mają wartość większą lub równą progowi kompletności, i nadpisz nimi
mags_pre. Zrób to samo dlamags_post. - Oblicz zaobserwowaną różnicę średnich magnitud, odejmując magnitudy trzęsień sprzed 2010 roku od tych z okresu po 2010 roku.
- Wygeneruj 10 000 replik permutacyjnych, używając
dcst.draw_perm_reps(). Jako argumentfuncpodajdcst.diff_of_means. - Oblicz i wyświetl wartość p, przyjmując, że „co najmniej tak ekstremalna jak" oznacza, że statystyka testowa jest mniejsza od zaobserwowanej.
Interaktywne ćwiczenie praktyczne
Spróbuj tego ćwiczenia, uzupełniając ten przykładowy kod.
# 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)