b-value की गणना
b-value किसी क्षेत्र की भूकंपीय सक्रियता का एक सामान्य मीट्रिक है. भूकंप डेटा के साथ काम करते समय आप इसे अक्सर निकालना चाहेंगे. ऐसे कार्यों के लिए जिन्हें आप बार-बार करेंगे, एक फंक्शन लिखना सबसे अच्छा होता है! तो, एक फंक्शन लिखिए जिसका सिग्नेचर b_value(mags, mt, perc=[2.5, 97.5], n_reps=None) हो, जो b-value और (वैकल्पिक रूप से, यदि n_reps None नहीं है) उसके confidence interval को magnitudes के सेट mags के लिए लौटाए. Completeness threshold mt से दिया गया है. perc कीवर्ड आर्ग्युमेंट confidence interval की निचली और ऊपरी सीमाओं के लिए percentiles देता है, और n_reps confidence interval की गणना में उपयोग होने वाले bootstrap replicates की संख्या है.
यह अभ्यास पाठ्यक्रम का हिस्सा है
Statistical Thinking के केस स्टडीज़
अभ्यास निर्देश
- एक फंक्शन परिभाषित करें जिसका सिग्नेचर
b_value(mags, mt, perc=[2.5, 97.5], n_reps=None)हो, जो निम्नलिखित करे:- Boolean इंडेक्सिंग का उपयोग करके
magsमें से completeness thresholdmtपर और उससे ऊपर की magnitudes को स्लाइस करें. परिणाम को वैरिएबलmमें स्टोर करें. - b-value का श्रेष्ठ आकलन निकालें. याद रखें, b-value का श्रेष्ठ आकलन है: b = (m - mt)·ln(10). परिणाम को वैरिएबल
bमें स्टोर करें. - यदि
n_repsNone नहीं है, तो निम्न करें.mके mean केn_repsbootstrap replicates ड्रॉ करें. परिणाम को वैरिएबलm_bs_repsमें स्टोर करें.mके mean के bootstrap replicates को b-value के replicates में बदलें. परिणामb_bs_repsमें स्टोर करें.- b-value के bootstrap replicates से confidence interval निकालें. परिणाम
conf_intमें स्टोर करें.
bऔरconf_intलौटाएँ, या केवलbलौटाएँ यदिn_repsNoneहै.
- Boolean इंडेक्सिंग का उपयोग करके
इंटरैक्टिव व्यावहारिक अभ्यास
इस अभ्यास को इस नमूना कोड को पूरा करके आज़माएँ।
____ ____:
"""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