---
title: "R Laboratuvarı 15: İki grup karşılaştırması"
format: html
execute:
  echo: true
  warning: false
---

İki bağımsız kaplama yöntemi varsayılan eşit olmayan varyanslı Welch yöntemiyle karşılaştırılır.

```{r}
a <- c(.27,.30,.34,.29,.36,.31,.25,.33,.28,.38,.32,.29)
b <- c(.29,.41,.35,.46,.31,.38,.27,.43,.36,.34)
fit <- t.test(a, b, var.equal = FALSE)
sonuc <- data.frame(
  fark = mean(a) - mean(b),
  standart_hata = unname(abs(diff(fit$conf.int)) / (2 * qt(.975, fit$parameter))),
  serbestlik_derecesi = unname(fit$parameter),
  iki_yonlu_p = fit$p.value,
  GA_alt = fit$conf.int[1],
  GA_ust = fit$conf.int[2]
)
knitr::kable(sonuc, digits = 4, caption = "Welch iki örneklem karşılaştırması")
```

Eşleştirilmiş deney, bağımsız gruplar gibi değil çift içi farklar üzerinden analiz edilir.

```{r}
d <- c(-1.2, -0.5, -0.9, -1.5, -0.4, -0.8, -1.1, -0.7)
paired_fit <- t.test(d, mu = 0)
paired <- data.frame(
  ortalama_fark_um = mean(d),
  standart_hata_um = sd(d) / sqrt(length(d)),
  t = unname(paired_fit$statistic),
  serbestlik_derecesi = unname(paired_fit$parameter),
  GA_alt_um = paired_fit$conf.int[1],
  GA_ust_um = paired_fit$conf.int[2]
)
stopifnot(isTRUE(all.equal(mean(d), -0.8875)))
knitr::kable(paired, digits = 4, caption = "Eşleştirilmiş pürüzlülük farkları")
```

Aynı iki kusur oranı, mutlak ve göreli ölçeklerde ayrı bilgiler verir.

```{r}
p1 <- 0.08
p2 <- 0.15
ikili_etki <- data.frame(
  olcu = c("risk farkı", "risk oranı", "olasılık oranı"),
  deger = c(p1 - p2, p1 / p2, (p1 / (1 - p1)) / (p2 / (1 - p2)))
)
knitr::kable(ikili_etki, digits = 4, caption = "İkili sonuç etki ölçüleri")
```