การประมาณค่าเฉลี่ยของช่วงเวลาระหว่างแผ่นดินไหว
จาก EDA เชิงกราฟิกในแบบฝึกหัดที่แล้ว จะเห็นได้ชัดว่าความถี่ของแผ่นดินไหวเปลี่ยนแปลงไปอย่างเห็นได้ชัดราวปี 2010 เพื่อเปรียบเทียบ ให้คำนวณค่าเฉลี่ยช่วงเวลาระหว่างแผ่นดินไหวขนาด 3 ขึ้นไปในช่วงปี 1980 ถึง 2009 และช่วงปี 2010 ถึงกลางปี 2017 พร้อมทั้งระบุช่วงความเชื่อมั่น 95% ของค่าเฉลี่ยด้วย โดยที่ตัวแปร dt_pre และ dt_post เก็บค่าช่วงเวลาระหว่างแผ่นดินไหวขนาดอย่างน้อย 3 ก่อนปี 2010 และหลังปี 2010 ตามลำดับ ในหน่วยวัน
แบบฝึกหัดนี้เป็นส่วนหนึ่งของหลักสูตร
กรณีศึกษาด้านการคิดเชิงสถิติ
คำแนะนำการฝึกหัด
- คำนวณค่าเฉลี่ยช่วงเวลาระหว่างแผ่นดินไหวสำหรับข้อมูลก่อนปี 2010 (
dt_pre) และหลังปี 2010 (dt_post) - สุ่ม bootstrap replicates ของค่าเฉลี่ยจำนวน 10,000 ครั้ง สำหรับชุดข้อมูลทั้งก่อนและหลังปี 2010
- ใช้
np.percentile()เพื่อคำนวณช่วงความเชื่อมั่น 95% ของค่าเฉลี่ยสำหรับชุดข้อมูลทั้งสอง - คลิก 'ส่งคำตอบ' เพื่อแสดงผลลัพธ์บนหน้าจอ
แบบฝึกหัดเชิงโต้ตอบแบบลงมือทำ
ลองทำแบบฝึกหัดนี้โดยเติมโค้ดตัวอย่างนี้ให้สมบูรณ์
# Compute mean interearthquake time
mean_dt_pre = ____
mean_dt_post = ____
# Draw 10,000 bootstrap replicates of the mean
bs_reps_pre = ____
bs_reps_post = ____
# Compute the confidence interval
conf_int_pre = ____
conf_int_post = ____
# Print the results
print("""1980 through 2009
mean time gap: {0:.2f} days
95% conf int: [{1:.2f}, {2:.2f}] days""".format(mean_dt_pre, *conf_int_pre))
print("""
2010 through mid-2017
mean time gap: {0:.2f} days
95% conf int: [{1:.2f}, {2:.2f}] days""".format(mean_dt_post, *conf_int_post))