3.2 Onafhankelijke t-toets

Als voorbeeld weer de sleepdataset. Nu is de vraag of de gemiddelde hoeveelheid extra slaap verschilt tussen de groep met slaapmiddel 1 en met slaapmiddel 2.

Standaard voert t.test() de Welch t-test uit, waarbij niet aangenomen wordt dat er gelijke variantie is tussen de groepen. Heb je geen duidelijke aanwijzingen dat variantie gelijk is, dan is de Welch t-test de beste optie.

data(sleep)
t.test(sleep$extra~sleep$group)

## 
##  Welch Two Sample t-test
## 
## data:  sleep$extra by sleep$group
## t = -1.8608, df = 17.776, p-value = 0.07939
## alternative hypothesis: true difference in means between group 1 and group 2 is not equal to 0
## 95 percent confidence interval:
##  -3.3654832  0.2054832
## sample estimates:
## mean in group 1 mean in group 2 
##            0.75            2.33

Heb je wel duidelijke aanwijzingen dat de variantie gelijk is tussen bij beide groepen, dan kan je dat aangeven met een argument:

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

## 
##  Two Sample t-test
## 
## data:  sleep$extra by sleep$group
## t = -1.8608, df = 18, p-value = 0.07919
## alternative hypothesis: true difference in means between group 1 and group 2 is not equal to 0
## 95 percent confidence interval:
##  -3.363874  0.203874
## sample estimates:
## mean in group 1 mean in group 2 
##            0.75            2.33