Kiểm định giả thuyết: b-value có khác nhau không?
Thực hiện kiểm định giả thuyết đã phác thảo ở bài trước. Các biến mags_pre và mags_post đã được nạp vào không gian làm việc của bạn, cùng với mt = 3.
Bài tập này là một phần của khóa học
Các nghiên cứu tình huống về tư duy thống kê
Hướng dẫn bài tập
- Lọc ra các độ lớn động đất trước năm 2010 có độ lớn lớn hơn (hoặc bằng) ngưỡng đầy đủ và ghi đè
mags_prebằng kết quả. Làm tương tự chomags_post. - Tính chênh lệch quan sát được giữa các độ lớn trung bình, bằng cách lấy trung bình sau 2010 trừ trung bình trước 2010.
- Tạo 10.000 mẫu lặp hoán vị bằng
dcst.draw_perm_reps(). Dùngdcst.diff_of_meanslàm đối số chofunc. - Tính và in p-value với cách hiểu "ít nhất là cực đoan như" nghĩa là thống kê kiểm định nhỏ hơn giá trị quan sát được.
Bài tập tương tác thực hành trực tiếp
Hãy thử làm bài tập này bằng cách hoàn thành đoạn mã mẫu này.
# 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)