Ước lượng thời gian giữa các trận động đất ở Parkfield
Trong bài tập này, bạn sẽ trước tiên tính các ước lượng tốt nhất cho tham số của hai mô hình Phân phối mũ (Exponential) và Chuẩn (Gaussian) cho thời gian giữa các trận động đất. Sau đó, bạn sẽ vẽ các CDF lý thuyết cho từng mô hình cùng với ECDF chính quy của dữ liệu thời gian giữa các trận động đất thực tế ở Parkfield.
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ời gian trung bình giữa các trận động đất và lưu vào
mean_time_gap. Các khoảng thời gian giữa những trận động đất lớn, theo năm, được lưu trongtime_gap. - Tính độ lệch chuẩn của các khoảng thời gian giữa các trận động đất và lưu vào
std_time_gap. - Dùng
np.random.exponential()để lấy 10.000 mẫu từ phân phối Mũ với trung bình phù hợp. Lưu vào biếntime_gap_exp. - Dùng
np.random.normal()để lấy 10.000 mẫu từ phân phối Chuẩn với trung bình và độ lệch chuẩn phù hợp. Lưu vào biếntime_gap_norm. - Vẽ các CDF lý thuyết, mỗi mô hình một đường, bằng cách dùng phương pháp
*dcst.ecdf()đã giới thiệu trước đó trong chương này. - Vẽ ECDF với các đối số từ khóa
formal=True,min_x=-10, vàmax_x=50.
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 mean time gap: mean_time_gap
mean_time_gap = ____
# Standard deviation of the time gap: std_time_gap
std_time_gap = ____
# Generate theoretical Exponential distribution of timings: time_gap_exp
time_gap_exp = ____
# Generate theoretical Normal distribution of timings: time_gap_norm
time_gap_norm = ____
# Plot theoretical CDFs
_ = plt.plot(*____)
_ = plt.plot(*____)
# Plot Parkfield ECDF
_ = plt.plot(*____(____, ____=____, ____=____, ____=____))
# Add legend
_ = plt.legend(('Exp.', 'Norm.'), loc='upper left')
# Label axes, set limits and show plot
_ = plt.xlabel('time gap (years)')
_ = plt.ylabel('ECDF')
_ = plt.xlim(-10, 50)
plt.show()