Začněte nyníZačněte zdarma

Výpočet hodnoty b

Hodnota b je běžnou metrikou seizmicity dané oblasti. Při práci se zemětřesnými daty ji budeš chtít počítat opakovaně. Pro takové opakující se úkoly je nejlepší napsat funkci! Napiš tedy funkci se signaturou b_value(mags, mt, perc=[2.5, 97.5], n_reps=None), která vrátí hodnotu b a (volitelně, pokud n_reps není None) její interval spolehlivosti pro sadu magnitud mags. Práh úplnosti je zadán hodnotou mt. Klíčový argument perc určuje percentily pro dolní a horní hranici intervalu spolehlivosti a n_reps udává počet bootstrapových replikátů použitých při jeho výpočtu.

Toto cvičení je součástí kurzu

Případové studie ve statistickém myšlení

Zobrazit kurz

Pokyny k cvičení

  • Definuj funkci se signaturou b_value(mags, mt, perc=[2.5, 97.5], n_reps=None), která provede následující:
    • Pomocí booleovského indexování vyber z mags magnitudy na úrovni prahu úplnosti mt a výše. Výsledek ulož do proměnné m.
    • Vypočítej nejlepší odhad hodnoty b. Pamatuj, že nejlepší odhad je b = (m - mt)·ln(10). Výsledek ulož do proměnné b.
    • Pokud n_reps není None, proveď následující:
      • Vytáhni n_reps bootstrapových replikátů průměru m. Výsledek ulož do proměnné m_bs_reps.
      • Převeď bootstrapové replikáty průměru m na replikáty hodnoty b. Výsledek ulož do b_bs_reps.
      • Z bootstrapových replikátů hodnoty b vypočítej interval spolehlivosti. Výsledek ulož do conf_int.
    • Vrať b a conf_int, nebo pouze b, pokud je n_reps rovno None.

Interaktivní cvičení na vyzkoušení si v praxi

Vyzkoušejte si toto cvičení dokončením tohoto ukázkového kódu.

____ ____:
    """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
Upravit a spustit kód