15  Poweranalyse

De huidige cursus geeft alleen een basisinleiding tot poweranalyse (de analyse van het onderscheidend vermogen van een statistische toets). Package pwr biedt poweranalyse voor onder andere t-toetsen, proportietoetsen en chikwadraattoetsen. Het gebruik van package pwr is geen examenstof.

Bereken de power voor een t-toets met onafhankelijke steekproeven met een gegeven steekproefgrootte (n = 30) en effectgrootte (Cohen’s d = 0,5):

library(pwr)
pwr.t.test(n = 30, d = 0.5, sig.level = .05, type = "two.sample",
           alternative = "two.sided")

     Two-sample t test power calculation 

              n = 30
              d = 0.5
      sig.level = 0.05
          power = 0.4778965
    alternative = two.sided

NOTE: n is number in *each* group

Als we willen weten welke steekproefgrootte nodig is om een effect van d = 0,8 met 80% power te detecteren, laten we n weg (maar specificeren wel p = 0.8):

pwr.t.test(d = 0.3, power = 0.8, sig.level = .05, type = "two.sample",
           alternative = "two.sided")

     Two-sample t test power calculation 

              n = 175.3847
              d = 0.3
      sig.level = 0.05
          power = 0.8
    alternative = two.sided

NOTE: n is number in *each* group

Meer informatie over pwr.t.test is te krijgen door ?pwr.t.test te typen in de R console.