Kiểm định hoán vị trên dữ liệu loài ếch
Lực đập trung bình của Ếch A là 0,71 Newton (N), và của Ếch B là 0,42 N, chênh lệch 0,29 N. Có thể hai con ếch thực ra đập với cùng một lực và chênh lệch quan sát được chỉ là do ngẫu nhiên. Bạn sẽ tính xác suất thu được chênh lệch ít nhất 0,29 N trong lực đập trung bình dưới giả thuyết rằng phân phối lực đập của hai con ếch là giống hệt nhau. Ta dùng kiểm định hoán vị với thống kê kiểm định là hiệu trung bình để kiểm tra giả thuyết này.
Để thuận tiện, dữ liệu đã được lưu trong các mảng force_a và force_b.
Bài tập này là một phần của khóa học
Tư duy Thống kê với Python (Phần 2)
Hướng dẫn bài tập
- Định nghĩa một hàm với chữ ký gọi
diff_of_means(data_1, data_2)trả về hiệu giữa trung bình của hai tập dữ liệu, trung bình củadata_1trừ trung bình củadata_2. - Dùng hàm này để tính chênh lệch trung bình thực nghiệm đã quan sát ở hai con ếch.
- Vẽ 10.000 bản lặp hoán vị của hiệu trung bình.
- Tính p-value.
- In p-value.
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.
def diff_of_means(data_1, data_2):
"""Difference in means of two arrays."""
# The difference of means of data_1, data_2: diff
diff = ____
return diff
# Compute difference of mean impact force from experiment: empirical_diff_means
empirical_diff_means = ____
# Draw 10,000 permutation replicates: perm_replicates
perm_replicates = draw_perm_reps(____, ____,
____, size=10000)
# Compute p-value: p
p = np.sum(____ >= ____) / len(____)
# Print the result
print('p-value =', p)