Le calcul de la puissance statistique est le pivot de toute planification expérimentale rigoureuse. Il permet de répondre à la question cruciale : « De combien de sujets ai-je besoin pour détecter un effet réel avec un niveau de confiance satisfaisant ? ». Ce script illustre l'utilisation de deux procédures majeures, PROC GLMPOWER pour les modèles structurés (ANOVA) et PROC POWER pour les tests d'hypothèses classiques.
Pourquoi cette démarche est-elle indispensable ?
Anticipation du Risque de Type II (β) : En visant des puissances de 0.7, 0.8 et 0.9, vous quantifiez le risque de ne pas détecter une différence qui existe pourtant réellement. Une puissance de 0.8 est le standard industriel, signifiant que vous avez 80 % de chances de rejeter l'hypothèse nulle si l'effet attendu est présent.
Adaptabilité au Modèle (GLM vs T-Test) :
PROC GLMPOWER est idéale lorsque vous avez des données pilotes (comme la table sleep). Elle utilise la structure de votre modèle linéaire pour extrapoler les besoins futurs.
PROC POWER (twosamplemeans) est l'outil de précision pour les tests de comparaison simples. En spécifiant une différence de moyennes (meandiff=3) et un écart-type (stddev=30), vous définissez la "taille de l'effet" que vous jugez cliniquement ou scientifiquement pertinente.
Visualisation Decisionnelle : L'instruction PLOT génère une courbe de puissance. C'est un outil de communication essentiel pour justifier un budget ou des ressources : elle montre visuellement le "coût" en effectifs pour chaque gain marginal de puissance.
Type : CREATION_INTERNE
Les données sont générées via une étape DATA avec des instructions DATALINES (valeurs 30 et 33).
| 1 | DATA sleep; |
| 2 | DO hours = 1 to 2; |
| 3 | INPUT recall; |
| 4 | OUTPUT; |
| 5 | END; |
| 6 | DATALINES; |
| 7 | 30 |
| 8 | 33 |
| 9 | ; |
| 10 | RUN; |
| 1 | PROC GLMPOWER DATA=sleep; |
| 2 | class hours; |
| 3 | model recall = hours; |
| 4 | power |
| 5 | stddev = 30 |
| 6 | ntotal = . |
| 7 | power = 0.7 0.8 0.9; |
| 8 | plot x=power min=0.7 max=0.95; |
| 9 | RUN; |
| 1 | PROC POWER; |
| 2 | twosamplemeans |
| 3 | meandiff=3 |
| 4 | alpha=0.05 |
| 5 | ntotal = . |
| 6 | power = 0.7 0.8 0.9 |
| 7 | stddev = 30 |
| 8 | ; |
| 9 | RUN; |