Le test d’égalité teste l’hypothèse H0 : µ1 = µ2. La puissance d’un test est la probabilité de rejeter l’hypothèse H0 sans commettre une erreur, c’est-à-dire lorsque µ1 est effectivement différente de µ2.
Nous voulons donc calculer la puissance du test avec un nombre n d’individus grâce à la fonction power.t.test().
Exemple : dans une expérience nous avons un écart type de 1.7, une moyenne de 1 et on prend le seuil classique α=5%. On calcule la puissance du test avec 15 individus par groupe.
power.t.test(n=15, delta=1, sd=1.7, sig.level=0.05)$power
[1] 0.3430475
Si l’on décide de n’utiliser que 15 individus dans l’expérimentation, alors on aura 34% de chances de détecter une différence entre les moyennes des différents groupes, ce qui est peu.