Poweranalyse Anteile in verbundenen Stichproben
Verfasst: Di Jun 13, 2017 12:39 pm
Hi!
Sagen wir, ich will Probanden mit und ohne Brille Bilder zeigen und beweisen, dass mit Brille mehr Bilder erkann werden, als ohne. Die Frage lautet, wieviele Probanden brauche ich, um das mit einer gegebenen Power (sagen wir 90%) berechnen zu können. Ich bin bereit anzunehmen, dass in der einen Messbedingung 40% der Bilder richtig erkannt werden, in der anderen 60%
Ich kenne die Funktion power.prop.test, die wie folgt aussieht:
Und die mir eine Fallzahl von 2 mal 129 Probanden vorschreibt.
Blöderweise ist das eine Power-Analyse für unverbundene Stichproben, ich will aber verbundene Stichproben untersuchen: Jeder Proband soll einmal mit und einmal ohne "Brille" Bilder erkennen. Ich erhoffe mir dadurch auch eine geringere Fallzahl.
Normalerweise stütze ich mich bei Poweranalysen gerne auf das pwr-Paket. Die Hilfe zu pwr::pwr.2p.test äußert sich aber nicht zur Frage der Verbundenheit von Stichproben und das Ergebnis enspricht dem von oben
Kann mir jemand mit einer Formel oder einer R Funktion für die Fallzahlanalyse helfen?
LG,
Bernhard
Sagen wir, ich will Probanden mit und ohne Brille Bilder zeigen und beweisen, dass mit Brille mehr Bilder erkann werden, als ohne. Die Frage lautet, wieviele Probanden brauche ich, um das mit einer gegebenen Power (sagen wir 90%) berechnen zu können. Ich bin bereit anzunehmen, dass in der einen Messbedingung 40% der Bilder richtig erkannt werden, in der anderen 60%
Ich kenne die Funktion power.prop.test, die wie folgt aussieht:
Code: Alles auswählen
power.prop.test( p1=.4, p2=.6, power=.9)
Blöderweise ist das eine Power-Analyse für unverbundene Stichproben, ich will aber verbundene Stichproben untersuchen: Jeder Proband soll einmal mit und einmal ohne "Brille" Bilder erkennen. Ich erhoffe mir dadurch auch eine geringere Fallzahl.
Normalerweise stütze ich mich bei Poweranalysen gerne auf das pwr-Paket. Die Hilfe zu pwr::pwr.2p.test äußert sich aber nicht zur Frage der Verbundenheit von Stichproben und das Ergebnis enspricht dem von oben
Code: Alles auswählen
pwr.2p.test(h=ES.h(0.40, 0.60), power=.9)
LG,
Bernhard