การวัดค่า b-value

จากกราฟที่สร้างในแบบฝึกหัดก่อนหน้า สามารถใช้ค่าขีดความสมบูรณ์ mt = 3 ได้อย่างปลอดภัย ให้ใช้ค่านี้คำนวณค่า b สำหรับช่วงปี 1980 ถึง 2009 และช่วงปี 2010 ถึงกลางปี 2017 ฟังก์ชัน b_value() ที่เขียนไว้ในบทที่แล้ว ซึ่งคำนวณค่า b และช่วงความเชื่อมั่นจากชุดขนาดแผ่นดินไหวและค่าขีดความสมบูรณ์ พร้อมใช้งานใน namespace แล้ว รวมถึง numpy array mags_pre และ mags_post จากแบบฝึกหัดก่อนหน้า และค่า mt

แบบฝึกหัดนี้เป็นส่วนหนึ่งของหลักสูตร

กรณีศึกษาด้านการคิดเชิงสถิติ

ดูคอร์ส

คำแนะนำการฝึกหัด

คำนวณค่า b และช่วงความเชื่อมั่น 95% สำหรับแผ่นดินไหวในช่วงปี 1980 ถึง 2009 โดยใช้ bootstrap 10,000 ครั้ง
คำนวณค่า b และช่วงความเชื่อมั่น 95% สำหรับแผ่นดินไหวในช่วงปี 2010 ถึงกลางปี 2017 โดยใช้ bootstrap 10,000 ครั้ง
คลิก 'ส่งคำตอบ' เพื่อแสดงผลลัพธ์บนหน้าจอ

แบบฝึกหัดเชิงโต้ตอบแบบลงมือทำ

ลองทำแบบฝึกหัดนี้โดยเติมโค้ดตัวอย่างนี้ให้สมบูรณ์

# 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))

แก้ไขและรันโค้ด