การทดสอบสมมติฐาน: ค่า b-value แตกต่างกันหรือไม่?

ดำเนินการทดสอบสมมติฐานตามที่ได้วางแผนไว้ในแบบฝึกหัดก่อนหน้า ตัวแปร mags_pre และ mags_post รวมถึง mt = 3 ถูกโหลดไว้ใน namespace ของคุณแล้ว

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

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

ดูคอร์ส

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

กรองค่าขนาดของแผ่นดินไหวก่อนปี 2010 ที่มีขนาดมากกว่าหรือเท่ากับ completeness threshold แล้วเขียนทับ mags_pre ด้วยผลลัพธ์ที่ได้ จากนั้นทำเช่นเดียวกันกับ mags_post
คำนวณความแตกต่างของค่าเฉลี่ยขนาดแผ่นดินไหวที่สังเกตได้ โดยลบค่าเฉลี่ยของแผ่นดินไหวก่อนปี 2010 ออกจากค่าเฉลี่ยของแผ่นดินไหวหลังปี 2010
สร้าง permutation replicates จำนวน 10,000 ครั้งโดยใช้ dcst.draw_perm_reps() และกำหนด dcst.diff_of_means เป็นอาร์กิวเมนต์ของ func
คำนวณและแสดงค่า p-value โดยนิยาม "สุดขีดอย่างน้อยเท่ากัน" ว่าค่าสถิติทดสอบมีค่า น้อยกว่า ที่สังเกตได้

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

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

# Only magnitudes above completeness threshold
mags_pre = ____[____ >= ____]
mags_post = ____[____ >= ____]

# Observed difference in mean magnitudes: diff_obs
diff_obs = ____ - ____

# Generate permutation replicates: perm_reps
perm_reps = ____

# Compute and print p-value
p_val = ____(____ < ____) / 10000
print('p =', p_val)

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