Planung des Stichprobenumfangs

Im Zuge eines Ethikantrags stellt sich immer wieder die Frage - Wie erstelle ich eine Fallzahlberechnung?

Ester Schritt ist die Formulierung der Forschungsfrage und die Klärung des Studientyps. Hier wird zwischen Forschung an Primärdaten und Forschung Sekundärdaten unterschieden. Wenn eine Forschung an Sekundärdaten (systematische Reviews und Metaanalysen) vorliegt brauchen wir keine Fallzahlberechnung. Es genügt die Zahlen aus Bortz (2006) Seite 672 verwenden. Mehr zum Thema finden sich im Artikel von Röhrig, Prel, and Blettner (2009).

Bei einer Forschung an Primärdaten (Experiment, klinischen Prüfungen oder epidemiologischen Forschung) ist er zweite Schritt die Definition der Zielgröße und die Bestimmung der Verteilungseigenschaft. Also ist es eine binäre Variable wie ZB. Krank/Gesund, sind es gezählte Daten wie die Anzahl an Tagen bis zum Eintreten eines Ereignisses oder sind es Messwerte aus einem Fragebogen oder einer Laboranalyse. Aus dieser Information lässt sich das geplante statistische Verfahren ableiten.

Dritter und wichtigster Punkt ist die Abschätzung des erwarteten Effekts also der Änderung durch die Intervention. Die Effektstärke ist dabei von entscheidender Bedeutung und sie wird entweder aus Literatur abgeleitet oder der Wissenschaftler schätzt sie aus seiner Erfahrung. Beim Abschätzen der Effektgröße genügt eine Verbale Beschreibung und zwar entweder mit kleiner Effekt, mittlerer Effekt oder großer Effekt. Die Einteilung in klein - mittel - groß stammt von Cohen (1988), nach seiner Daumenregel lassen sich die Verbalen Beschreibungen in Zahlen (Effektstärken) umwandeln.

Zur Berechnung der Fallzahl werden noch zwei Zahlen benötigt, die zwar theoretisch der Wissenschaftler festlegt aber in der Praxis vorgegeben sind den Alpha-Fehler und die Power Der Alpha-Fehler (Signifikanz Level oder Fehler Type I) wird meist mit 5% gewählt und als Power (Fehler Type II) wird 80% gewählt.

Nicht zu vernachlässigen ist auch die erwartete Drop-out Rate festzulegen. Die tatsächliche Stichprobe errechnet sich dann durch die Formel \(N=\frac{n}{1-d}\) (mit d Drop-out Rate).

# Berechnung mit R und dem Packet pwr
# library(pwr)
 
#Conventional effect size from Cohen (1982)
(efct<-cohen.ES(test=c("t"), "medium")$effect.size)

## [1] 0.5

(stchprb<-pwr.t.test(d=efct, 
                 power=0.80, 
                 sig.level=0.05, 
                 type="two.sample" ))

## 
##      Two-sample t test power calculation 
## 
##               n = 63.76561
##               d = 0.5
##       sig.level = 0.05
##           power = 0.8
##     alternative = two.sided
## 
## NOTE: n is number in *each* group