Mật độ ngẫu nhiên hóa

Sử dụng 100 lần lặp giúp bạn hiểu cơ chế tráo hoán. Tuy nhiên, 100 là chưa đủ để quan sát đầy đủ dải giá trị có khả năng xảy ra cho chênh lệch tỷ lệ dưới giả thuyết không.

Hãy nhớ 4 bước suy luận. Đây là 4 bước bạn sẽ dùng cho mọi bài tập suy luận trong khóa học này và các khóa suy luận thống kê tiếp theo. Hãy dùng tên hàm để gợi nhớ quy trình phân tích.

specify để chỉ định biến phản hồi và biến giải thích.
hypothesize để tuyên bố giả thuyết không.
generate để tạo mẫu lại, tráo hoán, hoặc mô phỏng.
calculate để tính các thống kê tóm tắt.

Trong bài này, bạn sẽ lặp quy trình 1000 lần để thấy toàn bộ phân phối của các chênh lệch tỷ lệ dưới giả thuyết không.

Bài tập này là một phần của khóa học

Nền tảng Suy luận trong R

Xem khóa học

Hướng dẫn bài tập

Các gói dplyr, ggplot2, NHANES và infer đã được nạp sẵn cho bạn.

Tạo 1000 chênh lệch tỷ lệ bằng cách xáo trộn biến HomeOwn theo cú pháp của infer. Nhắc lại cú pháp infer:
- specify mối quan hệ quan tâm là HomeOwn theo Gender và một lần thành công trong ngữ cảnh này là sở hữu nhà, success = "Own".
- hypothesize rằng giả thuyết không đúng với null = "independence" (nghĩa là giới tính và sở hữu nhà không liên quan).
- generate 1000 phép tráo hoán; đặt reps là 1000.
- calculate thống kê stat = "diff in props" với thứ tự c("male", "female").
Chạy đoạn mã vẽ biểu đồ mật độ để tạo biểu diễn trực quan làm mịn của phân phối chênh lệch. Đường cong có hình dạng như thế nào?

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.

# Perform 1000 permutations
homeown_perm <- homes %>%
  # Specify HomeOwn vs. Gender, with `"Own" as success
  ___(___ ~ ___, success = "___") %>%
  # Use a null hypothesis of independence
  ___(___) %>% 
  # Generate 1000 repetitions (by permutation)
  ___(reps = ___, type = "permute") %>% 
  # Calculate the difference in proportions (male then female)
  ___(___, order = ___))

# Density plot of 1000 permuted differences in proportions
ggplot(homeown_perm, aes(x = stat)) + 
  geom_density()

Chỉnh sửa và Chạy Mã