Farklar için aralıklar

R ile Kategorik Verilerde Çıkarım

Andrew Bray

Assistant Professor of Statistics at Reed College

İki değişkenli bir soru

Kadınlar ve erkekler farklı oranlarda mı inanıyor?

p, ölüm sonrası yaşama inananların oranı olsun.

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

Kadınlar ve erkeklerin ölüm sonrası yaşama dair görüşleri farklı mı?

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

Kadınlar ve erkeklerin ölüm sonrası yaşama dair görüşleri farklı mı?

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

Kadınlar ve erkeklerin ölüm sonrası yaşama dair görüşleri farklı mı?

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

d_hat <- diff(p_hats)
d_hat

0.1472851

H0’dan veri üretme

$H_{0} : p_{female} - p_{male} = 0$
Ahiret inancı ile deneğin cinsiyeti arasında ilişki yoktur.
postlife değişkeni sex değişkeninden bağımsızdır.

⇒ Permütasyonla veri üretin

Kadınlar ve erkeklerin ölüm sonrası yaşama dair görüşleri farklı mı?

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

Kadınlar ve erkeklerin ölüm sonrası yaşama dair görüşleri farklı mı?

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

Kadınlar ve erkeklerin ölüm sonrası yaşama dair görüşleri farklı mı?

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

Kadınlar ve erkeklerin ölüm sonrası yaşama dair görüşleri farklı mı?

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.

Kadınlar ve erkeklerin ölüm sonrası yaşama dair görüşleri farklı mı?

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

Bu veriler, ölüm sonrası inançta cinsiyetler arasında fark olduğunu gösterir.

Hadi pratik yapalım!

R ile Kategorik Verilerde Çıkarım