Khoảng tin cậy cho chênh lệch

Suy luận cho dữ liệu phân loại trong R

Andrew Bray

Assistant Professor of Statistics at Reed College

Câu hỏi với hai biến

Phụ nữ và nam giới có tỷ lệ tin khác nhau?

Gọi $p$ là tỷ lệ tin vào sự sống sau khi chết.

$H_{0} : p_{female} - p_{male} = 0$
$H_{A} : p_{female} - p_{male} \ne 0$

Phụ nữ và nam giới có ý kiến khác về sự sống sau khi chết?

ggplot(gss2016, aes(x = sex, fill = postlife)) +
  geom_bar()

Biểu đồ cột tỷ lệ postlife theo giới

Phụ nữ và nam giới có ý kiến khác về sự sống sau khi chết?

ggplot(gss2016, aes(x = sex, fill = postlife)) +
  geom_bar(position = "fill")

Biểu đồ cột chồng theo tỷ lệ postlife theo giới

Phụ nữ và nam giới có ý kiến khác về sự sống sau khi chết?

p_hats <- gss2016 %>%
  group_by(sex) %>%
  summarize(mean(postlife == "YES", na.rm = TRUE)) %>%
  pull()

d_hat <- diff(p_hats)
d_hat

0.1472851

Tạo dữ liệu từ H0

$H_{0} : p_{female} - p_{male} = 0$
Không có liên hệ giữa niềm tin vào thế giới bên kia và giới tính.
Biến postlife độc lập với biến sex.

⇒ Tạo dữ liệu bằng hoán vị

Phụ nữ và nam giới có ý kiến khác về sự sống sau khi chết?

gss2016 %>%
  specify(
    response = postlife, 
    explanatory = sex, 
    success = "YES"
  ) %>%
  hypothesize(null = "independence") %>%
  generate(reps = 1, type = "permute")

Phụ nữ và nam giới có ý kiến khác về sự sống sau khi chết?

gss2016 %>%
  specify(
    postlife ~ sex,  # this line is new
    success = "YES"
  ) %>%
  hypothesize(null = "independence") %>%
  generate(reps = 1, type = "permute")

Response: postlife (factor)
Explanatory: sex (factor)
Null Hypothesis:  independence 
# A tibble: 137 x 3
# Groups:   replicate [1]
   postlife sex    replicate
   <fct>    <fct>      <int>
 1 YES      FEMALE         1
 2 YES      MALE           1
 3 YES      FEMALE         1
 4 YES      MALE           1
 5 YES      MALE           1
 6 YES      FEMALE         1
 7 NO       FEMALE         1

Phụ nữ và nam giới có ý kiến khác về sự sống sau khi chết?

gss2016 %>%
  specify(
    postlife ~ sex, 
    success = "YES"
  ) %>%
  hypothesize(null = "independence") %>%
  generate(reps = 1, type = "permute")

Response: postlife (factor)
Explanatory: sex (factor)
Null Hypothesis:  independence 
# A tibble: 137 x 3
# Groups:   replicate [1]
   postlife sex    replicate
   <fct>    <fct>      <int>
 1 YES      FEMALE         1
 2 NO       MALE           1
 3 NO       FEMALE         1
 4 YES      MALE           1
 5 YES      MALE           1
 6 YES      FEMALE         1
 7 YES      FEMALE         1

Phụ nữ và nam giới có ý kiến khác về sự sống sau khi chết?

gss2016 %>%
  specify(postlife ~ sex, success = "YES") %>%
  hypothesize(null = "independence") %>%
  generate(reps = 500, type = "permute") %>%
  calculate(stat = "diff in props", order = c("FEMALE", "MALE"))

Warning message:
Removed 13 rows containing missing values.

Phụ nữ và nam giới có ý kiến khác về sự sống sau khi chết?

ggplot(null, aes(x = stat)) +
  geom_density() +
  geom_vline(xintercept = d_hat, color = "red")

Dữ liệu cho thấy có khác biệt theo giới trong niềm tin về sự sống sau khi chết.

Biểu đồ mật độ cho thống kê hoán vị, với đường đỏ tại d_hat

Ayo berlatih!

Suy luận cho dữ liệu phân loại trong R