1. Uczyć się
    2. /
  2. Courses
    3. /
  3. Studia przypadków w myśleniu statystycznym
Connected

Exercise

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.

Instrukcje

100 XP
  • 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.