b-values का मात्रात्मक आकलन
पिछले अभ्यास में आपने जो प्लॉट बनाया था, उसके आधार पर आप mt = 3 का completeness threshold निश्चिंत होकर इस्तेमाल कर सकते हैं। इस threshold का उपयोग करके 1980 से 2009 के बीच और 2010 से 2017 के मध्य तक की अवधि के लिए b-values निकालें। पिछले चैप्टर में आपने जो b_value() फंक्शन लिखा था, जो मैग्नीट्यूड्स और completeness threshold के आधार पर b-value और confidence interval निकालता है, वह आपके namespace में उपलब्ध है। पिछले अभ्यास से numpy के ऐरे mags_pre और mags_post, और mt भी उपलब्ध हैं।
यह अभ्यास पाठ्यक्रम का हिस्सा है
Statistical Thinking के केस स्टडीज़
अभ्यास निर्देश
- 1980 से 2009 तक के भूकंपों के लिए 10,000 bootstrap replicates का उपयोग करके b-value और 95% confidence interval निकालें.
- 2010 से 2017 के मध्य तक के भूकंपों के लिए 10,000 bootstrap replicates का उपयोग करके b-value और 95% confidence interval निकालें.
- परिणाम स्क्रीन पर प्रिंट करने के लिए 'उत्तर सबमिट करें' दबाएँ.
इंटरैक्टिव व्यावहारिक अभ्यास
इस अभ्यास को इस नमूना कोड को पूरा करके आज़माएँ।
# Compute b-value and confidence interval for pre-2010
b_pre, conf_int_pre = ____(____, ____, ____=[____, ____], ____=____)
# Compute b-value and confidence interval for post-2010
b_post, conf_int_post = ____
# Report the results
print("""
1980 through 2009
b-value: {0:.2f}
95% conf int: [{1:.2f}, {2:.2f}]
2010 through mid-2017
b-value: {3:.2f}
95% conf int: [{4:.2f}, {5:.2f}]
""".format(b_pre, *conf_int_pre, b_post, *conf_int_post))