ÎncepețiÎncepe gratuit

Calcularea valorii b

Valoarea b este o metrică frecvent utilizată pentru a caracteriza seismicitatea unei regiuni. Când lucrezi cu date seismice, vei dori să o calculezi des. Pentru sarcini repetitive ca aceasta, cel mai bine este să scrii o funcție! Scrie o funcție cu signatura b_value(mags, mt, perc=[2.5, 97.5], n_reps=None), care returnează valoarea b și (opțional, dacă n_reps nu este None) intervalul de încredere asociat, pentru un set de magnitudini mags. Pragul de completitudine este dat de mt. Argumentul cheie perc specifică percentilele pentru limitele inferioară și superioară ale intervalului de încredere, iar n_reps reprezintă numărul de replici bootstrap utilizate la calcularea acestuia.

Acest exercițiu face parte din cursul

Studii de caz în gândirea statistică

Vezi cursul

Instrucțiuni pentru exercițiu

  • Definește o funcție cu signatura b_value(mags, mt, perc=[2.5, 97.5], n_reps=None) care face următoarele:
    • Extrage din mags magnitudinile care ating sau depășesc pragul de completitudine mt, folosind indexare booleană. Stochează rezultatul în variabila m.
    • Calculează cea mai bună estimare a valorii b. Reamintește-ți că formula este b = (m - mt)·ln(10). Stochează rezultatul în variabila b.
    • Dacă n_reps nu este None, efectuează următorii pași:
      • Extrage n_reps replici bootstrap ale mediei lui m. Stochează rezultatul în variabila m_bs_reps.
      • Convertește replicile bootstrap ale mediei lui m în replici ale valorii b. Stochează rezultatul în b_bs_reps.
      • Calculează intervalul de încredere pe baza replicilor bootstrap ale valorii b. Stochează rezultatul în conf_int.
    • Returnează b și conf_int, sau doar b dacă n_reps este None.

Exercițiu interactiv practic

Încearcă acest exercițiu completând acest cod de exemplu.

____ ____:
    """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
Editează și rulează codul