Zacznij terazZacznij za darmo

Obliczanie wartości b

Wartość b to popularna miara sejsmiczności danego regionu. Pracując z danymi sejsmicznymi, zapewne będziesz obliczać ją wielokrotnie. Do zadań, które wykonujesz regularnie, warto napisać funkcję! Napisz więc funkcję o sygnaturze b_value(mags, mt, perc=[2.5, 97.5], n_reps=None), która zwraca wartość b oraz (opcjonalnie, jeśli n_reps nie jest None) jej przedział ufności dla zestawu magnitud mags. Próg kompletności jest podany przez mt. Argument kluczowy perc określa percentyle dla dolnej i górnej granicy przedziału ufności, a n_reps to liczba bootstrapowych replikacji używanych przy jego obliczaniu.

To ćwiczenie jest częścią kursu

Studia przypadków w myśleniu statystycznym

Zobacz kurs

Instrukcje do ćwiczenia

  • Zdefiniuj funkcję o sygnaturze b_value(mags, mt, perc=[2.5, 97.5], n_reps=None), która wykonuje następujące kroki:
    • Wyodrębnij z mags magnitudy równe lub wyższe od progu kompletności mt, używając indeksowania boolowskiego. Wynik zapisz w zmiennej m.
    • Oblicz najlepsze oszacowanie wartości b. Pamiętaj, że najlepsze oszacowanie to b = (m - mt)·ln(10). Wynik zapisz w zmiennej b.
    • Jeśli n_reps nie jest None, wykonaj poniższe kroki.
      • Wylosuj n_reps bootstrapowych replikacji średniej z m. Wynik zapisz w zmiennej m_bs_reps.
      • Przelicz bootstrapowe replikacje średniej z m na replikacje wartości b. Wynik zapisz w zmiennej b_bs_reps.
      • Oblicz przedział ufności na podstawie bootstrapowych replikacji wartości b. Wynik zapisz w zmiennej conf_int.
    • Zwróć b i conf_int, albo samo b, jeśli n_reps jest None.

Interaktywne ćwiczenie praktyczne

Spróbuj tego ćwiczenia, uzupełniając ten przykładowy kod.

____ ____:
    """Compute the b-value and optionally its confidence interval."""
    # Extract magnitudes above completeness threshold: m
    m = ____[____ >= ____]

    # Compute b-value: b
    b = ____

    # Draw bootstrap replicates
    if n_reps is None:
        return b
    else:
        m_bs_reps = ____

        # Compute b-value from replicates: b_bs_reps
        b_bs_reps = (____ - ____) * ____

        # Compute confidence interval: conf_int
        conf_int = ____
    
        return b, conf_int
Edytuj i uruchom kod