Analyse de puissance dans [R] pour Anova bidirectionnelle

Question

J'essaie de calculer la taille d'échantillon nécessaire pour un plan factoriel 2x2. J'ai deux questions.

1) J'utilise le paquet pwr et une ANOVA fonction pour calculer la taille de l'échantillon nécessaire en utilisant le code suivant

pwr.anova.test(k = , n = , f = , sig.level = , power =) 

Cependant, je voudrais examiner deux ANOVA façon, puisque c'est plus efficace à l'estimation des moyens de groupe que oneova anova. Il n'y a pas de fonction anova bidirectionnelle que j'ai pu trouver. Existe-t-il un paquet ou une routine dans [R] pour faire cela?

2) De plus, suis-je sûr de supposer que puisque j'utilise des calculs de puissance anova à sens unique, que la taille de l'échantillon sera plus conservatrice (c'est-à-dire plus grande)?

Answer 1

Dans une ANOVA 2x2 impliquant le facteur A, le facteur B et AxB, vous obtiendrez des estimations de puissance statistique distinctes pour chacun de ces trois effets.

G Puissance 3 fournit des logiciels libres et des tutoriels clairs pour estimer la puissance des effets dans les conceptions factorielles: http://www.psycho.uni-duesseldorf.de/abteilungen/aap/gpower3/user-guide-by-design

Answer 2

Après la recherche - je ne pouvais pas trouver de solution pour cette ligne.

Ce que je vous suggère de faire (si vous savez comment) est de programmer ceci en utilisant une simulation. Si vous ne savez pas comment le faire, alors écrivez une question SO sur "Comment puis-je écrire une simulation de deux anova, pour réaliser l'analyse de puissance" et voir ce que les gens pourraient vous aider :)

en outre, vous pouvez commencer par examiner le code ici:

http://www.rforge.net/doc/packages/NCStats/power.sim.html

Pour commencer le calcul de la puissance grâce à la simulation. Notez ce que Jeromy a écrit - cette analyse de puissance est pour de multiples résultats. Sujet intéressant - J'adorerais faire un suivi là-dessus.

Best,

Tal