BaşlayınÜcretsiz Başlayın

b-değerini hesaplama

b-değeri, bir bölgenin sismikliği için yaygın bir ölçüttür. Deprem verileriyle çalışırken bunu sık sık hesaplamak isteyeceğini düşünebilirsin. Bu tür sık yapacağın işler için en iyisi bir fonksiyon yazmaktır! Bu yüzden, b_value(mags, mt, perc=[2.5, 97.5], n_reps=None) imzasına sahip ve b-değerini ve (isteğe bağlı olarak, n_reps None değilse) bir büyüklük kümesi mags için güven aralığını döndüren bir fonksiyon yaz. Tamlık eşiği mt ile veriliyor. perc anahtar argümanı, güven aralığının alt ve üst sınırları için yüzdelikleri belirtir ve n_reps, güven aralığını hesaplarken kullanılacak bootstrap tekrarlarının sayısıdır.

Bu egzersiz

İstatistiksel Düşünmede Vaka Çalışmaları

kursunun bir parçasıdır
Kursu Görüntüle

Egzersiz talimatları

  • Aşağıdakileri yapan b_value(mags, mt, perc=[2.5, 97.5], n_reps=None) imzasına sahip bir fonksiyon tanımla:
    • Boolean indeksleme kullanarak mags içinden tamlık eşiği mt ve üzerindeki büyüklükleri dilimle. Sonucu m değişkeninde sakla.
    • b-değerinin en iyi kestirimini hesapla. Unutma, b-değeri için en iyi kestirim b = (m - mt)·ln(10). Sonucu b değişkeninde sakla.
    • n_reps None değilse, şunları yap:
      • m ortalamasının n_reps adet bootstrap tekrarını çiz. Sonucu m_bs_reps değişkeninde sakla.
      • m ortalamasının bootstrap tekrarlarını b-değeri tekrarlarına dönüştür. Sonucu b_bs_reps değişkeninde sakla.
      • b-değeri bootstrap tekrarlarından güven aralığını hesapla. Sonucu conf_int değişkeninde sakla.
    • n_reps None ise sadece byi; değilse b ve conf_inti döndür.

Uygulamalı interaktif egzersiz

Bu örnek kodu tamamlayarak bu egzersizi bitirin.

____ ____:
    """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
Kodu Düzenle ve Çalıştır