3.3 Gepaarde t-toets

Als je goed naar de data kijkt (bijv. via View(sleep)), zie je dat dezelfde proefpersonen (weergegeven met ID) allebei de slaapmiddelen hebben getest. De waarnemingen zijn dus niet onafhankelijk, maar gepaard. Voor zulke datasets kan je de gepaarde t-toets gebruiken. Met deze toets bekijk je of gemiddeld de ene behandeling een hogere of lagere respons geeft dan de andere behandeling.

In dit geval staan de ID’s van de proefpersonen op volgorde, en kunnen we veilig de volgende code uitvoeren voor een gepaarde t-toets:

data(sleep)
t.test(sleep$extra~sleep$group, paired = TRUE)

Veiliger is om de twee behandelingen als aparte variabelen te nemen:

library(tidyverse)
data(sleep)
sleep <- data %>% 
  pivot_wider(names_from = group, 
              values_from = extra) %>% 
t.test(sleep$1, sleep$2, paired = TRUE)

met de functie pivot_wider() uit de package tidyverse kan je de variabele extra verdelen over twee variabelen op basis van groep.