Kiểm định giả thuyết: tần suất động đất có thay đổi không?
Rõ ràng, tần suất động đất đã tăng vọt sau khi bắt đầu bơm nước thải. Dù vậy, bạn vẫn sẽ thực hiện một phép kiểm định giả thuyết để luyện tập. Bạn sẽ không kiểm định giả thuyết rằng thời gian giữa các trận động đất có cùng phân phối trước và sau năm 2010, vì việc bơm nước thải có thể ảnh hưởng đến phân phối. Thay vào đó, bạn sẽ giả định chúng có cùng giá trị trung bình. Vì vậy, hãy tính p-value cho giả thuyết rằng thời gian giữa các trận động đất trước và sau năm 2010 có cùng trung bình, sử dụng hiệu giữa trung bình khoảng thời gian trước 2010 trừ trung bình khoảng thời gian sau 2010 làm thống kê kiểm định.
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
- Tính thống kê kiểm định quan sát được. Các biến
mean_dt_prevàmean_dt_posttừ các bài trước đã có trong không gian tên của bạn. - Dịch chuyển dữ liệu sau năm 2010 để có cùng trung bình với dữ liệu trước năm 2010. Lưu kết quả vào
dt_post_shift. - Vẽ 10.000 bản sao bootstrap cho mỗi trung bình của
dt_prevàdt_post_shift. Lưu lần lượt vàobs_reps_prevàbs_reps_post. - Tính các bản sao của hiệu trung bình bằng cách trừ
bs_reps_postkhỏibs_reps_pre. - Tính và in ra p-value. Xem "ít nhất là cực trị như vậy" nghĩa là thống kê kiểm định lớn hơn hoặc bằng 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.
# Compute the observed test statistic
mean_dt_diff = ____ - ____
# Shift the post-2010 data to have the same mean as the pre-2010 data
dt_post_shift = ____ - ____ + ____
# Compute 10,000 bootstrap replicates from arrays
bs_reps_pre = ____
bs_reps_post = ____
# Get replicates of difference of means
bs_reps = ____ - ____
# Compute and print the p-value
p_val = ____(____ >= ____) / 10000
print('p =', p_val)