Kiểm định giả thuyết cho tỷ lệ

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

Andrew Bray

Assistant Professor of Statistics at Reed College

Một nửa người Mỹ ủng hộ tử hình?

gss2016 %>%
  ggplot(aes(x = cappun)) +
  geom_bar()

p_hat <- gss2016 %>%
  summarize(mean(cappun == "FAVOR")) %>%
  pull()
p_hat

0.5666667

biểu đồ mật độ ch2v1

Một nửa người Mỹ ủng hộ tử hình?

null <- gss2016 %>%
  specify(
    response = cappun, 
    success = "FAVOR"
  ) %>%

  hypothesize(
    null = "point", 
    p = 0.5
  ) %>%

  generate(
    reps = 500, 
    type = "simulate"
  ) %>%

  calculate(stat = "prop")

A tibble: 500 x 2
   replicate  stat
   <fct>     <dbl>
 1 1         0.48 
 2 2         0.447
 3 3         0.48 
 4 4         0.44 
 5 5         0.407
 6 6         0.52 
 7 7         0.413
 8 8         0.553
 9 9         0.52 
10 10        0.467
# … with 490 more rows

Một nửa người Mỹ ủng hộ tử hình?

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

null %>%
  summarize(mean(stat > p_hat)) %>%
  pull() * 2

biểu đồ mật độ ch2v1

Kiểm định giả thuyết

Giả thuyết không: mô tả về trạng thái thực tế.
Phân phối theo giả thuyết không: phân phối của thống kê kiểm định khi giả thuyết không đúng.
p-value: mức độ phù hợp giữa giả thuyết không và quan sát.
- p-value cao: phù hợp (p-val > alpha)
- p-value thấp: không phù hợp (p-val < alpha)

Ayo berlatih!

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