Kiểm định hoán vị: wild type so với heterozygote
Kiểm định giả thuyết rằng phân phối độ dài bout của heterozygote và wild type là giống hệt nhau bằng kiểm định hoán vị.
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 hiệu trung bình (heterozygote trừ wild type) của các bộ dữ liệu thực tế và lưu vào biến
diff_means_exp. Hai mảngnumpylàbout_lengths_wtvàbout_lengths_hetđã có sẵn trong không gian tên của bạn. - Rút 10.000 bản sao hoán vị của hiệu trung bình bằng
dcst.draw_perm_reps(). Bạn cũng có thể dùng hàmdcst.diff_of_means(); lưu kết quả vàoperm_reps. - Tính p-value, với định nghĩa "ít nhất cũng cực đoan như vậy" là hiệu trung bình dưới giả thuyết không lớn hơn hoặc bằng giá trị quan sát thực nghiệm.
- In p-value ra màn hình.
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 difference of means: diff_means_exp
diff_means_exp = ____ - ____
# Draw permutation replicates: perm_reps
perm_reps = ____(____, ____,
____, size=____)
# Compute the p-value: p-val
p_val = ____(____ >= ____) / len(____)
# Print the result
print('p =', p_val)